CN115472043B - 一种基于p-中值定理的机场群航路飞行时间场景生成方法 - Google Patents
一种基于p-中值定理的机场群航路飞行时间场景生成方法 Download PDFInfo
- Publication number
- CN115472043B CN115472043B CN202211059423.1A CN202211059423A CN115472043B CN 115472043 B CN115472043 B CN 115472043B CN 202211059423 A CN202211059423 A CN 202211059423A CN 115472043 B CN115472043 B CN 115472043B
- Authority
- CN
- China
- Prior art keywords
- flight
- time
- scene
- flight time
- scenes
- 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
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/91—Radar or analogous systems specially adapted for specific applications for traffic control
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/0043—Traffic management of multiple aircrafts from the ground
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Aviation & Aerospace Engineering (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于p‑中值定理的机场群航路飞行时间场景生成方法。选定机场和航路点后,我们计算通过历史数据计算每一次航路飞行时间,经过剔除异常值,进行飞行时间离散化,得到航路飞行时间集合。飞行时间集合以分钟为单位,若集合中的数据种类数大于等于拟生成场景数,则运行场景生成模型,否则补充若干出现概率为0的场景直到达到所需场景数。其中场景生成模型是一个线性规划模型,以最小化场景生成前后概率分布的康托洛维奇距离(Kantorovich distance)为目标,我们经过转化得到一个p‑中值模型,通过求解模型可以得到各个场景下的飞行时间和相应的概率。本发明旨在对机场群中任意航路生成指定场景个数下的场景,并得到其相应的飞行时间和概率。
Description
技术领域
本发明涉及航路场景生成方法,特别涉及一种基于p-中值定理的航路飞行时间场景生成方法。
背景技术
随着区域化经济的不断加强,机场之间的联合作用也越来越明显。一些地理位置相近的机场之间形成了机场群。十几年来,机场群中不同机场之间的制约和依赖变得更加突出,这些冲突主要体现在不同机场共享的进离场点和进离场航路之间。在进行机场群的航班时刻配置时,限制航路点容流平衡是机场之间协同优化的重点。然而,航路点的流量并不能直接得到,需要根据机场的航班时刻表和机场到航路点的飞行时间推算得到。而对于机场群内的机场和航路点,它们之间的飞行时间存在不确定性,主要基于以下原因:1)天气原因,2)军方活动,以及3)流量控制等。在不同因素的综合作用下,航路飞行被分为不同场景,每种场景对应的飞行时间不同。
不同的航路飞行时间场景,对应不同的航路飞行时间,进而导致了不同的航路点流量,而战术层面上实际发生的场景并不能在战术层面上被确定。因此,使用正确的方法生成合适的场景以及得到相应的发生概率,对战略航班时刻的配置、提高航班时刻表的稳定性具有重要意义。
发明内容
发明目的:本发明目的是提供一种基于p-中值定理的机场群航路飞行时间场景生成方法。
技术方案:本发明所述的基于p-中值定理的机场群航路飞行时间场景生成方法,包括如下步骤:
步骤1,对于给定航路点和机场,收集历史航班运行数据和雷达轨迹数据,得出机场到航路点每一次的航路飞行时间;
步骤2,对已收集到的数据集进行处理,剔除异常值,并进行离散化,得到飞行时间集合;
步骤3,判断飞行时间集合中数据种类的个数,选择场景生成模型或场景补充得到给定场景个数下的生成的场景;
步骤4,根据模型求解的结果得到各个场景对应的航路飞行时间和发生概率。
优选的,所述步骤1的具体过程为:
步骤11:通过历史运行数据,得到每个航班的进/离场时间;
步骤12:通过历史雷达轨迹数据,得到每个航班到达航路点的时间;
步骤13:筛选同时经过给定机场和航路点的航班,分别对他们的进/离场时间和过点时间作差,得到每天所有航班的飞行时间数据。
优选的,所述步骤2的具体过程为:
步骤21:分析飞行时间的频率分布,设定阈值为Q3+2·(Q3-Q1),其中Q3和Q1分别为数据的第3个4分位数和第1个4分位数;
步骤22:剔除原始数据中阈值以外数据;
步骤23:以1分钟为单位,对飞行时间进行离散化,得到飞行时间集合。
优选的,所述步骤3的具体过程为:
步骤31:判断飞行时间集合中数据的种类个数S′,即集合中飞行时间有几种不同的取值,若个数S′≥S,则进行场景生成模型,否则进程场景补充;
步骤32:建立场景生成模型,场景生成模型是一个混合整数规划模型,使用康托洛维奇距离来衡量场景生成前后飞行时间概率分布之间的距离,求解2个概率分布之间的康托洛维奇距离可以表述为如下的线性规划模型:
D表示2个概率分布之间的康托洛维奇距离,表示原场景集中飞行时间为i分钟出现的概率,pj表示生成的场景集中飞行时间为j分钟出现的概率;模型中的决策变量为它是一个连续变量,表示i到j之间转移的概率的数量;cij表示单位数量概率转移的代价;
步骤33:将历史数据中每一次航路飞行时间都视为场景集中的一个数据点,我们把原始飞行时间集合记为 表示原飞行时间集合中元素的数量;/>中包含了每一次航路飞行的飞行时间{τ′1,τ′2,...,τ′S′},分别对应其发生概率{p′1,p′2,...,p′S′},索引集合记为N,N={1,2,...,S′},每一次发生的概率相等,则/>生成的飞行时间集合记为/>包含的元素个数为S,/>其发生概率分别为{p1,p2,...,pS};将原飞行时间集合中每次的飞行时间都视为一个取值,而不是把相同的飞行时间进行加和来统计发生概率,那么原集合中飞行时间概率分布可以看作平均分布,每次飞行时间发生的概率都相等,而且不可再分,此时最小化康托洛维奇距离可以转化为一个求解p-中值模型的问题。p-中值模型是指在一个给定数量和位置的需求集合和一个候选设施位置的集合下,分别为p个设施找到合适的位置,并指派每个需求点到一个特定的设施,使之达到在工厂和需求点之间的运输费用最低;在本问题中,生成飞行时间集合中τ1,τ2,...,τS可以视为设施,原飞行时间集合中τ′1,τ′2,...,τ′S′可以视为需求点;
步骤34:将p-中值模型用公式表示,设置决策变量为λs′s和γs′,γs,γi,γj为0-1变量,s,s′,i,j∈N;λs′s=1表示原场景集中τ′s′被映射到生成场景集中的τs;γs′=1、γs=1、γi=1、γj=1表示原场景集中的τ′s′,τ′s,τ′i,τ′j被选为生成场景集中的场景;则目标函数为:
ds′s表示2个场景之间的距离,即飞行时间的差值,其计算方式为:
ds′s=|τ′s′-τs|,s,s′∈N
限制生成场景集中场景数为S:
限制生成场景集中的每个场景不同:
|τi-τj|≥γi·γj,i,j∈N
原场景集中的场景只能被映射到新场景集中的一个场景:
原场景集中场景的映射目标必须被选中成为生成场景集中的元素:
在以上约束和目标函数下求解模型;
步骤35:若飞行时间集合中数据的种类个数不足S,则补充若干出现概率为0的场景直到达到所需场景数。
优选的,所述步骤4具体过程为:
步骤41:若S′≥S,则生成的场景集中飞行时间分别为其对应的概率计算公式为:
步骤42:若S′<S,则生成的场景集中飞行时间分别为原场景集中的飞行时间,对应概率为在原场景集中出现的频率;和(S-S′)个M,对应的概率为0。
有益效果:机场群内机场到航路点的飞行时间存在不确定性,在战略层面的航路点流量计算中,需要考虑到这一不确定性,从而提高航班时刻表的鲁棒性,减少运行时的航班延误,解决不确定性最常见的方法就是随机规划,而随机规划的求解又依赖于场景生成。本发明解决了航路场景生成过程,不但将飞行时间进行离散化,还给出了给定场景数下的各个场景对应的航路飞行时间以及相应的发生概率,这对后续不确定性建模和随机规划的求解提供了基础;另外,场景生成过程可以对场景数进行缩减,得到发生概率较大的几个场景,从而忽略对极少数异常情况的考量,这种取舍可以大大简化模型的求解过程,保证在现有的计算能力下得到最优的解决方案。
附图说明
图1是本发明方法的流程图;
图2是实施例中机场到航路点的飞行时间分布箱线图。
具体实施方式
本实施例将本发明方法应用于粤港澳大湾区机场群,包括5个机场:广州白云国际机场(ZGGG)、深圳宝安国际机场(ZGSZ)、珠海金湾国际机场(ZGSD)、澳门国际机场(VMMC)和惠州平潭机场(ZGHZ)5个机场。我们选择机场群内5个机场共用的进场点GYA作为研究对象,分别对5个机场到GYA的航路飞行时间进行场景生成。本实施例中针对ZGGG、ZGSD、ZGSZ、ZGHZ、VMMC 5个机场以及GYA航路点进行航路飞行时间场景生成。如图1所示,一种基基于p-中值定理的机场群航路飞行时间场景生成方法,包括如下步骤:
步骤1,对于ZGGG、ZGSD、ZGSZ、ZGHZ、VMMC 5个机场和GYA航路点,收集历史航班运行数据和雷达轨迹数据,由2019年12月的历史航班运行和雷达轨迹数据得出机场到航路点每一次的航路飞行时间;
步骤11:通过历史运行数据,得到每个航班的进/离场时间;
步骤12:通过历史雷达轨迹数据,得到每个航班到达航路点的时间;
步骤13:筛选同时经过给定机场和航路点的航班,分别对他们的进/离场时间和过点时间作差,得到每天所有航班的飞行时间数据。
步骤2,对已收集到的数据集进行处理,剔除异常值,并进行离散化,得到飞行时间集合;
步骤21:分析飞行时间的频率分布,设定阈值为Q3+2·(Q3-Q1),其中Q3和Q1分别为数据的第3个4分位数和第1个4分位数;
步骤22:剔除原始数据中阈值以外数据;
步骤23:以1分钟为单位,对飞行时间进行离散化,得到飞行时间集合,每个机场到GYA的飞行时间分布箱线图如图2所示。
步骤3,判断飞行时间集合中数据种类的个数,选择场景生成模型或场景补充得到给定场景个数下的生成的场景;
步骤31:设定S=2,判断飞行时间集合中数据的种类个数S′,即集合中飞行时间有几种不同的取值,若个数S′≥S,则进行场景生成模型,否则进程场景补充,经过判断5个机场到GYA的飞行时间集合中数据的种类数均满足:S′≥S,可以进行场景生成;
步骤32:建立场景生成模型,场景生成模型是一个混合整数规划模型,我们使用康托洛维奇距离来衡量场景生成前后飞行时间概率分布之间的距离;
步骤33:将历史数据中每一次航路飞行时间都视为场景集中的一个数据点,把最小化康托洛维奇距离转化成一个p-中值模型;
步骤34:将p-中值模型用公式表示,把ZGGG、ZGSD、ZGSZ、ZGHZ、VMMC 5个机场到GYA航路点飞行时间集合带入模型中进行求解;
步骤4,根据模型求解的结果得到各个场景对应的航路飞行时间和发生概率。
经过判断5个机场到GYA的飞行时间集合中数据的种类数均满足:S′≥S,则生成的场景集中飞行时间分别为其对应的概率计算公式为:
生成的场景对应的飞行时间和出现的概率如表1所示
表1
机场 | 航路点 | 飞行时间1/分钟 | 概率1 | 飞行时间2/分钟 | 概率2 |
ZGGG | GYA | 16 | 0.779 | 24 | 0.221 |
ZGSD | GYA | 21 | 0.706 | 28 | 0.224 |
ZGSZ | GYA | 22 | 0.549 | 25 | 0.451 |
ZGHZ | GYA | 24 | 0.604 | 26 | 0.396 |
VMMC | GYA | 24 | 0.460 | 27 | 0.540 |
Claims (3)
1.一种基于p-中值定理的机场群航路飞行时间场景生成方法,其特征在于:包括如下步骤:
步骤1,对于给定航路点和机场,收集历史航班运行数据和雷达轨迹数据,得出机场到航路点每一次的航路飞行时间;
步骤2,对已收集到的数据集进行处理,剔除异常值,并进行离散化,得到飞行时间集合;
步骤3,判断飞行时间集合中数据种类的个数,选择场景生成模型或场景补充得到给定场景个数下的生成的场景;所述步骤3的具体过程为:
步骤31:判断飞行时间集合中数据的种类个数S′,即集合中飞行时间有几种不同的取值,若个数S′≥S,则进行场景生成模型,否则进程场景补充;
步骤32:建立场景生成模型,场景生成模型是一个混合整数规划模型,使用康托洛维奇距离来衡量场景生成前后飞行时间概率分布之间的距离,求解2个概率分布之间的康托洛维奇距离可以表述为如下的线性规划模型:
D表示2个概率分布之间的康托洛维奇距离,p′i表示原场景集中飞行时间为i分钟出现的概率,pj表示生成的场景集中飞行时间为j分钟出现的概率;模型中的决策变量为它是一个连续变量,表示i到j之间转移的概率的数量;cij表示单位数量概率转移的代价;
步骤33:将历史数据中每一次航路飞行时间都视为场景集中的一个数据点,把原始飞行时间集合记为 表示原飞行时间集合中元素的数量;/>中包含了每一次航路飞行的飞行时间{τ′1,τ′2,...,τ′S′},分别对应其发生概率{p′1,p′2,...,p′S′},索引集合记为N,N={1,2,...,S′},每一次发生的概率相等,则/>生成的飞行时间集合记为/>包含的元素个数为S,/>其发生概率分别为{p1,p2,...,pS};将原飞行时间集合中每次的飞行时间都视为一个取值,而不是把相同的飞行时间进行加和来统计发生概率,那么原集合中飞行时间概率分布可以看作平均分布,每次飞行时间发生的概率都相等,而且不可再分,此时最小化康托洛维奇距离可以转化为一个求解p-中值模型的问题,p-中值模型是指在一个给定数量和位置的需求集合和一个候选设施位置的集合下,分别为p个设施找到合适的位置,并指派每个需求点到一个特定的设施,使之达到在工厂和需求点之间的运输费用最低;在本问题中,生成飞行时间集合中τ1,τ2,...,τS可以视为设施,原飞行时间集合中τ′1,τ′2,...,τ′S′可以视为需求点;
步骤34:将p-中值模型用公式表示,设置决策变量为λs′s和γs′,γs,γi,γj为0-1变量,s,s′,i,j∈N;λs′s=1表示原场景集中τ′s′被映射到生成场景集中的τs;γs′=1、γs=1、γi=1、γj=1表示原场景集中的τ′s′,τ′s,τ′i,τ′j被选为生成场景集中的场景;则目标函数为:
ds′s表示2个场景之间的距离,即飞行时间的差值,其计算方式为:
ds′s=|τ′s′-τs|,s,s′∈N
限制生成场景集中场景数为S:
限制生成场景集中的每个场景不同:
|τi-τj|≥γi·γj,i,j∈N
原场景集中的场景只能被映射到新场景集中的一个场景:
原场景集中场景的映射目标必须被选中成为生成场景集中的元素:
在以上约束和目标函数下求解模型;
步骤35:若飞行时间集合中数据的种类个数不足S,则补充若干出现概率为0的场景直到达到所需场景数;
步骤4,根据模型求解的结果得到各个场景对应的航路飞行时间和发生概率;所述步骤4的具体过程为:
步骤41:若S′≥S,则生成的场景集中飞行时间分别为其对应的概率计算公式为:
步骤42:若S′<S,则生成的场景集中飞行时间分别为原场景集中的飞行时间,对应概率为在原场景集中出现的频率;和(S-S′)个M,对应的概率为0。
2.根据权利要求1所述基于p-中值定理的机场群航路飞行时间场景生成方法,其特征在于:所述步骤1的过程为:
步骤11:通过历史运行数据,得到每个航班的进/离场时间;
步骤12:通过历史雷达轨迹数据,得到每个航班到达航路点的时间;
步骤13:筛选同时经过给定机场和航路点的航班,分别对他们的进/离场时间和过点时间作差,得到每天所有航班的飞行时间数据。
3.根据权利要求2所述基于p-中值定理的机场群航路飞行时间场景生成方法,其特征在于:所述步骤2的具体过程为:
步骤21:分析飞行时间的频率分布,设定阈值为Q3+2·(Q3-Q1),其中Q3和Q1分别为数据的第3个4分位数和第1个4分位数;
步骤22:剔除原始数据中阈值以外数据;
步骤23:以1分钟为单位,对飞行时间进行离散化,得到飞行时间集合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211059423.1A CN115472043B (zh) | 2022-08-30 | 2022-08-30 | 一种基于p-中值定理的机场群航路飞行时间场景生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211059423.1A CN115472043B (zh) | 2022-08-30 | 2022-08-30 | 一种基于p-中值定理的机场群航路飞行时间场景生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115472043A CN115472043A (zh) | 2022-12-13 |
CN115472043B true CN115472043B (zh) | 2023-09-29 |
Family
ID=84368976
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211059423.1A Active CN115472043B (zh) | 2022-08-30 | 2022-08-30 | 一种基于p-中值定理的机场群航路飞行时间场景生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115472043B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104867184A (zh) * | 2015-06-12 | 2015-08-26 | 哈尔滨工业大学 | 应用于飞行场景仿真的场景生成模块及生成方法 |
CN105679102A (zh) * | 2016-03-03 | 2016-06-15 | 南京航空航天大学 | 一种全国飞行流量时空分布预测推演***及方法 |
CN113257042A (zh) * | 2021-05-18 | 2021-08-13 | 中国电子科技集团公司第二十八研究所 | 一种基于时空分布特征的多机场终端区航班流生成方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8781651B2 (en) * | 2012-09-21 | 2014-07-15 | Georgia Tech Research Corporation | Systems and methods providing a fuel-efficient RTA implementation with uncertain winds |
US10110881B2 (en) * | 2014-10-30 | 2018-10-23 | Microsoft Technology Licensing, Llc | Model fitting from raw time-of-flight images |
-
2022
- 2022-08-30 CN CN202211059423.1A patent/CN115472043B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104867184A (zh) * | 2015-06-12 | 2015-08-26 | 哈尔滨工业大学 | 应用于飞行场景仿真的场景生成模块及生成方法 |
CN105679102A (zh) * | 2016-03-03 | 2016-06-15 | 南京航空航天大学 | 一种全国飞行流量时空分布预测推演***及方法 |
CN113257042A (zh) * | 2021-05-18 | 2021-08-13 | 中国电子科技集团公司第二十八研究所 | 一种基于时空分布特征的多机场终端区航班流生成方法 |
Non-Patent Citations (1)
Title |
---|
面向适航审定的飞行场景研究;徐敏敏等;民用飞机设计与研究;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115472043A (zh) | 2022-12-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107529651B (zh) | 一种基于深度学习的城市交通客流预测方法和设备 | |
Chen et al. | Travel time prediction system based on data clustering for waste collection vehicles | |
Yu et al. | Headway-based bus bunching prediction using transit smart card data | |
Ma et al. | Short-term traffic flow forecasting by selecting appropriate predictions based on pattern matching | |
Wu et al. | DeepETA: a spatial-temporal sequential neural network model for estimating time of arrival in package delivery system | |
CN105678398A (zh) | 基于大数据技术的电力负荷预测方法及基于该方法的研究应用*** | |
CN115115278B (zh) | 一种离场航班推出速率动态自适应决策方法 | |
Wu et al. | An improved svm model for flight delay prediction | |
CN110751359B (zh) | 一种自动化航线网络评估方法、电子设备及存储介质 | |
CN110110339B (zh) | 一种日前水文预报误差校正方法及*** | |
CN112508261A (zh) | 一种基于神经网络的配电变压器负荷分层预测方法及装置 | |
Zhang et al. | Travel time prediction of urban public transportation based on detection of single routes | |
CN110490365B (zh) | 一种基于多源数据融合预测网约车订单量的方法 | |
CN111582603A (zh) | 基于多源信息融合的煤与瓦斯突出的智能预警方法 | |
CN115472043B (zh) | 一种基于p-中值定理的机场群航路飞行时间场景生成方法 | |
Luo et al. | A multi-output deep learning model based on Bayesian optimization for sequential train delays prediction | |
Chen et al. | A fuzzy collaborative forecasting approach based on XAI applications for cycle time range estimation | |
Xu et al. | Short-term passenger flow prediction during station closures in subway systems | |
Zjavka | Multi-site post-processing of numerical forecasts using a polynomial network substitution for the general differential equation based on operational calculus | |
CN116070385B (zh) | 一种境外矿产资源供应链风险自动识别方法及*** | |
Takacs | Predicting flight arrival times with a multistage model | |
Zhu et al. | Validating rail transit assignment models with cluster analysis and automatic fare collection data | |
US20230075475A1 (en) | Outlier Determination Among Models of Future Event Prediction | |
Grafeeva et al. | Passenger traffic analysis based on St. Petersburg public transport | |
Zhou et al. | Refined taxi demand prediction with ST-Vec |
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 |