WO2023005442A1 - 一种多维度航班放行效能评估方法 - Google Patents

一种多维度航班放行效能评估方法 Download PDF

Info

Publication number
WO2023005442A1
WO2023005442A1 PCT/CN2022/097769 CN2022097769W WO2023005442A1 WO 2023005442 A1 WO2023005442 A1 WO 2023005442A1 CN 2022097769 W CN2022097769 W CN 2022097769W WO 2023005442 A1 WO2023005442 A1 WO 2023005442A1
Authority
WO
WIPO (PCT)
Prior art keywords
flight
control area
flights
calculate
departure
Prior art date
Application number
PCT/CN2022/097769
Other languages
English (en)
French (fr)
Inventor
黄吉波
丁辉
田靖
童明
董斌
王煊
包旭
陈曦
Original Assignee
中国电子科技集团公司第二十八研究所
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 中国电子科技集团公司第二十八研究所 filed Critical 中国电子科技集团公司第二十八研究所
Priority to US17/822,812 priority Critical patent/US11710416B2/en
Publication of WO2023005442A1 publication Critical patent/WO2023005442A1/zh

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/0095—Aspects of air-traffic control not provided for in the other subgroups of this main group
    • 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
    • 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

Definitions

  • the invention belongs to the field of air traffic management, in particular to a multi-dimensional flight release efficiency evaluation method.
  • the technical problem to be solved by the present invention is to provide a multi-dimensional flight release performance evaluation method for the deficiency of the technical method for evaluating the flight release effect.
  • This method obtains civil aviation operation production data through the data exchange platform, starts from the perspectives of flight release stability, scope of influence, delay distribution, fairness, etc., covers macro-level overall operation effects, micro-regional efficiency evaluation methods, and constructs a set of multi-dimensional A collection of flight release performance evaluation indicators, calculate and analyze flight jump variables, delays, normal rate and other indicators, and use lists, histograms, line charts, radar charts and other forms to visualize the evaluation indicators.
  • the present invention discloses a multi-dimensional flight release performance evaluation method, including the following steps:
  • Step 1 analyze and process operational information such as airspace capacity and flight plan
  • Step 2 Based on the flight plan information, identify flight objects affected by flow control, flight objects restricted by flow control, and flight objects departing from various civil aviation control areas;
  • Step 3 compare and analyze the matching situation between the flight flow of the specified airspace object and the airspace capacity
  • Step 4 taking the civil aviation control area as the object, analyze the release time jump situation of the flights departing from each control area, including: jump sorties, jump levels;
  • Step 5 Taking the civil aviation control area as the object, analyze the distribution of release delays of flights departing from each control area, including: maximum delay, minimum delay, total delay, number of delayed flights, and average delay;
  • Step 6 Taking the civil aviation control area as the object, analyze the normal rate of departure flights from each control area, including: normal flight number, abnormal flight number, and normal rate;
  • Step 7 Analyze the fairness of release of flights departing from civil aviation control areas, including: delay fairness, jump fairness, and restriction fairness;
  • Step 8 analyze the overall performance of controlled flight release, including: delay category, stability category, impact category.
  • the delay category includes: total delays, total delays, and average delays
  • the stability category includes: total jumps, total jumps, and average jumps
  • the impact category includes: affected flights, restricted flights, and restricted proportions.
  • Step 1 includes the following steps:
  • Step 1-1 Analyze the basic flight plan information FlightInfo, including: flight number Callsign, departure airport DepAirport, landing airport DesAirport, departure area DepRegion, planned block removal time Sobt, planned landing time Sldt, and the calculation obtained from the latest release calculation Block-off time CurCobt, calculated departure time CurCtot from the latest release calculation, last block-off time LastCobt calculated from the last release calculation, last calculated departure time LastCtot, flight affected by flow control flag Affected, flight restricted by flow control Mark Restricted; the comprehensive information of the ith flight is recorded as FlightInfo i ;
  • Step 1-3 Analyze the capacity information Capacity of the airspace object to be evaluated, including: airspace object name Airsapce, airspace capacity Capacity.
  • Step 2 includes the following steps:
  • Step 2-1 Traverse the basic information FlightInfo of each flight plan
  • Step 2-2 If DepRegion is Xinjiang in the basic information of the flight plan traversed, then identify the flight as the departure flight object of the Xinjiang control area, and put the flight into the Xinjiang control area departure flight collection FlightSet 1 ; if DepRegion is Northeast, then Identify this flight as the departure flight object of the Northeast control area, and put this flight into the flight set FlightSet 2 of the Northeast control area; if DepRegion is northwest, then identify this flight as the departure flight object of the Northwest control area, and put this flight into Northwest Flight Set 3 of the departure flight collection of the control area; if DepRegion is southwest, then identify this flight as the departure flight object of the central-south control area, and put this flight into the collection FlightSet 4 of departure flight flights of the southwest control area; if DepRegion is East China, then identify this flight as The flight object departing from the East China control area, and put the flight into the FlightSet 5 of the departure flight set of the East China control area; if
  • Step 2-3 Traversing the flight set FlightSet i departing from the control area, where i takes a value from 1 to 7;
  • Step 2-4 Traverse each flight information FlightInfo j in FlightSet i , where j takes a value from 1 to N i , where N i represents the total number of flights departing from the corresponding control area;
  • Step 3 includes the following steps:
  • Step 3-1 Analyze the set evaluation parameters, including: evaluation airspace object name Airsapce, evaluation start time BgnTime, evaluation time span TimeSpan, evaluation period number TimeNum;
  • Step 3-3 Extract the period capacity information of the airspace object to be evaluated from the capacity information Capacity of the airspace object to be evaluated obtained by analysis, and obtain the airspace capacity C i of the i-th time slice;
  • Step 3-4 traverse the CrossAirsapce j in the four-dimensional trajectory information of the flight
  • i Flow i + 1, where Flow i represents the predicted flow of the i-th evaluation time slice;
  • Step 3-6 After obtaining the airspace forecast flow of each time slice, compare the matching situation between the forecast flow and the airspace capacity, and calculate the excess operation level OverFlowLv i of each time slice:
  • Step 4 includes the following steps:
  • Step 4-1 Traversing the flight set FlightSet i departing from each control area of civil aviation;
  • Step 4-2 traverse each flight information FlightInfo j in FlightSet i ;
  • ChangeLv1Num i indicates the flight number of the i-th control area jump level 1
  • ChangeLv2Num i indicates the i-th control area jump level 2 flight sorties
  • ChangeLv3Num i indicates the flight sorties of the i-th control area jump level 3
  • ChangeLv4Num i indicates the flight sorties of the i-th control area jump level 4
  • VSP 1 indicates the jump level interval threshold (unit: minute) ;
  • Step 4-5 Calculate the total jump variable ChangeTotal i of the flight departing from the corresponding control area
  • Step 5 includes the following steps:
  • Step 5-1 Traversing the flight set FlightSet i departing from each control area of civil aviation;
  • Step 5-2 traverse each flight information FlightInfo j in FlightSet i ;
  • DelayMax i MAX(DelayMax i ,CurCobt j -Sobt j ,0);
  • DelayMin i Min(DelayMin i ,Max(CurCobt j -Sobt j ,0));
  • Step 5-6 Calculate the total flight delay DelayTotal i in the corresponding control area, which represents the sum of flight delays in the control area:
  • Step 5-7 Calculate the average delay DelayAve i of the departure flight in the corresponding control area, indicating the average delay caused by the delayed flight in the control area:
  • Step 6 includes the following steps:
  • Step 6-1 Traversing the flight set FlightSet i departing from each control area of civil aviation;
  • Step 6-2 traverse each flight information FlightInfo j in FlightSet i ;
  • Step 7 includes the following steps:
  • Step 7-1 Calculate the average DelayAveAve of the average delay of departure flights in each control area:
  • Step 7-2 Calculate the delay fairness index DelayFairness, which represents the standard deviation of the average delay of departure flights in each control area:
  • Step 7-3 Calculate the average jump variable ChangeAve of flights departing from each control area:
  • Step 7-4 Calculate the jump fairness index ChangeFairness, which represents the standard deviation of the jump variable of flights departing from each control area:
  • Step 7-5 Calculate the average number of flights restricted by flow control RestrictedAve of flights departing from each control area:
  • Step 7-6 Calculate the controlled fairness index RestrictedFairness, which represents the standard deviation of the number of flights departing from each control area subject to flow control restrictions:
  • Step 8 includes the following steps:
  • Step 8-1 Calculate the total change number of flights ChangeNum, which means the number of flights with different off-block times obtained in two adjacent release calculations of all flights,
  • Step 8-2 Calculate the total jump variable ChangeTotal of the flight, which represents the cumulative absolute difference of the OBT calculated in two adjacent release calculations of all flights,
  • Step 8-3 Calculate the flight average jump variable ChangeTotalAve, which represents the average absolute difference of the OBT calculated in two adjacent release calculations of all flights:
  • Step 8-4 Calculate the total number of delayed flights DelayNum, which means the number of flights whose calculated off-block time is later than the planned off-block time,
  • Step 8-5 Calculate the total flight delay DelayTotal, which represents the cumulative difference between the calculated OBT of all flights and the planned OBT,
  • Step 8-6 Calculate the average flight delay DelayTotalAve, which means the average difference between the calculated OBT and the planned OBT for all flights:
  • Step 8-7 Calculate the total number of flights affected by AffectedNum, indicating the number of flights affected by flow control:
  • Step 8-8 Calculate the total number of controlled flights RestrictedNum, indicating the number of flights restricted by flow control:
  • Step 8-9 Calculate the proportion of flights controlled by RestrictedRato, indicating the proportion of flights restricted by flow control:
  • Steps 8-10 Change the total number of flight changes ChangeNum, total flight changeTotal, flight average changeTotalAve, total flight delay DelayNum, total flight delay DelayTotal, average flight delay DelayTotalAve, total affected flight number AffectedNum, flight
  • Normalization is a way of simplifying calculations, that is, transforming dimensional expressions into dimensionless expressions.
  • the above-mentioned evaluation values are converted into proportions.
  • the results of flight release performance evaluation indicators are displayed in a variety of graphical ways. Objects of flights subject to flow control restrictions are displayed in a list, the matching of airspace flow and capacity is displayed in a histogram and line graph overlay, the flight jump degree index is displayed in a circular histogram and a table overlay, and the flight delay distribution index is displayed in a histogram and Line segments are superimposed to display, flight fairness indicators are displayed by cards, flight regularity indicators are displayed by histogram and table superposition, and overall performance indicators are displayed by radar charts.
  • Figure 1 is a schematic flow chart of the airport capacity-demand balance forecasting analysis method.
  • Fig. 2 is a schematic diagram of a list of flight information restricted by flow control.
  • Figure 3 is a schematic diagram of the matching analysis of airspace flight flow and capacity.
  • Figure 4 is a schematic diagram of the analysis of flight release jump degree.
  • Figure 5 is a schematic diagram of the distribution analysis of flight release delays.
  • Figure 6 is a schematic diagram of flight release fairness analysis.
  • Figure 7 is a schematic diagram of flight release normality analysis.
  • Figure 8 is a schematic diagram of the overall performance analysis of flight release.
  • the dates of the relevant time information are all 2021/07/06.
  • the basic information of the flight plan is obtained and analyzed through the integrated operation information interaction platform, which mainly includes: flight call sign, departure airport, landing airport, calculation of off-block time, calculation of off-block time, calculation of departure time, Take-off area, marks affected by flow control, marks restricted by flow control, etc.
  • the basic information details of a flight plan CSN3999, ZGGG, ZBAA, 13:05, 13:05, 13:35, Zhongnan, 1, 1.
  • the acquisition and analysis of flight four-dimensional trajectory prediction information is a set of airspace objects that the flight flies through, including: the name of the airspace object that flies through, and the time of entering the airspace object.
  • the details of flight CSZ1156 four-dimensional trajectory prediction information ZBAA, 08:15; ZBACC02, 08:21; ...; ZSNJ, 10:45.
  • airspace capacity information is obtained and analyzed, including: airspace object name, capacity start time, capacity end time, capacity time granularity, and capacity value.
  • airspace ZBACC10 capacity information details: 08:00, 09:00, 60 minutes, 40; 09:00, 12:00, 60 minutes, 30; 12:00, 14:00, 60 minutes, 35; 14:00 , 16:00, 60 minutes, 40.
  • evaluation airspace object name Airsapce As shown in Figure 1, according to the needs of operation management, set the basic parameter information required for performance evaluation, including: evaluation airspace object name Airsapce, evaluation start time BgnTime, evaluation time span TimeSpan, evaluation period number TimeNum, jump level interval threshold VSP 1. Denotes the normal flight delay threshold VSP 2 .
  • the performance evaluation parameter details of a certain setting ZBACC10, 08:00, 60 minutes, 8, 15 minutes, 15 minutes.
  • the flight restricted by flow control is identified, the flight call sign, departure airport, landing airport, planned off-block time, calculated off-block time, Calculate the takeoff time, takeoff area, etc., and display them in a list.
  • display the details of a flight restricted by flow control CSN3999, ZGGG, ZBAA, 13:05, 13:05, 13:35, Zhongnan.
  • the attention period of the airspace object of interest is divided into several consecutive time slices, and the information of a single time slice includes: time slice start time, time slice end time.
  • time slice start time For example, a group of divided time slice details: 08:00-09:00, ..., 12:00-13:00, ..., 15:00-16:00.
  • the time slice capacity is extracted from the airspace capacity information, including: capacity start time, capacity end time, and capacity value, and displayed in the form of a broken line. For example, display the capacity information of ZBACC10 airspace time slice: 08:00-09:00, 40 sorties; ...; 12:00-13:00, 35 sorties; .; 15:00-16:00, 40 sorties .
  • the flight flow of the time slice flying through the airspace of interest is predicted, including: forecast start time, forecast end time, flow forecast value, and the histogram form display.
  • display a group of ZBACC10 airspace flight flow forecast information: 08:00-09:00, 42 sorties; 7-8; 12:00-13:00, 35 sorties; 7-8; 14:00-15:00, 42 sorties.
  • the excess level of each time slice is displayed, and the traffic forecast histogram is used to distinguish and display: the excess level 0 is displayed in gray, and the excess level 1 is displayed Green, yellow for excess level 2, orange for excess level 3, red for excess level 4.
  • display a set of accommodation matching information of ZBACC10 airspace 08:00-09:00 excess level is 2, flow histogram is yellow; ...; 12:00-13:00 excess level is 1, flow column The graph is green; certainly; 15:00-16:00 The excess level is 2, and the traffic histogram is yellow.
  • the jumping conditions of flights departing from each control area are displayed in tabular form, including: the control area to which they belong, the number of jumps, and the level of jumps.
  • display the jump situation of the North China control area jump level 1 is 3 sorties, jump level 2 is 1 sortie, jump level 3 is 14 sorties, and jump level 4 is 15 sorties.
  • the characteristic delay distribution of departure flights in each control area including: control area, maximum delay, minimum delay, and average delay.
  • the maximum delay and the minimum delay are displayed in the form of the top/bottom of the histogram, and the average delay is displayed in the form of a line segment.
  • display the characteristic delay distribution information of the central and southern control area the maximum delay is 50 minutes, the minimum delay is 0 minutes, and the average delay is 22 minutes.
  • the fairness of flight delays including: the mean value of the average delay in the control area, the standard deviation of the average delay in the control area, and display it in the form of a card. For example: displaying a delay fairness indicator: the mean value of the average delay in the control area is 18.5 minutes, and the standard deviation of the average delay in the control area is 7.6 minutes.
  • the fairness of flight restrictions including: the average number of flights restricted by flow control in the control area, and the standard deviation of the number of flights restricted by flow control in the control area, and display them in the form of cards.
  • Display a restricted fairness index the average number of sorties restricted by flow control in the control area is 106 sorties, and the standard deviation of the number of sorties restricted by flow control in the control area is 50.6 sorties.
  • each control area including: the control area to which it belongs, the number of normal flights, the number of abnormal flights, and the normal rate, and display it in the form of a table. For example, it shows the normality index of East China: 165 normal flights, 28 abnormal flights, and a normal rate of 85.49%.
  • the overall delay of flight release including: total number of delayed flights, total delays, and average delays, and display them in tabular form. For example, a set of overall performance indicators for delays is displayed: the total number of delayed sorties is 193, the total delayed amount is 3904 minutes, and the average delayed amount is 20.2 minutes.
  • the indicators of the above three categories and nine subcategories are normalized to generate a dimensionless value with a value range of [0, 100], which is displayed in the form of a radar chart.
  • a value range of [0, 100] which is displayed in the form of a radar chart.
  • stable indicators form a yellow quadrilateral area
  • delay indicators form a blue quadrilateral area
  • impact indicators form a red quadrilateral area (since the drawings in the manual can only be grayscale images, the The color cannot be seen accurately, hereby explain).
  • the invention realizes multi-dimensional evaluation and analysis of flight release efficiency such as calm flow matching, jump degree, delay distribution, fairness, and normality, and adopts forms such as tables, circular histograms, traditional histograms, line graphs, and radar graphs.
  • the visual graphic display of the evaluation indicators is comprehensive and intuitive, and provides method support for flight release performance analysis in the field of civil aviation flight flow management.
  • the present invention provides a multi-dimensional flight release performance evaluation method. There are many methods and approaches to realize this technical solution.
  • the above description is only a preferred embodiment of the present invention. , under the premise of not departing from the principle of the present invention, some improvements and modifications can also be made, and these improvements and modifications should also be regarded as the protection scope of the present invention. All components that are not specified in this embodiment can be realized by existing technologies.

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明提供了一种多维度航班放行效能评估方法，通过业务信息综合处理平台，获取空管生产运行数据，主要包括：空域容量信息、航班计划基本信息、航班四维轨迹信息等；通过对运行数据的处理，识别受流控影响航班对象及受流控限制航班对象、分析航班放行时间跳变程度，计算航班放行延误分布，评价受控航班放行公平性，预测受控航班放行正常率，对比空域容流匹配情况，建立一套多维度航班放行效能评估指标集，并采用列表、柱状图、折线图、雷达图等方式对评估指标进行可视化展示。

Description

一种多维度航班放行效能评估方法 技术领域
本发明属于空中交通管理领域，特别涉及一种多维度航班放行效能评估方法。
背景技术
针对民航空中交通运行有序化管理的实际需求，国际民航组织发布了航空系统组块升级计划ASBU，我国民航局发布了中国航空系统组块升级发展与实施策略，其中均提出飞行流量与空域容量相适应的核心管理概念，意义在于平衡飞行需求与服务能力，通过对航班进行放行调度，达到疏导空中交通拥堵，缓解航班延误的目的。因此，各级空中交通管理局建设了相应的流量管理系统、协同放行系统等，围绕如何科学合理地对航班进行放行管理，提供了对于空中交通运行各阶段的监视、预测、分析、决策、实施、评估等功能，保障航班安全、效率、经济运行。
由于民航运行环境具有复杂性高、时变性强的特点，航班的放行管理决策始终处于不断更新，不断适应的迭代过程中，如何全面客观、科学合理地对航班放行方案进行评估分析，成为当前运行管理的焦点问题。国内外已有相关的技术方法及应用系统主要偏重于对飞行流量以及航班放行延误的预测分析，评价方法单一，系统功能局限，目前尚缺少一种针对航班放行效能的多维度的评估分析方法及系统。
发明内容
发明目的：本发明所要解决的技术问题是针对航班放行效果的评估技术方法的不足，提供一种多维度航班放行效能评估方法。该方法通过数据交换平台获取民航运行生产数据，从航班放行的稳定程度、影响范围、延误分布、公平情况等视角入手，涵盖宏观总体运行效果，微观区域性效能的评估方法，构建一套多维度航班放行效能评估指标集合，计算分析航班跳变量、延误量、正常率等指标，并采用列表、柱状图、折线图、雷达图等形式对评估指标进行可视化展示。
技术方案：为实现上述发明目的，本发明公开了一种多维度航班放行效能评估方法，包括以下步骤：
步骤1，解析处理空域容量、航班计划等运行信息；
步骤2，根据航班计划信息，识别受流控影响航班对象、受流控限制航班对象、民航各管制区域起飞航班对象；
步骤3，对指定空域对象的飞行流量与空域容量的匹配情况进行对比分析；
步骤4，以民航管制区域为对象，分析各管制区域起飞航班的放行时间跳变情况，包括：跳变架次、跳变等级；
步骤5，以民航管制区域为对象，分析各管制区域起飞航班的放行延误的分布情况，包括：最大延误、最小延误、总延误、延误架次、平均延误；
步骤6，以民航管制区域为对象，分析各管制区域起飞航班的放行正常率情况，包括：正常航班架次、非正常航班架次、正常率；
步骤7，分析民航管制区域起飞航班的放行公平性情况，包括：延误公平性、跳变公平性、限制公平性；
步骤8，分析受控航班放行总体效能，包括：延误类、稳定类、影响类。其中，延误类包含：总延误架次、总延误、平均延误；稳定类包含：总跳变架次、总跳变量、平均跳变量；影响类包含：受影响架次、受限制架次、受限制占比。
步骤1包括以下步骤：
步骤1-1：解析航班计划基本信息FlightInfo，包括：航班号Callsign、起飞机场DepAirport、降落机场DesAirport、起飞区域DepRegion、计划撤轮档时间Sobt、计划降落时间Sldt、最新一次放行计算得到的计算撤轮档时间CurCobt、最新一次放行计算得到的计算起飞时间CurCtot、上一次放行计算得到的计算撤轮档时间LastCobt、上次计算起飞时间LastCtot、航班受流控影响标记Affected、航班受流控限制标记Restricted；第i个航班综合信息记为FlightInfo i；
步骤1-2：解析航班四维轨迹信息Flight4DT＝{CrossAirsapce 1，……，CrossAirsapce j}，其中CrossAirsapce j是航班飞行所经过的第j个空域对象，包含：飞经空域对象名称AirsapceName j、进入空域时间EntreTime j；第i个航班综合信息记为Flight4DT i；
步骤1-3：解析需评估空域对象的容量信息Capacity，包括：空域对象名称Airsapce、空域容 量Capacity。
步骤2包括以下步骤：
步骤2-1：遍历每一个航班计划基本信息FlightInfo；
步骤2-2：如果所遍历航班计划基本信息中DepRegion为新疆，则识别该航班为新疆管制区域起飞航班对象，并将该航班放入新疆管制区起飞航班集合FlightSet 1；如果DepRegion为东北，则识别该航班为东北管制区域起飞航班对象，并将该航班放入东北管制区起飞航班集合FlightSet 2；如果DepRegion为西北，则识别该航班为西北管制区域起飞航班对象，并将该航班放入西北管制区起飞航班集合FlightSet 3；如果DepRegion为西南，则识别该航班为中南管制区域起飞航班对象，并将该航班放入西南管制区起飞航班集合FlightSet 4；如果DepRegion为华东，则识别该航班为华东管制区域起飞航班对象，并将该航班放入华东管制区起飞航班集合FlightSet 5；如果DepRegion为华北，则识别该航班为华北管制区域起飞航班对象，并将该航班放入华北管制区起飞航班集合FlightSet 6；如果DepRegion为中南，则识别该航班为中南管制区域起飞航班对象，并将该航班放入中南管制区起飞航班集合FlightSet 7；
步骤2-3：遍历管制区起飞航班集合FlightSet i，i取值为1～7；
步骤2-4：遍历FlightSet i中的每一个航班信息FlightInfo j，j取值为1～N i，其中N i表示相应管制区起飞航班总架次；
步骤2-5：如果航班受流控影响标记Affected＝1，则识别该航班为受流控影响航班对象；如果航班受流控限制标记Restricted＝1，则识别该航班为受流控限制航班对象；
步骤2-6：计算相应管制区起飞航班的受流控影响架次AffectedNum i，如果航班受流控影响标记Affected＝1，则AffectedNum i＝AffectedNum i+1，并将航班信息放入受流控影响航班集合AffectedSet；
步骤2-7：计算相应管制区起飞航班的受流控限制架次RestrictedNum i，如果航班受流控限制标记Restricted＝1，则RestrictedNum i＝RestrictedNum i+1并将航班信息放入受流控限制航班集合RestrictedSet。
步骤3包括以下步骤：
步骤3-1：解析所设置的评估参数，包括：评估空域对象名称Airsapce、评估开始时间BgnTime、评估时间跨度TimeSpan、评估时段数量TimeNum；
步骤3-2：对评估时段进行划分，生成连续评估时间片，放入评估时间片集合TimeSpanSet＝{[BT 1,ET 1),..........,[BT i,ET i)},i∈[1,TimeNum]，其中，BT i＝BgnTime+TimeSpan*(i-1)，表示第i个评估时间片的开始时间；ET i＝BgnTime+TimeSpan*i，表示第i个评估时间片的结束时间；
步骤3-3：从解析得到的需评估空域对象的容量信息Capacity中，提取所要评估空域对象的时段容量信息，得到第i个时间片的空域容量C i；
步骤3-4：遍历航班四维轨迹信息中的CrossAirsapce j；
步骤3-5：如果AirsapceName＝Airsapce，EntreTime j∈[BT i,ET i)，表示航班飞经空域对象名称与所需评估空域对象相同，如果进入空域时间在时间片起止时段内时，则Flow i＝Flow i+1，其中Flow i表示第i个评估时间片的预测流量；
步骤3-6：得到每个时间片的空域预测流量后，对比预测流量与空域容量的匹配情况，计算每个时间片的超量运行等级OverFlowLv i：
Figure PCTCN2022097769-appb-000001
步骤4包括以下步骤：
步骤4-1：遍历民航各管制区起飞航班集合FlightSet i；
步骤4-2：遍历FlightSet i中的每一个航班信息FlightInfo j；
步骤4-3：如果航班CurCtot-LastCtot∈(-∞,-VSP 1]，则ChangeLv1Num i＝ChangeLv1Num i+1；如果航班CurCtot-LastCtot∈(-VSP 1,0)，则ChangeLv2Num i＝ChangeLv2Num i+1；如果航班CurCtot-LastCtot∈(0,VSP 1]，则ChangeLv3Num i＝ChangeLv3Num i+1；如果航班CurCtot-LastCtot∈(VSP 1,+∞)，则ChangeLv4Num i＝ChangeLv1Num i+1；其中，CurCtot-LastCtot表示最新两次放行计算得到的计算起飞时间的差值(即航班跳变量)，ChangeLv1Num i表示第i个管制区跳变等级1的航班架次，ChangeLv2Num i表示第i个管制区跳变等级2的航班架次，ChangeLv3Num i表示第i个管制区跳变等级3的航班架次，ChangeLv4Num i表示第i个管制区跳变等级4的航班架次，VSP 1表示跳变等级区间阈值(单位：分钟)；
步骤4-4：计算相应管制区起飞航班的跳变总架次ChangeNum i，ChangeNum i＝ChangeLv1Num i+ChangeLv2Num i+ChangeLv3Num i+ChangeLv4Num i；
步骤4-5：计算相应管制区起飞航班的总跳变量ChangeTotal i，
Figure PCTCN2022097769-appb-000002
步骤5包括以下步骤：
步骤5-1：遍历民航各管制区起飞航班集合FlightSet i；
步骤5-2：遍历FlightSet i中的每一个航班信息FlightInfo j；
步骤5-3：计算相应管制区的航班延误架次DelayNum i，如果CurCobt j-Sobt j>0，则DelayNum i＝DelayNum i+1，表示管制区起飞的航班计算撤轮挡时间大于计划撤轮挡时间，则为管制区延误航班，其中CurCobt j-Sobt j表示航班延误量；
步骤5-4：初始化DelayMax i＝0，计算相应管制区的航班最大延误DelayMax i，如果CurCobt j-Sobt j>DelayMax i，则DelayMax i＝CurCobt j-Sobt，也表示为：
DelayMax i＝MAX(DelayMax i,CurCobt j-Sobt j,0)；
步骤5-5：初始化DelayMin i＝0，计算相应管制区的航班最小延误DelayMin i，如果CurCobt j-Sobt j<DelayMin i，则DelayMin i＝CurCobt j-Sobt j，也表示为：
DelayMin i＝Min(DelayMin i,Max(CurCobt j-Sobt j,0))；
步骤5-6：计算相应管制区的航班总延误DelayTotal i，表示该管制区航班延误总和：
Figure PCTCN2022097769-appb-000003
步骤5-7：计算相应管制区起飞航班平均延误DelayAve i，表示管制区延误航班所产生的平均延误量：
Figure PCTCN2022097769-appb-000004
步骤6包括以下步骤：
步骤6-1：遍历民航各管制区起飞航班集合FlightSet i；
步骤6-2：遍历FlightSet i中的每一个航班信息FlightInfo j；
步骤6-3：计算第i个管制区起飞的正常航班架次NormalNum i，如果CurCobt j-Sobt j∈(-∞,VSP 2]，则NormalNum i＝NormalNum i+1，表示航班延误不大于VSP 2的航班为正常航班；
步骤6-4：计算非正常航班架次UnNormalNum i，如果CurCobt j-Sobt j∈(VSP 2,+∞)，则UnNormalNum i＝UnNormalNum i+1，表示航班延误大于VSP 2的航班为非正常航班；
步骤6-5：航班遍历结束后，计算相应管制区起飞航班正常率NormalRate i，表示管制区起飞航班中正常航班架次的比例NormalRate i＝NormalNum i/N i。
步骤7包括以下步骤：
步骤7-1：计算各管制区域起飞航班平均延误的平均值DelayAveAve：
Figure PCTCN2022097769-appb-000005
步骤7-2：计算延误公平性指标DelayFairness，表示各大管制区起飞航班平均延误的标准差：
Figure PCTCN2022097769-appb-000006
步骤7-3：计算各管制区域起飞航班的平均跳变量ChangeAve：
Figure PCTCN2022097769-appb-000007
步骤7-4：计算跳变公平性指标ChangeFairness，表示各大管制区起飞航班跳变量的标准差：
Figure PCTCN2022097769-appb-000008
步骤7-5：计算各各管制区域起飞航班的平均受流控限制架次RestrictedAve：
Figure PCTCN2022097769-appb-000009
步骤7-6：计算受控公平性指标RestrictedFairness，表示各各管制区域起飞航班受流控限制架次的标准差：
Figure PCTCN2022097769-appb-000010
步骤8包括以下步骤：
步骤8-1：计算航班总跳变架次ChangeNum，表示所有航班相邻两次放行计算中所得计算撤轮挡时间不同的架次，
Figure PCTCN2022097769-appb-000011
步骤8-2：计算航班总跳变量ChangeTotal，表示所有航班相邻两次放行计算中所得计算撤轮挡时间的累积绝对差值，
Figure PCTCN2022097769-appb-000012
步骤8-3：计算航班平均跳变量ChangeTotalAve，表示所有航班相邻两次放行计算中所得计算撤轮挡时间的平均绝对差值：
Figure PCTCN2022097769-appb-000013
步骤8-4：计算航班总延误架次DelayNum，表示所有航班计算撤轮挡时间晚于计划撤轮挡时间的架次，
Figure PCTCN2022097769-appb-000014
步骤8-5：计算航班总延误量DelayTotal，表示所有航班计算撤轮挡时间晚于计划撤轮挡时间 的累积差值，
Figure PCTCN2022097769-appb-000015
步骤8-6：计算航班平均延误量DelayTotalAve，表示所有航班计算撤轮挡时间晚于计划撤轮挡时间的平均差值：
Figure PCTCN2022097769-appb-000016
步骤8-7：计算航班总影响架次AffectedNum，表示受流控影响的航班架次：
Figure PCTCN2022097769-appb-000017
步骤8-8：计算航班总受控架次RestrictedNum，表示受流控限制的航班架次：
Figure PCTCN2022097769-appb-000018
步骤8-9：计算航班受控占比RestrictedRato，表示受流控限制的航班架次占比：
Figure PCTCN2022097769-appb-000019
步骤8-10：对航班总跳变架次ChangeNum、航班总跳变量ChangeTotal、航班平均跳变量ChangeTotalAve、航班总延误架次DelayNum、航班总延误量DelayTotal、航班平均延误量DelayTotalAve、航班总影响架次AffectedNum、航班总受控架次RestrictedNum、航班受控占比RestrictedRato进行归一化处理，采用雷达图对总体效能指标进行分析。归一化处理是一种简化计算的方式，即将有量纲的表达式，经过变换，化为无量纲的表达式，这里将上述各评估数值转化为比例。
对于航班放行效能评估指标结果采用多种图形化方式进行展示。受流控限制航班对象采用列表方式展示，空域流量与容量匹配情况采用柱状图及折线图叠加方式展示，航班跳变程度指标采用环形柱状图及表格叠加方式展示，航班延误分布指标采用柱状图及线段叠加方式展示，航班公平性指标采用卡片方式展示，航班正常性指标采用柱状图及表格叠加方式展示，总体效能指标采用雷达图方式展示。
有益效果：本发明显著的优点包括：
1、为保障空中交通运行提供了一套数字化评估指标集合；
2、为航班放行管理提供了一种多维度效能评估方法；
3、提升指标结果可视化显示能力，采用多元组合方式增加直观性、综合性；
4、为航班放行决策技术研究及系统研制提供了有力支撑。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明，本发明的上述和/或其他方面的优点将会变得更加清楚。
图1是机场容需平衡预测分析方法流程示意图。
图2是受流控限制航班信息列表示意图。
图3是空域飞行流量与容量匹配分析示意图。
图4是航班放行跳变程度分析示意图。
图5是航班放行延误分布分析示意图。
图6是航班放行公平性分析示意图。
图7是航班放行正常性分析示意图。
图8是航班放行总体效能分析示意图。
具体实施方式
实施例
本实施例中，相关时间信息的日期均为2021/07/06，为获取生产运行数据，首先构建综合运行 信息交互平台，解析航班飞行计划基本信息、航班四维轨迹预测信息，以及空域容量信息；根据运行管理需要，可以对所需关注的空域对象进行配置，设置关注的开始时间、时间颗粒度以及时段数量等；识别受流控限制的航班对象；预测所配置空域的飞行流量态势，分析飞行流量与容量的匹配情况，超量程度等；对于最近两次航班放行计算结果进行对比，分析航班放行结果的跳变情况，跳变架次、跳变等级等；分析航班放行的延误分布情况，最大/最小延误、延误架次、平均延误等；评估航班放行正常性情况，正常航班架次、非正常航班架次、正常航班占比等；评价航班放行公平性情况，分析跳变公平性、延误公平性、限制公平性等；从不同方面计算航班放行总体效能指标，并采用归一法进行综合评价。
下面结合图1～图8，详细说明基于本发明方法实现的一种航班放行效能评估方法处理流程步骤：
1、获取综合运行信息
1-1、处理航班计划基本信息
如图1所示，通过综合运行信息交互平台获取并解析航班飞行计划基本信息，主要包含：航班呼号、起飞机场、降落机场、计算撤轮挡时间、计算撤轮挡时间、计算起飞时间、起飞区域、受流控影响标记、受流控限制标记等。例如，某条航班飞行计划基本信息详情：CSN3999、ZGGG、ZBAA、13:05、13:05、13:35、中南、1、1。
1-2、处理航班四维轨迹预测信息
如图1所示，获取并解析航班四维轨迹预测信息，是一组航班飞经的空域对象集合，包括：飞经空域对象名称、进入空域对象时间。例如，航班CSZ1156四维轨迹预测信息详情：ZBAA，08:15；ZBACC02，08:21；……；ZSNJ，10:45。
1-3、处理空域容量信息
如图1所示，获取并解析空域容量信息，包括：空域对象名称、容量开始时间、容量结束时间、容量时间颗粒度、容量值。例如，空域ZBACC10容量信息详情：08:00、09:00、60分钟、40；09:00、12:00、60分钟、30；12:00、14:00、60分钟、35；14:00、16:00、60分钟、40。
2、设置效能评估参数
如图1所示，根据运行管理需要，设置效能评估所需的基本参数信息，包括：评估空域对象名称Airsapce、评估开始时间BgnTime、评估时间跨度TimeSpan、评估时段数量TimeNum、跳变等级区间阈值VSP 1、表示正常航班延误阈值VSP 2。例如，某次设置的效能评估参数详情：ZBACC10、08:00、60分钟、8、15分钟、15分钟。
3、受流控限制航班分析及图形化显示
如图2所示，根据航班飞行计划基本信息中的受流控限制标记，识别出受流控限制航班，航班呼号、起飞机场、降落机场、计划撤轮挡时间、计算撤轮挡时间、计算起飞时间、起飞区域等，并以列表形式进行展示。例如，显示一条受流控限制航班详情：CSN3999、ZGGG、ZBAA、13:05、13:05、13:35、中南。
4、容流匹配分析及图形化显示
4-1、时间片划分
如图3所示，根据所设置的效能评估参数，将关注空域对象的关注时段划分成若干个连续时间片，单个时间片信息包括：时间片开始时间、时间片结束时间。例如，一组划分后的时间片详情：08:00-09:00、…..、12:00-13:00、…..、15:00-16:00。
4-2、空域容量划分
如图3所示，根据划分好的时间片信息，从空域容量信息中提取时间片容量，包括：容量开始时间、容量结束时间、容量值，并以折线形式进行展示。例如，显示ZBACC10空域时间片容量信息：08:00-09:00、40架次；…..；12:00-13:00、35架次；…..；15:00-16:00、40架次。
4-3、飞行流量预测
如图3所示，根据划分好的时间片信息以及航班四维轨迹信息，对时间片飞经关注空域的飞行流量进行预测，包括：预测开始时间、预测结束时间、流量预测值，并以柱状图形式展示。例如，显示一组ZBACC10空域飞行流量预测信息：08:00-09:00、42架次；…..；12:00-13:00、35架次；…..；14:00-15:00、42架次。
4-4、容流匹配分析
如图3所示，根据时间片容量以及流量预测信息，对每个时间片的超量等级进行展示，以流量 预测柱状图填充色进行区分显示：超量等级0显示灰色，超量等级1显示绿色，超量等级2显示黄色，超量等级3显示橙色，超量等级4显示红色。例如，显示ZBACC10空域的一组容留匹配信息：08:00-09:00超量等级为2，流量柱状图为黄色；…..；12:00-13:00超量等级为1，流量柱状图为绿色；…..；15:00-16:00超量等级为2，流量柱状图为黄色。
5、航班跳变分析及图形化显示
5-1、跳变柱状图分析
如图4左侧区域所示，对比最新两次航班放行计算得到的计算撤轮挡时间，分析各大管制区起飞航班跳变情况，包括：所属管制区、跳变架次、跳变等级，并分别以环形柱状图形式进行展示，不同跳变等级的跳变架次采用柱状图填充色进行区分显示：跳变等级1显示绿色、跳变等级2显示褐色、跳变等级3显示黄色、跳变等级4显示蓝色。例如，显示华北管制区跳变情况：跳变等级1为最内圈柱状图，跳变等级4为最外圈柱状图。
5-2、跳变表格化分析
如图4右侧区域所示，以表格形式显示各大管制区起飞航班跳变情况，包括：所属管制区、跳变架次、跳变等级。例如，显示华北管制区跳变情况：跳变等级1为3架次、跳变等级2为1架次、跳变等级3为14架次、跳变等级4为15架次。
6、延误分布分析及图形化显示
6-1、特征延误分布分析
如图5上侧区域所示，分析各大管制区起飞航班特征延误分布情况，包括：所属管制区、最大延误、最小延误、平均延误。其中，最大延误及最小延误采用柱状图顶/底的形式展示，平均延误采用线段形式展示。例如，显示中南管制区特征延误分布信息：最大延误50分钟、最小延误0分钟、平均延误22分钟。
6-2、总体延误分布分析
如图5下侧区域所示，分析各大管制区起飞航班总体延误分布情况，包括：所属管制区、总延误、总延误架次，并采用左右双轴折线形式进行展示。例如：显示中南管制区总体延误分布信息：总延误40分钟、总延误架次880架次。
7、公平性分析及图形化显示
7-1、延误公平性分析
如图6左列所示，分析航班延误公平情况，包括：管制区平均延误的均值、管制区平均延误标准差，并采用卡片形式进行展示。例如：显示一份延误公平性指标：管制区平均延误的均值18.5分钟，管制区平均延误标准差7.6分钟。
7-2、跳变公平性分析
如图6中列所示，分析航班跳变公平情况，包括：管制区平均跳变、管制区跳变标准差，并采用卡片形式进行展示。例如：显示一份跳变公平性指标：管制区平均跳变295分钟，管制区跳变标准差153.7分钟。
7-3、限制公平性分析
如图6右列所示，分析航班受限制公平情况，包括：管制区平均受流控限制架次、管制区受流控限制架次标准差，并采用卡片形式进行展示。例如：显示一份限制公平性指标：管制区平均受流控限制架次106架次，管制区受流控限制架次标准差50.6架次。
8、正常性分析及图形化显示
8-1、正常性柱状图分析
如图7左侧区域所示，分析各大管制区起飞航班正常性情况，包括：所属管制区、正常率，并采用柱状图形式进行展示。例如，显示华东正常性指标：正常率85.49％。
8-2、正常性表格化分析
如图7右侧区域所示，分析各大管制区起飞航班正常性情况，包括：所属管制区、正常航班架次、非正常航班架次、正常率，并采用表格形式进行展示。例如，显示华东正常性指标：正常航班165架次、非正常航班28架次、正常率85.49％。
9、总体效能分析及图形化显示
9-1、稳定类总体效能分析
如图8右侧区域所示，分析航班放行总体稳定情况，包括：总跳变架次、总跳变量、平均跳变 量，并用表格形式进行展示。例如，显示一组稳定类总体效能指标：总跳变架次2065架次、总跳变量166分钟、平均跳变量12.4分钟。
9-2、延误类总体效能分析
如图8右侧区域所示，分析航班放行总体延误情况，包括：总延误架次、总延误量、平均延误量，并用表格形式进行展示。例如，显示一组延误类总体效能指标：总延误架次193架次、总延误量3904分钟、平均延误量20.2分钟。
9-3、影响类总体效能分析
如图8右侧区域所示，分析航班放行总体影响情况，包括：受流控影响总架次、受流控限制总架次、受流控限制占比，并用表格形式进行展示。例如，显示一组影响类总体效能指标：受流控影响总架次873架次、受流控限制总架次746架次、受流控限制占比85.45％。
9-4、归一分析
如图8左侧区域所示，将上述三大类九小类指标进行归一化处理，生成无量纲数值，数值区间为[0，100]，并采用雷达图方式进行展示。例如，显示一组总体效能指标：稳定类指标形成黄色四边形区域，延误类指标形成蓝色四边形区域，影响类指标形成红色四边形区域(由于说明书附图只能是灰度图，因此附图中的颜色不能准确看出来，特此说明)。
本发明实现了从容流匹配、跳变程度、延误分布、公平性、正常性等多维度对航班放行效能进行评估分析，并采用表格、环形柱状图、传统柱状图、折线图、雷达图等形式对评估指标进行可视化图形展示，全面直观，为民航飞行流量管理领域提供了航班放行效能分析的方法支撑。
本发明提供了一种多维度航班放行效能评估方法，具体实现该技术方案的方法和途径很多，以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。

Claims (9)

  1. 一种多维度航班放行效能评估方法，其特征在于，包括如下步骤：
    步骤1，解析处理运行信息；
    步骤2，根据航班计划信息，识别受流控影响航班对象、受流控限制航班对象、民航各管制区域起飞航班对象；
    步骤3，对指定空域对象的飞行流量与空域容量的匹配情况进行对比分析；
    步骤4，以民航管制区域为对象，分析各管制区域起飞航班的放行时间跳变情况，包括：跳变架次、跳变等级；
    步骤5，以民航管制区域为对象，分析各管制区域起飞航班的放行延误的分布情况，包括：最大延误、最小延误、总延误、延误架次、平均延误；
    步骤6，以民航管制区域为对象，分析各管制区域起飞航班的放行正常率情况，包括：正常航班架次、非正常航班架次、正常率；
    步骤7，分析民航管制区域起飞航班的放行公平性情况，包括：延误公平性、跳变公平性、受控公平性；
    步骤8，分析受控航班放行总体效能。
  2. 根据权利要求1所述的方法，其特征在于，步骤1包括以下步骤：
    步骤1-1：解析航班计划基本信息FlightInfo，包括：航班号Callsign、起飞机场DepAirport、降落机场DesAirport、起飞区域DepRegion、计划撤轮档时间Sobt、计划降落时间Sldt、最新一次放行计算得到的计算撤轮档时间CurCobt、最新一次放行计算得到的计算起飞时间CurCtot、上一次放行计算得到的计算撤轮档时间LastCobt、上次计算起飞时间LastCtot、航班受流控影响标记Affected、航班受流控限制标记Restricted；第i个航班综合信息记为FlightInfo i；
    步骤1-2：解析航班四维轨迹信息Flight4DT＝{CrossAirsapce 1，……，CrossAirsapce j}，其中CrossAirsapce j是航班飞行所经过的第j个空域对象，包含：飞经空域对象名称AirsapceName j、进入空域时间EntreTime j；第i个航班综合信息记为Flight4DT i；
    步骤1-3：解析需评估空域对象的容量信息Capacity，包括：空域对象名称Airsapce、空域容量Capacity。
  3. 根据权利要求2所述的方法，其特征在于，步骤2包括以下步骤：
    步骤2-1：遍历每一个航班计划基本信息FlightInfo；
    步骤2-2：如果所遍历航班计划基本信息中DepRegion为新疆，则识别该航班为新疆管制区域起飞航班对象，并将该航班放入新疆管制区起飞航班集合FlightSet 1；如果DepRegion为东北，则识别该航班为东北管制区域起飞航班对象，并将该航班放入东北管制区起飞航班集合FlightSet 2；如果DepRegion为西北，则识别该航班为西北管制区域起飞航班对象，并将该航班放入西北管制区起飞航班集合FlightSet 3；如果DepRegion为西南，则识别该航班为中南管制区域起飞航班对象，并将该航班放入西南管制区起飞航班集合FlightSet 4；如果DepRegion为华东，则识别该航班为华东管制区域起飞航班对象，并将该航班放入华东管制区起飞航班集合FlightSet 5；如果DepRegion为华北，则识别该航班为华北管制区域起飞航班对象，并将该航班放入华北管制区起飞航班集合FlightSet 6；如果DepRegion为中南，则识别该航班为中南管制区域起飞航班对象，并将该航班放入中南管制区起飞航班集合FlightSet 7；
    步骤2-3：遍历管制区起飞航班集合FlightSet i，i取值为1～7；
    步骤2-4：遍历FlightSet i中的每一个航班信息FlightInfo j，j取值为1～N i，其中N i表示相应管制区起飞航班总架次；
    步骤2-5：如果航班受流控影响标记Affected＝1，则识别该航班为受流控影响航班对象；如果航班受流控限制标记Restricted＝1，则识别该航班为受流控限制航班对象；
    步骤2-6：计算相应管制区起飞航班的受流控影响架次AffectedNum i，如果航班受流控影响标记Affected＝1，则AffectedNum i＝AffectedNum i+1，并将航班信息放入受流控影响航班集合AffectedSet；
    步骤2-7：计算相应管制区起飞航班的受流控限制架次RestrictedNum i，如果航班受流控限制标记Restricted＝1，则RestrictedNum i＝RestrictedNum i+1并将航班信息放入受流控限制航班集合RestrictedSet。
  4. 根据权利要求3所述的方法，其特征在于，步骤3包括以下步骤：
    步骤3-1：解析所设置的评估参数，包括：评估空域对象名称Airsapce、评估开始时间BgnTime、 评估时间跨度TimeSpan、评估时段数量TimeNum；
    步骤3-2：对评估时段进行划分，生成连续评估时间片，放入评估时间片集合TimeSpanSet＝{[BT 1,ET 1),..........,[BT i,ET i)},i∈[1,TimeNum]，其中，BT i＝BgnTime+TimeSpan*(i-1)，表示第i个评估时间片的开始时间；ET i＝BgnTime+TimeSpan*i，表示第i个评估时间片的结束时间；
    步骤3-3：从解析得到的需评估空域对象的容量信息Capacity中，提取所要评估空域对象的时段容量信息，得到第i个时间片的空域容量C i；
    步骤3-4：遍历航班四维轨迹信息中的CrossAirsapce j；
    步骤3-5：如果AirsapceName＝Airsapce，EntreTime j∈[BT i,ET i)，表示航班飞经空域对象名称与所需评估空域对象相同，如果进入空域时间在时间片起止时段内时，则Flow i＝Flow i+1，其中Flow i表示第i个评估时间片的预测流量；
    步骤3-6：得到每个时间片的空域预测流量后，对比预测流量与空域容量的匹配情况，计算每个时间片的超量运行等级OverFlowLv i：
    Figure PCTCN2022097769-appb-100001
  5. 根据权利要求4所述的方法，其特征在于，步骤4包括以下步骤：
    步骤4-1：遍历民航各管制区起飞航班集合FlightSet i；
    步骤4-2：遍历FlightSet i中的每一个航班信息FlightInfo j；
    步骤4-3：如果航班CurCtot-LastCtot∈(-∞,-VSP 1]，则ChangeLv1Num i＝ChangeLv1Num i+1；如果航班CurCtot-LastCtot∈(-VSP 1,0)，则ChangeLv2Num i＝ChangeLv2Num i+1；如果航班CurCtot-LastCtot∈(0,VSP 1]，则ChangeLv3Num i＝ChangeLv3Num i+1；如果航班CurCtot-LastCtot∈(VSP 1,+∞)，则ChangeLv4Num i＝ChangeLv1Num i+1；其中，CurCtot-LastCtot表示最新两次放行计算得到的计算起飞时间的差值，ChangeLv1Num i表示第i个管制区跳变等级1的航班架次，ChangeLv2Num i表示第i个管制区跳变等级2的航班架次，ChangeLv3Num i表示第i个管制区跳变等级3的航班架次，ChangeLv4Num i表示第i个管制区跳变等级4的航班架次，VSP 1表示跳变等级区间阈值；
    步骤4-4：计算相应管制区起飞航班的跳变总架次ChangeNum i，ChangeNum i＝ChangeLv1Num i+ChangeLv2Num i+ChangeLv3Num i+ChangeLv4Num i；
    步骤4-5：计算相应管制区起飞航班的总跳变量ChangeTotal i，
    Figure PCTCN2022097769-appb-100002
  6. 根据权利要求5所述的方法，其特征在于，步骤5包括以下步骤：
    步骤5-1：遍历民航各管制区起飞航班集合FlightSet i；
    步骤5-2：遍历FlightSet i中的每一个航班信息FlightInfo j；
    步骤5-3：计算相应管制区的航班延误架次DelayNum i，如果CurCobt j-Sobt j>0，则DelayNum i＝DelayNum i+1，表示管制区起飞的航班计算撤轮挡时间大于计划撤轮挡时间，则为管制区延误航班，其中CurCobt j-Sobt j表示航班延误量；
    步骤5-4：初始化DelayMax i＝0，计算相应管制区的航班最大延误DelayMax i，如果CurCobt j-Sobt j>DelayMax i，则DelayMax i＝CurCobt j-Sobt，也表示为：
    DelayMax i＝MAX(DelayMax i,CurCobt j-Sobt j,0)；
    步骤5-5：初始化DelayMin i＝0，计算相应管制区的航班最小延误DelayMin i，如果CurCobt j-Sobt j<DelayMin i，则DelayMin i＝CurCobt j-Sobt j，也表示为：
    DelayMin i＝Min(DelayMin i,Max(CurCobt j-Sobt j,0))；
    步骤5-6：计算相应管制区的航班总延误DelayTotal i，表示该管制区航班延误总和：
    Figure PCTCN2022097769-appb-100003
    步骤5-7：计算相应管制区起飞航班平均延误DelayAve i，表示管制区延误航班所产生的平均延误量：
    Figure PCTCN2022097769-appb-100004
  7. 根据权利要求6所述的方法，其特征在于，步骤6包括以下步骤：
    步骤6-1：遍历民航各管制区起飞航班集合FlightSet i；
    步骤6-2：遍历FlightSet i中的每一个航班信息FlightInfo j；
    步骤6-3：计算第i个管制区起飞的正常航班架次NormalNum i，如果CurCobt j-Sobt j∈(-∞,VSP 2]，则NormalNum i＝NormalNum i+1，表示航班延误不大于VSP 2的航班为正常航班；
    步骤6-4：计算非正常航班架次UnNormalNum i，如果CurCobt j-Sobt j∈(VSP 2,+∞)，则UnNormalNum i＝UnNormalNum i+1，表示航班延误大于VSP 2的航班为非正常航班；
    步骤6-5：航班遍历结束后，计算相应管制区起飞航班正常率NormalRate i，表示管制区起飞航班中正常航班架次的比例NormalRate i＝NormalNum i/N i。
  8. 根据权利要求7所述的方法，其特征在于，步骤7包括以下步骤：
    步骤7-1：计算各管制区域起飞航班平均延误的平均值DelayAveAve：
    Figure PCTCN2022097769-appb-100005
    步骤7-2：计算延误公平性指标DelayFairness，表示各大管制区起飞航班平均延误的标准差：
    Figure PCTCN2022097769-appb-100006
    步骤7-3：计算各管制区域起飞航班的平均跳变量ChangeAve：
    Figure PCTCN2022097769-appb-100007
    步骤7-4：计算跳变公平性指标ChangeFairness，表示各大管制区起飞航班跳变量的标准差：
    Figure PCTCN2022097769-appb-100008
    步骤7-5：计算各各管制区域起飞航班的平均受流控限制架次RestrictedAve：
    Figure PCTCN2022097769-appb-100009
    步骤7-6：计算受控公平性指标RestrictedFairness，表示各各管制区域起飞航班受流控限制架次的标准差：
    Figure PCTCN2022097769-appb-100010
  9. 根据权利要求8所述的方法，其特征在于，步骤8包括以下步骤：
    步骤8-1：计算航班总跳变架次ChangeNum，表示所有航班相邻两次放行计算中所得计算撤轮挡时间不同的架次，
    Figure PCTCN2022097769-appb-100011
    步骤8-2：计算航班总跳变量ChangeTotal，表示所有航班相邻两次放行计算中所得计算撤轮挡时间的累积绝对差值，
    Figure PCTCN2022097769-appb-100012
    步骤8-3：计算航班平均跳变量ChangeTotalAve，表示所有航班相邻两次放行计算中所得计算撤轮挡时间的平均绝对差值：
    Figure PCTCN2022097769-appb-100013
    步骤8-4：计算航班总延误架次DelayNum，表示所有航班计算撤轮挡时间晚于计划撤轮挡时间的架次，
    Figure PCTCN2022097769-appb-100014
    步骤8-5：计算航班总延误量DelayTotal，表示所有航班计算撤轮挡时间晚于计划撤轮挡时间的累积差值，
    Figure PCTCN2022097769-appb-100015
    步骤8-6：计算航班平均延误量DelayTotalAve，表示所有航班计算撤轮挡时间晚于计划撤轮挡时间的平均差值：
    Figure PCTCN2022097769-appb-100016
    步骤8-7：计算航班总影响架次AffectedNum，表示受流控影响的航班架次：
    Figure PCTCN2022097769-appb-100017
    步骤8-8：计算航班总受控架次RestrictedNum，表示受流控限制的航班架次：
    Figure PCTCN2022097769-appb-100018
    步骤8-9：计算航班受控占比RestrictedRato，表示受流控限制的航班架次占比：
    Figure PCTCN2022097769-appb-100019
    步骤8-10：对航班总跳变架次ChangeNum、航班总跳变量ChangeTotal、航班平均跳变量ChangeTotalAve、航班总延误架次DelayNum、航班总延误量DelayTotal、航班平均延误量DelayTotalAve、航班总影响架次AffectedNum、航班总受控架次RestrictedNum、航班受控占比RestrictedRato进行归一化处理，采用雷达图对总体效能指标进行分析。
PCT/CN2022/097769 2021-07-30 2022-06-09 一种多维度航班放行效能评估方法 WO2023005442A1 (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US17/822,812 US11710416B2 (en) 2021-07-30 2022-08-29 Multi-dimensional flight release efficiency evaluation method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202110869090.8A CN113554898B (zh) 2021-07-30 2021-07-30 一种多维度航班放行效能评估方法
CN202110869090.8 2021-07-30

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US17/822,812 Continuation US11710416B2 (en) 2021-07-30 2022-08-29 Multi-dimensional flight release efficiency evaluation method

Publications (1)

Publication Number Publication Date
WO2023005442A1 true WO2023005442A1 (zh) 2023-02-02

Family

ID=78104958

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/097769 WO2023005442A1 (zh) 2021-07-30 2022-06-09 一种多维度航班放行效能评估方法

Country Status (2)

Country Link
CN (1) CN113554898B (zh)
WO (1) WO2023005442A1 (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113554898B (zh) * 2021-07-30 2022-07-19 中国电子科技集团公司第二十八研究所 一种多维度航班放行效能评估方法
CN114360295B (zh) * 2021-11-08 2023-06-09 民航数据通信有限责任公司 一种空中交通容流平衡措施控制方法及装置
CN114627684B (zh) * 2022-01-24 2023-08-18 中国民用航空中南地区空中交通管理局海南分局 受多重流量管理策略影响下航班起飞时隙分配计算方法
CN114611975B (zh) * 2022-03-23 2022-10-04 中航材导航技术（北京）有限公司 一种自动评价航班放行质量的方法
CN115662197B (zh) * 2022-12-28 2023-03-17 中国电子科技集团公司第二十八研究所 基于信息差异赋权的空域灵活使用效能评估指标计算方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103093650A (zh) * 2013-01-08 2013-05-08 成都民航空管科技发展有限公司 管制运行质量诊断系统及方法
CN104933530A (zh) * 2015-06-29 2015-09-23 中国民用航空厦门空中交通管理站 空管实时管制效能评价系统
CN110660274A (zh) * 2019-09-30 2020-01-07 中国电子科技集团公司第二十八研究所 一种机场容需平衡预测方法
CN111462533A (zh) * 2020-03-13 2020-07-28 海南太美航空股份有限公司 航班起飞时刻管理方法及系统
CN113554898A (zh) * 2021-07-30 2021-10-26 中国电子科技集团公司第二十八研究所 一种多维度航班放行效能评估方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103473469B (zh) * 2013-09-25 2016-06-22 南京航空航天大学 一种基于客观指标的扇区交通态势多层次模糊评价方法
CN109741638B (zh) * 2019-01-23 2020-06-12 中国电子科技集团公司第二十八研究所 一种进离场管理系统一体化协同运行方法
CN113034980B (zh) * 2021-02-26 2022-08-19 中国电子科技集团公司第二十八研究所 一种基于时刻表的航班运行效能预先评估方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103093650A (zh) * 2013-01-08 2013-05-08 成都民航空管科技发展有限公司 管制运行质量诊断系统及方法
CN104933530A (zh) * 2015-06-29 2015-09-23 中国民用航空厦门空中交通管理站 空管实时管制效能评价系统
CN110660274A (zh) * 2019-09-30 2020-01-07 中国电子科技集团公司第二十八研究所 一种机场容需平衡预测方法
CN111462533A (zh) * 2020-03-13 2020-07-28 海南太美航空股份有限公司 航班起飞时刻管理方法及系统
CN113554898A (zh) * 2021-07-30 2021-10-26 中国电子科技集团公司第二十八研究所 一种多维度航班放行效能评估方法

Also Published As

Publication number Publication date
CN113554898A (zh) 2021-10-26
CN113554898B (zh) 2022-07-19

Similar Documents

Publication Publication Date Title
WO2023005442A1 (zh) 一种多维度航班放行效能评估方法
CN104484993B (zh) 用于交通小区划分的手机信令信息的处理方法
CN103530704B (zh) 一种终端空域空中交通动态容量预测系统及其方法
CN108152866B (zh) 基于天气影响航班量的航空气象预报质量评估方法
CN104021671A (zh) 一种svm与模糊判断相结合的道路实时路况的判断方法
WO2020083399A1 (zh) 基于交通流数据的协调干线线路规划方法及配置系统
CN111009155B (zh) 基于空域结构和航班流的空中交通流复杂性量化分析方法
CN205959030U (zh) 城乡规划监测管理系统
CN110060513A (zh) 基于历史轨迹数据的空中交通管制员工作负荷评估方法
CN106339358B (zh) 基于多元回归分析的航空器场面滑行时间预测方法
CN105225539B (zh) 基于主成分分析的扇区运行性能综合指数的方法及系统
CN110288177A (zh) 一种采用历史运行数据包络线评估机场公布容量的方法
CN104503406B (zh) 火电厂优化运行评价系统及其方法
CN110807917A (zh) 一种基于云计算的高速公路智能安防系统
CN106651088A (zh) 一种基于时态gis的飞行品质监控方法
CN110660274A (zh) 一种机场容需平衡预测方法
CN108053153A (zh) 一种生态旅游环境承载力预警系统
CN113643539A (zh) 一种基于大数据的智慧交通管控系统
CN107134136A (zh) 一种城市之间道路交通拥堵程度排序的方法和系统
CN111984878B (zh) 一种对流天气下机场终端区各进离港点流量控制方法
US20230035944A1 (en) Multi-dimensional flight release efficiency evaluation method
CN116384807A (zh) 基于模糊综合评价模型的航路网络拥挤态势评估方法
CN105096661A (zh) 基于垂直剖面视图的空中交通流管理系统与方法
CN114299763A (zh) 一种流量管理措施全生命周期建立方法
CN107038513A (zh) 一种指标的综合展示方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22848052

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE