Time-space conversion method of flight sequencing information

Abstract

A time-space conversion method of flight sequencing information aims to solve a problem that a flight sequencing suggestion generated by an arrival management system has poor effect in practical application due to abstract characteristics, comprising: generating a sequencing time and a delay suggestion of a flight in each key point in a terminal area and each runway by using a flight sequencing technology; predicting a flight status of the flight according to a current position and 4D trajectory information of the flight, and filtering a delay-consumed flight segment, and on that basis, generating a flight segment delay allocation strategy with reference to an aircraft performance and sequencing information, and obtaining a flight segment delay allocation result; and generating a visual spatial position reference target according to the flight segment delay allocation result and an operation deviation limit.

Claims

1. A time-space conversion method of flight sequencing information, comprising a computer readable medium operable on a computer with memory for the time-space conversion method, and comprising program instructions for executing the following steps of: step 1: sequencing and spacing a flight; and generating a sequencing time and a delay suggestion of the flight in each key point in a terminal area and each runway; step 2: allocating a flight segment delay; and predicting a flight status of the flight according to a current position and 4D trajectory information of the flight, and filtering a delay-consumed flight segment, and on that basis, generating a flight segment delay allocation strategy with reference to an aircraft performance and sequencing information, and obtaining a flight segment delay allocation result; wherein the step 2 comprises: predicting the flight status of the flight according to the current position and the 4D trajectory information of the flight, and filtering the delay-consumed flight segment, and on that basis, generating the flight segment delay allocation strategy with reference to the aircraft performance and the sequencing information, comprising the following steps of: step 2-1: defining variables; step 2-2: filtering the delay-consumed flight segment; wherein the step 2-2 comprises the following steps of: letting Pt.sub.i,bgn be an origin of the flight segments participating in delay consumption in the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i, and bgn denotes a serial number of the point Pt.sub.i,bgn in the queue PtList.sub.i; step 2-2-1: positioning an actual flight segment at which the flight is currently located: positioning the actual flight segment [Pt.sub.i,cur, Pt.sub.i,cur+1] at which the flight Flt.sub.i is currently located according to the estimated passing time ETO.sub.i,j of each point Pt.sub.i,j in the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i, and satisfying that SysTime∈[ETO.sub.i,cur, ETO.sub.i,cur+1]; wherein, Pt.sub.i,cur denotes an origin of the actual flight segment at which the flight Flt.sub.i is currently located, and Pt.sub.i,cur+1 denotes an end point of the actual flight segment at which the flight Flt.sub.i is currently located; cur denotes a serial number of the point Pt.sub.i,cur in the queue PtList.sub.i, and cur+1 denotes a serial number of the point Pt.sub.i,cur+1 in the queue PtList.sub.i; ETO.sub.i,cur denotes an estimated passing time of the flight Flt.sub.i at the point Pt.sub.i,cur, and ETO.sub.i,cur+1 denotes an estimated passing time of the flight Flt.sub.i at the point Pt.sub.i,cur+1; step 2-2-2: positioning a reference flight segment at which the flight is currently located: positioning the reference flight segment [Pt.sub.i,ref, Pt.sub.i,ref+1] at which the flight Flt.sub.i is currently located according to the sequencing passing time CTO.sub.i,j of each point Pt.sub.i,j in the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i, and satisfying that SysTime∈[CTO.sub.i,ref, CTO.sub.i,ref+1]; wherein, Pt.sub.i,ref denotes an origin of the reference flight segment at which the flight Flt.sub.i is currently located, and Pt.sub.i,ref+1 denotes an end point of the reference flight segment at which the flight Flt.sub.i is currently located; ref denotes a serial number of the point Pt.sub.i,ref in the queue PtList.sub.i, and ref+1 denotes a serial number of the point Pt.sub.i,ref+1 in the queue PtList.sub.i; CTO.sub.i,ref denotes a sequencing passing time of the flight Flt.sub.i at the point Pt.sub.i,ref, and CTO.sub.i,ref+1 denotes a sequencing passing time of the flight Flt.sub.i at the point Pt.sub.i,ref+1; step 2-2-3: searching for an adjacent sequencing key point in a preceding flight segment: letting Pt.sub.i,pre be the sequencing key point closest to a current position Pt.sub.i,cur in the preceding flight segment [Pt.sub.i,1, Pt.sub.i,cur] of the flight Flt.sub.i, wherein Pre denotes a serial number of the point Pt.sub.i,pre in the queue PtList.sub.i; letting Pt.sub.i,tmp be an intermediate variable of flight trajectory points of the flight Flt.sub.i in the calculation process of the method, wherein tmp denotes a serial number of the point Pt.sub.i,tmp in the queue PtList.sub.i; when satisfying that ∃Pt.sub.i,tmp∈[Pt.sub.i,1, Pt.sub.i,cur], and satisfying:
Min{(ETO.sub.i,cur−ETO.sub.i,tmp+1)*PtPro.sub.i,tmp}>0  (1) letting Pt.sub.i,pre=Pt.sub.i,tmp, and continuously executing step 2-2-4; otherwise, letting Pt.sub.i,pre=Ø, and executing step 2-2-5; step 2-2-4: judging whether a sequencing key point exists between the actual position of the flight and the reference flight segment: when satisfying that ETO.sub.i,cur>ETO.sub.i,ref and Pt.sub.i,pre∈[Pt.sub.i,ref+1, Pt.sub.i,cur], letting Pt.sub.i,bgn=Pt.sub.i,pre, and executing step 2-2-7; otherwise, continuously executing step 2-2-5; step 2-2-5: searching for an adjacent sequencing key point in a subsequent flight segment: letting Pt.sub.i,PtNum.sub.i be the last point in the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i; letting Pt.sub.i,aft be the sequencing key point closest to the current position Pt.sub.i,cur in the subsequent flight segment [Pt.sub.i,cur+1, Pt.sub.i,PtNum.sub.i] of the flight Flt.sub.i, wherein aft denotes a serial number of the point Pt.sub.i,aft in the queue PtList.sub.i; when satisfying that ∃Pt.sub.i,tmp∈[Pt.sub.i,cur+1, Pt.sub.i,PtNum.sub.i] and satisfying:
Min{(ETO.sub.i,tmp−ETO.sub.i,cur)*PtPro.sub.i,tmp}>0  (2) letting Pt.sub.i,aft=Pt.sub.i,tmp; otherwise, letting Pt.sub.i,aft=Pt.sub.i,PtNum.sub.i; step 2-2-6: judging whether the sequencing time of the subsequent sequencing key point has a certain fluctuation, wherein the method comprises: calculating a difference Div(CTO.sub.i,aft) between the currently allocated sequencing time and the allocated sequencing time during last calculation of the flight Flt.sub.i at the sequencing key point Pt.sub.i,aft, and filtering the origin Pt.sub.i,bgn of the delay-consumed flight segment according to Div(CTO.sub.i,aft), wherein the method is as follows: for the first operation, letting Div(CTO.sub.i,aft)=0; $\begin{array}{cc}{\mathrm{Pt}}_{i,\mathrm{bgn}}=\{\begin{array}{cc}{\mathrm{Pt}}_{i,\mathrm{ref}},& .Math.\mathrm{Div}\left({\mathrm{CTO}}_{i,\mathrm{aft}}\right).Math.<\mathrm{DivLimit}\\ {\mathrm{Pt}}_{i,\mathrm{cur}},& .Math.\mathrm{Div}\left({\mathrm{CTO}}_{i,\mathrm{aft}}\right).Math.\ge \mathrm{DivLimit}\&\&{\mathrm{Pt}}_{i,\mathrm{pre}}==\varnothing \\ {\mathrm{Pt}}_{i,\mathrm{pre}},& \mathrm{else}\end{array};& \left(3\right)\end{array}$ step 2-2-7: determining the delay-consumed flight segment: for the flight Flt.sub.i, the filtered flight segments participating in delay consumption are all the subsequent flight segment from Pt.sub.i,bgn in the flight trajectory point queue PtList.sub.i, which is [Pt.sub.i,bgn, Pt.sub.i,PtNum.sub.i]; and step 2-2-8: updating a sequencing time of the origin of the delay-consumed flight segment: $\begin{array}{cc}{\mathrm{CTO}}_{i,\mathrm{bgn}}=\{\begin{array}{cc}{\mathrm{ETO}}_{i,\mathrm{bgn}},& {\mathrm{Pt}}_{i,\mathrm{bgn}}\ne {\mathrm{Pt}}_{i,\mathrm{ref}}\&\&{\mathrm{ETO}}_{i,\mathrm{bgn}}\le \mathrm{SysTime}\\ {\mathrm{CTO}}_{i,\mathrm{bgn}},& \mathrm{else}\end{array}& \left(4\right)\end{array}$ for the first calculation, letting a sequencing time of all non-sequencing key points Pt.sub.i,j in PtList.sub.i be CTO.sub.i,j=ETO.sub.i,j; otherwise, using the last operation result for CTO.sub.i,j of all the non-sequencing key points step 2-3: dividing the flight segment according to a sequencing key point; step 2-4: dividing the flight segment according to the flight status; and step 2-5: allocating the flight segment delay to obtain the flight segment delay allocation result; step 3: generating a reference trajectory circle; and generating a visual spatial position reference target according to the flight segment delay allocation result and an operation deviation limit; and step 4: the aircraft is controlled by the time-space conversion method for taking off or landing.

2. The time-space conversion method of flight sequencing information according to claim 1, wherein the defining the variables in the step 2-1 comprises: SysTime: denoting a current time of a system; Flt.sub.i: denoting an i.sup.th flight; PtList.sub.i: denoting a flight trajectory point queue of the flight Flt.sub.i, wherein the queue is generated according to a flight plan of the flight and contains waypoints in a flight route and interpolation point information between waypoints; PtNum.sub.i: denoting a number of trajectory points in the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i; Pt.sub.i,j: denoting a j.sup.th trajectory point of the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i, wherein Pt.sub.i,j∈PtList.sub.i; ETO.sub.i,j: denoting an estimated passing time of the aircraft Flt.sub.i at Pt.sub.i,j; CTO.sub.i,j: denoting a sequencing passing time of the aircraft Flt.sub.i at Pt.sub.i,j, wherein an initial value is ETO.sub.i,j; Div(CTO.sub.i,j): denoting a difference of the sequencing passing time CTO.sub.i,j of the aircraft Flt.sub.i at Pt.sub.i,j allocated in two adjacent computation periods; DivLimit: denoting an upper limit of the difference of the sequencing passing time of the aircraft at the sequencing key point allocated in adjacent computation periods; Pt Pr o.sub.i,j: denoting a sequencing attribute of the point Pt.sub.i,j, wherein a value of 1 denotes that the point is a sequencing key point, and a value of 0 denotes others; PtDelay.sub.i,j: denoting an estimated passing time delay of the aircraft Flt.sub.i at Pt.sub.i,j; [Pt.sub.i,j, Pt.sub.i,k]: denoting a flight segment from the point Pt.sub.i,j to a point Pt.sub.i,k in PtList.sub.i, wherein Pt.sub.i,j denotes an origin of the flight segment, Pt.sub.i,k denotes an end point of the flight segment, and k denotes a serial number of the point Pt.sub.i,k in the queue PtList.sub.i; SegNum.sub.j,k: denoting a number of flight segments contained in the flight segment [Pt.sub.i,j, Pt.sub.i,k]; SegDis.sub.j,k: denoting a flight length of the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k]; SegMinSpeed.sub.j,k: denoting a minimum safe flight speed of the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k], wherein the parameter is acquired according to the aircraft performance or an airspace operation limit; SegMaxSpeed.sub.j,k: denoting a maximum safe flight speed of the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k], wherein the parameter is acquired according to the aircraft performance or the airspace operation limit; Seg Pr o.sub.j,k: denoting a flight status of the aircraft Flt.sub.i at the flight segment [Pt.sub.i,k], wherein a value of 1 denotes level flight, and a value of 0 denotes others; SegDelay.sub.j,k: denoting a delay value that needs to be consumed by the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k]; SegDelayTmp.sub.j,k: denoting a temporary variable of the delay that needs to be consumed by the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k]; SegMaxDelay.sub.j,k: denoting a maximum positive delay value capable of being consumed by the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k]; and SegMaxAcc.sub.j,k: denoting a maximum negative delay value capable of being consumed by the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k].

3. The time-space conversion method of flight sequencing information according to claim 1, wherein the dividing the flight segment according to the sequencing key point in the step 2-3 comprises: according to sequencing key point information in the flight trajectory point queue PtList.sub.i of the aircraft Flt.sub.i, dividing the delay-consumed flight segment [Pt.sub.i,bgn, Pt.sub.i,PtNum.sub.i] selected in the step 2-2 into a plurality of the flight segments with overlapped boundary points, wherein the method is as follows: step 2-3-1: setting a starting point of a sub-segment to be divided in the flight trajectory point queue PtList.sub.i of the flight, which is denoted with Pt.sub.i,m, and letting Pt.sub.i,m=Pt.sub.i,bgn, wherein m denotes a serial number of the point Pt.sub.i,m in the queue PtList.sub.i; step 2-3-2: searching for the closest sequencing key point after Pt.sub.i,m from PtList.sub.i, which is used as an end point of the sub-segment to be divided, and denoted with Pt.sub.i,n, wherein n denotes a serial number of the point Pt.sub.i,n in the queue PtList.sub.i, then the divided sub-flight segment is denoted with [Pt.sub.i,m, Pt.sub.i,n]; and step 2-3-3: letting Pt.sub.i,m=Pt.sub.i,n, and repeating the step 2-3-2 till the last point in PtList.sub.i; all the sub-flight segments divided in the step satisfy the following conditions: condition 1: [Pt.sub.i,m, Pt.sub.i,n]≠Ø condition 2: [Pt.sub.i,m, Pt.sub.i,n]⊂[Pt.sub.i,bgn, Pt.sub.i,PtNum.sub.i] and condition 3: ∪[Pt.sub.i,m, Pt.sub.i,n]=[Pt.sub.i,bgn, Pt.sub.i,PtNum.sub.i].

4. The time-space conversion method of flight sequencing information according to claim 3, wherein the method of dividing the flight segment according to the flight status in the step 2-4 comprises: subdividing each flight segment [Pt.sub.i,m, Pt.sub.i,n] divided in the step 2-3 according to the flight status of the aircraft in the flight segment, wherein the method comprises: dividing the flight segment [Pt.sub.i,m, Pt.sub.i,n] into the plurality of the flight segments [Pt.sub.i,o, Pt.sub.i,p] with overlapped boundary points according to the flight status of the aircraft, wherein o denotes a serial number of the point Pt.sub.i,o in the queue PtList.sub.i, p denotes a serial number of the point Pt.sub.i,p in the queue PtList.sub.i, and movement status of the flight Flt.sub.i in each small flight segment are consistent; and synchronously recording a flight status tag Seg Pr o.sub.o,p of the flight Flt.sub.i in each small flight segment [Pt.sub.i,o, Pt.sub.i,p]; and recording a number of the flight segments divided from the flight segment [Pt.sub.i,m, Pt.sub.i,n] in the step as SegNum.sub.m,n, wherein all the flight segments satisfy the following conditions: condition 1: [Pt.sub.i,o, Pt.sub.i,p]≠Ø condition 2: [Pt.sub.i,o, Pt.sub.i,p].Math.[Pt.sub.i,m, Pt.sub.i,n] and condition 3: U.sub.1.sup.SegNum.sup.m,n[Pt.sub.i,o, Pt.sub.i,p]=[Pt.sub.i,m, Pt.sub.i,n].

5. The time-space conversion method of flight sequencing information according to claim 4, wherein the method of allocating the flight segment delay in the step 2-5 comprises: allocating the sequencing delay of the aircraft at the sequencing key point to each flight segment according to the flight status and the flight segment length of the aircraft on each flight segment and the aircraft performance; performing the following steps in turn for each sub-flight segment [Pt.sub.i,m, Pt.sub.i,n] of the flight Flt.sub.i divided in the step 2-3: step 2-5-1: calculating a total delay of the flight segment: calculating the total delay that needs to be consumed by the flight Flt.sub.i at the flight segment [Pt.sub.i,m, Pt.sub.i,n] as:
SegDelay.sub.m,n=(CTO.sub.i,n−ETO.sub.i,n)−(CTO.sub.i,m−ETO.sub.i,m)  (5) step 2-5-2: calculating an upper limit of the flight segment delay: estimating the upper limit of the delay of the aircraft at the flight segment from a perspective of performance according to a level flight distance, a minimum level flight speed and a maximum level flight speed of the aircraft in the flight segment; referring to the step 2-4, and letting each small flight segment divided in the flight segment [Pt.sub.i,m, Pt.sub.i,n] be [Pt.sub.i,o, Pt.sub.i,p], then a maximum positive delay capable of being allocated by the flight Flt.sub.i at each small flight segment is: $\begin{array}{cc}{\mathrm{SegMaxDelay}}_{o,p}=\left(\frac{{\mathrm{SegDis}}_{o,p}}{{\mathrm{SegMinSpeed}}_{o,p}}-\left({\mathrm{ETO}}_{i,p}-{\mathrm{ETO}}_{i,0}\right)\right)*{\mathrm{SegPro}}_{o,p}& \left(6\right)\end{array}$ a maximum advance value capable of being allocated by the flight Flt.sub.i at the small flight segment [Pt.sub.i,o, Pt.sub.i,p] is: $\begin{array}{cc}{\mathrm{SegMaxAcc}}_{o,p}=\left(\left({\mathrm{ETO}}_{i,p}-{\mathrm{ETO}}_{i,0}\right)-\frac{{\mathrm{SegDis}}_{o,p}}{{\mathrm{SegMaxSpeed}}_{o,p}}\right)*{\mathrm{SegPro}}_{o,p}& \left(7\right)\end{array}$ then, the upper limit of the delay capable of being allocated by the flight Flt.sub.1 at the small flight segment [Pt.sub.i,m, Pt.sub.i,n] is:
SegMaxDelay.sub.m,n=Σ.sub.1.sup.SegNum.sup.m,nSegMaxDelay.sub.o,p  (8)
SegMaxAcc.sub.m,n=Σ.sub.1.sup.SegNum.sup.m,nSegMaxAcc.sub.o,p  (9) step 2-5-3: judging whether the flight segment is delayed: when any of the following conditions are satisfied, it is indicated that the delay of the flight Flt.sub.i at the flight segment [Pt.sub.i,m, Pt.sub.i,n] is excessive, prompting that the flight needs to be manually guided for delay absorption, and returning to step 2-5-1 to continuously processing next sub-flight segment; otherwise, indicating that the delay is within an acceptable range, and continuing subsequent steps; excessive flight segment delay conditions are as follows: condition 1: SegDelay.sub.m,n>0&&SegDelay.sub.m,n>SegMaxDelay.sub.m,n and condition 2: SegDelay.sub.m,n<0&&|SegDelay.sub.m,n|>SegMaxAcc.sub.m,n step 2-5-4: allocating the flight segment delay: allocating a delay for each level flight segment according to a distance and a delay absorption capacity of each level flight segment, wherein the method comprises: letting SegDelayTmp.sub.m,n be the delay needing to be received by the flight Flt.sub.i at the flight segment [Pt.sub.i,m, Pt.sub.i,n] in the step, and SegDelayTmp.sub.m,n=SegDelay.sub.m,n; and initializing the delay of each small flight segment [Pt.sub.i,o, Pt.sub.i,p] in the flight segment [Pt.sub.i,m, Pt.sub.i,n] and letting SegDelay.sub.o,p=0; and step 2-5-5: allocating a waypoint delay: calculating a sequencing time and a delay of the flight Flt.sub.i in a non-sequencing key point of the flight segment [Pt.sub.i,m, Pt.sub.i,n] according to the delay SegDelay.sub.o,p allocated for the flight Flt.sub.i at each small flight segment [Pt.sub.i,o, Pt.sub.i,p] in the flight segment [Pt.sub.i,m, Pt.sub.i,n] in the step 2-5-4, wherein the sequencing time and the delay of the sequencing key point are generated by the step 1; starting from an origin of the flight segment [Pt.sub.i,m, Pt.sub.i,n], for any point Pt.sub.i,q in the flight segment, q denoting a serial number of the point Pt.sub.i,q in the queue PtList.sub.i, positioning a small flight segment that the point belongs to, i.e. Pt.sub.i,q∈[Pt.sub.i,o, Pt.sub.i,p], then the delay needing to be allocated to the point is as follows: $\begin{array}{cc}{\mathrm{PrDelay}}_{i,q}=\{\begin{array}{cc}\begin{array}{c}\frac{{\mathrm{SegDis}}_{q-1,q}}{{\mathrm{SegDis}}_{o,p}}*{\mathrm{SegDelay}}_{o,p}+\\ {\mathrm{PtDelay}}_{i,q-1},\end{array}& \begin{array}{c}{\mathrm{Pt}}_{q-1}\in \left[{\mathrm{Pt}}_{i,o},{\mathrm{Pt}}_{i,p}\right]\&\&\\ {\mathrm{PtPro}}_{i,q}=0\end{array}\\ {\mathrm{CTO}}_{i,q}-{\mathrm{ETO}}_{i,q},& \mathrm{else}\end{array}& \left(13\right)\end{array}$ the sequencing time of the point Pt.sub.i,q is:
CTO.sub.i,q=ETO.sub.i,q+PtDelay.sub.i,q  (14).

6. The time-space conversion method of flight sequencing information according to claim 5, wherein the step 2-5-4 comprises the following steps of: step 2-5-4-1: allocating the flight segment delay: for each small flight segment [Pt.sub.i,o, Pt.sub.i,p] in the flight segment [Pt.sub.i,m, Pt.sub.i,n], the calculating method for allocating the delay according to the level flight distance of each small flight segment is as follows: $\begin{array}{cc}{\mathrm{SegDelay}}_{o,p}={\mathrm{SegDelay}}_{o,p}+{\mathrm{SegDelayTmp}}_{m,n}*\frac{{\mathrm{SegDis}}_{o,p}*{\mathrm{SegPro}}_{o,p}}{{.Math.}_1^{{\mathrm{SegNum}}_{m,n}}{\mathrm{SegDis}}_{o,p}*{\mathrm{SegPro}}_{o,p}}& \left(10\right)\end{array}$ correcting according to a flight segment delay absorption capacity, wherein the method is as follows: $\begin{array}{cc}{\mathrm{SegDelay}}_{o,p}=\{\begin{array}{cc}{\mathrm{SegMaxDelay}}_{o,p},& {\mathrm{SegDelay}}_{o,p}>{\mathrm{SegMaxDelay}}_{o,p}\\ -1*{\mathrm{SegMaxAcc}}_{o,p},& {\mathrm{SegDelay}}_{o,p}<\left(-1*{\mathrm{SegMaxAcc}}_{o,p}\right)\\ {\mathrm{SegDelay}}_{o,p},& \mathrm{else}\end{array}& \left(11\right)\end{array}$ step 2-5-4-2: detecting a delay allocation result: detecting whether any delay which is not allocated exists in the flight segment [Pt.sub.i,m, Pt.sub.i,n], wherein the method is as follows:
SegDelayTmp.sub.m,n=SegDelay.sub.m,n−Σ.sub.1.sup.SegNum.sup.m,nSegDelay.sub.o,p  (12) when satisfying that SegDelayTmp.sub.m,n=0, indicating that the delay of the flight segment [Pt.sub.i,m, Pt.sub.i,n] is completely allocated, and continuing step 2-5-5; otherwise, returning to step 2-5-4-1, and updating the flight segment delay allocation result.

7. The time-space conversion method of flight sequencing information according to claim 6, wherein the step 3 comprises the following steps of: step 3-1: defining variables; step 3-2: positioning the reference flight segment; step 3-3: generating a reference position point; and step 3-4: generating the reference trajectory circle; wherein the defining the variables in the step 3-1 comprises: Lat(Pt.sub.i,j): denoting a latitude of the flight Flt.sub.i at the waypoint Pt.sub.i,j; Lon(Pt.sub.i,j): denoting a longitude of the flight Flt.sub.i at the waypoint Pt.sub.i,j; Velo(Pt.sub.i,j): denoting a speed suggestion of the flight Flt.sub.i passing through the waypoint Pt.sub.i,j; Res Div: denoting a deviation limit between an actual passing time at the sequencing key point and the sequencing passing time of the flight; ATO.sub.i,j: denoting an actual passing time of the flight Flt.sub.i at the waypoint Pt.sub.i,j; and R: denoting a radius of the reference trajectory circle; step 3-2: positioning the reference flight segment: re-positioning the reference flight segment [Pt.sub.i,ref, Pt.sub.i,ref+1] at which the flight is Flt.sub.i currently located according to the sequencing passing time CTO.sub.i,j of each point Pt.sub.i,j in the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i calculated in the step 2-5-5, and satisfying that SysTime∈[CTO.sub.i,ref, CTO.sub.i,ref+1]; and step 3-3: generating the reference position point: approximately solving a real-time reference position point of the aircraft at the current time of the system by using a difference method, wherein the method comprises: letting Pt.sub.i,x be the current reference position point of the flight Flt.sub.i, wherein x denotes a serial number of a point Pt.sub.i,x in the queue PtList.sub.i, and letting the current reference flight segment at which the flight Flt.sub.i is located be [Pt.sub.i,ref, Pt.sub.i,ref+1] with reference to the step 3-2, then the method for calculating related information of the reference position point is as follows: $\begin{array}{cc}\mathrm{Lat}\left({\mathrm{Pt}}_{i,x}\right)=\mathrm{Lat}\left({\mathrm{Pt}}_{i,\mathrm{ref}}\right)+\frac{\mathrm{SysTime}-{\mathrm{CTO}}_{i,\mathrm{ref}}}{{\mathrm{CTO}}_{i,\mathrm{ref}+1}-{\mathrm{CTO}}_{i,\mathrm{ref}}}*\left(\mathrm{Lat}\left({\mathrm{Pt}}_{i,\mathrm{ref}+1}\right)-\mathrm{Lat}\left({\mathrm{Pt}}_{i,\mathrm{ref}}\right)\right)& \left(15\right)\end{array}$ $\begin{array}{cc}\mathrm{Lon}\left({\mathrm{Pt}}_{i,x}\right)=\mathrm{Lon}\left({\mathrm{Pt}}_{i,\mathrm{ref}}\right)+\frac{\mathrm{SysTime}-{\mathrm{CTO}}_{i,\mathrm{ref}}}{{\mathrm{CTO}}_{i,\mathrm{ref}+1}-{\mathrm{CTO}}_{i,\mathrm{ref}}}*\left(\mathrm{Lon}\left({\mathrm{Pt}}_{i,\mathrm{ref}+1}\right)-\mathrm{Lon}\left({\mathrm{Pt}}_{i,\mathrm{ref}}\right)\right)& \left(16\right)\end{array}$ $\begin{array}{cc}\mathrm{Velo}\left({\mathrm{Pt}}_{i,x}\right)=\frac{{\mathrm{SegDis}}_{\mathrm{ref},\mathrm{ref}+1}}{{\mathrm{CTO}}_{i,\mathrm{ref}+1}-{\mathrm{CTO}}_{i,\mathrm{ref}}}.& \left(17\right)\end{array}$

8. The time-space conversion method of flight sequencing information according to claim 7, wherein the step 3-4 comprises the following step of: when the actual passing time ATO.sub.i,j of the flight at the sequencing key point satisfies that:
ATO.sub.i,j∈[CTO.sub.i,j−Res Div,CTO.sub.i,j+Res Div] indicating that the flight runs according to the sequencing time; introducing the reference trajectory circle on the basis of the reference position point in the step 3-3, which is a circle taking the reference position point as a center point and R as a radius; wherein the calculation formula of the radius of the reference trajectory circle is:
R=Velo(Pt.sub.i,x)*Res Div  (18).

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The advantages of the above and/or other aspects of the present invention will become more apparent by further explaining the present invention with reference to the following drawings and detailed description.

(2) FIG. 1 is a general flow chart of the present invention.

(3) FIG. 2 is a schematic diagram showing sequencing key point information display of the present invention.

(4) FIG. 3 is a schematic diagram a flight segment delay allocation processing flow of the present invention.

(5) FIG. 4 is a schematic diagram of the first example of a delay-consumed flight segment filtering scenario of the present invention.

(6) FIG. 5 is a schematic diagram of the second example of the delay-consumed flight segment filtering scenario of the present invention.

(7) FIG. 6 is a schematic diagram of the third example of the delay-consumed flight segment filtering scenario of the present invention.

(8) FIG. 7 is a schematic diagram of the fourth example of the delay-consumed flight segment filtering scenario of the present invention.

(9) FIG. 8 is a schematic diagram of dividing a flight segment according to a sequencing key point of the present invention.

(10) FIG. 9 is a schematic diagram of dividing a flight segment according to a flight status of the present invention.

(11) FIG. 10 is a schematic diagram showing a processing flow of generating a reference trajectory circle of the present invention.

(12) FIG. 11 is an application example diagram of a reference trajectory circle of an aircraft of the present invention.

DETAILED DESCRIPTION

(13) The present invention is further explained with reference to the drawings and embodiments hereinafter.

(14) As shown in FIG. 1, a method of the present invention comprises the following steps of:

(15) step 1: sequencing and spacing a flight; and generating a sequencing time and a delay suggestion of the flight in each key point in a terminal area and each runway by using a flight sequencing technology;

(16) step 2: allocating a flight segment delay; and predicting a flight status of the flight according to a current position and 4D trajectory information of the flight, and filtering a delay-consumed flight segment, and on that basis, and generating a flight segment delay allocation strategy with reference to an aircraft performance and sequencing information; and

(17) step 3: generating a reference trajectory circle; and generating a visual spatial position reference target according to the flight segment delay allocation result and an operation deviation limit, which provides a visual reference for an air traffic controller to guide the aircraft according to the sequencing time.

(18) Step 1: sequencing and spacing the flight.

(19) The function of this step is that: a flight sequencing algorithm can comprehensively consider various operation limits in a terminal area, and generate sequencing time and delay suggestions for the flight in each key point in the terminal area and airport runway, so as to ensure the safe, sequential and efficient circulation of a traffic flow. For the specific methods, please refer to a previous patent “Multi-efficiency Sequencing and Spacing Method for Arrival Flights in Airport with Multi-runways”.

(20) Note 1: the key point in the terminal area usually refers to a corridor entrance point, an internal convergence point or an attention point designated by the user, and the like, and an aircraft passing through this point may be allocated with a sequencing time. In this method, this kind of point in the flight trajectory (i.e., a 4D trajectory of the flight) and the airport runway are called sequencing key points. As shown by the circled points in FIG. 2, the airplane in the figure may pass through three sequencing key points from the current position to airport landing.

(21) Note 2: the 4D trajectory information of the flight is generated by a 4D trajectory predicting technology, wherein the 4D trajectory predicting technology is a general technology used by the civil aviation industry to predict the flight trajectory.

(22) Step 2: allocating the flight segment delay.

(23) The function of this step is: to distribute the sequencing delay suggestions of the flight in each key point in the terminal area and each runway to each flight segment of the flight for absorption, so as to more accurately execute the sequencing time of the flight at the sequencing key point. Therefore, the processing flow of the predicting the flight status of the flight according to the current position and the 4D trajectory information of the flight, and filtering the delay-consumed flight segment, and on that basis, generating the flight segment delay allocation strategy with reference to the aircraft performance and the sequencing information is as shown in FIG. 3.

(24) The following steps are comprised:

(25) step 2-1: defining variables;

(26) step 2-2: filtering the delay-consumed flight segment;

(27) step 2-3: dividing the flight segment according to a sequencing key point;

(28) step 2-4: dividing the flight segment according to the flight status; and

(29) step 2-5: allocating the flight segment delay.

(30) Step 2-1: defining the variables:

(31) SysTime: denoting a current time of a system;

(32) Flt.sub.i: denoting an i.sup.th flight;

(33) PtList.sub.i: denoting a flight trajectory point queue of the flight Flt.sub.i, wherein the queue is generated according to a flight plan of the flight by using a 4D trajectory predicting technology and contains waypoints in a flight route and interpolation point information between waypoints;

(34) PtNum.sub.i: denoting a number of trajectory points in the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i;

(35) Pt.sub.i,j: denoting a j.sup.th trajectory point of the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i, wherein Pt.sub.i,j∈PtList.sub.i;

(36) ETO.sub.i,j: denoting an estimated passing time of the aircraft Flt.sub.i at Pt.sub.i,j and generated by the 4D trajectory predicting technology;

(37) CTO.sub.i,j: denoting a sequencing passing time of the aircraft Flt.sub.i at Pt.sub.i,j, wherein an initial value is ETO.sub.i,j; the sequencing time of the sequencing key point in the method is provided by the step 1; and sequencing time of other types of points is updated by the step;

(38) Div(CTO.sub.i,j): denoting a difference of the sequencing passing time CTO.sub.i,j of the aircraft Flt.sub.i at Pt.sub.i,j allocated in two adjacent computation periods, in a unit of second;

(39) DivLimit: denoting an upper limit of the difference of the sequencing passing time of the aircraft at the sequencing key point allocated in adjacent computation periods, in a unit of second, which is set as 60 seconds in the method, and may be set by a user according to the needs thereof;

(40) Pt Pr o.sub.i,j: denoting a sequencing attribute of the point Pt.sub.i,j, wherein a value of 1 denotes that the point is a sequencing key point, and a value of 0 denotes others;

(41) PtDelay.sub.i,j: denoting an estimated passing time delay of the aircraft Flt.sub.i at Pt.sub.i,j, in a unit of second;

(42) [Pt.sub.i,j, Pt.sub.i,k]: denoting a flight segment from the point Pt.sub.i,j to a point Pt.sub.i,k in PtList.sub.i, wherein Pt.sub.i,j denotes an origin of the flight segment, and Pt.sub.i,k denotes an end point of the flight segment;

(43) SegNum.sub.j,k: denoting a number of small flight segments contained in the flight segment [Pt.sub.i,j, Pt.sub.i,k];

(44) SegDis.sub.j,k: denoting a flight length of the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k];

(45) SegMinSpeed.sub.j,k: denoting a minimum safe flight speed of the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k], wherein the parameter is acquired according to the aircraft performance or an airspace operation limit;

(46) SegMaxSpeed.sub.j,k: denoting a maximum safe flight speed of the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k], wherein the parameter is acquired according to the aircraft performance or the airspace operation limit;

(47) Seg Pr o.sub.j,k: denoting a flight status of the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k], wherein a value of 1 denotes level flight, and a value of 0 denotes others;

(48) SegDelay.sub.j,k: denoting a delay value that needs to be consumed by the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k], in a unit of second;

(49) SegDelayTmp.sub.j,k: denoting a temporary variable of the delay that needs to be consumed by the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k] used in the calculating process of the method, and in a unit of second;

(50) SegMaxDelay.sub.j,k: denoting a maximum positive delay value capable of being consumed by the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k], in a unit of second; and

(51) SegMaxAcc.sub.j,k: denoting a maximum negative delay value capable of being consumed by the aircraft Flt.sub.i at the flight segment [Pt.sub.i,j, Pt.sub.i,k], in a unit of second.

(52) Step 2-2: filtering the delay-consumed flight segment:

(53) In order to improve a stability of the reference trajectory circle generated by this method in periodic calculation (to reduce the problem that users cannot use the reference trajectory circle normally due to frequent hopping in the position of the reference trajectory circle changing with the time), and meanwhile, take into account the need of timely correction when there is a big deviation between the reference trajectory circle and the actual position of the aircraft, or when the sequencing time of the sequencing key point fluctuates greatly, the following ways are adopted in this method to filter the delay-consumed flight segment. The concept of the reference trajectory circle is described in the step 3.

(54) Letting Pt.sub.i,bgn be an origin of the flight segments participating in delay consumption in the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i.

(55) Step 2-2-1: positioning an actual flight segment at which the flight is currently located:

(56) positioning the actual flight segment [Pt.sub.i,cur, Pt.sub.i,cur+1] at which the flight Flt.sub.i is currently located according to the estimated passing time ETO.sub.i,j of each point Pt.sub.i,j in the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i, and satisfying that SysTime∈[ETO.sub.i,cur, ETO.sub.i,cur+1].

(57) Step 2-2-2: positioning a reference flight segment at which the flight is currently located:

(58) positioning the reference flight segment [Pt.sub.i,ref, Pt.sub.i,ref+1] at which the flight Flt.sub.i is currently located according to the sequencing passing time CTO.sub.i,j of each point Pt.sub.i,j in the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i, and satisfying that SysTime∈[CTO.sub.i,ref, CTO.sub.i,ref+1];

(59) Step 2-2-5: searching for an adjacent sequencing key point in a preceding flight segment:

(60) letting Pt.sub.i,pre be the sequencing key point closest to a current position Pt.sub.i,cur in the preceding flight segment [Pt.sub.i,1, Pt.sub.i,cur] of the flight Flt.sub.i;

(61) letting Pt.sub.i,tmp be an intermediate variable of flight trajectory points of the flight Flt.sub.i in the calculation process of the method;

(62) when satisfying that ∃Pt.sub.i,tmp∈[Pt.sub.i,1, Pt.sub.i,cur], such that
Min{(ETO.sub.i,cur−ETO.sub.i,tmp+1)*PtPro.sub.i,tmp}>0  (1)

(63) letting Pt.sub.i,pre=Pt.sub.i,tmp, and continuously executing step 2-2-4; otherwise, letting Pt.sub.i,pre=Ø, and skipping to step 2-2-5.

(64) step 2-2-4: judging whether a sequencing key point exists between the actual position of the flight and the reference flight segment:

(65) when satisfying that ETO.sub.i,cur>ETO.sub.i,ref and Pt.sub.i,pre∈[Pt.sub.i,ref+1, Pt.sub.i,cur], letting Pt.sub.i,bgn=Pt.sub.i,pre, and skipping to step 2-2-7; otherwise, continuously executing step 2-2-5;

(66) Step 2-2-5: searching for an adjacent sequencing key point in a subsequent flight segment:

(67) letting Pt.sub.i,PtNum.sub.i be the last point in the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i

(68) letting Pt.sub.i,aft be the sequencing key point closest to the current position Pt.sub.i,cur in the subsequent flight segment [Pt.sub.i,cur+1, Pt.sub.i,PtNum.sub.i] of the flight Flt.sub.i;

(69) when satisfying that ∃Pt.sub.i,tmp∈[Pt.sub.i,cur+1, Pt.sub.i,PtNum.sub.i], such that
Min{(ETO.sub.i,tmp−ETO.sub.i,cur)*PtPro.sub.i,tmp}>0  (2)

(70) letting Pt.sub.i,aft=Pt.sub.i,tmp; otherwise, letting Pt.sub.i,aft=Pt.sub.i,PtNum.sub.i, wherein Pt.sub.i,PtNum.sub.i is the last point in PtList.sub.i, i.e., an airport runway, and the airport runway is also set as the sequencing key point in the method.

(71) Step 2-2-6: judging whether the sequencing time of the subsequent sequencing key point has a larger fluctuation:

(72) calculating a difference Div(CTO.sub.i,aft) between the currently allocated sequencing time and the allocated sequencing time during last calculation of the flight Flt.sub.i at the sequencing key point Pt.sub.i,aft, and filtering the origin Pt.sub.i,bgn of the delay-consumed flight segment according to Div(CTO.sub.i,aft); for the first operation, letting Div(CTO.sub.i,aft)=0:

(73) $\begin{array}{cc}{\mathrm{Pt}}_{i,\mathrm{bgn}}=\{\begin{array}{cc}{\mathrm{Pt}}_{i,\mathrm{ref}},& .Math.\mathrm{Div}\left({\mathrm{CTO}}_{i,\mathrm{aft}}\right).Math.<\mathrm{DivLimit}\\ {\mathrm{Pt}}_{i,\mathrm{cur}},& .Math.\mathrm{Div}\left({\mathrm{CTO}}_{i,\mathrm{aft}}\right).Math.\ge \mathrm{DivLimit}\&\&{\mathrm{Pt}}_{i,\mathrm{pre}}==\varnothing \\ {\mathrm{Pt}}_{i,\mathrm{pre}},& \mathrm{else}\end{array}& \left(3\right)\end{array}$

(74) Step 2-2-7: determining the delay-consumed flight segment:

(75) for the flight Flt.sub.i, the flight segments participating in delay consumption filtered in the method are all the subsequent flight segment from Pt.sub.i,bgn in the flight trajectory point queue PtList.sub.i, which is [Pt.sub.i,bgn, Pt.sub.i,PtNum.sub.i].

(76) Step 2-2-8: updating a sequencing time of the origin of the delay-consumed flight segment:

(77) 0$\begin{array}{cc}{\mathrm{CTO}}_{i,\mathrm{bgn}}=\{\begin{array}{cc}{\mathrm{ETO}}_{i,\mathrm{bgn}},& {\mathrm{Pt}}_{i,\mathrm{bgn}}\ne {\mathrm{Pt}}_{i,\mathrm{ref}}\&\&{\mathrm{ETO}}_{i,\mathrm{bgn}}\le \mathrm{SysTime}\\ {\mathrm{CTO}}_{i,\mathrm{bgn}},& \mathrm{else}\end{array}& \left(4\right)\end{array}$

(78) As shown in FIG. 4, the origin of the reference flight segment of the flight in the figure is Pt.sub.i,ref, and the end point of the reference flight segment is Pt.sub.i,ref+1, the origin of the actual flight segment is Pt.sub.i,cur, and the end point of the actual flight segment is Pt.sub.i,cur+1. The actual flight segment of the aircraft in the figure is ahead of (that is, closer to a destination airport) the reference flight segment, and the sequencing key point Pt.sub.i,pre exists between the origin Pt.sub.i,cur of the actual flight segment of the aircraft and the end point Pt.sub.i,ref+1 of the reference flight segment, then the sequencing key point is used as the origin Pt.sub.i,bgn of the delay-consumed flight segment.

(79) As shown in FIG. 5, if the sequencing key point does not exist between the actual flight segment of the aircraft in the figure and the reference flight segment, and the subsequent adjacent sequencing key point Pt.sub.i,aft satisfies that |Div(CTO.sub.i,aft)|<DivLimit, then the origin Pt.sub.i,ref of the reference flight segment in FIG. 5 is used as the origin Pt.sub.i,bgn of the delay-consumed flight segment.

(80) As shown in FIG. 6, the sequencing key point does not exist in the preceding flight segment of the current position of the aircraft, i.e., Pt.sub.i,pre==Ø is satisfied, and the subsequent adjacent sequencing key point Pt.sub.i,aft satisfies that |Div(CTO.sub.i,aft)|≥DivLimit, then the origin Pt.sub.i,cur of the actual flight segment in FIG. 6 is used as the origin Pt.sub.i,bgn of the delay-consumed flight segment.

(81) As shown in FIG. 7, the adjacent sequencing key point Pt.sub.i,pre exists in the preceding flight segment of the current position of the aircraft in the figure, and the subsequent adjacent sequencing key point Pt.sub.i,aft satisfies that |Div(CTO.sub.i,aft)|≥DivLimit, then the sequencing key point Pt.sub.i,pre in FIG. 7 is used as the origin Pt.sub.i,bgn of the delay-consumed flight segment.

(82) Note 3: for the first calculation, letting a sequencing time of all the non-sequencing key points Pt.sub.i,j in PtList.sub.i be CTO.sub.i,j=ETO.sub.i,j; otherwise, using the last operation result for CTO.sub.i,j of all the non-sequencing key points herein.

(83) Step 2-3: dividing the flight segment according to a sequencing key point:

(84) according to sequencing key point information in the flight trajectory point queue PtList.sub.i of the aircraft Flt.sub.i, dividing the delay-consumed flight segment [Pt.sub.i,bgn, Pt.sub.i,PtNum.sub.i] selected in the step 2-2 into a plurality of small flight segments with overlapped boundary points.

(85) The specific process is as follows:

(86) 1): setting a starting point of a sub-segment to be divided in the flight trajectory point queue PtList.sub.i of the flight;

(87) 2): searching for the closest sequencing key point after Pt.sub.i,m from PtList.sub.i, which is used as an end point of the sub-segment to be divided, and denoted with Pt.sub.i,n, then the divided sub-flight segment is [Pt.sub.i,m, Pt.sub.i,n]; and

(88) 3): letting Pt.sub.i,m=Pt.sub.i,n, and repeating the step 2-3-2 till the last point in PtList.sub.i.

(89) Then, all the sub-flight segments divided in the step need to satisfy the following conditions: 1) [Pt.sub.i,m, Pt.sub.i,n]≠Ø 2) [Pt.sub.i,m, Pt.sub.i,n]⊂[Pt.sub.i,bgn, Pt.sub.i,PtNum.sub.i] and 3) ∪[Pt.sub.i,m, Pt.sub.i,n]=[Pt.sub.i,bgn, Pt.sub.i,PtNum.sub.i].

(90) Taking FIG. 8 for example, the delay-consumed flight segment of the flight Flt.sub.i is divided into three segments according to the key point information participating in sequencing, which are respectively denoted with Seg1, Seg2, and Seg3.

(91) Step 2-4: dividing the flight segment according to the flight status:

(92) In practice, in view of the operation complexity and safety of the aircraft in a maneuvering stage such as climbing or descending, this method allocates the delay absorption process to the level flight stage of the aircraft. Therefore, each flight segment [Pt.sub.i,m, Pt.sub.i,n] divided in the step 2-3 of the method is subdivided according to the flight status of the aircraft in the flight segment, wherein the method specific is as follows:

(93) dividing the flight segment [Pt.sub.i,m, Pt.sub.i,n] into a plurality of small flight segments [Pt.sub.i,o, Pt.sub.i,p] with overlapped boundary points according to the flight status of the aircraft, wherein movement status of the flight Flt.sub.i in each small flight segment are consistent; and synchronously recording a flight status tag Seg Pr o.sub.o,p of the flight Flt.sub.i in each small flight segment [Pt.sub.i,o, Pt.sub.i,p]; and

(94) recording a number of the small flight segments divided from the flight segment [Pt.sub.i,m, Pt.sub.i,n] in the step as SegNum.sub.m,n, wherein all the small flight segments satisfy the following conditions: 1) [Pt.sub.i,o, Pt.sub.i,p]≠Ø 2) [Pt.sub.i,o, Pt.sub.i,p]⊂[Pt.sub.i,m, Pt.sub.i,n] and 3) ∪.sub.1.sup.SegNum.sup.m,n[Pt.sub.i,o, Pt.sub.i,p]=[Pt.sub.i,m, Pt.sub.i,n].

(95) As shown in FIG. 9, this figure shows a vertical section of the flight trajectory of the aircraft in a certain flight segment. According to the movement status of the aircraft in the flight segment, the flight segment is divided into three segments, wherein two level flight segments X1 and X3 are denoted with solid lines and one descending flight segment X2 is denoted with dotted lines. This method only selects the level flight segment in the figure for delay allocation.

(96) Step 2-5: allocating the flight segment delay:

(97) In this step, the sequencing delay of the aircraft at the sequencing key point is allocated to each flight segment according to the flight status and the flight segment length of the aircraft on each flight segment and the aircraft performance, which is convenient for the delay absorption process to be safer and more efficient.

(98) The following processing is performed in turn for each sub-flight segment [Pt.sub.i,m, Pt.sub.i,n] of the flight Flt.sub.i divided in the step 2-3.

(99) Step 2-5-1: calculating a total delay of the flight segment:

(100) calculating the total delay that needs to be consumed by the flight Flt.sub.i at the flight segment [Pt.sub.i,m, Pt.sub.i,n] as:
SegDelay.sub.m,n=(CTO.sub.i,n−ETO.sub.i,n)−(CTO.sub.i,m−ETO.sub.i,m)  (5)

(101) Step 2-5-2: calculating an upper limit of the flight segment delay:

(102) The method estimates the upper limit of the delay of the aircraft at the flight segment from a perspective of performance according to a level flight distance, a minimum level flight speed and a maximum level flight speed of the aircraft in the flight segment, and a user may adjust as needed.

(103) Referring to the step 2-4, and letting each small flight segment divided in the flight segment [Pt.sub.i,m, Pt.sub.i,n] be [Pt.sub.i,o, Pt.sub.i,p], then a maximum positive delay capable of being allocated by the flight Flt.sub.i at each small flight segment is:

(104) $\begin{array}{cc}{\mathrm{SegMaxDelay}}_{o,p}=\left(\frac{{\mathrm{SegDis}}_{o,p}}{{\mathrm{SegMinSpeed}}_{o,p}}-\left({\mathrm{ETO}}_{i,p}-{\mathrm{ETO}}_{i,0}\right)\right)*{\mathrm{SegPro}}_{o,p}& \left(6\right)\end{array}$

(105) a maximum advance value capable of being allocated by the flight Flt.sub.i at the small flight segment [Pt.sub.i,o, Pt.sub.i,p] is:

(106) $\begin{array}{cc}{\mathrm{SegMaxAcc}}_{o,p}=\left(\left({\mathrm{ETO}}_{i,p}-{\mathrm{ETO}}_{i,0}\right)-\frac{{\mathrm{SegDis}}_{o,p}}{{\mathrm{SegMaxSpeed}}_{o,p}}\right)*{\mathrm{SegPro}}_{o,p}& \left(7\right)\end{array}$

(107) then, the upper limit of the delay capable of being allocated by the flight Flt.sub.i at the small flight segment [Pt.sub.i,m, Pt.sub.i,n] is:
SegMaxDelay.sub.m,n=Σ.sub.1.sup.SegNum.sup.m,nSegMaxDelay.sub.o,p  (8)
SegMaxAcc.sub.m,n=Σ.sub.1.sup.SegNum.sup.m,nSegMaxAcc.sub.o,p  (9)

(108) Step 2-5-3: judging whether the flight segment delay is excessive:

(109) When any of the following conditions are satisfied, it is indicated that the delay of the flight Fits at the flight segment [Pt.sub.i,m, Pt.sub.i,n] is excessive, needing to prompt an air traffic controller that the flight needs to be manually guided for delay absorption, and returning to step 2-5-1 to continuously processing next sub-flight segment. Otherwise, it is indicated that the delay is within an acceptable range, and subsequent steps are continued.

(110) Excessive flight segment delay conditions are as follows: 1) SegDelay.sub.m,n>0&&SegDelay.sub.m,n>SegMaxDelay.sub.m,n and 2) SegDelay.sub.m,n<0&&|SegDelay.sub.m,n|>SegMaxAcc.sub.m,n.

(111) Step 2-5-4: allocating the flight segment delay:

(112) In order to reduce instructions issued by the air traffic controller and actions taken by a pilot in the delay absorption process, the flight segment delay allocation result should be as uniform as possible. This method allocates a delay for each level flight segment according to a distance and a delay absorption capacity of each level flight segment.

(113) Letting SegDelayTmp.sub.m,n be the delay needing to be received by the flight Flt.sub.i at the flight segment [Pt.sub.i,m, Pt.sub.i,n] in the step, and SegDelayTmp.sub.m,n=SegDelay.sub.m,n; and

(114) initializing the delay of each small flight segment [Pt.sub.i,o, Pt.sub.i,p] in the flight segment [Pt.sub.i,m, Pt.sub.i,n], and letting SegDelay.sub.o,p=0.

(115) Step 2-5-4-1: allocating the flight segment delay:

(116) for each small flight segment [Pt.sub.i,o, Pt.sub.i,p] in the flight segment [Pt.sub.i,j, Pt.sub.i,n], the calculating method for allocating the delay according to the level flight distance of each small flight segment is as follows:

(117) $\begin{array}{cc}{\mathrm{SegDelay}}_{o,p}={\mathrm{SegDelay}}_{o,p}+{\mathrm{SegDelayTmp}}_{m,n}*\frac{{\mathrm{SegDis}}_{o,p}*{\mathrm{SegPro}}_{o,p}}{{.Math.}_1^{{\mathrm{SegNum}}_{m,n}}{\mathrm{SegDis}}_{o,p}*{\mathrm{SegPro}}_{o,p}}& \left(10\right)\end{array}$

(118) In order to ensure the feasibility of the flight segment delay allocation result, correcting is made according to a flight segment delay absorption capacity, wherein the method is as follows:

(119) $\begin{array}{cc}{\mathrm{SegDelay}}_{o,p}=\{\begin{array}{cc}{\mathrm{SegMaxDelay}}_{o,p},& {\mathrm{SegDelay}}_{o,p}>{\mathrm{SegMaxDelay}}_{o,p}\\ -1*{\mathrm{SegMaxAcc}}_{o,p},& {\mathrm{SegDelay}}_{o,p}<\left(-1*{\mathrm{SegMaxAcc}}_{o,p}\right)\\ {\mathrm{SegDelay}}_{o,p},& \mathrm{else}\end{array}& \left(11\right)\end{array}$

(120) Step 2-5-4-2: detecting a delay allocation result:

(121) detecting whether any delay which is not allocated exists in the flight segment [Pt.sub.i,m, Pt.sub.i,n], wherein the method is as follows:
SegDelayTmp.sub.m,n=SegDelay.sub.m,n−Σ.sub.1.sup.SegNum.sup.m,nSegDelay.sub.o,p  (12)

(122) when satisfying that SegDelayTmp.sub.m,n=0, indicating that the delay of the flight segment [Pt.sub.i,m, Pt.sub.i,n] is completely allocated, and continuing step 2-5-5; otherwise, returning to step 2-5-4-1, and updating the flight segment delay allocation result.

(123) Step 2-5-5: allocating a waypoint delay:

(124) calculating a sequencing time and a delay of the flight Flt.sub.i in a non-sequencing key point of the flight segment [Pt.sub.i,m, Pt.sub.i,n] according to the delay SegDelay.sub.o,p allocated for the flight Flt.sub.i at each small flight segment [Pt.sub.i,o, Pt.sub.i,p] in the flight segment [Pt.sub.i,m, Pt.sub.i,n] in the step 2-5-4, wherein the sequencing time and the delay of the sequencing key point are generated by the step 1; and

(125) starting from an origin of the flight segment [Pt.sub.i,m, Pt.sub.i,n], for any point Pt.sub.i,q in the flight segment, positioning a small flight segment that the point belongs to, i.e. Pt.sub.i,q∈[Pt.sub.i,o, Pt.sub.i,p], then the delay needing to be allocated to the point is as follows:

(126) $\begin{array}{cc}{\mathrm{PtDelay}}_{i,q}=\{\begin{array}{cc}\begin{array}{c}\frac{{\mathrm{SegDis}}_{q-1,q}}{{\mathrm{SegDis}}_{o,p}}*{\mathrm{SegDelay}}_{o,p}+\\ {\mathrm{PtDelay}}_{i,q-1},\end{array}& \begin{array}{c}{\mathrm{Pt}}_{q-1}\in \left[{\mathrm{Pt}}_{i,o},{\mathrm{Pt}}_{i,p}\right]\&\&\\ {\mathrm{PtPro}}_{i,q}=0\end{array}\\ {\mathrm{CTO}}_{i,q}-{\mathrm{ETO}}_{i,q},& \mathrm{else}\end{array}& \left(13\right)\end{array}$

(127) the sequencing time of the point Pt.sub.i,q is:
CTO.sub.i,q=ETO.sub.i,q+PtDelay.sub.i,q  (14).

(128) Step 3: generating a reference trajectory circle:

(129) The function of this step is: to generate a visual spatial position reference target according to the flight segment delay allocation result and an operation deviation limit, solve an abstract problem of the flight sequencing time, and enhance a time-based operation ability of the air traffic controller. The processing flow of generating the reference trajectory circle is shown in FIG. 10.

(130) The following steps are comprised:

(131) step 3-1: defining variables;

(132) step 3-2: positioning the reference flight segment;

(133) step 3-3: generating a reference position point; and

(134) step 3-4: generating the reference trajectory circle;

(135) Step 3-1: defining the variables:

(136) Lat(Pt.sub.i,j): denoting a latitude of the flight Flt.sub.i at the waypoint Pt.sub.i,j;

(137) Lon(Pt.sub.i,j): denoting a longitude of the flight Flt.sub.i at the waypoint Pt.sub.i,j;

(138) Velo(Pt.sub.i,j): denoting a speed suggestion of the flight Flt.sub.i passing through the waypoint Pt.sub.i,j;

(139) Res Div: denoting a deviation limit between an actual passing time at the sequencing key point and the sequencing passing time of the flight, which may be set by the user as needed;

(140) ATO.sub.i,j: denoting an actual passing time of the flight Flt.sub.i at the waypoint Pt.sub.i,j; and

(141) R: denoting a radius of the reference trajectory circle.

(142) Step 3-2: positioning the reference flight segment:

(143) re-positioning the reference flight segment [Pt.sub.i,ref, Pt.sub.i,ref+1] at which the flight is Flt.sub.i currently located according to the sequencing passing time CTO.sub.i,j of each point Pt.sub.i,j in the flight trajectory point queue PtList.sub.i of the flight Flt.sub.i calculated in the step 2-5-5, and satisfying that SysTime∈[CTO.sub.i,ref, CTO.sub.i,ref+1].

(144) Step 3-3: generating the reference position point:

(145) An interval between two adjacent points in the flight trajectory point queue PtList.sub.i of the flight generated by 4D trajectory predicting technology is usually controlled at a second level (for example, 8 seconds). Although the interval is not large, it is not enough to support a real-time update demand. On the basis of this reference flight segment, this method uses a difference method to approximately solve the real-time reference position point of the aircraft at the current time of the system, so as to meet practical industrial application requirements.

(146) Letting Pt be the current reference position point of the flight Flt.sub.i, and letting the current reference flight segment at which the flight Flt.sub.i is located be [Pt.sub.i,ref, Pt.sub.i,ref+1] with reference to the step 3-2, then the method for calculating related information of the reference position point is as follows:

(147) $\begin{array}{cc}\mathrm{Lat}\left({\mathrm{Pt}}_{i,x}\right)=\mathrm{Lat}\left({\mathrm{Pt}}_{i,\mathrm{ref}}\right)+\frac{\mathrm{SysTime}-{\mathrm{CTO}}_{i,\mathrm{ref}}}{{\mathrm{CTO}}_{i,\mathrm{ref}+1}-{\mathrm{CTO}}_{i,\mathrm{ref}}}*\left(\mathrm{Lat}\left({\mathrm{Pt}}_{i,\mathrm{ref}+1}\right)-\mathrm{Lat}\left({\mathrm{Pr}}_{i,\mathrm{ref}}\right)\right)& \left(15\right)\end{array}$$\begin{array}{cc}\mathrm{Lon}\left({\mathrm{Pt}}_{i,x}\right)=\mathrm{Lon}\left({\mathrm{Pt}}_{i,\mathrm{ref}}\right)+\frac{\mathrm{SysTime}-{\mathrm{CTO}}_{i,\mathrm{ref}}}{{\mathrm{CTO}}_{i,\mathrm{ref}+1}-{\mathrm{CTO}}_{t,\mathrm{ref}}}*\left(\mathrm{Lon}\left({\mathrm{Pt}}_{i,\mathrm{ref}+1}\right)-\mathrm{Lon}\left({\mathrm{Pt}}_{i,\mathrm{ref}}\right)\right)& \left(16\right)\end{array}$$\begin{array}{cc}\mathrm{Velo}\left({\mathrm{Pt}}_{i,x}\right)=\frac{{\mathrm{SegDis}}_{\mathrm{ref},\mathrm{ref}+1}}{{\mathrm{CTO}}_{i,\mathrm{ref}+1}-{\mathrm{CTO}}_{i,\mathrm{ref}}}.& \left(17\right)\end{array}$

(148) Step 3-4: generating the reference trajectory circle:

(149) From the perspective of practical application, it is usually very difficult for the air traffic controller to guide the aircraft to pass through the sequencing key point precisely according to the sequencing time of the sequencing key point generated in the step 1. Therefore, a deviation limit Res Div (for example, 10 seconds) is usually set manually at the sequencing key point. As long as the actual passing time ATO.sub.i,j of the flight at the sequencing key point satisfies that ATO.sub.i,j∈[CTO.sub.i,j−Res Div,CTO.sub.i,j+Res Div], it means that the flight is running normally according to the sequencing time.

(150) Therefore, this method introduces the concept of the reference trajectory circle based on the reference position point in the step 3-3, that is, a circle with the reference position point as a center point and R as a radius. As shown in FIG. 11, an aircraft icon in the figure indicates a current actual position of flight CSC9376, while a circle with a flag of the flight number in the figure indicates a reference trajectory circle of the flight. The reference trajectory circle of the flight in the figure lags behind the actual position of the flight, indicating that the flight is flying too fast at present and needs to be adjusted according to the sequencing time and deviation limit of the flight at the key sequencing point.

(151) The calculation formula of the radius of the reference trajectory circle is:
R=Velo(Pt.sub.i,x)*Res Div  (18).

(152) In actual operation, if the air traffic controller can guide the aircraft into the corresponding reference trajectory circle, the airplane can be guided to pass through the sequencing key point on the premise of satisfying the deviation limit. This method can provide a time-based planning and guidance ability for air traffic controller, and provide technical support for future TBO operation.

(153) The aircraft (airplane) is planned and guided according to the result of the reference trajectory circle generated in the step 3.

(154) The time-space conversion method of flight sequencing information of this embodiment is loaded and operated in a processing server of an air traffic control system (ATC system).

(155) In a specific implementation, the present application provides a computer storage medium and a corresponding data processing unit, wherein the computer storage medium is capable of storing a computer program, and the computer program, when executed by the data processing unit, can run the inventive contents of the time-space conversion method of flight sequencing information provided by the present invention and some or all steps in various embodiments. The storage medium may be a magnetic disk, an optical disk, a Read Only Storage (ROM) or a Random Access Storage (RAM), and the like.

(156) Those skilled in the art can clearly understand that the technical solutions in the embodiments of the present invention can be realized by means of a computer program and a corresponding general hardware platform thereof. Based on such understanding, the essence of the technical solutions in the embodiments of the present invention or the part contributing to the prior art, may be embodied in the form of a computer program, i.e., a software product. The computer program, i.e., the software product is stored in a storage medium comprising a number of instructions such that a device (which may be a personal computer, a server, a singlechip, a MUU or a network device, and the like) comprising the data processing unit executes the methods described in various embodiments or some parts of the embodiments of the present invention.

(157) The present invention provides an idea and a process for the time-space conversion method of flight sequencing information. There are many methods and ways to realize the technical solutions. The above is only the preferred embodiments of the present invention. It should be pointed out that those of ordinary skills in the art can make some improvements and embellishments without departing from the principle of the present invention, and these improvements and embellishments should also be regarded as falling with the scope of protection of the present invention. All the unspecified components in the embodiments can be realized by the prior art.