CN116882590B - 基于时序网络模型的航班计划优化方法、***及储存介质 - Google Patents
基于时序网络模型的航班计划优化方法、***及储存介质 Download PDFInfo
- Publication number
- CN116882590B CN116882590B CN202311127516.8A CN202311127516A CN116882590B CN 116882590 B CN116882590 B CN 116882590B CN 202311127516 A CN202311127516 A CN 202311127516A CN 116882590 B CN116882590 B CN 116882590B
- Authority
- CN
- China
- Prior art keywords
- airport
- time sequence
- network
- flight
- node
- 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
- 238000000034 method Methods 0.000 title claims abstract description 70
- 238000005457 optimization Methods 0.000 title claims abstract description 38
- 238000003860 storage Methods 0.000 title claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims description 56
- 238000004891 communication Methods 0.000 claims description 18
- 238000011156 evaluation Methods 0.000 claims description 16
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000005315 distribution function Methods 0.000 claims description 8
- 238000012216 screening Methods 0.000 claims description 8
- 238000010183 spectrum analysis Methods 0.000 claims description 8
- 238000004088 simulation Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 4
- 230000010006 flight Effects 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 2
- 238000007781 pre-processing Methods 0.000 claims description 2
- 238000010276 construction Methods 0.000 abstract description 5
- 230000003068 static effect Effects 0.000 description 6
- 230000002123 temporal effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012916 structural analysis Methods 0.000 description 2
- 101100219189 Caenorhabditis elegans byn-1 gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06314—Calendaring for a resource
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/067—Enterprise or organisation modelling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Quality & Reliability (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Game Theory and Decision Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Educational Administration (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于时序网络模型的航班计划优化方法、***及储存介质，首先准确的对中国机场网络的网络结构和抗毁性进行评估，为提高机场网络的抗毁性提供依据，为机场网络的规划建设和航班计划的确定提供数据支持；其次通过输入全国预先航班计划，运用时序网络模型进行预先推演，在航路点和机场设置相应阀值，对按航班计划预先推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警；最后在时序网络模型内对航班计划的航班时刻进行优化，达到对航路和机场的拥堵情况发生预发现、预优化，从而提高航路和机场资源的利用率。解决了现在航班计划编制时无法做到对航路拥堵和机场超负荷运行的预先识别和预先优化问题。
Description
技术领域
本发明涉及航班计划优化技术领域，具体涉及一种基于时序网络模型的航班计划优化方法、***及储存介质。
背景技术
机场网、道路网以及电力网实际上并不是静态的，会随着多种原因而进行改变，例如突发事件、人类作息时间等，并且增加或修改航线（道路）均会改变网络的拓扑结构。时序网络方法将关于事件的时间信息从动态***移动到网络，即动态发生的底层结构，可以更加准确的描述真实***的特性。
在航空网络的研究中，以往仅是将航空网络抽象成静态网络进行研究，在静态网络中，节点的连接忽略了接触的持续时间、接触间隔时间、重复接触和接触的时间顺序，且无法对航班计划进行预演和优化。航空业的快速发展对航空机场网络效率和安全提出了更高的要求，迫切需要一种加有效的方法对航空机场网络的进行评估和解决现在航班计划编制时无法做到对航路拥堵和机场超负荷运行的预先识别和预先优化问题。
发明内容
为解决现有技术中存在的问题，本发明提供了一种基于时序网络模型的航班计划优化方法、***及储存介质，能更准确的对机场网络的网络结构和抗毁性进行评估，为提高机场网络的抗毁性提供依据，为机场网络的规划建设和航班计划的确定提供数据支持，弥补了目前运用时序网络模型对机场网络进行分析的空白，解决了现在航班计划编制时无法做到对航路拥堵和机场超负荷运行的预先识别和预先优化问题，做到预发现、预优化，有较强的创新性和实际意义。
为实现上述目的，本发明提供如下技术方案：一种基于时序网络模型的航班计划优化方法，包括如下步骤：
S101、采集航班数据，并运用功率谱分析方法分析航班数据的周期性，确定航班数据固有周期；
S102、基于时序网络模型构建机场多层时序网络模型；
S103、构建机场时序网络模型拓扑结构分析指标，包括时序度、平均时序距离和节点中心性指标，并利用拓扑结构分析指标对机场时序网络模型结构进行分析；
S104、构建机场时序网络模型抗毁性分析指标，包括时序网络效率和时序最大联通子图，并利用抗毁性分析指标对机场时序网络模型进行抗毁性分析；
S105、机场网络抗毁性评估：采用蓄意攻击、随机攻击和非完全信息攻击的策略对机场网络进行模拟攻击，得到机场时序网络模型网络效率损失曲线，对得到的网络效率损失曲线进行分析，识别其拐点，完成对机场网络抗毁性的评估；
S106、输入预先航班计划，运用时序网络模型进行预先推演，在航路点和机场设置相应阀值，对按航班计划预先推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警，并在时序网络模型内对航班计划时刻进行优化，输出优化后航班计划。
优选的，在步骤S101中，所述的航班数据包含民航航班信息***中的民航航班数据库文档数据以及航班计划文档数据。
优选的，在步骤S102中，所述基于时序网络模型构建机场多层时序网络模型具体包括如下：
S201、构建机场时序网络邻接矩阵，构建表示机场时序网络，其中/>为网络中所包含的机场节点集，/>，/>为机场间的航班连边集合，，其中i表示起始节点即起飞机场（以机场编号表示）j表示目标节点即到达机场（以机场编号表示），/>是航班运行时间集合，,其中/>表示从节点i到节点j航班起飞时刻，/>表示从节点i到节点j航班运行时间；
S202、将时序网络中航班数据的总时间 按其固有周期划分为相应的层数，在每层中将对应的时序数据增加到网络中，构建机场多层时序网络模型。
优选的，在步骤S103中，利用时序度分析指标对网络结构中的时序度进行分析，包括时序入度，时序出度，时序度，以及幂律性的分析，使用幂律分布函数为，其中/>为时序入度幂律分布函数待求的三个参数，/>为节点时序入度值，/>为时序入度值为k的机场节点的概率值，/>为时序出度幂律分布函数待求的三个参数，/>为节点时序出度值，/>为时序出度值为k的机场节点的概率值，并对数据进行拟合求出待求参数；利用节点中心性指标对网络结构进行分析，所述节点中心性指标包括时序度中心性，时序介数中心性以及时序接近度中心性。
优选的，在步骤S105中，具体包括如下步骤：
S501、随机攻击策略为在机场节点集合中随机选择攻击顺序并形成攻击顺序集合其中，n为受攻击机场节点编号，v为攻击次数， 1attack代表第1次攻击，2attack代表第2次攻击，vattack代表第v次攻击，/>为攻击类型，也表示随机攻击，并按集合顺序依次删除网络中对应节点并计算网络效率和最大联通子图；
S502、蓄意攻击策略分为按时序度强度排序攻击、按时序介数中心性排序攻击、按时序接近度中心性强度攻击，分别形成时序度攻击顺序集合其中，n为受攻击机场节点编号，v为攻击次数， 1attack代表第1次攻击，2attack代表第2次攻击，vattack代表第v次攻击，/>表示时序度强度排序攻击，时序介数攻击顺序集合，其中，n为受攻击机场节点编号，v为攻击次数， 1attack代表第1次攻击，2attack代表第2次攻击，vattack代表第v次攻击，p表示时序介数值排序攻击，时序接近度攻击顺序集合/>，其中，n为受攻击机场节点编号，v为攻击次数， 1attack代表第1次攻击，2attack代表第2次攻击，vattack代表第v次攻击，其中b表示时序接近度值排序攻击，并按顺序集合依次删除节点并计算网络效率和最大联通子图；
S503、非完全信息攻击的攻击方式为在给定已知信息的情况下对已知信息的节点集合运用蓄意攻击策略进行攻击，对未知信息的节点运用随机攻击的策略，并按先已知后未知节点的顺序形成非完全信息攻击顺序集合，其中，n为受攻击机场节点编号，v为攻击次数， 1attack代表第1次攻击，2attack代表第2次攻击，vattack代表第v次攻击，m表示非完全信息攻击,最后按顺序集合依次删除节点并计算网络效率和最大联通子图，通过步骤S501-S503，得到机场时序网络模型网络效率损失曲线；
S504、对得到的网络效率损失曲线进行分析，识别其拐点，完成对机场网络抗毁性的评估。
优选的，在步骤S106中，具体包括如下步骤：
S601、将航班计划数据进行导入，对数据进行筛选，筛除空数据并筛选航班编号、航班起飞时间、飞行时间、起飞机场和目的地机场关键数据，基于处理过后的数据将带有时间属性的数据添加入邻接矩阵，矩阵中每个元素不再用简单的0，1表示是否存在连边，而是用列表代替，并将节点对间的航班时间数据添加入列表中进行储存形成航班计划时序邻接矩阵，导入航路信息和机场信息，并设定航路点和机场阀值；
S602、运用时序网络模型进行预先推演，对按航班计划推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警；
S603、在时序网络模型内对航班计划的航班时刻进行优化，输出优化后航班计划。
另一方面，为实现上述目的，本发明还提供了如下技术方案：一种基于时序网络模型的航班计划优化***，所述***包括：
航班数据采集模块：用于采集航班数据，并运用功率谱分析方法分析航班数据的周期性，确定航班数据固有周期；
多层时序网络模型构建模块：用于基于时序网络模型构建机场多层时序网络模型；
网络模型拓扑结构分析模块：用于构建机场时序网络模型拓扑结构分析指标，包括时序度、平均时序距离和节点中心性指标，并利用拓扑结构分析指标对机场时序网络模型结构进行分析；
网络模型抗毁性分析模块：用于构建机场时序网络模型抗毁性分析指标，包括时序网络效率和时序最大联通子图，并利用抗毁性分析指标对机场时序网络模型进行抗毁性分析；
抗毁性评估模块：用于对机场网络进行抗毁性评估，采用蓄意攻击、随机攻击和非完全信息攻击的策略对机场网络进行模拟攻击，得到机场时序网络模型网络效率损失曲线，对得到的网络效率损失曲线进行分析，识别其拐点，完成对机场网络抗毁性的评估；
航班计划优化模块：用于输入预先航班计划，运用时序网络模型进行预先推演，在航路点和机场设置相应阀值，对按航班计划预先推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警，并在时序网络模型内对航班计划时刻进行优化，输出优化后航班计划。
另一方面，为实现上述目的，本发明还提供了如下技术方案：一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理执行时实现所述的航班计划优化方法。
本发明的有益效果是：
1）首先本发明能更准确的对机场网络的网络结构和抗毁性进行评估，为提高机场网络的抗毁性提供依据，为机场网络的规划建设和航班计划的确定提供数据支持。其次通过输入预先航班计划，运用时序网络模型进行预先推演，在航路点和机场设置相应阀值，对按航班计划预先推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警，最后在时序网络模拟内对航班计划的航班时刻进行优化，达到对航路和机场的拥堵情况发生预发现、预优化，从而提高航路和机场资源的利用率。弥补目前运用时序网络模型对机场网络进行分析的空白，解决现在航班计划编制时无法做到对航路拥堵和机场超负荷运行的预先识别和预先优化问题，做到预发现、预优化，有较强的创新性和实际意义。
2）本发明主要基于时序网络模型，将航班运行数据或航班计划数据加入网络进行分析，数据库数据集相比静态网络包含数据量大，数据特征多，结构分析更加精确。通过对航班运行数据进行分析，可以更加精确的分析机场网络结构特征并对机场网络的抗毁性进行评价。通过对航班计划数据进行分析，运用时序网络模型进行预先推演，在航路和机场设置相应阀值，对按航班计划推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警，并可在时序网络模型内对航班计划时刻进行优化，可以缓解航路和机场的拥堵发生，具有很高的创新和实际意义。解决了现在解决航路、机场拥堵只能被动调整，而不能做到预发现、预优化的问题。
附图说明
图1为航班数据处理和功率谱分析方法的流程图；
图2为多层时序网络模型结构示意图；
图3为机场网络结构分析和抗毁性分析流程图；
图4为基于时序网络模型的航班计划预演及航路、机场拥堵识别和预警流程图；
图5为本发明航班计划优化***模块示意图；
图中，110-航班数据采集模块；120-多层时序网络模型构建模块；130-网络模型拓扑结构分析模块；140-网络模型抗毁性分析模块；150-抗毁性评估模块；160-航班计划优化模块。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参阅图1-图4，本发明提供一种技术方案：一种基于时序网络模型的航班计划优化方法，具体流程包括：
（1）采集航班数据，如图1所示，并运用功率谱分析方法分析航班数据的周期性，确定航班数据固有周期；
（2）构建机场多层时序网络模型，多层时序网络结构如图2所示；
（3）构建机场时序网络模型拓扑结构分析指标，包括时序入度，时序出度，时序度，平均时序距离和机场节点中心性指标：时序度中心性，时序介数中心性，时序接近度中心性；
（4）构建机场时序网络模型抗毁性分析指标，并选取时序网络效率，和时序最大联通子图作为机场时序网络模型抗毁性分析指标；
（5）采用蓄意攻击，随机攻击和非完全信息攻击对机场网络进行模拟攻击，对得到的机场时序网络模型网络效率损失曲线进行分析，识别其拐点，以此对机场网络抗毁性进行评价。
（6）通过输入全国预先航班计划，运用时序网络模型进行预先推演，在航路点和机场设置相应阀值，对按航班计划预先推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警，并在时序网络模型内对航班计划时刻进行优化。
进一步的，步骤（1）中航班数据包含民航航班信息***中的民航航班数据库文档以及航班计划文档，以及对数据进行周期性分析。
进一步的，步骤（2）中构建机场多层时序网络模型是基于时序网络模型构建，包括以下步骤：
（2.1）构建机场时序网络邻接矩阵，构建表示的机场时序网络，其中（nodes）为网络中所包含的机场节点集，/>，/>（edge）为机场间的航班连边集合，/>，/>（time）是航班运行时间集合，,其中/>表示从节点i到节点j航班起飞时刻，/>表示从节点i到节点j航班运行时间。
其中为从节点i到节点j所包含的全部航班运行时间集合。
其中表示在观察期T中从节点i到节点j的第m次执飞航班的起飞时间，/>表示在观察期T中从节点i到节点j的第m次执飞航班的飞行时间。
（2.2）将时间序列网络中航班数据的总时间按其固有周期划分为相应的层数，在每层中将对应的时序数据增加到网络中，构建机场多层时序网络模型。
进一步的，所述步骤（3）是对机场网络拓扑结构和抗毁性进行分析，如图3所示，具体包括以下步骤：
（3.1）对网络结构中的时序度进行分析，包括时序入度，时序出度，时序度，以及其幂律性的分析，并使用幂律分布函数为，并对结果进行拟合，确定参数为时序入度和时序出度条件下的常数，其值表征在模型数据取值的情况下，用幂律分布函数能很好的进行拟合。
其中为时序入度，T为数据固有周期长度，/>为时序多层网络的网络层编号，为第/>层网络的入度值，/>为时序出度，/>为第/>层网络的出度值，/>为时序网络总度值，/>表示第/>层网络的入度值和出度值之和。
（3.2）对网络结构中的平均时序距离进行分析，平均时序距离是表征网络上各节点以消耗时间作为指标的平均各节点间远近的参数。
其中为节点i到节点j的平均时序距离，T为网络固有周期长度，/>为第n次航班的起飞时刻，/>为第n次航班的飞行时间。
（3.3）对网络结构中的节点中心性研究包括时序度中心性，时序介数中心性，时序接近度中心性等指标进行。
时序度中心性（temporal degree centrality, TD）
一个节点的时序度越大就意味着该节点的时序度中心性越高，此节点在时序网络中占有更加重要的地位。网络规模越大，时序度的最大可能值就越高。为了消除网络规模变化对度中心性的影响，需要除以N-1进行标准化。
其中为第/>层网络的时序度中心值，N为网络节点数，/>表示第/>层网络的入度值和出度值之和。
时序介数中心性（temporal betweenness centrality, TB）
计算网络中任意两个节点的所有最短路径，如果这些最短路径中很多条都经过了某个节点，那么就认为这个节点的中介中心性高,则节点的时序介数中心性计算公式如下。
其中表示编号为k节点第/>层网络的介数中心性值，k取值为1到N，/>是在观察时间/>中最短路径中经过/>节点的数量，/>是在观察时间/>中所有节点对的最短路径总数。
时序接近度中心性（temporal closeness centrality，TC）:
节点到图中其他节点的最短时序距离都很小，那么它的接近中心性就很高。相比中介中心性，接近中心性更接近几何上的中心位置，则节点的时序接近度中心性计算公式如下。
其中为编号为k节点第/>层网络的接近度中心性值，/>是节点个数，/>表示/>到所有其他节点的平均时序距离之和。
进一步的，步骤（4）中构建机场时序网络模型抗毁性分析指标，指标包括时序网络效率和时序最大联通子图，并将结果保存为时序网络效率数据库。
时序网络效率：
其中是节点个数，/>表示节点i到节点j平均时序距离。
进一步的，步骤（5）中的对机场网络使用的蓄意攻击，随机攻击和非完全信息攻击的攻击方式，对机场网络进行模拟攻击，对得到的机场时序网络模型网络效率和最大联通子图损失曲线进行分析，识别其拐点，以此对机场网络抗毁性进行评估。时序网络抗毁性分析指标包括时序网络效率和时序网络最大联通子图。
（5.1）随机攻击策略为在机场节点集合中随机选择攻击顺序并形成攻击顺序集合中右上标r为攻击类型上标，这里r也表示随机攻击，1attack代表第1次攻击，2attack代表第2次攻击，vattack代表第v次攻击， n为受攻击机场节点编号，并按顺序集合依次删除节点并计算网络效率和最大联通子图。
（5.2）蓄意攻击策略分为按时序度强度排序攻击、按时序介数中心性排序攻击、按时序接近度中心性强度攻击，分别形成时序度攻击顺序集合，中右上标k为攻击类型上标，其中k表示时序度强度排序攻击，1attack代表第1次攻击，2attack代表第2次攻击，vattack代表第v次攻击， n为受攻击机场节点编号，时序介数攻击顺序集合/>中右上标p为攻击类型上标，其中p表示时序介数值排序攻击，右下标为攻击次数下标vattack代表第v次攻击， n为受攻击机场节点编号，时序接近度攻击顺序集合/>，中右上标b为攻击类型上标，其中b表示时序接近度值排序攻击，右下标为攻击次数下标vattack代表第v次攻击，n为受攻击机场节点编号，并按顺序集合依次删除节点并计算网络效率和最大联通子图。
（5.3）非完全信息攻击的攻击方式为在给定已知信息的情况下对已知信息的节点集合运用蓄意攻击策略进行攻击，对未知信息的节点运用随机攻击的策略，并按先已知后未知节点的顺序形成非完全信息攻击顺序集合中右上标m为攻击类型上标，其中m表示非完全信息攻击，1attack代表第1次攻击，2attack代表第2次攻击，vattack代表第v次攻击， n为受攻击机场节点编号,最后按顺序集合依次删除节点并计算网络效率和最大联通子图，通过上述（5.1）-（5.3）步骤得到机场时序网络模型网络效率损失曲线。
（5.4）对得到的机场时序网络模型网络效率和最大联通子图损失曲线进行分析，识别其拐点，完成对机场网络抗毁性的评估。
进一步的，如图4 所示，步骤（6）通过输入预先航班计划，运用时序网络模型进行预先推演，在航路点和机场设置相应阀值，对按航班计划预先推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警，并在时序网络模型内对航班计划时刻进行优化：
（6.1）将航班计划数据进行导入，进行数据预处理，即对数据进行筛选，筛除空数据并筛选航班编号、航班起飞时间、飞行时间、起飞机场和目的地机场关键数据，基于处理过后的数据将带有时间属性的数据添加入邻接矩阵，矩阵中每个元素不再用简单的0，1表示是否存在连边，而是用列表代替，并将节点对间的航班时间数据添加入列表中进行储存形成航班计划时序邻接矩阵，导入航路信息和机场信息，并设定航路和机场阀值。
（6.2）运用时序网络模型进行预先推演，对按航班计划推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警。
（6.3）在时序网络模型内对航班计划的航班时刻进行优化，达到对航路和机场的拥堵情况发生预发现、预优化，从而提高航路和机场资源的利用率。
本发明构思合理，过程清晰，能提高对动态***结构分析的准确性，对提高机场网络的抗毁性提供依据，并为机场网络的规划建设提供数据支持，因为时间属性的加入也使得在网络模型内直接进行各种视角下的航班时刻优化成为可能。弥补目前运用时序网络模型对机场网络进行分析的空白，相比与静态网络分析其能克服现在常用的静态网络对动态***的不适应性，提高分析的准确性。
基于与上述方法实施例相同地发明构思，本申请实施例还提供了一种基于时序网络模型的航班计划优化***，用于实现上面实施例所述的基于时序网络模型的航班计划优化方法，如图5所示，所述***具体包括：
航班数据采集模块110：用于采集航班数据，并运用功率谱分析方法分析航班数据的周期性，确定航班数据固有周期；
多层时序网络模型构建模块120：用于基于时序网络模型构建机场多层时序网络模型；
网络模型拓扑结构分析模块130：用于构建机场时序网络模型拓扑结构分析指标，包括时序度、平均时序距离和节点中心性指标，并利用拓扑结构分析指标对机场时序网络模型结构进行分析；
网络模型抗毁性分析模块140：用于构建机场时序网络模型抗毁性分析指标，包括时序网络效率和时序最大联通子图，并利用抗毁性分析指标对机场时序网络模型进行抗毁性分析；
抗毁性评估模块150：用于对机场网络进行抗毁性评估，采用蓄意攻击、随机攻击和非完全信息攻击的策略对机场网络进行模拟攻击，得到机场时序网络模型网络效率损失曲线，对得到的网络效率损失曲线进行分析，识别其拐点，完成对机场网络抗毁性的评估；
航班计划优化模块160：用于输入预先航班计划，运用时序网络模型进行预先推演，在航路点和机场设置相应阀值，对按航班计划预先推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警，并在时序网络模型内对航班计划时刻进行优化，输出优化后航班计划。
基于与上述方法实施例相同地发明构思，本申请实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现所述的航班计划优化方法。
所述航班计划优化方法包括如下步骤：
采集航班数据，并运用功率谱分析方法分析航班数据的周期性，确定航班数据固有周期；
基于时序网络模型构建机场多层时序网络模型；
构建机场时序网络模型拓扑结构分析指标，包括时序度、平均时序距离和节点中心性指标，并利用拓扑结构分析指标对机场时序网络模型结构进行分析；
构建机场时序网络模型抗毁性分析指标，包括时序网络效率和时序最大联通子图，并利用抗毁性分析指标对机场时序网络模型进行抗毁性分析；
机场网络抗毁性评估：采用蓄意攻击、随机攻击和非完全信息攻击的策略对机场网络进行模拟攻击，得到机场时序网络模型网络效率损失曲线，对得到的网络效率损失曲线进行分析，识别其拐点，完成对机场网络抗毁性的评估；
输入预先航班计划，运用时序网络模型进行预先推演，在航路点和机场设置相应阀值，对按航班计划预先推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警，并在时序网络模型内对航班计划时刻进行优化，输出优化后航班计划。
通过本发明弥补了目前运用时序网络模型对机场网络进行分析的空白，解决现在航班计划编制时无法做到对航路拥堵和机场超负荷运行的预先识别和预先优化问题，做到预发现、预优化，有较强的创新性和实际意义。
尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换,凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于时序网络模型的航班计划优化方法，其特征在于，包括如下步骤：
S101、采集航班数据，并运用功率谱分析方法分析航班数据的周期性，确定航班数据固有周期；
S102、基于时序网络模型构建机场多层时序网络模型，具体包括如下：
S201、构建机场时序网络邻接矩阵，构建G＝{N，E，T}表示机场时序网络，其中N为网络中所包含的机场节点集，N(G)＝1,2,3,4,5...N，E为机场间的航班连边集合，E(G)＝(i1,j1),(i2,j2),(i3,j3),...,(iE,jE)，其中i表示起始节点即起飞机场，j表示目标节点即到达机场，T是航班运行时间集合，T(G)＝{(t12,σ12),(t13,σ13)…(tij,σij)},其中tij表示从节点i到节点j航班起飞时刻，σij表示从节点i到节点j航班运行时间；机场时序网络邻接矩阵如下：
其中T[i,j]为从节点i到节点j所包含的全部航班运行时间集合；
其中表示在观察期T中从节点i到节点j的第m次执飞航班的起飞时间，/>表示在观察期T中从节点i到节点j的第m次执飞航班的飞行时间；
S202、将时序网络中航班数据的总时间按其固有周期划分为相应的层数，在每层中将对应的时序数据增加到网络中，构建机场多层时序网络模型；
S103、构建机场时序网络模型拓扑结构分析指标，包括时序度、平均时序距离和节点中心性指标，并利用拓扑结构分析指标对机场时序网络模型结构进行分析；
在步骤S103中，利用时序度分析指标对网络结构中的时序度进行分析，包括时序入度，时序出度，时序度，以及幂律性的分析，使用幂律分布函数为和其中Ain,Bin,λin为时序入度幂律分布函数待求的三个参数，kin为节点时序入度值，P(kin)为时序入度值为k的机场节点的概率值，Aout,Bout,λout为时序出度幂律分布函数待求的三个参数，kout为节点时序出度值，P(kout)为时序出度值为k的机场节点的概率值，并对数据进行拟合求出待求参数；
对网络结构中的平均时序距离进行分析，平均时序距离是表征网络上各节点以消耗时间作为指标的平均各节点间远近的参数；
其中为节点i到节点j的平均时序距离，T为网络固有周期长度，tn为第n次航班的起飞时刻，Δtn为第n次航班的飞行时间；
利用节点中心性指标对网络结构进行分析，所述节点中心性指标包括时序度中心性，时序介数中心性以及时序接近度中心性；
时序度中心性：
一个节点的时序度越大就意味着该节点的时序度中心性越高，此节点在时序网络中占有更加重要的地位；网络规模越大，时序度的最大可能值就越高；为了消除网络规模变化对度中心性的影响，需要除以N-1进行标准化，表达如下：
其中为第Ti层网络的时序度中心值，N为网络节点数，/>表示第Ti层网络的入度值和出度值之和；
时序介数中心性：
计算网络中任意两个节点的所有最短路径，如果这些最短路径中很多条都经过了某个节点，那么就认为这个节点的中介中心性高,则节点i的时序介数中心性计算公式如下：
其中表示编号为k节点第Ti层网络的介数中心性值，k取值为1到N，/>是在观察时间T中最短路径中经过k节点的数量，di，j是在观察时间T中所有节点对的最短路径总数；
时序接近度中心性:
节点到图中其他节点的最短时序距离都很小，那么它的接近中心性就很高，相比中介中心性，接近中心性更接近几何上的中心位置，则节点i的时序接近度中心性计算公式如下：
其中为编号为k节点第Ti层网络的接近度中心性值，N是节点个数，/>表示i到所有其他节点的平均时序距离之和；
S104、构建机场时序网络模型抗毁性分析指标，包括时序网络效率和时序最大联通子图，并利用抗毁性分析指标对机场时序网络模型进行抗毁性分析；
S105、机场网络抗毁性评估：采用蓄意攻击、随机攻击和非完全信息攻击的策略对机场网络进行模拟攻击，得到机场时序网络模型网络效率损失曲线，对得到的网络效率损失曲线进行分析，识别其拐点，完成对机场网络抗毁性的评估；
S106、输入预先航班计划，运用时序网络模型进行预先推演，在航路点和机场设置相应阀值，对按航班计划预先推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警，并在时序网络模型内对航班计划时刻进行优化，输出优化后航班计划；
在步骤S106中，具体包括如下步骤：
S601、将航班计划数据进行导入，并对航班数据进行预处理，即对数据进行筛选，筛除空数据并筛选航班编号、航班起飞时间、飞行时间、起飞机场和目的地机场关键数据，基于处理过后的数据将带有时间属性的数据添加入邻接矩阵，矩阵中每个元素不再用简单的0，1表示是否存在连边，而是用列表代替，并将节点对间的航班时间数据添加入列表中进行储存形成航班计划时序邻接矩阵，导入航路信息和机场信息，并设定航路点和机场阀值；
S602、运用时序网络模型进行预先推演，对按航班计划推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警；
S603、在时序网络模型内对航班计划的航班时刻进行优化，输出优化后航班计划。
2.根据权利要求1所述的基于时序网络模型的航班计划优化方法，其特征在于：在步骤S101中，所述的航班数据包含民航航班信息***中的民航航班数据库文档数据以及航班计划文档数据。
3.根据权利要求1所述的基于时序网络模型的航班计划优化方法，其特征在于：在步骤S105中，具体包括如下步骤：
S501、随机攻击策略为在机场节点集合中随机选择攻击顺序并形成攻击顺序集合其中，n为受攻击机场节点编号，v为攻击次数，1attack代表第1次攻击，vattack代表第v次攻击，并按集合顺序依次删除网络中对应节点并计算网络效率和最大联通子图；
S502、蓄意攻击策略分为按时序度强度排序攻击、按时序介数中心性排序攻击、按时序接近度中心性强度攻击，分别形成时序度攻击顺序集合k表示时序度强度排序攻击，时序介数攻击顺序集合/>p表示时序介数值排序攻击，时序接近度攻击顺序集合/>b表示时序接近度值排序攻击，并按顺序集合依次删除节点并计算网络效率和最大联通子图；
S503、非完全信息攻击的攻击方式为在给定已知信息的情况下对已知信息的节点集合运用蓄意攻击策略进行攻击，对未知信息的节点运用随机攻击的策略，并按先已知后未知节点的顺序形成非完全信息攻击顺序集合m表示非完全信息攻击，最后按顺序集合依次删除节点并计算网络效率和最大联通子图，通过步骤S501-S503，得到机场时序网络模型网络效率损失曲线；
S504、对得到的网络效率损失曲线进行分析，识别其拐点，完成对机场网络抗毁性的评估。
4.一种根据权利要求1-3中任一项所述基于时序网络模型的航班计划优化方法的***，其特征在于：所述***包括：
航班数据采集模块：用于采集航班数据，并运用功率谱分析方法分析航班数据的周期性，确定航班数据固有周期；
多层时序网络模型构建模块：用于基于时序网络模型构建机场多层时序网络模型；
网络模型拓扑结构分析模块：用于构建机场时序网络模型拓扑结构分析指标，包括时序度、平均时序距离和节点中心性指标，并利用拓扑结构分析指标对机场时序网络模型结构进行分析；
网络模型抗毁性分析模块：用于构建机场时序网络模型抗毁性分析指标，包括时序网络效率和时序最大联通子图，并利用抗毁性分析指标对机场时序网络模型进行抗毁性分析；
抗毁性评估模块：用于对机场网络进行抗毁性评估，采用蓄意攻击、随机攻击和非完全信息攻击的策略对机场网络进行模拟攻击，得到机场时序网络模型网络效率损失曲线，对得到的网络效率损失曲线进行分析，识别其拐点，完成对机场网络抗毁性的评估；
航班计划优化模块：用于输入预先航班计划，运用时序网络模型进行预先推演，在航路点和机场设置相应阀值，对按航班计划预先推演中出现的航路拥堵点和超过枢纽机场高峰服务能力的情况进行识别和预警，并在时序网络模型内对航班计划时刻进行优化，输出优化后航班计划。
5.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理执行时实现如权利要求1-3中任意一项所述的航班计划优化方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311127516.8A CN116882590B (zh) | 2023-09-04 | 2023-09-04 | 基于时序网络模型的航班计划优化方法、***及储存介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311127516.8A CN116882590B (zh) | 2023-09-04 | 2023-09-04 | 基于时序网络模型的航班计划优化方法、***及储存介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116882590A CN116882590A (zh) | 2023-10-13 |
| CN116882590B true CN116882590B (zh) | 2023-11-24 |
Family
ID=88264730
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311127516.8A Active CN116882590B (zh) | 2023-09-04 | 2023-09-04 | 基于时序网络模型的航班计划优化方法、***及储存介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116882590B (zh) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104751683A (zh) * | 2015-03-27 | 2015-07-01 | 中国民用航空飞行学院 | 机场航班滑行调度指挥***及方法 |
| CN111191843A (zh) * | 2019-12-30 | 2020-05-22 | 南京航空航天大学 | 一种基于时序网络传播动力学方程的机场延误预测方法 |
| CN111738425A (zh) * | 2020-08-07 | 2020-10-02 | 北京航空航天大学 | 一种基于图挖掘的机场延误原因分析方法 |
| CN112489498A (zh) * | 2020-11-25 | 2021-03-12 | 中国电子科技集团公司第二十八研究所 | 一种航路交通精细改航规划方法 |
| CN113759894A (zh) * | 2020-05-22 | 2021-12-07 | 株式会社东芝 | 信息处理装置、信息处理方法、信息处理***以及计算机程序 |
| CN114464017A (zh) * | 2022-01-27 | 2022-05-10 | 南京航空航天大学 | 一种基于排队论的机场群航班延误估算方法 |
| CN115311902A (zh) * | 2022-03-25 | 2022-11-08 | 中国航空无线电电子研究所 | 一种基于多层时序网络的实时航路规划及优化方法 |
| CN115564188A (zh) * | 2022-09-16 | 2023-01-03 | 中国民航大学 | 一种航班计划鲁棒性评价***及方法 |
| CN115713872A (zh) * | 2022-11-11 | 2023-02-24 | 中国航空无线电电子研究所 | 一种基于soa的环境自适应航路规划方法 |
| CN116246496A (zh) * | 2023-05-10 | 2023-06-09 | 北京航空航天大学 | 一种基于航班与流量协同环境的飞行计划规划方法 |
| CN116523138A (zh) * | 2023-05-08 | 2023-08-01 | 南京航空航天大学 | 一种考虑旅客中转出行选择的枢纽机场航班时刻优化方法 |
| CN116659504A (zh) * | 2023-05-31 | 2023-08-29 | 中国民用航空飞行学院 | 一种基于卡尔曼滤波法的实时飞行数据估计方法 |
-
2023
- 2023-09-04 CN CN202311127516.8A patent/CN116882590B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104751683A (zh) * | 2015-03-27 | 2015-07-01 | 中国民用航空飞行学院 | 机场航班滑行调度指挥***及方法 |
| CN111191843A (zh) * | 2019-12-30 | 2020-05-22 | 南京航空航天大学 | 一种基于时序网络传播动力学方程的机场延误预测方法 |
| CN113759894A (zh) * | 2020-05-22 | 2021-12-07 | 株式会社东芝 | 信息处理装置、信息处理方法、信息处理***以及计算机程序 |
| CN111738425A (zh) * | 2020-08-07 | 2020-10-02 | 北京航空航天大学 | 一种基于图挖掘的机场延误原因分析方法 |
| CN112489498A (zh) * | 2020-11-25 | 2021-03-12 | 中国电子科技集团公司第二十八研究所 | 一种航路交通精细改航规划方法 |
| CN114464017A (zh) * | 2022-01-27 | 2022-05-10 | 南京航空航天大学 | 一种基于排队论的机场群航班延误估算方法 |
| CN115311902A (zh) * | 2022-03-25 | 2022-11-08 | 中国航空无线电电子研究所 | 一种基于多层时序网络的实时航路规划及优化方法 |
| CN115564188A (zh) * | 2022-09-16 | 2023-01-03 | 中国民航大学 | 一种航班计划鲁棒性评价***及方法 |
| CN115713872A (zh) * | 2022-11-11 | 2023-02-24 | 中国航空无线电电子研究所 | 一种基于soa的环境自适应航路规划方法 |
| CN116523138A (zh) * | 2023-05-08 | 2023-08-01 | 南京航空航天大学 | 一种考虑旅客中转出行选择的枢纽机场航班时刻优化方法 |
| CN116246496A (zh) * | 2023-05-10 | 2023-06-09 | 北京航空航天大学 | 一种基于航班与流量协同环境的飞行计划规划方法 |
| CN116659504A (zh) * | 2023-05-31 | 2023-08-29 | 中国民用航空飞行学院 | 一种基于卡尔曼滤波法的实时飞行数据估计方法 |
Non-Patent Citations (7)
| Title |
|---|
| Robustness Analysis for China’s Airport Network Based on Multi-Layer Temporal Complex Network Model;Jiuxia Guo 等;《International Conference on Transportation and Development 2023》;第14-24页 * |
| 中国航空网络时序特征分析;牟建红;黄格;吕欣;;电子科技大学学报(第03期);第144-150页 * |
| 基于冗余度的复杂网络抗毁性及节点重要度评估模型;王梓行;姜大立;漆磊;陈星;赵禹博;;复杂***与复杂性科学(第03期);第81-88页 * |
| 基于大数据的作战体系协同时序网络分析;伍文峰;郭圣明;贺筱媛;胡晓峰;;指挥与控制学报(第02期);第150-159页 * |
| 基于网络模体特征攻击的网络抗毁性研究;贾承丰;韩华;完颜娟;吕亚楠;;复杂***与复杂性科学(第04期);第46-53页 * |
| 航空公司机队的鲁棒性规划启发式算法;汪瑜;孙宏;;***工程理论与实践(第04期);第963-970页 * |
| 航空公司机队集中调度理论研究;孙宏;张翔;徐杰;;中国管理科学(第01期);第86-89页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116882590A (zh) | 2023-10-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ren et al. | Characterizing air traffic networks via large-scale aircraft tracking data: A comparison between China and the US networks | |
| CN109784636A (zh) | 欺诈用户识别方法、装置、计算机设备及存储介质 | |
| CN104102745B (zh) | 基于局部最小边的复杂网络社团挖掘方法 | |
| CN108683448B (zh) | 适用于航空网络的影响力节点识别方法及*** | |
| US9310204B2 (en) | Route modeler for improving desired environmental and economic flight characteristics | |
| Akbari et al. | An online optimization approach to post-disaster road restoration | |
| Bianco et al. | Locating sensors to observe network arc flows: exact and heuristic approaches | |
| CN116468186B (zh) | 一种航班链延误时间预测方法、电子设备及存储介质 | |
| Zhang et al. | Research on a model of node and path selection for traffic network congestion evacuation based on complex network theory | |
| CN102571431A (zh) | 针对复杂网络的基于群思想改进的Fast-Newman聚类方法 | |
| Ai et al. | A deep learning approach to predict the spatial and temporal distribution of flight delay in network | |
| CN112735188B (zh) | 基于复杂网络理论的空中交通网络脆弱性分析*** | |
| CN112270445B (zh) | 基于统计分析和分类预测的航班延误波及综合评估方法 | |
| de Souza et al. | MOCHA: A tool for mobility characterization | |
| CN112530206A (zh) | 空中交通网络脆弱性分析方法 | |
| Ren et al. | A measure of identifying influential waypoints in air route networks | |
| CN109887280B (zh) | 一种交通网络节点关键度评估方法 | |
| CN114781704A (zh) | 一种基于过站航班保障流程的航班延误预测方法 | |
| CN110991913A (zh) | 繁忙机场高峰时间拥堵风险分析方法 | |
| Yang et al. | Air traffic network optimization via Laplacian energy maximization | |
| Nasiboglu | Dijkstra solution algorithm considering fuzzy accessibility degree for patch optimization problem | |
| CN116882590B (zh) | 基于时序网络模型的航班计划优化方法、***及储存介质 | |
| CN117391653A (zh) | 一种机场航班延误传播预测方法及*** | |
| KR102081377B1 (ko) | 중소교량 네트워크의 유지관리 최적성능목표 결정 방법 및 시스템 | |
| Chunyu et al. | Discrete bat algorithm and application in community detection |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |