Sage Journals: Discover world-class research

Abstract

At present, airspace congestion and flight delays have become widespread concerns. This study aims to optimize the sequencing of arrival flights in the terminal area of multirunway airports. Considering the constraints of multiple runways, slant intervals and moving flight positions, this article establishes an optimization model for arrival flight sequencing in a multirunway airport terminal area. Accordingly, an improved sparrow search algorithm (ISSA) is proposed based on Chebyshev chaotic mapping, the golden sine strategy, and the variable neighborhood strategy. Through six basic test functions, the ISSA is compared with particle swarm optimization, the whale optimization algorithm, the genetic algorithm, and other algorithms to verify its superiority. Finally, two sets of instance data from Kunming Changshui Airport were used for experiments. The results show that the total delay times (TDTs) of small-scale flights (number of aircraft: 29) and large-scale flights (number of aircraft: 147) are 55.3% and 20.5% lower, respectively, than those of the first-come-first-served algorithm. The superiority of the ISSA designed in this article is verified, and it can significantly reduce the TDTs of arrival flights. It is suitable for optimizing arrival flights during peak hours at most airports. This approach provides theoretical support for optimizing the sorting of flights in terminal areas.

Keywords

Flight sequence optimization algorithm sparrow algorithm air traffic management

Introduction

Considering the post-COVID-19 era, China's civil aviation transportation market experienced explosive growth in 2023. Subsequently, there were many challenges in air traffic flow, such as the low utilization rate of runways in many busy airports and the poor efficiency of air traffic control systems. Therefore, effectively and efficiently scheduling and sorting arriving flights is one of the main problems in improving airport safety, capacity, and efficiency.

Many operational studies and management-related methods have been applied to solve the modeling problem of arrival flight sequencing and scheduling. With respect to flight optimization models, Wieland et al.¹ transformed the multiairport flight scheduling problem into a single-airport problem. Then, multiple single airports were combined to form a multiairport scheduling module. Wang et al.² established a multiairport collaborative decision-making arrival and departure flight scheduling model for the terminal area based on a multiconstraint situation. Zhou et al.³ established an optimal scheduling model for approach and departure flights with maximum throughput, minimum flight delay loss and robustness as the objective functions. Tian et al.⁴ constructed a multiobjective terminal area resource allocation problem by using the arrival fix allocation problem, taking into account the environmental constraints related to emissions and noise. Jiang et al.⁵ proposed a comprehensive mixed-integer programming model for arrival sequencing that considered multiple airports, runways, and routes.

There are many current research results on approach sequencing, but they still have shortcomings, such as the following: (1) model construction methods seldom consider the runway configuration and the current operation status, and this is mostly done for single-runway flight optimization; (2) the tail flow interval standard under RECAT-CN is rarely considered, so the needs of actual operation cannot be met; (3) most scholars adopt the first-come-first-served (FCFS) algorithm for flight scheduling, which is simple and easy to implement but is also prone to causing a large number of flight delays for a busy airport or an uneven distribution of multiple runways; and (4) most studies have not considered the fact that in the flight sequencing problem, the delays brought about by different aircraft types and different flight modes are also different.

In terms of solution algorithms, Bae and Clarke⁶ introduced the delay saturation weighting factor and proposed an improved mixed integer planning algorithm. In addition, the scatter search algorithm,⁷ branch and bound method,⁸ genetic algorithm (GA),^4.9,10,11 greedy algorithm,¹² and simulated annealing algorithm¹³ have been used to solve the approach scheduling problem. Most scholars use traditional optimization algorithms, but developing a more effective and accurate solution algorithm in addition to the commonly used traditional algorithms is an urgent problem; for example, the swarm intelligence algorithm has a strong global search ability and high convergence speed, so the introduction of swarm intelligence algorithms into the study of the flight scheduling problem may be a helpful tactic. The sparrow search algorithm has the advantages of a global search and few adjustable parameters, and it has been applied by some scholars in image processing,¹⁴ path planning¹⁵ and other practical engineering optimization problems. However, the SSA has the problems of insufficient search accuracy and easily falling into a local optimum solution. In this regard, Lv et al.¹⁶ used tent chaotic mapping to perturb the population to reduce local searches. Tang et al.¹⁷ drew on the idea of GAs and introduced crossover and mutation into the SSA to protect the diversity of populations while simultaneously increasing the probability of individuals jumping out of a local optimum. In addition, other scholars have conducted research on improvement methods, including k-means clustering,¹⁸ Cauchy's mutation,¹⁹ and other strategies. However, there is still a lack of search ability in the late stage of the optimization search, and premature convergence easily occurs. Further research is needed to explore better improvement strategies to enhance the convergence accuracy and global optimization performance of the algorithm.

This study focuses on dual-runway airports and considers constraints such as multiple runways, flight fairness, and terminal area arrival flow. The latest RECAT technology is used for wake turbulence separation with the objective of minimizing the total delay time (TDT). A multirunway arrival flight sorting optimization model is established, and a new improved sparrow search algorithm (ISSA) is proposed. Based on the standard Sparrow Search Algorithm, the ISSA initializes the population using the Chebyshev chaotic mapping strategy²⁰ to ensure population diversity and avoid becoming trapped in local optima. The ISSA incorporates the golden sine strategy²¹ to update the discoverers, addressing the issues of local optima and a large search space. The ISSA also incorporates the Lévy flight strategy²² to update the followers, enhancing the algorithm's exploration capability in the solution space and increasing the likelihood of escaping local extreme values. Finally, the ISSA incorporates a variable neighborhood search strategy to escape local optima.

Arrival flight sequencing model establishment

Parameter setting

The variable symbols used and their meanings are provided in Table 1.

Table 1.

Variable symbols and meanings.

Parameter	Description
$F$	Flight collection, $F = {f_{1}, f_{2}, \dots, f_{n}}$
$R$	Runway collection, $R = {r_{1}, \dots r_{m}}$
$i, j$	Arrival flight sequence
$E s t_{i r}$	Expected arrival time of flight i at runway r
$A c t_{i r}$	Actual arrival time of flight i at runway r
$T_{d e l a y}$	Total flight delay time
$W_{r i j}$	The minimum wake interval of aircraft landing on the same runway
$M_{a}$	Runway 1 can accommodate the wake separation matrix of the aircraft.
$M_{b}$	Runway 2 can accommodate the wake separation matrix of the aircraft.
$P_{i 0}$	The location of flight i in the entire queue without algorithm processing
$P_{i 1}$	Location of optimized flight i in the group
$M P S$	Maximum number of movements per aircraft
$W_{a i}$	According to the controller’s command ability, the terminal area waits for the incoming flight capacity.
$F_{a}$	The actual number of flights waiting for arrival
$O_{i j}$	Oblique distance interval of aircraft approaching successively on adjacent runways in correlated approach mode
$X_{i r}$	$X_{i r} = {\begin{matrix} 1, Flight i landed on runway r \\ 0, otherwise \end{matrix}$

Model assumptions

When there is an airport in the terminal area, only the situation of arriving flights is considered.

To maintain fairness and impartiality, this paper does not consider VIP flights and does not classify flights except for the purpose of dividing the wake interval.

The scheduled arrival times and arrival routes of the flights are known.

Objective function

It is assumed that there are n flights waiting to land in the terminal area during a certain period of time and that there are m runways for flight landing. The overall optimization goal is to minimize the total flight delay time:

min T_{d e l a y} = min \sum_{f \in F} \sum_{r \in R} | E s t_{i r} - A c t_{i r} |

(1)

Constraint condition

(1) RUNWAY CONSTRAINTS:

\sum_{r = 1}^{m} X_{i r} = 1, i = 1, 2, \dots, n

(2)

The above formula indicates that each flight can only land on a unique runway according to its own aircraft type.

(2) FLIGHT WAKE INTERVAL CONSTRAINT:

During the approach and landing phases of a flight, the aircraft can be affected by wake turbulence. To facilitate modeling and calculations, the more refined RECAT-CN wake turbulence separation standards²³ based on categorization can improve the efficiency of flight sequencing. The required time intervals can be determined based on the aircraft's category and wake turbulence intensity, as shown in Table 2.

A c t_{j r} - A c t_{i r} \geq W_{i j} \forall i, j \in f, r \in R

(3)

Table 2.

RECAT wake turbulence separation standards (unit: seconds).

	Next flight
Previous flight	J	B	C	M	L
J	60	113	135	160	180
B	60	60	90	113	160
C	60	60	60	75	135
M	60	60	60	60	113
L	60	60	60	60	60

The above formula indicates that the flight will be affected by the wake when it lands, and the difference between the actual landing times of the next flight j and the previous flight i is greater than the specified time interval standard. When r = 1, $W_{r i j} \in M_{a}$ ; when $r = 2$ , $W_{r i j} \in M_{b}$ , where

M_{a} = (\begin{matrix} 60 113 135 160 180 \\ 60 60 90 113 160 \\ 60 60 60 75 135 \\ 60 60 60 60 113 \\ 60 60 60 60 60 \end{matrix}), M_{b} = (\begin{matrix} 60 60 90 113 160 \\ 60 60 60 75 135 \\ 60 60 60 60 113 \\ 60 60 60 60 60 \end{matrix}) .

(3) OBLIQUE SPACING CONSTRAINT:

A c t_{j r} - A c t_{i (r + 1)} \geq O_{i j} \forall i, j \in f, r \in R

(4)

A c t_{j r} - A c t_{i (r - 1)} \geq O_{i j} \forall i, j \in f, r \in R

(5)

The above formula is intended to ensure that the oblique distance between two adjacent flights is greater than the required oblique distance interval in the parallel double-runway configuration. Here,

O_{i j}

is 4 km in the correlated approach mode. Otherwise, the interval constraint is not considered when approaching independently.

(4) ARRIVAL FLOW IN THE TERMINAL AREA:

F_{a} < W_{a i}

(6)

Considering the communication pressure of the controller, the arrival flow in the terminal area is controlled. The number of aircraft in the sorting area of the terminal area should not exceed the command capacity of the controller.

(5) NUMBER OF FLIGHT MOVING POSITIONS:

| P_{i 1} - P_{i 0} | < M P S

(7)

Considering the fairness of flights and the workload of controllers, the order of aircraft moving forward or backward during scheduling cannot exceed the maximum number of flight movement positions. This article takes

M P S = 5

Improved sparrow algorithm

Algorithm design

In this paper, four optimization strategies are used in the improved algorithm: (1) the Chebyshev chaotic mapping strategy, (2) the golden sine strategy, (3) the Lévy flight strategy, and (4) the variable neighborhood search strategy. A schematic of the improved algorithm strategy is shown in Figure 1.

(1) Chebyshev chaotic mapping population initialization

Figure 1.

Schematic diagram of the code.

In the sparrow algorithm,²⁴ randomly generating the initial population will increase the difference in the population and make the distribution of the search space uneven. The introduction of Chebyshev chaotic mapping can increase the difference and diversity of the population and effectively prevent premature convergence of the sparrow algorithm. The definition of this mapping is as follows:

x_{n + 1} = \cos (k \arccos x_{n}), x_{n} \in [- 1, 1], k \in [0, 4]

(8)

In the above equation,

x_{n}

is the nth generation population, and k is the order of chaos. The sequence produced for

k \geq 2

is chaotic and ergodic; in this paper, k = 4.

(2) Golden sine strategy update discoverer

In the sparrow algorithm, the discoverer's extensive global search for food may easily fall into a local optimal solution at the beginning of the search. To avoid premature convergence, the golden sine algorithm is introduced into the process of the discoverer searching for food. The improved position update formula is as follows:

X_{i, d}^{t + 1} = {\begin{matrix} X_{i, d}^{t} * | \sin (r_{1}) | + r_{2} * \sin (r_{1}) * | x_{1} * X_{b e s t}^{t} - x_{2} * X_{i, d}^{t} |, R_{2} < S T \\ X_{i, d}^{t} + Q * L, R_{2} \geq S T \end{matrix}

(9)

In the formula,

X_{b e s t}^{t}

is the current global optimal position;

r_{1}

and

r_{2}

are random numbers;

r_{1} \in [0, 2 π]

and

r_{2} \in [0, π]

determine the distance and direction of the individual position iteration, respectively;

x_{1}

and

x_{2}

are the golden section coefficients, calculated as

x_{1} = - π + (1 - τ) * 2 π

and

x_{2} = - π + τ * 2 π

, respectively;

τ

is the golden mean, taking the value of

(\sqrt{5} - 1) / 2

; Q is a random number obeying a normal distribution; L is a matrix whose elements are all 1;

R_{2} \in [0, 1]

is the warning value; and

S T \in [0.5, 1]

is the safety threshold.

(3) Lévy flight strategy update follower

In this paper, the Lévy flight strategy is introduced to generate a random step size in each iteration of the follower's position update, which can improve the global search ability. After introducing this strategy, the follower position update formula for the sparrow population can be expressed as:

X_{i, d}^{t + 1} = {\begin{matrix} Q * \exp (\frac{X_{w}^{t} - X_{i, d}^{t}}{α * i t e r_{max}}), i > n / 2 \\ X_{p}^{t + 1} + | X_{i, d}^{t} - X_{p}^{t + 1} | * A^{+} * L * s t e p, o t h e r w i s e \end{matrix}

(10)

In the formula, A is a matrix with elements −1 and 1, and

\overset{+}{A} = \overset{T}{A} (A \overset{T}{A} \overset{- 1}{)}

; when

i > n / 2

, the sparrow is in a hungry state, making it search for food.

s t e p

is the random step size.

(4) Position update of the vigilant observer

The vigilant observer is based on the discoverer and the joiner, and its position is updated as follows:

X_{i, d}^{t + 1} = {\begin{matrix} X_{b e s t}^{t} + δ * | X_{i, d}^{t} - X_{b e s t}^{t} |, f_{i} > f_{g} \\ X_{i, d}^{t} + K * (\frac{| X_{i, d}^{t} - X_{w o r s t}^{t} |}{(f_{i} - f_{w}) + ε}), f_{i} = f_{g} \end{matrix}

(11)

In the formula,

X_{b e s t}^{t}

represents the current global optimal position;

δ

is a random number that conforms to a normal distribution and has a variance of 1;

K \in [- 1, 1]

is a random number; and

f_{i}

represents the fitness value of the ith sparrow.

f_{g}

represents the current best global fitness value;

f_{w}

represents the current worst global fitness value; and the value of ε is any small constant.

(5) Variable neighborhood search strategy

In this paper, the variable neighborhood search strategy²⁵ is used to make the algorithm jump out of local optima by changing the neighborhood structure. In the jitter stage, the method of moving the current solution toward the individual historical optimal solution and the population historical optimal solution is adopted. The neighborhood search is implemented with three operators, exchange, insertion and 2-opt, and the neighborhood search is nested within the jitter process.

Algorithm performance test

To verify the effectiveness of the ISSA, this paper selects six sets of standard functions and compares the results of the ISSA with those of the traditional GA, the widely used particle swarm optimization algorithm (PSO), the whale optimization algorithm (WOA), and the sparrow search algorithm (SSA). In the test function, F1, F2, F4, and F7 are unimodal test functions, and F8 and F11 are multimodal test functions; the specific functions are shown in Table 3 . The experimental environment is a Windows 10 system, the software used is MATLAB2021A, the population number N is 100, and the maximum number of iterations is 1000. The hyperparameter information of the algorithm is as follows: the learning factor of the PSO algorithm is $c_{1} = c_{2} = 2$ , the maximum speed is $v_{max} = 6$ , the crossover probability of the GA algorithm is pc = 0.8, and the mutation probability is pm = 0.05. For the SSA and ISSA, the number of discoverers is 20, the number of participants is 80, the number of warnings is 10% of the entire population, and the safety threshold is 0.6.

Table 3.

Test functions.

Expression	Function name	Dim	Definition
$F_{1} = \sum_{i = 1}^{n} x_{i}^{2}$	Sphere	30	[−100,100]
$F_{2} = \sum_{i = 1}^{n} \| x_{i} \| + \prod_{i = 1}^{n} \| x_{i} \|$	Schwefel 2.22	30	[−10,10]
$F_{4} = max_{i} {\| x_{i} \|, 1 \leq i \leq n}$	Schwefel 2.21	30	[−100,100]
$F_{7} = \sum_{i = 1}^{n} i x_{i}^{4} + r a n d o m [0, 1)$	Quartic	30	[−1.28,1.28]
$F_{8} = \sum_{i = 1}^{n} - x_{i} \sin (\sqrt{\| x_{i} \|})$	Schwefel 2.26	30	[−500,500]
$F_{11} = \frac{1}{4000} \sum_{i = 1}^{n} x_{i}^{2} - \prod_{i = 1}^{n} \cos (\frac{x_{i}}{\sqrt{i}}) + 1$	Griewank	30	[−600,600]

The iterative results obtained are shown in Figure 2.

Figure 2.

Comparison of algorithm results.

The six test functions in the table from top to bottom correspond to (a)–(f) in Figure 2. The above diagram clearly demonstrates that in the case of the four unimodal test functions (a), (b), (c), and (d), the ISSA exhibits the highest accuracy and the fastest iteration speed compared to the other four algorithms. Compared with the traditional SSA, which can find the optimal solution in 600 iterations, the ISSA requires only 300 iterations. Additionally, in the multipeak test function F5 shown in Figure 2(e), although the solution result of the ISSA is comparable to those of the other algorithms, its optimal value remains lower and more stable. In conclusion, the ISSA demonstrates faster convergence speed and superior optimization ability, providing evidence for its feasibility.

Optimization of arrival flight sorting based on the ISSA

Data encoding and decoding

In this paper, the flight is encoded based on the flight number, the number of encoding layers is an integer input, and the encoding type is a one-dimensional array input. Assume that there are currently six flights sorted as follows: $f_{1}, f_{2}, \dots, f_{6}$ . The corresponding flight numbers are 1–6. When the population is initialized, random number coding is used to determine the value range of each layer according to the maximum number of moving positions of the flight. The time of the flight is converted from hours and minutes to seconds; that is, $x (h) = x \times 3600 (s)$ , $x (min) = x \times 60 (s)$ .

Fitness function

This paper studies the sequencing of arrival flights in the terminal area. To minimize the TDT of flights, the fitness function is set as the reciprocal of the objective function of the model. The fitness function is defined as $f (x) = \frac{1}{T_{d e l a y}} = \frac{1}{\sum_{f \in F} \sum_{r \in R} | E s t_{i r} - A c t_{i r} |}$ .

Algorithm flow

The specific algorithm process is shown in Figure 3.

Figure 3.

Algorithm flow chart.

Example verification

Algorithm design

To conduct accurate research on delayed flights, in this paper, we select a certain period of time on a certain day at Kunming Changshui Airport to conduct targeted research on the sequencing of arrival flights: the arrival flight data during the peak period of a typical festival are chosen, and the total number of arrival flights is 457. The specific flight data are shown in Tables 4 and 5. The data from small-scale flights (29) and large-scale flights (147) are selected for the simulation experiments, and the algorithm proposed in this paper is verified on two dimensions.

Table 4.

Arrival delay flight data (small-scale flight data).

No.	Flight number	RFC	Aircraft type	ETA
1	DR6522	Chengdu-Kunming	B738	0001
2	8L9802	Lijiang-Kunming	A320	0003
…	…	…	…	…
28	8L9928	Baoshan-Kunming	A320	1111
29	8L9504	Zhengzhou-Kunming	B738	1117

Table 5.

Large-scale flight data.

No.	Flight number	RFC	Aircraft type	ETA
1	DR6522	Chengdu-Kunming	B738	0001
2	A67118	Lijiang-Kunming	A320	0251
3	8L9866	Nanchang-Kunming	A320	0013
4	8L9872	Ningbo-Yichang-Kunming	A320	0018
…	…	…	…	…
143	MU7892	Shanghai-Kunming	B738	1217
145	8M756	Yangguang-Kunming	A320	1327
146	A67106	Mangshi-Kunming	A320	1251
147	MU5864	Chongqing-Kunming	B738	1210

Analysis of effectiveness

(1) Small-scale flight delays

A total of 29 delayed flights from 0:00–12:00 on a certain day were selected for simulation. The relationship between the value of the adaptive function and the number of iterations after 100 iterations using the five algorithms is shown in Figure 4(a). A comparison of the cumulative TDTs of the 29 delayed flights is shown in Figure 4(b).

Figure 4.

(a) Convergence curves of the objective function. (b) Cumulative total delay time comparison chart.

In Figure 4(a), it is evident that the flight delay time progressively decreases with increasing number of iterations. Both PSO and the GA converge after approximately 40 iterations, while the ISSA proposed in this paper achieves the optimal value by the 10th iteration. Despite the high convergence speed of the WOA, it yields the poorest results. Therefore, the results of the ISSA presented in this paper are optimal, and it exhibits faster convergence speed than the other algorithms. Figure 4(b) shows a comparison chart of the cumulative TDT of all the delayed flights. Currently, the widely used method for sorting aircraft takeoff and landing on runways is the FCFS algorithm, which determines the order in which aircraft use the runway for landing according to the flight's estimated time of arrival (ET) on the runway. The advantage of this algorithm is that it is simple and easy to implement and is widely used in many studies.^2–4,10,11 Therefore, this paper uses it as one of the comparative algorithms. Upon comparing the flight delay results obtained by the ISSA with those obtained by the FCFS algorithm, it is apparent that from a global perspective, the ISSA incurs the least total flight delay time while ensuring fairness.

(2) Large-scale flight delays

A total of 147 flights from 0:00–12:00 on a certain day were selected for simulation. The relationship between the value of the post-adaptation function and the number of iterations after 100 iterations of the five algorithms is shown in Figure 5(a). A comparison of the cumulative TDTs of the 29 delayed flights is shown in Figure 5(b). The data corresponding to each algorithm are shown in Table 6 .

Figure 5.

(a) Convergence curves of the objective function. (b) Cumulative total delay time comparison chart.

Table 6.

Comparison of the results of different algorithms.

Algorithms	Small-scale flight delay time/s	Algorithm optimization comparison	Large-scale flight delay time/s	Algorithm optimization comparison
FCFS	11,224	50.33%	460,021	20.47%
GA	8037	30.63%	424,821	13.88%
WOA	7760	28.16%	412,892	11.39%
PSO	7223	22.82%	404,640	9.58%
SSA	7194	22.5%	398,662	8.22%
ISSA	5575	—	365,876	—

It can be clearly seen from Figure 5(a) that the ISSA operation results are optimal. Figure 5(b) shows that the total flight delay time of the optimized algorithm is lower than that of the FCFS algorithm, and from a global perspective, the total flight delay time of the ISSA is the smallest under the large number of flights. Combined with the results of Table 5 , the TDT of all flights optimized in the case of small-scale flights is 5575 s, which is 50% lower than that of the FCFS algorithm and 30.36%, 28.16%, 22.82%, and 22.5% lower than that of the GA, WOS, PSO, and SSA methods, respectively. In the case of large-scale flights, the TDT of all flights after optimization is 365,876 s, which is 20.47% lower than that of the FCFS algorithm and 13.88%, 11.39%, 9.58%, and 8.22% lower than that of the GA, WOS, PSO, and SSA methods, respectively. The ISSA proposed in this paper is generalizable and can effectively reduce the total flight delay time, especially for small-scale flight scheduling problems. The use of airport runways before and after the optimization of small-scale flight sequencing is shown in Figure 6(a) and (b), respectively. A comparison of the order of arriving flights is shown in Figure 7.

Figure 6.

(a) Runway usage before optimization. (b) Runway usage after optimization.

Figure 7.

Comparison before and after flight sequence optimization.

In Figure 6(a) and (b), 1 and 2, respectively, represent the east runway (R1) and the west runway (R2). The optimized runway usage is more evenly distributed, which solves the problem of local runway usage being too dense and effectively alleviates runway congestion. Figure 7 shows a comparison of the optimized flight sequences before and after optimization. For example, the first flight sequence is adjusted to the second one after improvement, and the third and twentieth flights have the largest movement, moving five positions. However, the maximum number of movements of all flights after improvement does not exceed five, which corresponds to the constraint conditions mentioned earlier and ensures fairness in flight delays and coordination of the interests of various airlines.

This paper has been compared with other studies to highlight the strengths of the research presented here. The results of the comparison are shown in Table 7.

Table 7.

Comparison with other studies.

The compared studies	Feature	Advantages of this article
^9,10,11,13	Running time of algorithm >10 min	Running time of algorithm <3 min
^2,3,4,10,11	The algorithm is only compared with the FCFS algorithm.	A horizontal comparison is conducted with the FCFS algorithm and other optimization algorithms.
^5,11	Modeling is conducted without considering the impact of the model.	The latest RECAT-CN wake spacing standards are adopted for different aircraft models.
¹⁶	The optimal results of the improved algorithm are as follows: F1: 6.186 × 10⁻⁷⁸, F2:1.745 × 10⁻⁴⁰ F4: 5.188 × 10⁻³⁹, F7:0: F8: 8.882 × 10⁻¹⁶, F11: −3.306	The optimal results of the improved algorithm are as follows: F1: 5.7651 × 10⁻¹²⁵, F2:2.01736 × 10⁻⁴⁷ F4: 4.00208 × 10⁻⁸⁰, F7: −0.0565 F8: 1.2569 × 10⁻⁴, F11: −3.409 Therefore, the results in this paper are better.

When there are many congested flights at an airport during a peak flight period, the method in this paper can be used to optimize the sequencing of incoming flights and effectively alleviate airport congestion. However, this paper has several shortcomings. Weather factors are not considered in the modeling. In fact, weather will also affect the sequencing of airport flights. In addition, there will be departing flights and arriving flights participating in the sequencing of flights in the terminal area of the same airport.

Conclusion

To alleviate the pressure on terminal areas and air traffic control personnel and reduce airline operating costs, this article focuses on inbound delayed flights. The main research areas are as follows:

A multirunway airport arrival flight sorting optimization model is established to address airport arrival flight delays, reduce passenger delay times, and address various factors, such as airlines, airports, and air traffic control personnel. This model considers runway, wake, and terminal area inflow factors to achieve comprehensive optimization.

Based on RECAT technology, the improved RECAT-CN wake interval standard is used to convert the interval time to its minimum value, which is convenient for model calculations. This approach makes the interval time more compact to improve the operation efficiency of the terminal area.

We improved upon the traditional SSA. The improved algorithm was tested on six basic test functions and compared with other algorithms, such as PSO, the WOA, and the GA. The results showed that the designed algorithm has certain advantages over other algorithms.

Using two sets of actual flight data at small and large scales for the simulation experiments, the algorithm obtained good results, verifying its practicality and effectiveness in the study of flight scheduling optimization problems. This paper provides some theoretical support for reducing flight delays in busy airports.

When constructing the model, the article did not consider the impact of weather. The proposed method can only optimize flights in a certain period, which is challenging for optimizing many real-time flights. Therefore, future research should further consider departure flight factors, weather factors, and sequencing optimization for a large number of flights to optimize the problem of arrival and departure flight delays.

Footnotes

Declaration of conflicting interests

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

This work is supported by Yunnan Fundamental Research Projects (grant No. 202101AU070041).

ORCID iD

Weixi Zhao

Author biographies

Weixi Zhao is a master's student. Her main research areas are transportation planning and management, swarm intelligence, and deep learning.

Te Liang is currently pursuing a PhD degree with the School of instrumentation and Optoelectronic Engineering from Beihang University, Beijing, China.

References

Wieland

Tyagi

Kumar

, et al. METROSIM: A metroplex-wide route planning and airport scheduling tool. In: 14th AIAA Aviation Technology, Integration, and Operations Conference. 2014.

Wang

Zhang

. Modeling and algorithm of arrival and departing aircraft sequencing in multi-airport terminal area collaborative decision. Flight Dynamics 2016; 34: 90–94.

Zhou

Jiang

. Research on arrival/departure scheduling of flights on multirunways based on genetic algorithm. Math Probl Eng 2014; 2014: 1–13.

Tian

Wan

Han

. Optimization of terminal airspace operation with environmental considerations. Transport Res D-TR E 2018; 63: 872–889.

Jiang

Zhang

. Optimal sequencing of arrival flights at metroplex airports: a study on shared waypoints based on path selection and rolling horizon control. Aerospace 2023; 10: 81.

Bae S

Park

Clarke J

. Modified mixed integer linear program for airport departure scheduling. In: AIAA Guidance, Navigation, and Control (GNC) Conference. 2013. doi:https://doi.org/10.2514/6.2013-4885.

Pinol

Beasley J

. Scatter search and bionomic algorithms for the aircraft landing problem[J]. Eur J Oper Res 2006; 171: 439–462.

Wang

Zhang

. Dynamic departure sequencing based on branch and bound algorithm. J Nanjing Univ Aeronaut Astronaut 2015; 47: 547–552.

Xiao-Bing

Di Paolo

. Binary-representation-based genetic algorithm for aircraft arrival sequencing and scheduling. IEEE T Intell Transp 2008; 9: 301–310.

10.

Zhang

, et al. Optimized method for collaborative arrival sequencing and scheduling in metroplex terminal area. Acta Aeronaut 2015; 36: 2279–2290.

11.

Zhang

Zhong

. Research on the problem of sorting parallel runway aircraft in terminal area based on genetic algorithm. IOP Conf Ser: Mater Sci Eng 2019; 612: 042040.

12.

Hancerliogullari

Rabadi

Al-Salem

, et al. Greedy algorithms and metaheuristics for a multiple runway combined arrival-departure aircraft sequencing problem. J Air Transp Manag 2013; 32: 39–48.

13.

Salehipour

Modarres

Naeni

. An efficient hybrid meta-heuristic for aircraft landing problem. Comput Oper Res 2013; 40: 207–213.

14.

MU X

LVX

. Multi-threshold image segmentation based on improved sparrow search algorithm. Sys Eng Electron 2021; 43: 318–327.

15.

Liu

Shu

Liang

, et al. A modified sparrow search algorithm with application in 3d route planning for UAV. Sensors 2021; 21: 1–16.

16.

MU X

LVX

, et al. Chaotic sparrow search optimization algorithm. J Beijing Univ Aeronaut Astronaut 2021; 47: 1712–1720.

17.

Tang

, et al. A fusion crossover mutation sparrow search algorithm. Math Probl Eng 2021; 21: 1224.

18.

Kusy

Kowalski

. Architecture reduction of a probabilistic neural network by merging k-means and k-nearest neighbour algorithms. Appl Soft Comput 2022; 128: 109387.

19.

Zeng

Wei

. Multi-objective particle swarm optimization algorithm using Cauchy mutation and improved crowding distance. Int J Intell Comput 2023; 16: 250–276.

20.

Wang

, et al. Restoration of high-frequency vibration blurred hyperspectral image based on dynamic chaos disturbance genetic algorithm. Spectrosc Spect Anal 2021; 41: 2219–2225.

21.

Tanyildizi

Demir

. Golden sine algorithm: a novel math-inspired algorithm. Adv Electr Comput En 2017; 17: 71.

22.

Yang

Ling

. Parameter identification of proton exchange membrane fuel cells model based on improved chicken swarm optimization algorithm. Acta Energiae Solaris Sinica 2023; 44: 269–278.

23.

Zhao

, et al. Safety assessment method of the aircraft wake vortex re-categorizing standard. J Saf Environ 2020; 20: 1277–1283.

24.

Xue

Shen

. A novel swarm intelligence optimization approach: sparrow search algorithm. Syst Sci Control Eng 2020; 8: 22–34.

25.

Hansen

Mladenović

. Variable neighborhood search: principles and applications. Eur J Oper Res 2001; 130: 449–467.

Optimization of terminal area arrival flight sorting based on an improved sparrow search algorithm

Abstract

Keywords

Introduction

Arrival flight sequencing model establishment

Parameter setting

Model assumptions

Objective function

Constraint condition

Improved sparrow algorithm

Algorithm design

Algorithm performance test

Optimization of arrival flight sorting based on the ISSA

Data encoding and decoding

Fitness function

Algorithm flow

Example verification

Algorithm design

Analysis of effectiveness

Conclusion

Footnotes

Declaration of conflicting interests

Funding

ORCID iD

Author biographies

References