CN117351786A - Flight integrated scheduling method under multi-element constraint - Google Patents

Abstract

The invention provides a flight integrated scheduling method under multiple constraints, which comprises the following steps: flight plan basic information, traffic management measure information, etc.; initializing space-time resources of a high-level flow management measure according to different characteristics of the flow management measure; constructing a multidimensional flight priority evaluation model, calculating the flight priority according to flight plan information, generating a flight scheduling queue, and calculating flights in sequence in the subsequent scheduling process; in the process of flight scheduling calculation, calculating target passing time, target departure time, target passing height and total delay of flights at flow management measure limiting points and sub delays caused by various flow management measures; and updating the available time-space resources and the unavailable time-space resources of the flow management measures according to the flight scheduling calculation result.

Description

Flight integrated scheduling method under multi-element constraint

Technical Field

The invention belongs to the field of air traffic management, and particularly relates to an integrated scheduling method for flights under multiple constraints.

Background

The decision of flow management is mainly based on the management measures of space restriction point release interval management, total amount control, height restriction and the like, namely the time and space restriction, and then the flight distribution is adjusted by carrying out release scheduling on flights, so that the traffic jam is relieved. Because the actual running environment is complex, the flow management target is often required to be realized through combined measures, and how to solve the scheduling management problem of flights under multiple running constraint through a scientific and efficient method becomes a hot spot problem of the current running management. In the aspect of flight scheduling management, the existing related technical methods at home and abroad are mainly focused on solving by adopting a probability search algorithm, so that the calculation time is long, the interpretation is poor, and the actual operation requirements are difficult to meet; the existing application system mainly adopts an iteration test method to solve, has poor controllability of an iteration process and has certain randomness. At present, a stable, efficient and highly interpretable flight scheduling method under multiple constraints is not available.

Disclosure of Invention

The invention aims to: the invention aims to solve the technical problem of providing a scientific and effective integrated flight scheduling method aiming at the defects of the flight scheduling method under the multi-element constraint. The method comprises the steps of acquiring civil aviation production operation data through a comprehensive service information platform, analyzing and preprocessing flight plan information and flow management measure information, describing information in a structured form, and initializing flow management measure high-level space-time resources; constructing a multidimensional flight priority evaluation model, generating a flight scheduling queue order, and scheduling and calculating flights according to the order; carrying out space-time resource analysis on the flow management measures related to the flight, and calculating to obtain the target passing point time, the target passing point height, the target departure time, the total delay of the flight and the sub delay of the flow management measures of the flight; and finally updating the space-time resources of the flow management measure high-level according to the flight calculation result. The invention can provide fine management capability for flight operation, and can quickly form a flight scheduling scheme in actual operation, relieve air traffic jam and ensure orderly operation.

The technical scheme is as follows: in order to achieve the aim of the invention, the invention discloses a flight integrated scheduling method under the multi-element constraint, which comprises the following steps:

Step 1, analyzing and processing comprehensive operation information, wherein the comprehensive operation information comprises: flight plan information and traffic management measure information;

step 2, initializing space-time resources of a high-level flow management measure according to the flow management measure information;

step 3, constructing a multidimensional flight priority evaluation model, calculating the flight priority according to the comprehensive operation information, and generating a flight scheduling queue;

step 4, calculating target passing point time, target passing point height, target departure time and total delay of the flight at the flow management measure limiting point;

step 5, decomposing the total delay of the flights according to the space-time resources of the high-level layers of the flow management measures to obtain the sub delays of the flights caused by the flow management measures;

step 6, updating the space-time resources of the flow management measure height layer according to the calculated target passing point time of the flow management measure limiting point;

step 1 comprises the following steps:

step 1-1: acquiring and resolving empty pipe height layer sets from comprehensive operation platformHeightSet={Height ₁ ，…，Height _i }, whereinHeight _i Representing the height value of the ith height layer used in the normal operation of the empty pipe,HtNumfor all the number of layers of the height,，Height _i >Height _i-1 ；

acquiring the total set of all the time of the operation day as TimeSet=[00:00，23:59]；

Step 1-2: acquiring and resolving flight plan information from an integrated operation platformFlightInfoFlight plan informationFlightInfoComprising the following steps: flight call signCallsignUnique number of flightPlanIDStatus of flightStatusExemption attribute of flight associated traffic management measureExemptScheduled gear withdrawal timeSOBTPredicted take-off timeETOTAirport for taking offDepAirportAirport for landingDesAirportFlow management measures associated with flights limit the number of pointsResPtNumAnd flight restriction informationPtResourceSet；

The flight status includes: not being executedNoExecutedGround executionGroundHas taken offDEPHas fallen downARR；

Exemption attribute for traffic management measuresExempt=1, indicating that flights are exempted from all flow management measures; if the flight exempts from one flow management measure, the flight is indicated to be free from considering the limiting condition of the flow management measure in the calculation process;

PtResourceSet={PtResource ₁ ，…，PtResource _i }，PtResource _i the ith restriction point information indicating the arrival of the flight includes: limiting point nameRestrictPt _i Number of traffic management measures associated with limiting pointMeasureNum _i Time of flight arrival at limit pointRestrictPTO _i Flight arrival limit point altitude layer sequence numberRestrictPHt _i Flight arrival limit point altitude rangeRestrictHtRange _i Actual passing point indicationPassMark _i Set of limit point traffic management measuresPtMsResourceSet _i ，RestrictPTO _i >RestrictPTO _i-1 ；

If it isPassMark _i =1, indicating that the flight has flown past the ith limit point,RestrictPTO _i the actual passing point time of the ith limiting point; if it is PassMark _i =0, indicating that the flight did not fly past the ith limit point,RestrictPTO _i the predicted passing time for the ith limit point;

PtMsResourceSet _i ={PtMsResource _i,1 ，…，PtMsResource _i,j }，PtMsResource _i,j the j-th associated traffic management measure information indicating that the flight is at the i-th restriction point, comprising: associated traffic management measure unique numberPtMeasureID _i,j Associated traffic management measure classPtRank _i,j Exemption indication for associated flow management measuresExpMark _i,j ；

If it isExpMark _i,j =1, indicating that the flight is exempted from the jth associated traffic management measure at the ith limit point; if it isExpMark _i,j =0, indicating that the flight is restricted for the jth associated traffic management measure at the ith restriction point;

RestrictHtRange _i ={RestrictHt _i,1 ，…，RestrictHt _i,m }，RestrictHt _i,m indicating that the flight is at the highest level possible at the ith limit pointHeightSetThe number of which is a sequence number,RestrictHt _i,m >RestrictHt _i,m ；

step 1-3: the following information is extended:

at the position ofFlightInfoIn extension information: flight target departure timeCTOT；

At the position ofPtResource _i In extension information: target transit time for flight to reach limit pointRestrictCTO _i Limiting point target overpass heightRestrictCHt _i Limiting point availabilitySource time periodPtUsPeriod _i Limiting point high-level available resource setPtHtUsPeriodSet _i ={PtHtUsPeriod _i,1 ，…，PtHtUsPeriod _i,j }，PtHtUsPeriod _i,j Representing the availability period of the jth altitude layer of the flight at the ith limit point;

at the position ofPtMsResource _i,j The extension definition in (a): flight delay caused by traffic management measuresPartDelay _i,j High-level availability period set for traffic management measuresPtMsHtUsPeriodSet _i,j ={PtMsHtUsPeriod _i,j,1 ，…，PtMsHtUsPeriod _i,j,k }，PtMsHtUsPeriod _i,j,k Indicating the period of availability of the flight at the kth altitude layer of the jth traffic management measure at the ith limit point, ；

Step 1-4: creating a set of traffic management measuresMeasurInfoSet={MeasurInfo ₁ ，…，MeasurInfo _i Acquiring and analyzing flow management measure information from the comprehensive operation platformMeasureInfo _i Traffic management measure informationMeasureInfo _i Comprising the following steps: unique number of flow management measureMeasureID _i Start timeBgnTime _i End timeEndTime _i Duration of effectivenessDuration _i Grade (grade)Rank _i Total amount limitationCountRes _i Time granularity of total limitCountReTimesSize _i Spacing limitIntervalRes _i Limiting the heightHeightResSet _i Height limiting typeHeightType _i ；

Duration _i =EndTime _i -BgnTime _i ，Duration _i /CountReTimesSize _i =PeriodNum _i ，Duration _i %CountReTimesSize _i =0，PeriodNum _i Indicating that the ith traffic management action is in effect period of timeCountReTimesSize _i The total amount obtained by dividing limits the number of the monitoring sub-periods;

Rank _i comprising three levels 1,2, 3;

HeightResSet _i ={ResHtIndex _i,1 ，…，ResHtIndex _i,j }，ResHtIndex _i,j limited height layer representing the ith traffic management measureHeightSetThe number of which is a sequence number,ResHtIndex _{i, j} >ResHtIndex _i,j-1 ；

if it isHeightType _i =0Indicating that flights pass through flow management measure limiting points at any altitude, the flights are limited by total amount and interval;

if it isHeightType _i =1If the altitude layer sequence number indicating that the flight passes the flow management measure limiting point is contained inHeightResSet _i If the flight is not limited by the total amount, the interval is limited by the total amount;

if it isHeightType _i =2Indicating that flights can only be included in sequence numbersHeightResSet _i The internal altitude layers pass through the flow management measure limiting points and are limited by total amount and interval, and flights cannot pass through other altitude layers in the measure period [ BgnTime _i ，EndTime _i ]The flow management measure limiting point is passed through;

step 1-5: the following information is extended:

at the position ofMeasureInfo _i The extension definition in (a): traffic management measure height layerSpace-time resourceHtResourceSet _i ={HtResource _i,1 ，…，HtResource _i,j }，HtResource _i,j Resources representing the ith traffic management measure at the jth level include: resource management indicatorHtResourceMgt _i,j Period of high-level available resourcesHtUsPeriod _i,j Period of high-level unavailable resourceHtFdPeriod _i,j Aggregate of highly-layer aggregate resource monitoring periodsHtMonitor _i,j Resource counting in a high-level aggregate resource monitoring periodHtCount _i,j ，；

If it isHtResourceMgt _i,j =0, meaning that the ith traffic management measure is not limited by total amount, interval at the jth high level group, thenThe method comprises the steps of carrying out a first treatment on the surface of the If it isHtResourceMgt _i,j =1, indicating that the ith traffic management measure is limited by the total amount, spacing, at the jth level group;

HtMonitor _i,j ={HtMonPeriod _i,j,1 ，…，HtMonPeriod _i,j,m }，HtMonPeriod _i,j,m an mth monitoring sub-period representing an jth level of an ith traffic management measure,HtMonPeriod _{i, j,m} =[BgnTime _i +(m-1)* CountReTimesSize _i ，BgnTime _i +m*CountReTimesSize _i )]，；

HtCount _i,j ={HtMonCount _i,j,1 ，…，HtMonCount _i,j,m }，HtMonCount _i,j,m the number of remaining flight resources for the mth monitoring sub-period of the jth elevation layer representing the ith traffic management measure.

Step 2 comprises the following steps:

step 2-1: traversing traffic management measure informationMeasureInfo _i ；

Step 2-2; initializing space-time resources of a high-level layer of traffic management measures:

；

if it isHeightType _i =0，

；

；

If it isHeightType _i =1，

；

；

If it isHeightType _i =2，

；

；

，

HtFdPeriod _i,j =[BgnTime _i ,EndTime _i ],HtUsPeriod _i,j =TimeSet-HtFdPeriod _i,j 。

Step 3 comprises the following steps:

step 3-1: constructing a multidimensional flight priority assessment model comprising 13 packets:

PriorityGroup ₁ : the flight group has fallen, grouping conditions:Status=ARR；

PriorityGroup ₂ : group of flights taken off, grouping conditions: does not satisfyPriorityGroup ₁ Grouping conditions, andStatus=DEP；

PriorityGroup ₃ : airlines have been executed on the ground, grouping conditions: does not satisfyPriorityGroup ₁ ~PriorityGroup ₂ Grouping conditions, andStatus=Ground；

PriorityGroup ₄ : full-exemption airlines, grouping conditions: does not satisfyPriorityGroup ₁ ~PriorityGroup ₃ The grouping condition is that,Exempt=1；

PriorityGroup ₅ : and (3) associating class 1 traffic management measures with the full-exemption airlines, and grouping conditions: does not satisfyPriorityGroup ₁ ~PriorityGroup ₄ The grouping condition is that,PtRank _i,j maximum value is 1, and whenPtRank _i,j When=1, satisfyExpMark _i,j =1；

PriorityGroup ₆ : and associating class 1 traffic management measures with semi-exemption airlines, and grouping conditions: does not satisfyPriorityGroup ₁ ~PriorityGroup ₅ The grouping condition is that,PtRank _i,j maximum value is 1, and whenPtRank _i,j When=1, there isExpMark _i,j =0；

PriorityGroup ₇ : class 1-associated traffic management measure full-limited flight group, grouping conditions: does not satisfyPriorityGroup ₁ ~PriorityGroup ₆ The grouping condition is that,PtRank _i,j maximum value is 1, and whenPtRank _i,j When=1, satisfyExpMark _i,j =0；

PriorityGroup ₈ : and (3) associating class 2 traffic management measures with the full-exemption airlines, and grouping conditions: does not satisfyPriorityGroup ₁ ~PriorityGroup ₇ The grouping condition is that,PtRank _i,j maximum value is 2, and whenPtRank _i,j When=2, satisfyExpMark _i,j =1；

PriorityGroup ₉ : and associating class 2 traffic management measures with semi-exemption airlines, and grouping conditions: does not satisfyPriorityGroup ₁ ~PriorityGroup ₈ The grouping condition is that,PtRank _i,j maximum value is 2, and whenPtRank _i,j When=2, there isExpMark _i,j =0；

PriorityGroup ₁₀ : traffic management measures of association class 2 are all restricted airlines, grouping conditions: does not satisfyPriorityGroup ₁ ~PriorityGroup ₉ The grouping condition is that, PtRank _i,j Maximum value is 2, and whenPtRank _i,j When=2, satisfyExpMark _i,j =0；

PriorityGroup ₁₁ : and (3) associating class 3 traffic management measures with the full-exemption airlines, and grouping conditions: does not satisfyPriorityGroup ₁ ~PriorityGroup ₁₀ The grouping condition is that,PtRank _i,j maximum value is 3, and whenPtRank _i,j When=3, satisfyExpMark _i,j =1；

PriorityGroup ₁₂ : and associating class 3 traffic management measures with semi-exemption airlines, and grouping conditions: does not satisfyPriorityGroup ₁ ~PriorityGroup ₁₁ The grouping condition is that,PtRank _i,j maximum value is 3, and whenPtRank _i,j When=3, there isExpMark _i,j =0；

PriorityGroup ₁₃ : traffic management measures of association level 3 are all restricted airlines, grouping conditions: does not satisfyPriorityGroup ₁ ~PriorityGroup ₁₂ The grouping condition is that,PtRank _i,j maximum value is 3, and whenPtRank _i,j When=3, satisfyExpMark _i,j =0；

Step 3-2: traversing each flight;

step 3-3: judging whether the flight meets the condition of priority grouping or not, if soPriorityGroup _i For the packet condition of (a), then the flight priority packetsPriorityIndex=i；

Step 3-4: after traversing the flights, judging every two flightsFlight _i AndFlight _j is of the flight priority of (2)FlightPriority _i AndFlightPriority _j is a size relationship of (a):

if it isPriorityIndex _i >PriorityIndex _j ThenFlightPriority _i >FlightPriority _j ；

If it isPriorityIndex _i =PriorityIndex _j And (2) andSOBT _i <SOBT _j thenFlightPriority _i > FlightPriority _j ；

If it isPriorityIndex _i =PriorityIndex _j And (2) andSOBT _i =SOBT _j and (2) andPlanID _i >PlanID _j thenFlightPriority _i >FlightPriority _j ；

The flight unique numbers will not be the same, i.e. ifThen->；

Step 3-5: ordering the flights according to the ascending order of the priority of the flights to obtain a flight scheduling queueFlightArray，FlightArray={Flight ₁ ，…，Flight _i }, whereinFlightPriority _i > FlightPriority _i-1 。

Step 4 comprises the steps of:

step 4-1: at the position ofFlightArraySequentially calculates for each flight, and records the currently calculated flight as Fight；

Step 4-2: traversing restricted information for flightsPtResourceSetCalculating target passing time of flight at limiting pointRestrictCTO _i ；

Step 4-3: calculating a target passing point altitude of a flight at a limiting pointRestrictCHt _i ；

Step 4-4: calculating total delay of flights due to flow management measuresMeasureTotalDelayTakeoff time of flight targetCTOT。

Step 4-2 includes: if it isPassMark _i =1, thenRestrictCTO _i =RelatedPTO _i The target passing time of the limiting point indicating that the flight has flown is equal to the actual passing time;

if it isPassMark _i =0, then:

when (when)If at the timeExpMark _i,j =0，The available resource period of the kth altitude layer is equal to the value of the (j) th associated traffic management measure representing the (j) th limit point of the (i) th limit pointMeasurInfoSetM-th stream of (b)The volume management measure is at the time of the available resources of the kth level, wherein the kth level must be the available level of the flight at the ith limit point;

if it isExpMark _i,j =1，When the j-th associated flow management measure of the flight at the i-th limiting point is exempted, the available resource period of the k-th height layer is a full period, wherein the k-th height layer is the available height layer of the flight at the i-th limiting point;

the available resource period of the unavailable altitude layer of the j-th associated traffic management measure of the flight at the i-th limiting point is represented as empty;

representing that the available resource period of the kth high-level of the flight at the ith limit point is equal to the intersection of the available resource periods of the high-level of the respective traffic management measures of the flight at the ith limit point;

Representing that the available resource period of the flight at the ith limit point is equal to the union of the available resource periods of the flights at each altitude layer of the ith limit point;

calculating a limit point reference expected transit time；

Calculating the difference between the predicted passing time of each limiting point and the reference predicted passing timePTODvalue _i ，PTODvalue _i =RestrictPTO _i -BasePTO；

Performing translation processing on the available resource time period of the limiting point to obtain the updated available resource time period of the limiting point；

Taking intersection of available resource time periods updated by all limiting points to obtain available resource time periods of flights；

Obtaining a reference target passing point time by taking the minimum value of the reference predicted passing point time of the limiting point in the available resource time period of the flightBaseCTO，；

Calculating target passing time of flight at each limiting pointRestrictCTO _i =RestrictPTO _i +(BaseCTO- BasePTO)。

Step 4-3 includes: if it isPassMark _i =1, thenRestrictCHt _i =RestrictPHt _i The target passing point height of the limiting point which indicates that the flight has flown is equal to the actual passing point height;

if it isPassMark _i =0, thenRestrictCHt _i =HSo thatThe target passing point height layer of the flight at the ith limiting point is the H-th height layer, the H-th height layer is closest to the predicted passing point height layer of the limiting point, and the available resource period of the H-th height layer comprises the target passing point time of the limiting point.

Step 4-4 includes: if presentPassMark _i =1, thenIndicating that the taken-off flight does not need to calculate the target take-off time;

If all ofPassMark _i =1, thenMeasureTotalDelay=0Indicating that the flight has passed all restriction points, the traffic management measures do not cause delays;

if presentPassMark _i =0, thenMeasureTotalDelay=BaseCTO-BasePTO；

If all ofPassMark _i =0, then:

MeasureTotalDelay=BaseCTO- BasePTO；

CTOT=ETOT+ MeasureTotalDelay。

step 5 comprises the steps of:

step 5-1: calculating a child delay of the flight caused by the j-th associated traffic management measure of the i-th restriction pointPartDelay _i,j ；

Step 5-2: if it isMeasureTotalDelay=0, thenPartDelay _i,j =0, indicating that the total delay of the flight is 0, and the child delay is 0;

step 5-3: if it isThen the following determination is made:

if it isExpMark _i,j =0, thenPartDelay _i,j =0, indicating that the associated exemption traffic management measure does not cause delays to flights;

if it isThen:

memory collection；

Memory collection，FdPeriod _i,j A period of unavailable resources representing a jth associated traffic management measure for the flight at the ith limit point;

recording device，PartFdLen _i,j Inedibility of the target transit point elevation layer of the jth associated traffic management measure representing a flight at the ith limit pointBy resource periodsTotal _i The number of time elements of the intersection of (a);

calculated to obtain。

Step 6 comprises the steps of:

step 6-1: after calculation for each flight, the traffic management measure high-level space-time resources are updatedHtResourceSet _i ；

Step 6-2: spacing restrictions based on traffic management measuresIntervalRes _i Updating:

if it isThe flow management measures are not limited at intervals, so that space-time resources of the flow management measures do not need to be updated;

If it isAnd (2) andHtResourceMgt _i,j =0, which indicates that the traffic management measure is not limited at the jth level, so that the level space-time resource does not need to be updated;

if it isAnd (2) andHtResourceMgt _i,j when the traffic management measure satisfies =1MeasureID _i =PtMeasureID _m,n When the method is performed, the following determination is made:

if it isExpMark _m,n =0, then generate the interval unavailable period of the ith traffic management measure at the jth levelIntervalFdPeriod _i,j =[RestrictCTO _m -IntervalRes _i ，RestrictCTO _m +IntervalRes _i ]；

If it isExpMark _m,n =1, obtainHtUsPeriod _i,j Middle greater than or equal toRestrictCTO _m Is set to be a minimum value of (c),recorded as the earliest time of availabilityThe method comprises the steps of carrying out a first treatment on the surface of the Generating a temporary unavailable period of traffic management measure intervalIntervalFdPeriod _i,j =[EarliestUsTime _i,j -IntervalRes _i ，EarliestUsTime _i,j +IntervalRes _i ]；

Updating the available period of the ith traffic management action at the jth level；

Updating the ith traffic management measure at the unavailable period of the jth level；

Step 6-3: total amount restriction based on traffic management measuresCountRes _i Updating:

if it isThe flow management measures are not limited in total amount, so that space-time resources of the flow management measures do not need to be updated;

if it isAnd (2) andHtResourceMgt _i,j =0, which indicates that the traffic management measure is not limited at the jth level, so that the level space-time resource does not need to be updated;

if it isAnd (2) andHtResourceMgt _i,j when the traffic management measure satisfies =1MeasureID _i =PtMeasureID _m,n When the method is performed, the following determination is made:

if it isExpMark _m,n =0, andthen updateHtMonCount _i,j,k =HtMonCount _i,j,k -1；

If it isExpMark _m,n =1, obtainHtUsPeriod _i,j Middle greater than or equal toRestrictCTO _m Is recorded as the earliest time of availability The method comprises the steps of carrying out a first treatment on the surface of the If it isThen updateHtMonCount _i,j,k =HtMonCount _i,j,k -1；

If after updatingGenerating a total unavailable period of the ith traffic management measure at the jth levelCountFdPeriod _i,j =HtMonPeriod _i,j,k ；

Updating the available period of the ith traffic management measure at the jth level；

Updating the ith traffic management measure at the unavailable period of the jth level。

Supplementary explanation:

setting presence setsA，BFormulas, representing the setAMedium element removal setBA new set is obtained by the elements in (a);card(A) Representing a collectionAThe number of elements in the matrix.Min(A) Representing a collectionAFor a time set, represent the setAIs the earliest time in time.

The time periods and time intervals mentioned in the present invention all represent time sets with accuracy of minutes, for example: [9:00,9:05] = { 9:00,9:01,9:02,9:03,9:04,9:05}.

The beneficial effects are that: the invention has the remarkable advantages that:

1. the flight scheduling method based on space-time resource management is provided for the air traffic control operation;

2. the multi-dimensional flight priority evaluation method is provided for the air traffic control operation;

3. the flight delay decomposition method facing the flow management measures is provided for the air traffic control operation;

4. provides powerful support for the research of the flight scheduling management technology and the system development.

Drawings

The foregoing and/or other advantages of the invention will become more apparent from the following detailed description of the invention when taken in conjunction with the accompanying drawings and detailed description.

Fig. 1 is a schematic diagram of an integrated service information processing flow.

FIG. 2 is a schematic diagram of a high-level space-time resource flow for initializing traffic management measures.

FIG. 3 is a schematic illustration of a process for calculating a flight restriction point target transit time.

FIG. 4 is a schematic diagram of a high-level space-time resource flow for updating traffic management measures.

Detailed Description

In one embodiment of the invention, the date of the related time information is 2021/08/17, firstly, a comprehensive business information platform is constructed, simulation data are generated, and flight plan information and flow management measure information are analyzed; selecting flow management measures with different restriction types as examples, and initializing the space-time resource of each flow management measure in the high layer; selecting a flight plan related to the flow management measures, and evaluating priority according to flight information to obtain a flight scheduling queue; in the flight scheduling queue, sequentially generating available resource time periods with various scales for flights, and finally calculating to obtain target passing point time, target passing point height, target departure time, total delay of the flights caused by flow management measures and sub-delay of the flights caused by the flow management measures of the flights at all limiting points; and after the calculation of the single flight is finished, updating and maintaining the space-time resources of the associated flow management measure high-level.

The following details the processing flow steps of a flight integrated scheduling method under the multicomponent constraint based on the method of the invention with reference to fig. 1 to fig. 4:

1. acquiring comprehensive service information

1-1, resolving air traffic control height layer information

As shown in fig. 1, the reading of the configuration file needs to be specifically described as follows: for facilitating the example analysis, a part of the height layer range is selected as the available height layer set for the empty pipe operation to obtainHeightSet={4200，4800，5400，7500，7800，8100，8400，8900，9200，9500}。

1-2, resolving flight plan information

As shown in fig. 1, acquiring and resolving flight plan information includes: flight call signCallsignUnique number of flightPlanIDStatus of flightStatusExemption attribute of flight associated traffic management measureExemptScheduled gear withdrawal timeSOBTPredicted take-off timeETOTAirport for taking offDepAirportAirport for landingDesAirportFlow management measures associated with flights limit the number of pointsResPtNumFlight restriction informationPtResourceSet. Four flights are selected for illustration in this example, and the flight schedule information is as follows:

Callsignthe method comprises the following steps of:DKH1287、CSN5566、CSC8867、CHH8654；

PlanIDthe method comprises the following steps of:NO128755、NO345785、NO357877、NO786543；

Statusthe method comprises the following steps of:DEP、NoExecuted、NoExecuted、NoExecuted；

Exemptthe method comprises the following steps of: 0. 0, 1;

SOBTthe method comprises the following steps of: 09:00, 09:10, 11:10, 10:58;

ETOTthe method comprises the following steps of: 09:30, 09:40, 11:40, 11:28;

DepAirportthe method comprises the following steps of:ZGGG、ZGGG、ZGGG、ZGGG；

DesAirportthe method comprises the following steps of:ZSSS、ZSSS、ZSSS、ZSSS；

ResPtNumthe method comprises the following steps of: 3. 3, 1;

PtResource ₁ the following are provided:

RestrictPt ₁ The method comprises the following steps of:NXD、NXD、PLT、PLT；

MeasureNum ₁ the method comprises the following steps of: 1. 1, 1;

RestrictPTO ₁ the method comprises the following steps of: 09:45, 09:55, 12:50, 12:30;

RestrictPHt ₁ the method comprises the following steps of: 2. 2, 6;

RestrictHtRange ₁ the method comprises the following steps of: {1,2,3}, {1,2}, {4,5,6,7,8,9};

PassMark ₁ the method comprises the following steps of: 1. 0, 0;

PtMsResource _1,1 the following are provided:

PtMeasureID _1,1 the method comprises the following steps of:R1、R1、R3、R3；

PtRank _1,1 the method comprises the following steps of: 1. 1, 3;

ExpMark _1,1 the method comprises the following steps of: 0. 0, 0;

PtResource ₂ the following are provided:

RestrictPt ₂ the method comprises the following steps of:PLT、PLT、without any means for、The method is free;

MeasureNum ₂ the method comprises the following steps of: 2. 2, have nothingThe method is free;

RestrictPTO ₂ the method comprises the following steps of: 10:40, 11:00, none;

RestrictPHt ₂ the method comprises the following steps of: 5. 5, none;

RestrictHtRange ₂ the method comprises the following steps of: {4,5,6,7,8}, none;

PassMark ₂ the method comprises the following steps of: 0. 0, none;

PtMsResource _2,1 the following are provided:

PtMeasureID _2,1 the method comprises the following steps of:R2、R2、without any means for、The method is free;

PtRank _2,1 the method comprises the following steps of: 2. 2, none;

ExpMark _2,1 the method comprises the following steps of: 1. 0, none;

PtMsResource _2,2 the following are provided:

PtMeasureID _2,2 the method comprises the following steps of:R3、R3、none, none;

PtRank _2,2 the method comprises the following steps of: 3. 3, none;

ExpMark _2,2 the method comprises the following steps of: 0. 0, none.

1-2, analyzing flow management measure information

As shown in fig. 1, acquiring and analyzing flow management measure information includes: unique number of flow management measureMeasureID、Measure start timeBgnTimeTime of measure endEndTimeDuration of measure validationDurationGrade (grade)RankTotal amount limitationCountRes(unit: installment), total amount limiting timeParticle sizeCountReTimesSize(units: minutes), interval limitation IntervalRes(unit: minutes), limit heightHeightResSetHeight limiting typeHeightType. A flow management measure set is generated, and three flow management measures are selected for illustration in this example, flow management measure information is as follows:

MeasureIDthe method comprises the following steps of:R ₁ ，R ₂ ，R ₃ ；

BgnTimethe method comprises the following steps of: 09:00, 10:30, 10:00;

EndTimethe method comprises the following steps of: 11:00, 11:30, 13:00;

Durationthe units of minutes are respectively as follows: 120, 60, 180;

Rankthe method comprises the following steps of: 1,2, 3;

CountResthe method comprises the following steps of: none, 1;

CountReTimesSizethe method comprises the following steps of: none, 30;

IntervalResthe method comprises the following steps of: 15, 20, none;

HeightResSetthe method comprises the following steps of: {1,2,3,4,5,6,7,8, 9,10}, {5,6,7,8, 9,10}, 4;

HeightTypethe method comprises the following steps of: 0, 1, 2;

MeasurInfoSet

the method comprises the following steps: { R ₁ ，R ₂ ，R ₃ }。

2. Initializing space-time resources for traffic management measures

2-1, initializing to generate an available resource period and an unavailable resource period

As shown in fig. 2, the available resource period, the unavailable resource period, are initialized on each empty pipe height layer for each traffic management measure according to the traffic management measure information. The results of this example process are shown below:

R ₁ (i=1)，, HtUsPeriod _i,j for [00:00,23:59 ]]，HtFdPeriod _i,j Is->；

R ₂ (i=2)，, HtUsPeriod _i,j For [00:00,23:59 ]] ，HtFdPeriod _i,j Is->；

,HtResourceMgt _i,j =1；

R ₃ (i=3), j=4, HtUsPeriod _i,j For [00:00,23:59 ]], HtFdPeriod _i,j Is that；

, HtUsPeriod _i,j Is->, HtFdPeriod _i,j Is->；

j=4,HtResourceMgt _i,j =1。

2-2, initializing the height layer total amount monitoring information

As shown in fig. 2, the high-level total resource monitoring period is initialized for each traffic management measure containing a total limit, and the high-level total resource monitoring period resource count information is initialized according to the traffic management measure information. The results of this example process are shown below:

R ₁ (i=1)、R ₂ (i=2)，HtMonitor _i,j Is that；HtCount _i,j Is->；

R3(i=3)，HtMonitor _i,j Is { [10:00, 10:30), [10:30, 11:00), [11:00, 11:30), [11:30, 12:00), [12:00, 12:30), [12:30, 13:00) };

HtCount _i,j is {1, 1}.

3. Generating a flight schedule queue

3-1, judging the priority grouping of the flight

And judging which priority group the flight belongs to according to the flight priority evaluation model. The results of this example process are shown below:

DKH1287the flight status isDEPThe related measures areF1, F2, F3The method comprises the steps of carrying out a first treatment on the surface of the Whether the gap is avoided is respectively as follows: no, highest ranking, 1,2,3, priority grouping ofPriorityGroup ₁ ；

CSN5566The flight status isNoExecuted,The association measure is thatF1, F2, F3The method comprises the steps of carrying out a first treatment on the surface of the Whether the gap is avoided is respectively as follows: no, highest ranking 1,2,3, priority grouping ofPriorityGroup ₇ ；

CSC8867The flight status isNoExecuted,The association measure is thatF3, Whether to exempt from: if not, the highest level is 3, the priority group isPriorityGroup ₁₃ ；

CHH8654The flight status isNoExecuted, The association measure is thatF3, Whether to exempt from: if not, the highest level is 3, the priority group isPriorityGroup ₁₃ 。

3-2, judging the priority order of flights

And judging the priority order of the flights according to the flight schedule information, and generating a flight scheduling queue. The results of this example process are shown below:

DKH1287,PlanID：NO128755,SOBT:09:00, priority grouping:PriorityGroup ₁ priority order: 1, a step of;

CSN5566,PlanID：NO345785,SOBT:09:10, priority grouping:PriorityGroup ₇ priority order: 2;

CSC8867,PlanID：NO357877,SOBT:11:10, priority packets:PriorityGroup ₁₃ priority order: 4, a step of;

CHH8654，PlanID：NO786543,SOBT:10:58, priority packets:PriorityGroup ₁₃ priority order: 3, a step of;

PlanID,NO128755,NO345785,NO357877,NO786543；

FlightArrayis { as }DKH1287,CSN5566,CSC8867,CHH8654}。

4. Calculating flight scheduling decision information

4-1, calculating the target passing time of the limiting point

As shown in fig. 3, each traffic management measure of generating a flight at a limit point is at an available resource period, an altitude layer available resource period, a limit point available resource period, an updated limit point available resource period, and a limit point target transit time of each altitude layer. FlightDKH1287The analysis results are shown below:

the limiting point is NXD, i=1,PassMark _i =1, height layer (k): the method is free; flow management measures: r is R ₁ ，j=1；ExpMarki, jIs 0; corresponding toHtUsPeriod _m,k : the method is free;PtMsHtUsPeriod _i,j,k : the method is free;PtHtUsPeriod _i,k : the method is free;PtUsPeriod _i : the method is free;BasePTO：09:45；NewPtUsPeriod _i : the method is free;UsPeriod: the method is free;BaseCTO: the method is free;RestrictCTO _i ：09:45；

the limiting point is PLT, i=2,PassMark _i =0, height layer (k): 1, a step of; flow management measures: r is R ₂ ，j=2；ExpMarki,j1 is shown in the specification; corresponding toHtUsPeriod _m,k ：[00:00,23:59]；PtMsHtUsPeriod _i,j,k ：；PtHtUsPeriod _i,k ：/>；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 1, a step of; flow management measures: r is R ₃ ，j=3；ExpMarki,jIs 0; corresponding toHtUsPeriod _m,k ：[00:00,10:00) (13:00,23:59]；PtMsHtUsPeriod _i,j,k ：；PtHtUsPeriod _i,k ：/>；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, highDegree layer (k): 2; flow management measures: r is R ₂ ，j=2；ExpMarki,j1 is shown in the specification; corresponding toHtUsPeriod _m,k ：[00:00,23:59]；PtMsHtUsPeriod _i,j,k ：；PtHtUsPeriod _i,k ：/>；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 2; flow management measures: r is R ₃ ，j=3；ExpMarki,jIs 0; corresponding toHtUsPeriod _m,k ：[00:00,10:00) (13:00,23:59]；PtMsHtUsPeriod _i,j,k ：；PtHtUsPeriod _i,k ：/>；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 3, a step of; flow management measures: r is R ₂ ，j=2；ExpMarki,j1 is shown in the specification; corresponding toHtUsPeriod _m,k ：[00:00,23:59]；PtMsHtUsPeriod _i,j,k ：；PtHtUsPeriod _i,k ：/>；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 3, a step of; flow management measures: r is R ₃ ，j=3；ExpMarki,jIs 0; corresponding toHtUsPeriod _m,k ：[00:00,10:00) (13:00,23:59]；PtMsHtUsPeriod _i,j,k ：；PtHtUsPeriod _i,k ：/>；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 4, a step of; flow management measures: r is R ₂ ，j=2；ExpMarki,j1 is shown in the specification; corresponding toHtUsPeriod _m,k ：[00:00,23:59]；PtMsHtUsPeriod _i,j,k ：[00:00,23:59]；PtHtUsPeriod _i,k ：[00:00,23:59]；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 4, a step of; flow management measures: r is R ₃ ，j=3；ExpMarki,jIs 0; corresponding toHtUsPeriod _m,k ：[00:00,23:59]；PtMsHtUsPeriod _i,j,k ：[00:00,23:59]；PtHtUsPeriod _i,k ：[00:00,23:59]；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 5, a step of; flow management measures: r is R ₂ ，j=2；ExpMarki,j1 is shown in the specification; corresponding toHtUsPeriod _m,k ：[00:00,23:59]；PtMsHtUsPeriod _i,j,k ：[00:00,23:59]；PtHtUsPeriod _i,k ：[00:00,10:00) (13:00,23:59]；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 5, a step of; flow management measures: r is R ₃ ，j=3；ExpMarki,jIs 0; corresponding toHtUsPeriod _m,k ：[00:00,10:00) (13:00,23:59]；PtMsHtUsPeriod _i,j,k ：[00:00,13:00) (13:00,23:59]；PtHtUsPeriod _i,k ：[00:00,10:00) (13:00,23:59]；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 6, preparing a base material; flow management measures: r is R ₂ ，j=2；ExpMarki,j1 is shown in the specification; corresponding toHtUsPeriod _m,k ：[00:00,23:59]；PtMsHtUsPeriod _i,j,k ：[00:00,23:59]；PtHtUsPeriod _i,k ：[00:00,10:00) (13:00,23:59]；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 6, preparing a base material; flow management measures: r is R ₃ ，j=3；ExpMarki,jIs 0; corresponding toHtUsPeriod _m,k ：[00:00,10:00) (13:00,23:59]；PtMsHtUsPeriod _i,j,k ：[00:00,13:00) (13:00,23:59]；PtHtUsPeriod _i,k ：[00:00,10:00) (13:00,23:59]；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 7, preparing a base material; flow management measures: r is R ₂ ，j=2；ExpMarki,j1 is shown in the specification; corresponding toHtUsPeriod _m,k ：[00:00,23:59]；PtMsHtUsPeriod _i,j,k ：[00:00,23:59]；PtHtUsPeriod _i,k ：[00:00,10:00) (13:00,23:59]；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 7, preparing a base material; flow management measures: r is R ₃ ，j=3；ExpMarki,jIs 0; corresponding toHtUsPeriod _m,k ：[00:00,10:00) (13:00,23:59]；PtMsHtUsPeriod _i,j,k ：[00:00,13:00) (13:00,23:59]；PtHtUsPeriod _i,k ：[00:00,10:00) (13:00,23:59]；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 8, 8; flow management measures: r is R ₂ ，j=2；ExpMarki,j1 is shown in the specification; corresponding toHtUsPeriod _m,k ：[00:00,23:59]；PtMsHtUsPeriod _i,j,k ：[00:00,23:59]；PtHtUsPeriod _i,k ：[00:00,10:00) (13:00,23:59]；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 8, 8; flow management measures: r is R ₃ ，j=3；ExpMarki,jIs 0; corresponding toHtUsPeriod _m,k ：[00:00,10:00) (13:00,23:59]；PtMsHtUsPeriod _i,j,k ：[00:00,13:00) (13:00,23:59]；PtHtUsPeriod _i,k ：[00:00,10:00) (13:00,23:59]；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 9, a step of performing the process; flow management measures: r is R ₂ ，j=2；ExpMarki,j1 is shown in the specification; corresponding toHtUsPeriod _m,k ：[00:00,23:59]；PtMsHtUsPeriod _i,j,k ：；PtHtUsPeriod _i,k ：/>；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 9, a step of performing the process; flow management measures: r is R ₃ ，j=3；ExpMarki,jIs 0; corresponding toHtUsPeriod _m,k ：[00:00,10:00) (13:00,23:59]；PtMsHtUsPeriod _i,j,k ：；PtHtUsPeriod _i,k ：/>；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 10; flow management measures: r is R ₂ ，j=2；ExpMarki, j1 is shown in the specification; corresponding toHtUsPeriod _m,k ：[00:00,23:59]；PtMsHtUsPeriod _i,j,k ：；PtHtUsPeriod _i,k ：/>；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40；

The limiting point is PLT, i=2,PassMark _i =0, height layer (k): 10; traffic managementMeasures are as follows: r is R ₃ ，j=3；ExpMarki, jIs 0; corresponding toHtUsPeriod _m,k ：[00:00,10:00) (13:00,23:59]；PtMsHtUsPeriod _i,j,k ：；PtHtUsPeriod _i,k ：/>；PtUsPeriod _i ：[00:00,23:59]；BasePTO：PLT(10:40)；NewPtUsPeriod _i ：[00:00,23:59]；UsPeriod：[00:00,23:59]；BaseCTO：10:40；RestrictCTO _i ：10:40。

4-2, calculating the target passing point height of the limiting point

And calculating the target passing point height of the flight at the limiting point according to the target passing point time of the flight at the limiting point and the available resource time period of each altitude layer. FlightDKH1287The analysis results are shown below:

the limiting point is NXD, i=1,PassMark _i =1, height layer (j): the method is free;PtHtUsPeriod _i,j : the method is free; whether or not to includeRestrictCTO _i : the method is free;RestrictCTO _i : the method is free;RestrictPHt _i ：2；RestrictCHt _i ：2；

the limiting point is PLT, i=2,PassMark _i =0, height layer (j): 1, a step of;PtHtUsPeriod _i,j ：the method comprises the steps of carrying out a first treatment on the surface of the Whether or not to includeRestrictCTO _i : if not, then judging whether the current is equal to or greater than the preset threshold;RestrictCTO _i ：10:40；RestrictPHt _i ：5；RestrictCHt _i ：4；

the limiting point is PLT, i=2,PassMark _i =0, height layer (j): 2;PtHtUsPeriod _i,j ：the method comprises the steps of carrying out a first treatment on the surface of the Whether or not to includeRestrictCTO _i : if not, then judging whether the current is equal to or greater than the preset threshold;RestrictCTO _i ：10:40；RestrictPHt _i ：5；RestrictCHt _i ：4；

the limiting point is PLT, i=2,PassMark _i =0, height layer (j): 3, a step of;PtHtUsPeriod _i,j ：the method comprises the steps of carrying out a first treatment on the surface of the Whether or not to includeRestrictCTO _i : if not, then judging whether the current is equal to or greater than the preset threshold;RestrictCTO _i ：10:40；RestrictPHt _i ：5；RestrictCHt _i ：4；

the limiting point is PLT, i=2,PassMark _i =0, height layer (j): 4, a step of;PtHtUsPeriod _i,j ：[00:00,23:59]the method comprises the steps of carrying out a first treatment on the surface of the Whether or not to includeRestrictCTO _i : is;RestrictCTO _i ：10:40；RestrictPHt _i ：5；RestrictCHt _i ：4；

the limiting point is PLT, i=2,PassMark _i =0, height layer (j): 5, a step of;PtHtUsPeriod _i,j ：[00:00,10:00) (13:00,23:59]the method comprises the steps of carrying out a first treatment on the surface of the Whether or not to includeRestrictCTO _i : if not, then judging whether the current is equal to or greater than the preset threshold;RestrictCTO _i ：10:40；RestrictPHt _i ：5；RestrictCHt _i ：4；

the limiting point is PLT, i=2,PassMark _i =0, height layer (j): 6, preparing a base material;PtHtUsPeriod _i,j ：[00:00,10:00) (13:00,23:59]the method comprises the steps of carrying out a first treatment on the surface of the Whether or not to include RestrictCTO _i : if not, then judging whether the current is equal to or greater than the preset threshold;RestrictCTO _i ：10:40；RestrictPHt _i ：5；RestrictCHt _i ：4；

the limiting point is PLT, i=2,PassMark _i =0, height layer (j): 7, preparing a base material;PtHtUsPeriod _i,j ：[00:00,10:00) (13:00,23:59]the method comprises the steps of carrying out a first treatment on the surface of the Is thatWhether or not to includeRestrictCTO _i : if not, then judging whether the current is equal to or greater than the preset threshold;RestrictCTO _i ：10:40；RestrictPHt _i ：5；RestrictCHt _i ：4；

the limiting point is PLT, i=2,PassMark _i =0, height layer (j): 8, 8;PtHtUsPeriod _i,j ：[00:00,10:00) (13:00,23:59]the method comprises the steps of carrying out a first treatment on the surface of the Whether or not to includeRestrictCTO _i : if not, then judging whether the current is equal to or greater than the preset threshold;RestrictCTO _i ：10:40；RestrictPHt _i ：5；RestrictCHt _i ：4；

the limiting point is PLT, i=2,PassMark _i =0, height layer (j): 9, a step of performing the process;PtHtUsPeriod _i,j ：the method comprises the steps of carrying out a first treatment on the surface of the Whether or not to includeRestrictCTO _i : if not, then judging whether the current is equal to or greater than the preset threshold;RestrictCTO _i ：10:40；RestrictPHt _i ：5；RestrictCHt _i ：4；

the limiting point is PLT, i=2,PassMark _i =0, height layer (j): 10;PtHtUsPeriod _i,j ：the method comprises the steps of carrying out a first treatment on the surface of the Whether or not to includeRestrictCTO _i : if not, then judging whether the current is equal to or greater than the preset threshold;RestrictCTO _i ：10:40；RestrictPHt _i ：5；RestrictCHt _i ：4。

4-3, calculating the target take-off time and the total delay of the flights

And calculating the departure time of the flight target and the total delay of the flight according to the information of the flight state, the reference predicted passing time, the reference target passing time and the like. FlightDKH1287The analysis results are shown below:

status:DEP；BasePTO：10:40；BaseCTO：10:40；MeasureTotalDelay：0；CTOT：。

5. calculating child delays of flights caused by traffic management measures

And decomposing the flight delay according to the unavailable resource time period of each flow management measure of each height layer of the flight at the limiting point, and calculating the sub delay of the flight caused by each flow management measure. FlightDKH1287The analysis results are shown below:

limiting pointNXD，i=1，PassMark _i =1，RestrictCHt _i :2; flow management measures:R ₁ ，j=1；PtMsHtUsPeriod _{i,j,RestrictCHti} ：[00:00,23:59]；FdPeriod _i,j ：；Total _i ：[09:45, 09:45]；FdPeriod _i,j Total _i ：/>；PartFdLen _i,j ：0；MeasureTotalDelay：0；PartDelay _i,j ：0；

limiting pointPLT，i=2，PassMark _i =1，RestrictCHt _i :4, a step of; flow management measures:R ₂ ，j=2；PtMsHtUsPeriod _{i,j,RestrictCHti} ：[00:00,23:59]；FdPeriod _i,j ：；Total _i ：[10:40, 10:40]；FdPeriod _i,j Total _i ：/>；PartFdLen _i,j ：0；MeasureTotalDelay：0；PartDelay _i,j ：0；

limiting pointPLT，i=2，PassMark _i =1，RestrictCHt _i :4, a step of; flow management measures:R ₃ ，j=3；PtMsHtUsPeriod _{i,j,RestrictCHti} ：[00:00,23:59]；FdPeriod _i,j ：；Total _i ：[10:40, 10:40]；FdPeriod _i,j Total _i ：/>；PartFdLen _i,j ：0；MeasureTotalDelay：0；PartDelay _i,j ：0。

6. space-time resource for updating each flow management measure

As shown in fig. 4, the available resource period and the unavailable resource period of each traffic management measure at different levels are updated according to the traffic management measure interval limit and the total amount limit. Traffic management measures R ₁ The analysis results are shown below:

height layer (j): 1, a step of;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：09:45；EarliestUsTime _i,j ：15；IntervalFdPeriod _i,j ：[09:30,10:00]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtUsPeriod _i,j ：[00:00,09:30) ( 10:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[09:30,10:00]；

Height layer (j): 2;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：09:45；EarliestUsTime _i,j ：15；IntervalFdPeriod _i,j ：[09:30,10:00]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtUsPeriod _i,j ：[00:00,09:30) ( 10:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[09:30,10:00]；

Height layer (j): 3, a step of;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：09:45；EarliestUsTime _i,j ：15；IntervalFdPeriod _i,j ：[09:30,10:00]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtUsPeriod _i,j ：[00:00,09:30) ( 10:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[09:30,10:00]；

Height layer (j): 4, a step of;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：09:45；EarliestUsTime _i,j ：15；IntervalFdPeriod _i,j ：[09:30,10:00]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtUsPeriod _i,j ：[00:00,09:30) ( 10:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[09:30,10:00]；

Height layer (j): 5, a step of;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：09:45；EarliestUsTime _i,j ：15；IntervalFdPeriod _i,j ：[09:30,10:00]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtUsPeriod _i,j ：[00:00,09:30) ( 10:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[09:30,10:00]；

Height layer (j): 6, preparing a base material;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：09:45；EarliestUsTime _i,j ：15；IntervalFdPeriod _i,j ：[09:30,10:00]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtUsPeriod _i,j ：[00:00,09:30) ( 10:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[09:30,10:00]；

Height layer (j): 7, preparing a base material;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：09:45；EarliestUsTime _i,j ：15；IntervalFdPeriod _i,j ：[09:30,10:00]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtUsPeriod _i,j ：[00:00,09:30) ( 10:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[09:30,10:00]；

Height layer (j): 8, 8;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：09:45；EarliestUsTime _i,j ：15；IntervalFdPeriod _i,j ：[09:30,10:00]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtUsPeriod _i,j ：[00:00,09:30) ( 10:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[09:30,10:00]；

Height layer (j): 9, a step of performing the process;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：09:45；EarliestUsTime _i,j ：15；IntervalFdPeriod _i,j ：[09:30,10:00]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtUsPeriod _i,j ：[00:00,09:30) ( 10:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[09:30,10:00]；

Height layer (j): 10;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：09:45；EarliestUsTime _i,j ：15；IntervalFdPeriod _i,j ：[09:30,10:00]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtUsPeriod _i,j ：[00:00,09:30) ( 10:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[09:30,10:00]。

Flow management measure R ₂ The analysis results are shown below:

height layer (j)）：1；HtResoureMgt _i,j :0; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；EarliestUsTime _i,j ：10:40；IntervalRes _i ：20；IntervalFdPeriod _i,j : none, updated HtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：。

Height layer (j): 2;HtResoureMgt _i,j :0; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；EarliestUsTime _i,j ：10:40；IntervalRes _i ：20；IntervalFdPeriod _i,j : none, updatedHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：。

Height layer (j): 3, a step of;HtResoureMgt _i,j :0; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；EarliestUsTime _i,j ：10:40；IntervalRes _i ：20；IntervalFdPeriod _i,j : none, updatedHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：。

Height layer (j): 4, a step of;HtResoureMgt _i,j :0; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；EarliestUsTime _i,j ：10:40；IntervalRes _i ：20；IntervalFdPeriod _i,j : none, updatedHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：。

Height layer (j): 5, a step of;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；EarliestUsTime _i,j ：10:40；IntervalRes _i ：20；IntervalFdPeriod _i,j ：[10:20,11:00]Updated, updatedHtUsPeriod _i,j ：[00:00, 10:20) (11:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:20,11:00]。

Height layer (j): 6, preparing a base material;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；EarliestUsTime _i,j ：10:40；IntervalRes _i ：20；IntervalFdPeriod _i,j ：[10:20,11:00]Updated, updatedHtUsPeriod _i,j ：[00:00, 10:20) (11:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:20,11:00]。

Height layer (j): 7, preparing a base material;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；EarliestUsTime _i,j ：10:40；IntervalRes _i ：20；IntervalFdPeriod _i,j ：[10:20,11:00]Updated, updatedHtUsPeriod _i,j ：[00:00, 10:20) (11:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:20,11:00]。

Height layer (j): 8, 8;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；EarliestUsTime _i,j ：10:40；IntervalRes _i ：20；IntervalFdPeriod _i,j ：[10:20,11:00]Updated, updatedHtUsPeriod _i,j ：[00:00, 10:20) (11:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:20,11:00]。

Height layer (j): 9, a step of performing the process;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；EarliestUsTime _i,j ：10:40；IntervalRes _i ：20；IntervalFdPeriod _i,j ：[10:20,11:00]Updated, updatedHtUsPeriod _i,j ：[00:00, 10:20) (11:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:20,11:00]。

Height layer (j): 10;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00,23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；EarliestUsTime _i,j ：10:40；IntervalRes _i ：20；IntervalFdPeriod _i,j ：[10:20,11:00]Updated, updatedHtUsPeriod _i,j ：[00:00, 10:20) (11:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:20,11:00]。

Flow management measure R ₃ The analysis results are shown below:

height layer (j): 1, a step of;HtResoureMgt _i,j :0; pre-updateHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；HtMonCount _i,j,k : the method is free;CountFdPeriod _i,j : the method is free; updatedHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:00,13:00]；

Height layer (j): 2;HtResoureMgt _i,j :0; pre-update HtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；HtMonCount _i,j,k : the method is free;CountFdPeriod _i,j : the method is free; updatedHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:00,13:00]；

Height layer (j): 3, a step of;HtResoureMgt _i,j :0; pre-updateHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；HtMonCount _i,j,k : the method is free;CountFdPeriod _i,j : the method is free; updatedHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:00,13:00]；

Height layer (j): 4, a step of;HtResoureMgt _i,j :1, a step of; pre-updateHtUsPeriod _i,j ：[00:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；HtMonCount _i,j,k Comprising 6:

[10:00,10:30)，1；

[10:30,11:00)，0；

[11:00,11:30)，1；

[11:30,12:00)，1；

[12:00,12:30)，1；

[12:30,13:00)，1；

CountFdPeriod _i,j : [10:30, 11:00); updatedHtUsPeriod _i,j ：[00:00, 10:30) [13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j Comprising 6: [10:30,11:00]；[10:30,11:00]；[10:30,11:00]；[10:30,11:00]；[10:30,11:00]；[10:30,11:00]；

Height layer (j): 5, a step of;HtResoureMgt _i,j :0; pre-updateHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；HtMonCount _i,j,k : the method is free;CountFdPeriod _i,j : the method is free; updatedHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:00,13:00]；

Height layer (j): 6, preparing a base material;HtResoureMgt _i,j :0; pre-updateHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；HtMonCount _i,j,k : the method is free;CountFdPeriod _i,j : the method is free; updatedHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:00,13:00]；

Height layer (j): 7, preparing a base material;HtResoureMgt _i,j :0; pre-updateHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；HtMonCount _i,j,k : the method is free;CountFdPeriod _i,j : the method is free; updatedHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:00,13:00]；

Height layer (j): 8, 8;HtResoureMgt _i,j :0; pre-updateHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；HtMonCount _i,j,k : the method is free;CountFdPeriod _i,j : the method is free; updatedHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:00,13:00]；

Height layer (j): 9, a step of performing the process;HtResoureMgt _i,j :0; pre-updateHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；HtMonCount _i,j,k : the method is free;CountFdPeriod _i,j : the method is free; updatedHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:00,13:00]；

Height layer (j): 10;HtResoureMgt _i,j :0; pre-updateHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the Corresponding toRestrictCTO _m ：10:40；HtMonCount _i,j,k : the method is free;CountFdPeriod _i,j : the method is free; updatedHtUsPeriod _i,j ：[00:00, 10:00) (13:00, 23:59]The method comprises the steps of carrying out a first treatment on the surface of the UpdatedHtFdPeriod _i,j ：[10:00,13:00]；

The invention realizes multidimensional evaluation of the flight priority, evaluates the flight scheduling sequence through the flight basic information and the flight limited information, generates the flight scheduling queue, and provides a fairness foundation facing the operation rule for subsequent flight sequencing calculation. In the process of scheduling the flights, a series of decision results such as target passing time and target passing height of the flights at the limiting points, target departure time and total delay of the flights are obtained by carrying out logic operation on available resource time periods of each altitude layer of the flow management measures related to the flights. Meanwhile, the total flight delay is decomposed by analyzing the unavailable resource time period of each flow management measure, and the flight sub delay caused by each flow management measure is obtained. And finally, updating and maintaining the available resource time period of each flow management measure according to the calculation result of the single flight. The calculation process has strong interpretability, stable calculation result and good calculation timeliness, and can meet the flight integrated scheduling requirement of the air traffic control operation under the complex scene for the constraint of the multiple flow management measures.

The invention provides a method for integrally dispatching flights under multiple constraints, and the method and the way for realizing the technical scheme are numerous, the above description is only a preferred embodiment of the invention, and it should be pointed out that, for a person skilled in the art, several improvements and modifications can be made without departing from the principle of the invention, and the improvements and modifications should be regarded as the protection scope of the invention. The components not explicitly described in this embodiment can be implemented by using the prior art.

Claims (10)

1. The integrated scheduling method for the flights under the multi-element constraint is characterized by comprising the following steps of:

step 1, analyzing and processing comprehensive operation information, wherein the comprehensive operation information comprises: flight plan information and traffic management measure information;

step 2, initializing space-time resources of a high-level flow management measure according to the flow management measure information;

step 3, constructing a multidimensional flight priority evaluation model, calculating the flight priority according to the comprehensive operation information, and generating a flight scheduling queue;

step 4, calculating target passing point time, target passing point height, target departure time and total delay of the flight at the flow management measure limiting point;

Step 5, decomposing the total delay of the flights according to the space-time resources of the high-level layers of the flow management measures to obtain the sub delays of the flights caused by the flow management measures;

and step 6, updating the space-time resources of the flow management measure height layer according to the calculated target point passing time of the flow management measure limiting point.

2. The method according to claim 1, wherein step 1 comprises the steps of:

step 1-1: acquiring and resolving empty pipe height layer sets from comprehensive operation platformHeightSet={Height ₁ ，…，Height _i }, whereinHeight _i Representing what is used by the empty pipe to run normallyThe height value of the ith height layer used,HtNumfor all the number of layers of the height,，Height _i > Height _i-1 ；

acquiring the total set of all the time of the operation day asTimeSet=[00:00，23:59]；

Step 1-2: acquiring and resolving flight plan information from an integrated operation platformFlightInfoFlight plan informationFlightInfoComprising the following steps: flight call signCallsignUnique number of flightPlanIDStatus of flightStatusExemption attribute of flight associated traffic management measureExemptScheduled gear withdrawal timeSOBTPredicted take-off timeETOTAirport for taking offDepAirportAirport for landingDesAirportFlow management measures associated with flights limit the number of pointsResPtNumAnd flight restriction informationPtResourceSet；

The flight status includes: not being executedNoExecutedGround executionGroundHas taken offDEPHas fallen downARR；

Exemption attribute for traffic management measures Exempt=1, indicating that flights are exempted from all flow management measures; if the flight exempts from one flow management measure, the flight is indicated to be free from considering the limiting condition of the flow management measure in the calculation process;

PtResourceSet={PtResource ₁ ，…， PtResource _i }， PtResource _i the ith restriction point information indicating the arrival of the flight includes: limiting point nameRestrictPt _i Number of traffic management measures associated with limiting pointMeasureNum _i Time of flight arrival at limit pointRestrictPTO _i Flight arrival limit point altitude layer sequence numberRestrictPHt _i Flight arrival limit point altitude rangeRestrictHtRange _i Actual passing point indicationPassMark _i Flow rate of limiting pointManagement measure setPtMsResourceSet _i ，RestrictPTO _i >RestrictPTO _i-1 ；

If it isPassMark _i =1, indicating that the flight has flown past the ith limit point,RestrictPTO _i the actual passing point time of the ith limiting point; if it isPassMark _i =0, indicating that the flight did not fly past the ith limit point,RestrictPTO _i the predicted passing time for the ith limit point;

PtMsResourceSet _i ={PtMsResource _i,1 ，…，PtMsResource _i,j }，PtMsResource _i,j the j-th associated traffic management measure information indicating that the flight is at the i-th restriction point, comprising: associated traffic management measure unique numberPtMeasureID _i,j Associated traffic management measure classPtRank _i,j Exemption indication for associated flow management measuresExpMark _i,j ；

If it isExpMark _i,j =1, indicating that the flight is exempted from the jth associated traffic management measure at the ith limit point; if it isExpMark _i,j =0, indicating that the flight is restricted for the jth associated traffic management measure at the ith restriction point;

RestrictHtRange _i ={RestrictHt _i,1 ，…，RestrictHt _i,m }，RestrictHt _i,m Indicating that the flight is at the highest level possible at the ith limit pointHeightSetThe number of which is a sequence number,RestrictHt _i,m >RestrictHt _i,m ；

step 1-3: the following information is extended:

at the position ofFlightInfoIn extension information: flight target departure timeCTOT；

At the position ofPtResource _i In extension information: flightTarget passing time to reach limit pointRestrictCTO _i Limiting point target overpass heightRestrictCHt _i Limiting the time period of available point resourcesPtUsPeriod _i Limiting point high-level available resource setPtHtUsPeriodSet _i ={PtHtUsPeriod _i,1 ，…，PtHtUsPeriod _i,j }，PtHtUsPeriod _i,j Representing the availability period of the jth altitude layer of the flight at the ith limit point;

at the position ofPtMsResource _i,j The extension definition in (a): flight delay caused by traffic management measuresPartDelay _i,j High-level availability period set for traffic management measuresPtMsHtUsPeriodSet _i,j ={PtMsHtUsPeriod _i,j,1 ，…，PtMsHtUsPeriod _i,j,k }，PtMsHtUsPeriod _i,j,k Indicating the period of availability of the flight at the kth altitude layer of the jth traffic management measure at the ith limit point,；

step 1-4: creating a set of traffic management measuresMeasurInfoSet={ MeasurInfo ₁ ，…，MeasurInfo _i Acquiring and analyzing flow management measure information from the comprehensive operation platformMeasureInfo _i Traffic management measure informationMeasureInfo _i Comprising the following steps: unique number of flow management measureMeasureID _i Start timeBgnTime _i End timeEndTime _i Duration of effectivenessDuration _i Grade (grade)Rank _i Total amount limitationCountRes _i Time granularity of total limitCountReTimesSize _i Spacing limitIntervalRes _i Limiting the heightHeightResSet _i Height limiting typeHeightType _i ；

Duration _i =EndTime _i -BgnTime _i ，Duration _i /CountReTimesSize _i =PeriodNum _i ，Duration _i %CountReTimesSize _i =0，PeriodNum _i Indicating that the ith traffic management action is in effect period of timeCountReTimesSize _i The total amount obtained by dividing limits the number of the monitoring sub-periods;

Rank _i Comprising three levels 1,2, 3;

HeightResSet _i ={ResHtIndex _i,1 ，…，ResHtIndex _i,j }，ResHtIndex _i,j limited height layer representing the ith traffic management measureHeightSetThe number of which is a sequence number,ResHtIndex _{i, j} > ResHtIndex _i,j-1 ；

if it isHeightType _i =0Indicating that flights pass through flow management measure limiting points at any altitude, the flights are limited by total amount and interval;

if it isHeightType _i =1If the altitude layer sequence number indicating that the flight passes the flow management measure limiting point is contained inHeightResSet _i If the flight is not limited by the total amount, the interval is limited by the total amount;

if it isHeightType _i =2Indicating that flights can only be included in sequence numbersHeightResSet _i The internal altitude layers pass through the flow management measure limiting points and are limited by total amount and interval, and flights cannot pass through other altitude layers in the measure period [BgnTime _i ，EndTime _i ]The flow management measure limiting point is passed through;

step 1-5: the following information is extended:

at the position ofMeasureInfo _i The extension definition in (a): traffic management measures high-level space-time resourcesHtResourceSet _i ={HtResource _i,1 ，…，HtResource _i,j }，HtResource _i,j Resources representing the ith traffic management measure at the jth level include: resource management indicatorHtResourceMgt _i,j Period of high-level available resourcesHtUsPeriod _i,j Period of high-level unavailable resourceHtFdPeriod _i,j Aggregate of highly-layer aggregate resource monitoring periodsHtMonitor _i,j Resource counting in a high-level aggregate resource monitoring periodHtCount _i,j ，；

If it isHtResourceMgt _i,j =0, meaning that the ith traffic management measure is not limited by total amount, interval at the jth high level group, then The method comprises the steps of carrying out a first treatment on the surface of the If it isHtResourceMgt _i,j =1, indicating that the ith traffic management measure is limited by the total amount, spacing, at the jth level group;

HtMonitor _i,j ={ HtMonPeriod _i,j,1 ，…，HtMonPeriod _i,j,m }，HtMonPeriod _i,j,m an mth monitoring sub-period representing an jth level of an ith traffic management measure,HtMonPeriod _{i, j,m} =[BgnTime _i +(m-1)* CountReTimesSize _i ，BgnTime _i +m*CountReTimesSize _i )]，；

HtCount _i,j ={ HtMonCount _i,j,1 ，…，HtMonCount _i,j,m }，HtMonCount _i,j,m the number of remaining flight resources for the mth monitoring sub-period of the jth elevation layer representing the ith traffic management measure.

3. The method according to claim 2, wherein step 2 comprises the steps of:

step 2-1: traversing traffic management measure informationMeasureInfo _i ；

Step 2-2; initializing space-time resources of a high-level layer of traffic management measures:

；

if it isHeightType _i =0，

；

；

If it isHeightType _i =1，

；

；

If it isHeightType _i =2，

；

；

，

HtFdPeriod _i,j =[BgnTime _i ,EndTime _i ], HtUsPeriod _i,j =TimeSet-HtFdPeriod _i,j 。

4. A method according to claim 3, wherein step 3 comprises the steps of:

step 3-1: constructing a multidimensional flight priority assessment model comprising 13 packets:

PriorityGroup ₁ : the flight group has fallen, grouping conditions:Status=ARR；

PriorityGroup ₂ : group of flights taken off, grouping conditions: does not satisfyPriorityGroup ₁ Grouping conditions, andStatus=DEP；

PriorityGroup ₃ : airlines have been executed on the ground, grouping conditions: does not satisfyPriorityGroup ₁ ~PriorityGroup ₂ Grouping conditions, andStatus=Ground；

PriorityGroup ₄ : full-exemption airlines, grouping conditions: does not satisfyPriorityGroup ₁ ~ PriorityGroup ₃ The grouping condition is that,Exempt=1；

PriorityGroup ₅ : and (3) associating class 1 traffic management measures with the full-exemption airlines, and grouping conditions: does not satisfyPriorityGroup ₁ ~ PriorityGroup ₄ The grouping condition is that,PtRank _i,j maximum value is 1, and when PtRank _i,j When=1, satisfyExpMark _i,j =1；

PriorityGroup ₆ : and associating class 1 traffic management measures with semi-exemption airlines, and grouping conditions: does not satisfyPriorityGroup ₁ ~ PriorityGroup ₅ The grouping condition is that,PtRank _i,j maximum value is 1, and whenPtRank _i,j When=1, there isExpMark _i,j =0；

PriorityGroup ₇ : traffic management measures of association level 1 are all restricted airlines, grouping conditions: does not satisfyPriorityGroup ₁ ~ PriorityGroup ₆ The grouping condition is that,PtRank _i,j maximum value is 1, and whenPtRank _i,j When=1, satisfyExpMark _i,j =0；

PriorityGroup ₈ : and (3) associating class 2 traffic management measures with the full-exemption airlines, and grouping conditions: does not satisfyPriorityGroup ₁ ~ PriorityGroup ₇ The grouping condition is that,PtRank _i,j maximum value is 2, and whenPtRank _i,j When=2, satisfyExpMark _i,j =1；

PriorityGroup ₉ : and associating class 2 traffic management measures with semi-exemption airlines, and grouping conditions: does not satisfyPriorityGroup ₁ ~ PriorityGroup ₈ The grouping condition is that,PtRank _i,j maximum value is 2, and whenPtRank _i,j When=2, there isExpMark _i,j =0；

PriorityGroup ₁₀ : traffic management measures of association class 2 are all restricted airlines, grouping conditions: does not satisfyPriorityGroup ₁ ~ PriorityGroup ₉ The grouping condition is that,PtRank _i,j maximum value is 2, and whenPtRank _i,j When=2, satisfyExpMark _i,j =0；

PriorityGroup ₁₁ : and (3) associating class 3 traffic management measures with the full-exemption airlines, and grouping conditions: does not satisfyPriorityGroup ₁ ~ PriorityGroup ₁₀ The grouping condition is that,PtRank _i,j maximum value is 3, and whenPtRank _i,j When=3, satisfyExpMark _i,j =1；

PriorityGroup ₁₂ : and associating class 3 traffic management measures with semi-exemption airlines, and grouping conditions: does not satisfyPriorityGroup ₁ ~ PriorityGroup ₁₁ The grouping condition is that,PtRank _i,j maximum value is 3, and whenPtRank _i,j When=3, there isExpMark _i,j =0；

PriorityGroup ₁₃ : traffic management measures of association level 3 are all restricted airlines, grouping conditions: does not satisfy PriorityGroup ₁ ~ PriorityGroup ₁₂ The grouping condition is that,PtRank _i,j maximum value is 3, and whenPtRank _i,j When=3, satisfyExpMark _i,j =0；

Step 3-2: traversing each flight;

step 3-3: judging whether the flight meets the condition of priority grouping or not, if soPriorityGroup _i For the packet condition of (a), then the flight priority packetsPriorityIndex=i；

Step 3-4: after traversing the flights, judging every two flightsFlight _i AndFlight _j is of the flight priority of (2)FlightPriority _i AndFlightPriority _j is a size relationship of (a):

if it isPriorityIndex _i >PriorityIndex _j ThenFlightPriority _i >FlightPriority _j ；

If it isPriorityIndex _i =PriorityIndex _j And (2) andSOBT _i <SOBT _j thenFlightPriority _i > FlightPriority _j ；

If it isPriorityIndex _i =PriorityIndex _j And (2) andSOBT _i =SOBT _j and (2) andPlanID _i >PlanID _j thenFlightPriority _i >FlightPriority _j ；

The flight unique numbers will not be the same, i.e. ifThen->；

Step 3-5: ordering the flights according to the ascending order of the priority of the flights to obtain a flight scheduling queueFlightArray，FlightArray={Flight ₁ ，…，Flight _i }, whereinFlightPriority _i > FlightPriority _i-1 。

5. The method of claim 4, wherein step 4 comprises the steps of:

step 4-1: at the position ofFlightArraySequentially calculates for each flight, and records the currently calculated flight asFight；

Step 4-2: traversing restricted information for flightsPtResourceSetCalculating target passing time of flight at limiting pointRestrictCTO _i ；

Step 4-3: calculating a target passing point altitude of a flight at a limiting pointRestrictCHt _i ；

Step (a)4-4: calculating total delay of flights due to flow management measuresMeasureTotalDelayTakeoff time of flight targetCTOT。

6. The method of claim 5, wherein step 4-2 comprises: if it is PassMark _i =1, thenRestrictCTO _i =RelatedPTO _i The target passing time of the limiting point indicating that the flight has flown is equal to the actual passing time;

if it isPassMark _i =0, then:

when (when)If at the timeExpMark _i,j =0，The available resource period of the kth altitude layer is equal to the value of the (j) th associated traffic management measure representing the (j) th limit point of the (i) th limit pointMeasurInfoSetThe mth traffic management measure at the kth level of available resources period, wherein the kth level of available levels must be the level of available levels for flights at the ith limit point;

if it isExpMark _i,j =1，When the j-th associated flow management measure of the flight at the i-th limiting point is exempted, the available resource period of the k-th height layer is a full period, wherein the k-th height layer is the available height layer of the flight at the i-th limiting point;

the available resource period of the unavailable altitude layer of the j-th associated traffic management measure of the flight at the i-th limiting point is represented as empty;

representing that the available resource period of the kth high-level of the flight at the ith limit point is equal to the intersection of the available resource periods of the high-level of the respective traffic management measures of the flight at the ith limit point;

representing that the available resource period of the flight at the ith limit point is equal to the union of the available resource periods of the flights at each altitude layer of the ith limit point;

Calculating a limit point reference expected transit time；

Calculating the difference between the predicted passing time of each limiting point and the reference predicted passing timePTODvalue _i ， PTODvalue _i = RestrictPTO _i -BasePTO；

Performing translation processing on the available resource time period of the limiting point to obtain the updated available resource time period of the limiting point；

Taking intersection of available resource time periods updated by all limiting points to obtain available resource time periods of flights；

Obtaining a reference target passing point time by taking the minimum value of the reference predicted passing point time of the limiting point in the available resource time period of the flightBaseCTO，；

Calculating target passing time of flight at each limiting pointRestrictCTO _i =RestrictPTO _i +(BaseCTO- BasePTO)。

7. The method of claim 6, wherein step 4-3 comprises: if it isPassMark _i =1, thenRestrictCHt _i =RestrictPHt _i The target passing point height of the limiting point which indicates that the flight has flown is equal to the actual passing point height;

if it isPassMark _i =0, thenRestrictCHt _i =HSo thatThe target passing point height layer of the flight at the ith limiting point is the H-th height layer, the H-th height layer is closest to the predicted passing point height layer of the limiting point, and the available resource period of the H-th height layer comprises the target passing point time of the limiting point.

8. The method of claim 7, wherein step 4-4 comprises: if presentPassMark _i =1, thenIndicating that the taken-off flight does not need to calculate the target take-off time;

If all ofPassMark _i =1, thenMeasureTotalDelay=0Indicating that the flight has passed all restriction points, the traffic management measures do not cause delays;

if presentPassMark _i =0, thenMeasureTotalDelay=BaseCTO-BasePTO；

If all ofPassMark _i =0, then:

MeasureTotalDelay=BaseCTO- BasePTO；

CTOT=ETOT+ MeasureTotalDelay。

9. the method of claim 8, wherein step 5 comprises the steps of:

step 5-1: calculating a child delay of the flight caused by the j-th associated traffic management measure of the i-th restriction pointPartDelay _i,j ；

Step 5-2: if it isMeasureTotalDelay=0, thenPartDelay _i,j =0, indicating that the total delay of the flight is 0, and the child delay is 0;

step 5-3: if it isThen the following determination is made:

if it isExpMark _i,j =0, thenPartDelay _i,j =0, indicating that the associated exemption traffic management measure does not cause delays to flights;

if it isThen:

memory collection；

Memory collection，FdPeriod _i,j A period of unavailable resources representing a jth associated traffic management measure for the flight at the ith limit point;

recording device，PartFdLen _i,j Time period of unavailable resources and target overpass for jth associated traffic management measure representing flights at ith limit pointTotal _i The number of time elements of the intersection of (a);

calculated to obtain。

10. The method according to claim 9, wherein step 6 comprises the steps of:

step 6-1: after calculation for each flight, the traffic management measure high-level space-time resources are updated HtResourceSet _i ；

Step 6-2: spacing restrictions based on traffic management measuresIntervalRes _i Updating:

if it isThe flow management measures are not limited at intervals, so that space-time resources of the flow management measures do not need to be updated;

if it isAnd (2) andHtResourceMgt _i,j =0, which indicates that the traffic management measure is not limited at the jth level, so that the level space-time resource does not need to be updated;

if it isAnd (2) andHtResourceMgt _i,j when the traffic management measure satisfies =1MeasureID _i =PtMeasureID _m,n When the method is performed, the following determination is made:

if it isExpMark _m,n =0, then generate the interval unavailable period of the ith traffic management measure at the jth levelIntervalFdPeriod _i,j =[RestrictCTO _m -IntervalRes _i ，RestrictCTO _m +IntervalRes _i ]；

If it isExpMark _m,n =1, obtainHtUsPeriod _i,j Middle greater than or equal toRestrictCTO _m Is recorded as the earliest time of availabilityThe method comprises the steps of carrying out a first treatment on the surface of the Generating a temporary unavailable period of traffic management measure intervalIntervalFdPeriod _i,j =[EarliestUsTime _i,j -IntervalRes _i ，EarliestUsTime _i,j +IntervalRes _i ]；

Updating the available period of the ith traffic management action at the jth level；

Updating the ith traffic management measure at the unavailable period of the jth level；

Step 6-3: total amount restriction based on traffic management measuresCountRes _i Updating:

if it isThe flow management measures are not limited in total amount, so that space-time resources of the flow management measures do not need to be updated;

if it isAnd (2) andHtResourceMgt _i,j =0, which indicates that the traffic management measure is not limited at the jth level, so that the level space-time resource does not need to be updated;

if it isAnd (2) and HtResourceMgt _i,j When the traffic management measure satisfies =1MeasureID _i =PtMeasureID _m,n When the method is performed, the following determination is made:

if it isExpMark _m,n =0, andthen updateHtMonCount _i,j,k =HtMonCount _i,j,k -1；

If it isExpMark _m,n =1, obtainHtUsPeriod _i,j Middle greater than or equal toRestrictCTO _m Is recorded as the earliest time of availabilityThe method comprises the steps of carrying out a first treatment on the surface of the If it isThen updateHtMonCount _i,j,k =HtMonCount _i,j,k -1；

If after updatingGenerating a total unavailable period of the ith traffic management measure at the jth levelCountFdPeriod _i,j = HtMonPeriod _i,j,k ；

Updating the available period of the ith traffic management measure at the jth level；

Updating the ith traffic management measure at the unavailable period of the jth level。