Research Article  Open Access
A User Equilibrium Assignment Flow Model for Multiairport Open Network System
Abstract
To reduce flight delays and promote fairness in air traffic management, we study the imbalance problem between supply and demand in airport network system from the view of both the system and the users. First, we establish an open multiairport oriented network flow system with the correlation between the arrival and departure in capacityconstrained airports, as well as the relevance between multiairports united flights. Then, based on the efficiency rationing principle, we propose an optimization model to reassign flow with user equilibrium constraints. These constraints include Gini coefficient, system capacity, and united flights. The model minimizes the total flight delay cost of capacityconstrained airports in the network system. We also introduce some evaluation indexes to quantitatively analyze fairness among airlines. Finally, with an open multiairport network system and its actual flights data in China, the model is verified. Test results show that the model can be used to coordinate and optimize matching the flow and capacity in the multiairport, make full use of the capacity of airports, and minimize the system delays. The findings can provide a powerful reference of developing scientific and rational assignment strategy for air traffic controllers.
1. Introduction
It is clear that the phenomenon of growing traffic demand should be met by a concomitant improvement in airport capacity. The longterm programs for increasing airport capacity include construction of new airports and expansion of runway systems at existing airports; the shortterm programs consider new operational methods in air traffic flow management (ATFM) and capacity utilization as potentially effective measures for improving the existing capacity resources [1, 2]. Analysis shows that optimization of the present airport system by the operational and technological measures might result in increasing current traffic flow by up to 50% [3]. ATFM aims to optimize flow distribution, improve resource utilization, and promote safety in air traffic by rationalizing the traffic demand and adjusting traffic flow, whose essential purpose is to adjust the air traffic demand and supply with coordinated manner [4]. With the rapid development of air transport, the contradiction between the dramatic increase of air traffic demand and the relative decrease of capacity supply is becoming increasingly prominent, and the air traffic congestion is getting worse. Airports are air traffic management bottlenecks, studying the imbalances between supply and demand of airport traffic has an important realization meaning in solving air traffic congestion, reducing flight delays, and promoting the safety of air traffic. The issue of imbalance between airports traffic demand and supply is mainly that the airports traffic demand exceeds the supply of airports capacity. It belongs to strategic and pretactical air traffic flow management issue. The essence is the coordination configuration between flow and capacity.
Air traffic flow allocation problem has been a growing concern on scholars home and abroad. In order to assign capacity and allocate flow, Gilbo, Ma, Chen, and others have, respectively, established a kind of optimization model to minimize delayed flights quantity or losses in arrival and departure queue. Gilbo presented a model which first considers the runways and arrival and departure fixes jointly as a single system resource and secondly considers arrivals and departures simultaneously as two interdependent processes [4]. The model takes into account the interaction between runway capacity and capacities of fixes to optimize the traffic flow through the airport system. The effects are achieved by dynamic timedependent allocation of airport capacity and flows between arrivals and departures coordinated with the operational constraints at runways and arrival and departure fixes as well as dynamics of traffic demand and weather. Ma et al. proposed a model for traffic flow allocation and optimal scheduling in the terminal area [5]. The model provides more efficient utilization of existing terminal capacity in each time interval to optimize flow distributions and alleviate consequences of flight delays. The model treats the arrivals and departures as two interdependent processes to optimize the airport flow distribution. The model also considers the dynamic characteristics of the timedependent allocation of airport capacity and the tradeoff between arrivals and departures along with the dynamic characteristics of traffic demand and weather. Chen et al. set up a multiobjective model to optimize the airside capacity utilization and the traffic flow distribution, which maximize flow through runway and arrival/departure fix by introducing utilization satisfaction level [6]. The result indicates that the proposed model is more effective to decrease overall flight demand queues than Gilbo model results. Within the considered time period, all flight demands are satisfied without any delay at the end of the last 15 min interval by optimal allocation of flow for each arrival/demand fix of airside. Thus, it guarantees the capacity utilization satisfaction level and gives an effective utilization of airport airside capacity and reduces flight delays. Dell’Olmo and Lulli studied the problem of finding the optimal tradeoff between the number of arrivals and the number of departures in order to reduce a delay function of all the flights [7], using a more realistic representation of the airport capacity, that is, the capacity envelope. Under the assumption of piecewise linear convex capacity envelopes and of the exact interpolation of all the paretooptimal operational points, they show that the problem can be formulated as a linear programming model. For general airport capacity envelopes, they proposed a dynamic programming formulation with a corresponding backward solution algorithm, which was robust and easy to implement and had a linear computational complexity.
With the development of the collaborative decisionmaking (CDM) in ATFM, collaborative optimization problem emerged. Under collaborative optimization, arrivals and departures are considered jointly and are treated as interdependent operations. When arrival capacity can be traded for departure capacity, collaborative optimization provides the best allocation of capacity to arrivals and departures to maximize the airport throughput and minimize delays. Furthermore, collaborative optimization uses flight priorities from airlines and other users in the optimization model. Dell’Olmo and Lulli presented a collaborative optimization model taking into account user priorities (preserving integrity of banks at hub airports), which can be used as a decision support tool for solving ATFM problems at airports [7]. During years 2008 and 2010, Zhang and Hu took into account airport capacity dynamic constraints, tradeoff between arrival and departure, or/and the interaction between multiple airports and, respectively, established some collaborative optimization models of capacity utilization and flow assignment for single and multiple airports terminal area [8–11]. To find optimal strategy, these models minimized the total flights’ weighted delay costs and introduced delay cost coefficients for each flight as airline’s decisionpreferential information. The results of some experiment have shown that these proposed models can not only have a balance between traffic flows and capabilities and make best use of system resource but also minimize delay cost with airline’s preferences.
Researches above are achieved by optimizing by flight lineup or with flight delay coefficient taking into account users’ expection or preference information. They make the issue more reasonable to some degrees but have not considered the fairness among users, that is, the delay degree of influence on the airlines, which is caused from the imbalance between supply and demand. What is more, they have not considered the spread phenomenon of delay from the system’s view. In this view, from both system optimization and user equilibrium perspectives, we research the issue of imbalance between demand and supply in an open airport network system, so as to find out a tradeoff among the balancing users’ needs, optimizing the system’s flow distribution and making full use of the system’s capacity.
2. Open Multiairport Network System
An open multiairport network system is a kind of air traffic network system in which some airports are known but the others are unknown, referring to Figure 1. The system considers not only the flights connecting among the airports, but also the flights turnover within a specific airport. These two elements make up the air traffic flow within the network system. An airport in the system is shown by a model which consists of points A, B, C, D, and E. Under normal circumstances, the arrival flights from other airports enter into point A; then, some flights enter into point B ending their journeys; the other continuous flights turn over after point C and turn into departure flights with the flights from point D. The selfloop arc at point C indicates the turnover; the selfloop arc at point D stands for the delay in Figure 1.
An open multiairport network system flow is made up by these airports’ flow. An airport flow consists of the arrival and the departure flow. The capacity of the network system is mainly reflected in the capacity of each airport, including arrival and departure capacity. Since the arrival and departure of the flights are closely related and interacting with each other, airport capacity also significantly changes with the weather conditions and is even directly related to the type of aircrafts, the ratio of combination, and the order of takeoff and landing. Researches have shown that the arrival and departure capacity can usually be displayed as the airport capacity curve figure as Figure 2 shows, which varies under different weather conditions [12, 13].
The vertical axis indicates the number of arrival flights every 15 minutes and the horizontal axis is the number of departure flights every 15 minutes. If the traffic demand exceeds the capacity limits, such as point P in Figure 2, delay is inevitable. In order to minimize the delay, a Pareto optimistic capacity configuration point (or flow distribution point) on or inside the curve should be found to rearrange the traffic flow and make the best use of capacity. As a result, in the airport network system, if a certain airport flow exceeds its capacity (as the grey airport shown in Figure 1), some airports in the system will be affected and the related flights will be delayed.
3. A User Equilibrium Assignment Model
3.1. Problem Description
The user equilibrium assignment problem for the open airport network system can be described as follows: given certain air traffic demand and capacity supply at each airport within this system in a certain period of time, how to reasonably rearranging these airport’s demand and capacities is the key to minimize the total delay loss among the interests of users fairly when some airports exceed their limits. The basic clue of modeling is to divide the time period into several time slots and deploy the best arriving/departure periods for the flights in capacityconstrained airports according to the principles of users equally delayed to rearrange the demands, so as to make these airports’ traffic flow and capacity harmonious in all these time slots. The model requires four hypotheses.(1)During a given period, the capacity curves of the capacityconstrained airports under different weather conditions are known in the network system.(2)The air traffic demand not accommodated in a given slot can be always arranged in the first slot outside of the period whose capacity has no upper limits, which assures that the problem can be solved.(3)During the given period, the capacities of the other airports in the network system have no upper limits.(4)All the flights cannot take off or arrive before their scheduling time and cannot be delayed indefinitely, which is rational in real life.
3.2. Symbol Definition
, the time interval of preassignment, is composed by a number of consecutive time segments whose length is Δ (usually 15 min), whereDefine as the next time segment which is out of the time interval and its capacity is infinite.
is the set of airlines involved in the capacityconstrained airports in the interval , .
is the set of the capacityconstrained airports which are in the multiairport open network, .
is the set of flights departing from airport and arriving at airport which are both in the network system, where , .
is the set of flights in the capacityconstrained airport , where , . is composed of the set of departure flights and the set of arrival flights , wherewhere is the set of flights that depart from the capacityconstrained airport and arrive at noncapacityconstrained airports. is the set of flights that depart from the capacityconstrained airport and arrive at the capacityconstrained airport :where is the set of flights that depart from the capacityconstrained airport and arrive at noncapacityconstrained airports and is the set of flights that depart from the capacityconstrained airport and arrive at the capacityconstrained airport .
is the set of flights for airline , where
It means the matching relationship between flights and airline, where is the flights set of airline at capacityconstrained airport, where is the set of flights of airline that depart from the capacityconstrained airport and arrive at the noncapacityconstrained airport and is the set of flights of airline that depart from the capacityconstrained airport and arrive at the capacityconstrained airport . is the set of flights of airline that depart from the noncapacityconstrained airport and arrive at the capacityconstrained airport , and is the set of flights of airline that depart from the capacityconstrained airport and arrive at the capacityconstrained airport .
is the set of connecting flights that have continuous range in , where is proceeding flight, and is subsequent flight.
is the minimum time interval between and in individual capacityconstrained airports, where
is the estimated time of arrival or departure of flight , where is the estimated takeoff time and is the estimated landing time.
is the set of probable departure time segment or land time segment for the flight f, where
is the arrival and departure capacity curve of airport in segment , where , .
is the maximum arrival capacity of airport in segment .
is the maximum departure capacity of airport in segment .
is the arrival flow of airport in segment , where
is the departure flow of airport in segment , where
The decision variables are
3.3. Mathematical Model
The objective function of the model is to minimize the total delay costs of flights in capacityconstrained airports of the open network system, including the delay costs of flights to take off from and land at capacityconstrained airports, flights to take off from capacityconstrained airport and land at capacityunconstrained airports, and flights to take off from capacityunconstrained airport and land at capacityconstrained airports. In order to avoid double counting the delay in the multiairport network system, the delay costs of flights to take off from and land at capacityconstrained airports are represented in departure delay; the delay costs of flights to take off from capacityconstrained airport and land at capacityunconstrained airport are represented in departure delay; the delay costs of flights to take off from capacityunconstrained airport and land at capacityconstrained airport are represented in arrival delay, as shown in the following formulation:
Then,and is subject to the following constraints.
Segment assignment constraint: it is used to ensure that each flight has only one departure/arrival time segment in the set of takeoff/landing time segments. Among them, formulation (15) is the constraint of departure time segment assignment and formulation (16) is the constraint of arrival time segment:
Capacityconstrained constraint: it constitutes the flow distribution point field of capacityconstrained airport (within or on the capacity curve), finding optimal distribution point in the point region to ensure that arrival and departure flows match its capacity, which is a tradeoff between arrival and departure. Among them, , , and are the coefficients of the capacity curve . As shown in Figure 2, , , , and constitute the flow distribution point field. Among them, formulation (17) is the largest capacity constraints; formulation (18) is the capacity balance constraints:
Continuous flight constraint: it is used to ensure connecting flights have enough time till turnaround:
Airborne hour constraint: it is used to meet estimated airborne hours of flights between capacityconstrained airports unchanged after optimization. As the goal of optimization is to get the flow of each time segment, the constraints are shown as in this equation to ensure the model has a feasible solution. However, the optimization times of arrival and departure and estimated airborne of these flights still meet equivalent relationships:
User equilibrium constraint: it is used to ensure the equal distribution of delay among different airlines. This paper uses the Gini coefficient which reflects the fairness to quantify the equilibrium and establish the constraint. Gini coefficient reflects the extent of unfairness in resources distribution problem; its value varies from 0 to 1. It indicates that the lower the number, the more equal the distribution. For the problem of this paper, the Gini coefficient can be expressed as (21), and its practical significance is the deviation of average delay between airlines. The lower the number, the smaller the loss deviation, and it indicates that the distribution tends toward more balance and fairness:where is the delay loss of airline , defined as the following formulation:
is flight proportion of airline . It indicates that flow is allotted according to airlines’ market share, but not absolute egalitarianism. So it will be more reasonable:
Then, set does not exceed a certain small parameter value to limit deviations, order , and then establish fairness constraints:
Define
Then,
Then, constraints (24) can be turned into the following linear form:
4. User Equilibrium Evaluation Indexes
The essence of user equilibrium assignment is that airlines share the total delay loss according to the proportion of their flights. Therefore, the pros and cons of equilibrium assignment are mainly a response to total delay loss (i.e., effectiveness) and fairness among airlines. So there are two evaluation criteria. One is effectiveness criteria: minimize the total delay loss of all flights; the assignment with less total loss is better; otherwise it is more undesirable. The other is fairness criteria: balance the total loss delays among airlines as much as possible. Overall evaluation index can be reflected by means of flight delay deviation, that is, Gini coefficient. The smaller the delay deviation between airlines, the more balanced the delay distribution; otherwise, it is more uneven. To quantitatively evaluate the fairness of a specific airline, an effectivenessfair degree is used, which is defined as the ratio of airline’s total loss proportion and its flight proportion. According to the principle of distribution, the loss proportion airline bears should be equal to its flight proportion. So the effectiveness fair degree of the airline which is close to 1 indicates a fair loss degree borne by the flight proportion; that is, the loss of the airline is more fair; on the contrary, the loss is more unfair. The effectivenessfair degree of airline is expressed as the following formulation:
5. Sample Simulation Analysis
5.1. Simulation Case Study
An open multiairport network system is constituted with some capacityconstrained airports, that is, Beijing Capital Airport, Shanghai PuDong Airport, and Guangzhou BaiYun Airport. Using a typical time interval (8:00–12:00, , min) of the actual flight data on a certain day in these three airports, the initial traffic demand distribution in three airports is shown in Table 1. The capacity curves supposed for the three airports and their match charts between demand and capacity in different periods are shown in Figure 3. Apparently, due to capacity constraints, supply and demand are imbalanced in the three major airports during some time segments. The demand of arrival and departure exceeds the capacity limit in some periods (the demand point is beyond the capacity curve in the figure), and an imbalanced demand and supply will naturally result in flight delay.

(a) Beijing Capital Airport
(b) Shanghai PuDong Airport
(c) Guangzhou BaiYun Airport
5.2. Model Optimization Analysis
Airport flow assignment problem is essential to reassign times of arrival or departure of affected flights with some constraints; it belongs to combinatorial optimization problem. The model is 01 integer programming model. There are many optimization methods to solve it. Because the model can deal with a lot of flights and airports in an open network, and the actual decisionmaking needs to address the balance between supply and demand during several hours, involving a large amount of data, a method suitable for solving largescale integer programming problem should be used; we use LINGO software for the model simulation, and set . Flow distribution results of each period and the matching conditions of flow and capacity in the three airports after optimization are shown in Table 2 and Figure 4; the total delays time is 3060 minutes. Apparently, Figure 4 shows that the optimized points in all traffic are distributed within or on the capacity curve, which indicates that traffic flow has been coordinated with the capacity. Figure 5 compares the traffic flow before and after optimization in the three airports and shows the optimized traffic distribution scheme, which not only effectively eliminates the traffic peak, but also balances the flow distribution among those times, takes full advantage of the limited capacities of these airports, and gives reasonable arrangements for the flights to take off or land.

(a) Beijing Capital Airport
(b) Shanghai PuDong Airport
(c) Guangzhou BaiYun Airport
(a) Beijing Capital Airport
(b) Shanghai PuDong Airport
(c) Guangzhou BaiYun Airport
Figure 6 shows that airlines share delay loss caused by allocation. We can see that the average flight delays of airlines are close to 3.3 min, in which is less than 0.3 min, so the delay is relatively balanced; those effectivenessfair degrees of airlines are between 0.98 and 1.06 and are closed to 1, so the delay is relatively fair. For example, the effectivenessfair degrees of airlines 1 and 2 are both equal to 1, which is shown to be absolutely fair to them.
5.3. Parameter Affects Analysis
In order to analyze the impact of the parameter , we selected 10 cases to optimize. The performance and the optimization results are shown in Table 3; the analysis results of average delay and effectivenessfair degree of airlines are shown in Figure 7. It can be seen from Table 3 that, with the increase of , the total delay time is declining and the effectiveness is gradually increasing, but when it increases to a certain extent, it tends to a constant. The increasing of standard deviation of delay reflects the reduction of fairness. The optimization duration tends to decrease, indicating that a more relaxed constraint makes it easy to obtain the optimal solution. While the value of is smaller (), the fairness and effectiveness are a tradeoff, basically conflicting. In contrast, while the value of is larger (), the total delay tends to a constant, which no longer affects the effectiveness.

In Figure 7, with the decrease of , the change of average delay of airlines is leveling off and effectivenessfair degrees trend to be 1. It means that the deviation of the average delay between airlines is continuing to decline and the injustice input on the single airline has also been reduced; that is to say, fairness continues to increase. Therefore, the selection of value has an important influence on the optimization results. In accordance with the provisions of United Nations organizations, the Gini coefficient which is between 0.2 and 0.3 shows balance. Considering effectiveness, fairness, and efficiency factors, we found that the values between 0.002 and 0.01 are better.
6. Conclusions
A user equilibrium assignment flow model for multiairport open network system is proposed in this paper; it can be used to balance supply and demand in airport network system with the view of both the system and users.
The test results show that the model can optimize the matching relationship of flow and capacity in a multiairport open network system and make full use of the capacity of airports and minimize the system delays, so that it can provide a powerful reference of developing scientific and rational assignment strategy for air traffic controllers.
If the model is extended to the whole air traffic network, which consists of airports, sectors, enroute, and so forth, further studies on the problem of demand and supply imbalance in more complicated largescale regional air traffic can be carried on.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This study is funded by the Fundamental Research Funds for the Central University (NS2014068).
References
 FAA, “1993 aviation system capacity plan,” Tech. Rep. DOT/FAA/ASC931, Office of System Capacity and Requirements, Washington, DC, USA, 1993. View at: Google Scholar
 S. G. Hamzawi, “Lack of airport capacity: exploration of alternative solutions,” Transportation Research Part A, vol. 26, no. 1, pp. 47–58, 1992. View at: Google Scholar
 Airport Association, Airport/Airway System Capacity Needs, Airports Council International—North America, Washington, DC, USA, 1992.
 H. H. Zhang, The key technologies of collaborative flow management in airport terminal area [Ph.D. thesis], Nanjing University of Aeronautics and Astronautics, Nanjing, China, 2009, (Chinese).
 Z. P. Ma, D. G. Cui, and Y. L. Xie, “Flow distribution and optimal scheduling in airport terminal area,” Journal of Tsinghua University (Science and Technology), vol. 43, no. 7, pp. 876–879, 2003 (Chinese). View at: Google Scholar
 X. Chen, X. Lu, and J. Zhu, “Optimization model of hubairport airside capacity utilization and flow allocation,” Journal of Nanjing University of Aeronautics and Astronautics, vol. 39, no. 5, pp. 680–684, 2007. View at: Google Scholar
 P. Dell'Olmo and G. Lulli, “A dynamic programming approach for the airport capacity allocation problem,” IMA Journal of Management Mathematics, vol. 14, no. 3, pp. 235–249, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 H. Zhang, M. Hu, and S. Chen, “Collaborative distribution strategy of airport arrival and departure traffic flow,” Journal of Nanjing University of Aeronautics and Astronautics, vol. 40, no. 5, pp. 641–645, 2008 (Chinese). View at: Google Scholar
 H. Zhang, M. Hu, and S. Chen, “Collaborative optimization of capacity utilization and flow assignment in airport terminal area,” Journal of Southwest Jiaotong University, vol. 44, no. 1, pp. 128–134, 2009. View at: Google Scholar
 H.H. Zhang and M.H. Hu, “Modeling and simulation of capacity utilization and flow assignment in multiairport terminal area,” Journal of System Simulation, vol. 21, no. 18, pp. 5639–5644, 2009 (Chinese). View at: Google Scholar
 Y. Tian, K. Song, and Y. H. Gu, “The research on rerouting problem in air traffic flow management,” Mathematics in Practice and Theory, vol. 38, no. 10, pp. 70–76, 2008 (Chinese). View at: Google Scholar
 E. P. Gilbo, “Optimizing airport capacity utilization in air traffic flow management subject to constraints at arrival and departure fixes,” IEEE Transactions on Control Systems Technology, vol. 5, no. 5, pp. 490–503, 1997. View at: Publisher Site  Google Scholar
 E. P. Gilbo, “Airport capacity: representation, estimation, optimization,” IEEE Transactions on Control Systems Technology, vol. 1, no. 3, pp. 144–154, 1993. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Honghai Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.