Research Article  Open Access
Haiwen Xu, Songchen Han, "Weighted TimeBand Approximation Model for Flight Operations Recovery considering Simplex Group Cycle Approaches in China", Mathematical Problems in Engineering, vol. 2016, Article ID 3201490, 17 pages, 2016. https://doi.org/10.1155/2016/3201490
Weighted TimeBand Approximation Model for Flight Operations Recovery considering Simplex Group Cycle Approaches in China
Abstract
The timeband approximation model for flight operations recovery following disruption (Bard, Yu, Arguello, IIE Transactions, 33, 931–947, 2001) is constructed by partitioning the recovery period into time bands and by approximating the delay costs associated with the possible flight connections. However, for disruptions occurring in a hubandspoke network, a large number of possible flight connections are constructed throughout the entire flight schedule, so as to obtain the approximate optimal. In this paper, we show the application of the simplex group cycle approach to hubandspoke airlines in China, along with the related weighted threshold necessary for controlling the computation time and the flight disruption scope and depth. Subsequently, we present the weighted timeband approximation model for flight operations recovery, which incorporates the simplex group cycle approach. Simple numerical experiments using actual data from Air China show that the weighted timeband approximation model is feasible, and the results of stochastic experiments using actual data from Sichuan Airlines show that the flight disruption and computation time are controlled by the airline operations control center, which aims to achieve a balance between the flight disruption scope and depth, computation time, and recovery value.
1. Introduction
Flight plans are very important for all airlines, as the efficiency of their construction can significantly affect profit margins and because they involve aircraft, crew, airports, and other important factors. Flight plans contain a series of flight schedules, which consist of the originating city, the departure time, the destination city, and the arrival time for flights the airline intends to operate. However, aircraft mechanical problems, crew unavailability due to illness, and inclement weather conditions often lead to disruptions and delay flight plan propagation [1]. With the rapid increase in scheduled flights and airlines in China, a greater number of airline schedule disruptions are occurring. Moreover, the aircraft, crew, airports, airspace ports, and airspace capacities are connected and restricted; thus, disruptions can lead to a domino effect that generates further flight disruptions. Such disruptions often result in irregular operations such as delays, cancellations, and aircraft reassignment. The majority of irregular operations can be attributed to four major factors, namely, shortage of airline resources (such as crew, aircraft, ground personnel, and gates), airport and airspace capacity shortages, inclement weather conditions, and managementlevel discontent that affects the related resources. In 2013, the economic losses due to airline schedule disruptions exceeded RMB43 billion, according to the Civil Aviation Resource Net of China. However, the Civil Aviation Development Statistical Bulletin reports that, at most, 50.49% of the schedule disruptions in 2014 can be attributed to airlines, with 26.41% being the average result. Therefore, the speed with which flight operations are recovered is very important in order to reduce economic losses due to schedule disruptions.
A large amount of research has been devoted to this topic. For example, Teodorović and Guberinić [2] have presented a network (representing an airline network) in which the nodes represent the flights and the arcs are the total time losses on individual flights. Hence, a branch and bound procedure was developed to minimize the total passenger delay when an aircraft is grounded. Further, Teodorović and Stojković [3] have introduced a lexicographic dynamic programming approach and a heuristic algorithm to minimize the total passenger delays and the number of cancelled flights. Two minimumcost network flow models to address aircraft shortages have also been introduced by Jarrah et al. [4]. The first is a delay model that allows for flights to be delayed or cancelled, aircraft to be alternated among scheduled flights, or the usage of surplus aircraft to be requested. The second model minimizes the total flight leg cancellation costs by similarly selecting multiple cancellations, employing aircraft substitution, and utilizing spare aircraft. However, accurate flight delay costs can not yet be quickly and easily determined, and the scope and depth of the flight disruption are not considered in the above models. In contrast, delay costs can be approximated with accuracy using the timeband approximation model introduced by Yu et al., which was developed in order to retain the superiority of the resource assignment model and the multicommodity flow model. Further, Yu [5] has introduced a timebased twocommodity network optimization model for aircraft routing. One commodity in the network model represents aircraft, while the other represents flight cancellations. In addition, Yan and Yang [6] have published a single model to initially combine flight cancellations and delays and to ferry flights. Using a basic timespace network representation, Yan and Lin [7] have developed network flow models to handle station closures, while Yan and Tu [8] have proposed a multiplefleet network flow model. Finally, Cao and Kanafani [9, 10] have developed a minimumcost network flow aircraft recovery model, along with network flow algorithms to determine solutions.
However, the mathematical solution derived from a timebased network does not represent the delay costs accurately. In order to develop an optimization model that can determine delay costs with greater accuracy than the multicommodity flow model, the timeband optimization model for aircraft routing reconstruction in response to groundings and delays experienced over the course of a day has been introduced [11, 12]. A mathematical model that captures the delay costs for each flight based on a timebased network has been developed; the motivation for this approach is to exploit the underlying network structure associated with aircraft routings, with the mathematical model being a singlecommodity integral minimumcost flow problem with side constraints. The concept behind this model is to partition the time horizon into time bands with a fixed length. The nodes in the resulting network represent possible activity at a station. Thus, in the established network, some nodes are positioned at fixed intervals on the time axis. The aggregation of activities for that station can then be attained in operations, where is the number of stations, is the number of flights, and All aircraft that arrive at a station within the time period represented by a node may then be connected to any flights scheduled to depart from that station. The delay time can be approximated with an accuracy that is dependent on the length of the time segment; therefore, the delay costs can be determined quickly.
When disruptions occur, the airline controllers in the airline operations control centers (AOCC) adjust the scheduled operations by delaying flight departures, cancelling flight legs, rerouting or swapping aircraft, summoning new crew (or reassigning existing crew), and so forth. In order to reassign resources and adjust flight schedules in the most economic manner, certain factors must be acknowledged: the flight delay costs, flight cancellation costs, aircraft substitution costs, aircraft assignment, crew schedules, maintenance requirements, aircraft positions, and so forth. However, it is difficult to determine the delay costs exactly. Further, flight cancellation costs involve both passenger opportunity costs and direct costs, and the delay time cannot be determined in advance. Therefore, the ability to compute both delay and flight cancellation costs is very important as regards the solution of flight operations recovery problems. Such solutions determine adjusted flight departure times, cancellations, and revised routing for affected aircraft.
Many researchers have conducted significant work with regard to flight operations recovery considering passenger factors. For example, Bratu and Barnhart [13] have proposed two flight operations recovery optimization models, the Disrupted Passenger Metric (DPM) and the Passenger Delay Metric (PDM) models, the objective of which is to find the optimal tradeoff between airline operation costs and passenger delay costs. The resultant recovery plans work to select flight departure times and flight cancellations, to assign aircraft, to reserve crews for each flight leg as necessary, and to ensure the constraints regarding crew regulations and aircraft maintenance requirements are considered. The objective of the DPM model is to minimize the sum of the operating and disrupted passenger costs. In contrast, the PDM model is constructed to more accurately compute delay costs. Thus, a special procedure to compute the costs due to disrupted passengers has been developed. In addition, Hu et al. [14] have proposed the integrated integer programming model, which is based on an approximate reduced timeband network and a passenger transiting relationship. The objective is to minimize the total cost associated with reassigning aircraft and passengers to flights. A feasibility analysis for the problem is conducted in order to obtain the conditions under which aircraft and passenger recovery is possible. Further, Zhao et al. [15] have introduced an approach in which the flight delay costs are subdivided, proposing the concept of passenger disappointment spillover cost, and demonstrated a mathematical model of an irregular flight schedule recovery problem that reschedules flight legs and reroutes aircraft by minimizing an objective function involving the flight delay costs or the delay time. The airline can choose one of selections of objective functions according to their needs. However, determining the passenger disappointment spillover cost function remains a problem, and the compensation provided by the airline is a very important factor influencing this cost function. For further research on recovery considering passenger factors, see [16–19]. Castro et al. [19] give a new approach for disruption management in airline operations control. They provide a Multiagent System Architecture that represents the AOCC, and introduce a method to calculate the quality costs from the passenger point of view and a system classification to classify disruption management systems. They have done much original work in the recovery operational costs; especially they first propose the quality operational costs (the cost that each minute of delay and flight cancellation has from a passengers point of view). They try to quantify in a monetary unit the cost that each minute of delay and flight cancellation has from a passengers point of view and, at the same time, use that value together with the direct operational costs to improve the decision making during the disruption management process. The quality costs of a flight are given out under the passenger goodwill; what is more, they propose the airline company case study of TAP Portugal, regarding the delay cost from the point of view of the passengers. However, we can not directly use a similar method to calculate the passenger recovery costs in China, in that the preconception of passenger’s compensation costs comes from Confucius’s idea of the average distribution two thousand years ago and some exemption rules from Civil Aviation Administration of China, such as the weather and other objective reasons, and the passengers can not get any compensation.
When one or several aircraft are grounded, a heuristic approach based on a greedy randomized adaptive search procedure (GRASP) for the reconstruction of aircraft routes, which has been developed by Argüello et al. [20], can be applied. (Further optimization methods are provided in [21–25].) However, a large proportion of the feasible optimization methods copy flights to the entire flight schedule, while other resources are reassigned in order to repair the disrupted schedules in the best possible manner; thus, the disruption scope or depth may be wide or great. Further, the computation time increases with increases in the number of potentially disrupted flights, when such disruptions occur in hubandspoke airline networks.
Petersen et al. [26] have introduced a disruptable flight set in an attempt to control the recovery search scope of an optimization approach to airline integrated recovery. This approach yields the first computed results for the fully integrated problem. In fact, the disruptable flight set may be comprised of all flights and the actual scope and depth of the disruption may not be controlled. In addition, Lan et al. [27] have introduced aircraft routing and flight departure time adjustment in order to reduce the number of disrupted passengers and to minimize the number of passenger misconnections. However, some passengers may have a long delay time or the flight scope may be wide. The computation time and the disruption scope and depth are not yet fully controllable.
At the 2012 Informs meeting, we first proposed the simplex group cycle approach as a means of controlling the potential disruption scope [28]. The approximated delay costs considering the random flying time around the planned flying time is first introduced in the literature [29], and the timeband approximation model on flight operations recovery problems is constructed based on the timeband network and approximated delay costs. Further, a series of numerical experiments show that the timeband approximation model on flight operations recovery with the random flying time may be more efficient than the timeband approximation model on flight operations recovery with the planned flying time. Using the thought of stochastic delay costs, the approximated cancellation costs formula is first introduced in this paper. In addition, the formula takes into account interval random number, the number of passengers, and approximate city accommodation costs. However, we did not introduce the feasible algorithm for aircraft and flight partitioning, the feasible algorithm to compute weighted value and simulation tests, and so forth, in our previous work.
Previously, Clausen et al. [30] have shown that robust planning or schedule robustness is an interesting research topic closely related to disruption management. The goal of robust planning is to reduce the sensitivity of flight and crew schedules and aircraft rotations to disruptions. Rosenberger et al. [31] have also proposed a robust fleetassignment model with hub isolation and short cycles and proven that fleet assignments with limited hub connectivity have a greater number of short cycles. A cycle is a sequence of flights that begins and ends at the same airport; further when a cycle contains fewer flights, this cycle is called short cycle. And it is possible to reduce the impact of a disruption at one hub on other hubs by isolating the hubs insofar as possible; hub isolation is zero hub connectivity which is the number of legs in a rotation that are in a route that begins at a hub and ends at a different hub. In fact, the aim of hub isolation and short cycles is to control the difficulty and scope of schedule recovery. However, limited hub isolation corresponds to a lower number of feasible schedule recovery solutions, while short cycles are associated with flight cancellations or longer delay times within the short cycles. The scope of schedule recovery is controlled within the short cycle, rather than by the AOCC controllers.
Moreover, airline network structure has an important effect on flight operations recovery. In America, the AOCC controllers formulating flight operations recovery plans must consider passenger itineraries within a hubandspoke network structure. A hubandspoke network [19] means a system of connections arranged like a chariot wheel, in which all traffic moves along spokes connected to the hub at the center. However, passenger itineraries of domestic flights in China containing a single flight leg are in the majority, and the international flights accounted for about 9% of all flights in 2015 according to flights data from the Civil Aviation Administration of China. Further, international flights and domestic flights are different in aircraft type and crews, so there is little implication between international flights and domestic flights. Thus, we do not consider passenger itineraries in this study on domestic flights in China. What is more, multihubandspoke airlines are associated with multibase airlines in China. Therefore, the hubandspoke airlines have some advantages as regards options for schedule recovery from irregular operations.
This article makes three new contributions compared to previous results. First, we describe more specifically the simplex group cycle approach and first present a feasible algorithm for aircraft and flight partitioning and a feasible algorithm to compute weighted value. Second, we use a weighted threshold to control the disruption scope and depth, computation time and recovery value for flight operations recovery for hubandspoke airlines. Hence, we describe the weighted timeband approximation model (WTBAM) for flight operations recovery using the simplex group cycle approach, referring to the concept of the timeband approximation model for flight operations recovery [12]. Third, we design some practical numerical experiments and the results demonstrate that use of the WTBAM incorporating the simplex group cycle is reasonable and feasible, and we provide first the feasible method to choose the weighted threshold . Furthermore, the results verify that the gap between the relaxed weighted timeband approximation model and the weighted timeband approximation model is zero under certain conditions. Finally, the results demonstrate that we can obtain a satisfactory recovery plan considering the flight disruption scope and depth, computation time and recovery value under some appropriate rules.
The remainder of this paper is organized as follows. In Section 2, we introduce some fundamental concepts, including important definitions and key aspects of the timeband approximation model. In Section 3, we present the mathematical formulations of the WTBAM. The model is then applied to actual data from two airlines and the computational results are presented in Section 4. In Section 5, the conclusions are presented with a brief discussion.
2. Fundamental Concepts
To construct the WTBAM, we first introduce several fundamental concepts.
2.1. Objective and Constraints
The objective of flight operations recovery is to make decisions that correspond to the lowestcost response to restoring the disrupted schedule to the planned schedule. In this paper, we firstly consider flight departure delays or flight leg cancellations to repair a flight schedule after a disruption and subsequently assign individual aircraft to a new flight schedule by the method [32]. Moreover, considering the characteristics of multihubandspoke airlines in China, we have no considering maintenance base for the overnight restrictions in the process of assigning individual aircraft to a new flight schedule. In fact, the AOCC controllers customarily give priority to assign the same type aircraft, and in general there are a lot of the same type aircraft in the multihubandspoke airlines. Further, the different type of aircraft will bring about a series of problems, such as crew, airport, and maintenance; the AOCC controller has to make sure the aircraft routing is feasible. Therefore, the relaxation of maintenance base restrictions is one of the advantages for multihubandspoke airlines in China.
Note that minimization of weighted delay costs and flight cancellation costs is a common objective. The delay costs include additional fuel expenses, overtime payments, further expenses associated with overtime payments when an aircraft is delayed in an airport, and passenger food expenses. Thus, delay costs are the direct costs incurred by flight delays, rather than associated costs that arise by chance. We consider the different weighted values of these costs in different simplex group cycles. In contrast, flight cancellation costs include loss of profit, passenger food and accommodation expenses, and passenger compensation. And flight cancellation costs were generated by the approximate estimation formula in this research, is the number of passengers, is a random number between 80 and 100 and generally includes meals and passenger compensation, and is the approximate city accommodation cost. Further, the recovery plan must also satisfy the following constraints in order to be considered feasible. First, the turnaround time must be more than the minimum threshold, that is, 40 min. Second, every flight on each aircraft route must depart from the station at which the immediately preceding flight arrived. Finally, the aircraft balance necessary for the flight schedule to be resumed the next day must be achieved.
2.2. Simplex Group Cycle
The network of a multihubandspoke airline is comprised of two or more hubandspoke structures. Thus, multihubandspoke airlines have more maintenance bases and a greater number of feasible recovery plans. In fact, the group of flights between one or more hubandspoke systems in the same airline network is defined as belonging to a simplex group cycle and labelled , where denotes the specific simplex group cycle. Those flights are assigned to certain fixed aircraft, and the originating and destination cities contain one or more hubandspoke systems. Thus, each flight belongs to only one . According to the definition of a simplex group cycle, a simplex group cycle partition for all aircraft can be determined, with each aircraft being assigned to one of various groupings. The simplex group cycle approach and its superiority are shown in the following example.
The Air China hubandspoke systems shown in the following flight schedule are Beijing (PEK), Shanghai (SHA), and Chengdou (CTU). The network comprised of the PEK, SHA, Xian (XIY), and CTU hubs is shown in Figure 1.
The flight schedule and aircraft assignments are shown in Table 1. Thus, we can determine a partition based on the natural aircraft sequence. The first simplex group cycle is the set comprised of CA4193, CA4102, CA1216, CA1215, CA1945, CA1946, CA425, and CA1950, and the related aircraft are , , and . Thus, the second simplex group cycle is the set comprised of CA1203, CA1204, CA4107, CA4108, CA4201, and CA4202, with the related aircraft being and .
 
Flight : flight number; Tail : tail number; Dep. : departure time; Arr. : arrival time; Ass. A: assigned aircraft; C. cost: cancellation cost. 
The partitioning of the flights or aircraft is usually determined by the AOCC controllers, according to their empirical data and objectives. However, partitioning the flights or aircraft in a manner that establishes the optimal groupings is a difficult problem. Therefore, we first present a feasible algorithm for aircraft and flight partitioning in Algorithm 1.

2.3. Weighted Set and Potential Disruption Set
When disruptions occur in a hubandspoke network, they may lead to a large amount of disruptions overall. Thus, schedule recovery should be accomplished in the most economic manner. However, if we only consider the flights themselves, the disruptions may have a negative influence on passengers. Therefore, in order to control the potential disruption scope in a hubandspoke network and the associated negative influence on passengers, can be determined based on the flight disruption scope and depth, the computation time, and the recovery value. Therefore, we employ the concept of a weighted set, as explained below.
Let be the set containing the simplex group cycle which is at least one schedule disruption in . And let be the aircraft partition related to , the set of all simplex group cycles. Then is the aircraft partition related to .
Definition 1. When a disruption occurs, the weighted value is denoted
When the flights of the th simplex group cycle influence those of the th simplex group cycle, the value of can usually be assigned by the AOCC controllers in accordance with their empirical data and objectives. Naturally, we define , , and , Here, is the amplification coefficient, and and are the average available flight times of each and , respectively. Then, the weighted set .
The simplex group cycle approach can exploit the advantages provided by a multihubandspoke network. For example, in Table 1, when aircraft (CA4193) departed from CTU arrived in PEK, the maintenance personnel detected mechanical problems. It was estimated that a period of 4 h was required to complete the repairs. Thus, using the simplex group cycle approach, we utilize to execute flight CA4102, followed by flight CA1950, with being utilized for flights CA1945 and CA1946. At 13:30, then executes flights CA1215 and CA425. The total delay time is 100 min, and we save 60 min. Note that only the disruption scope is controlled by the simplex group cycle approach. In addition, the developed recovery plan involved cooperation between China Airlines () and its branch company, China Southwest Airlines ().
Based on a weighted threshold , AOCC airline controllers can obtain a potential disrupted set , which is defined as follows.
Definition 2. The potential disrupted flights set is the set of flights in the simplex group cycle obtained when the condition , , is satisfied.
Hence, the AOCC airline controllers can set the potential disruption flights conveniently by considering the value of . However, it is difficult to compute the values for flights or aircraft. Therefore, a feasible algorithm to compute and is first given in Algorithm 2. The timeband network incorporating is constructed using the timeband network transformation procedure found in [11, 12]. The timeband network mainly consists of stationtime nodes, stationsink nodes, and flight arcs, so we refer to some concepts of the literature [11, 12] in the following Flight arcs: arcs representing specific scheduled flights. Stationtime nodes: nodes representing aggregate activity for a specific time segment at a station. Stationsink nodes: nodes representing the end of the recovery period at stations.

Subsequently all the possible delay costs for the flight arcs can be computed by considering the difference between the real takeoff time and the scheduled plan in the timeband network. Naturally, the WTBAM is constructed in accordance with the optimization model by considering a series of flight arcs with some potential flight delays. In the next section, we present the mathematical formulation of the timeband approximation model.
3. TimeBand Approximation Model Mathematical Formulation
In order to construct the mathematical formulation of the optimization model, some indices, sets, parameters, and variables are employed, which are given in Nomenclature.
Using and considering the recovery plan for only, we can express the WTBA approximation model as follows:The objective function (3) minimizes the delay and cancellation costs for . Constraint (4) shows that all flights in the simplex group cycle must be operated or cancelled. Further, constraint (5) satisfies the aircraft utilization requirement at the stationtime node. The stationsink node flow constraint (6) enforces aircraft balance. The binary and ensure that a flight can only be operated or cancelled; that is, it must not be partially serviced and then cancelled. The possible arcs must be integrated, which is enforced by constraint (9).
Remark 3. When , the solution of WTBAM is obtained after searching the flights in , and the passengers affected by this solution may experience long flights, or some flights in the simplex group cycle may be cancelled. In other words, the disruption scope can be controlled in the simplex group cycle , but the disruption depth may be very large.
In contrast to the optimization model with , the objective function and constraints in the model can be modified by . In this optimization model, is restricted to the disrupted flight set only. We can obtain the optimization model with , as follows:The objective function (10) minimizes the delay costs and the cancellation costs in . Further, constraint (11) shows that all flights in must be operated or cancelled. Constraint (12) satisfies the aircraft utilization requirement for at the stationtime node. The stationsink node flow constraint (13) enforces aircraft balance in . The binary and ensure that a flight can only be operated or cancelled; that is, it must not be partially serviced and then cancelled. Finally, the copy arcs must be integrated, which is enforced by constraint (16).
Remark 4. When has the maximum weighted value , corresponds to all of the flights. In that case, it may be very difficult for the AOCC to quickly obtain the optimal solution or the disrupted flight scope may be very broad. It may be very difficult to implement such a recovery plan for actual flight production, because it does not satisfy the requirements of timeliness and balance between the disruption scope, disruption depth, and recovery value.
For , it is difficult for the AOCC to determine the optimal that satisfies the balance between the passenger and airline requirements. However, the AOCC can choose the best satisfactory value under those conditions.
Remark 5. When a disruption occurs, the weighted value is given out by the Algorithm, Then
Proof. Because all the flights are contained within , when , , Therefore, we can obtain In fact, must be less than or equal to the total number of aircraft . And then the linear programming relaxation of the weighted timeband approximation model (LPRWTBAM) is expressed as follows:Naturally, and , , , . We can set Further, we assume that and satisfy the following properties.
WellOrdered Property. The nonzero value coefficient set (NVCS)is well ordered according to .
In fact, the weighted value, delayed time, and cancellation cost are often different in the NVCS; therefore, the following sequence properties may be satisfied in most casesNaturally, the NVCS frequently exhibits a wellordered property.
Let be the number of elements in the set. Then, we can obtain the following property in the majority of cases.
Possible Zero Property. or 1, , . In fact, when there is only a small number of available aircraft in , the aircraft utilization is generally reasonably high and flights cannot take off early; then or 1 often occurs.
Remark 6. If the feasible solution of WTBAM is nonempty, and are integers, and NVCS exhibits a wellordered property, then has the possible zero property, , . Then, the WTBAM and LPRWTBAM solutions are identical.
Proof. Let be the LPRWTBAM solution and be the WTBAM solution. Then, because the LPRWTBAM is the relaxed WTBAM, If is an integer, thenThe analysis of is as follows.
When or 1, , we adopt the assumption Then is well ordered. According to , , , , (19), and (20), we can obtainFurthermore, according to (21) and (22) and as and are integers, we know that , where is an integer. Finally, we conclude that is an integer.
4. Computational Results
Two sets of computational results are presented in this section: (i) the computational results obtained by application of the WTBAM to the sample data given in Table 1, which demonstrate that use of the WTBAM incorporating the simplex group cycle is reasonable and feasible; (ii) results of a test using data from actual aircraft belonging to Sichuan Airlines in China, which involves a hubandspoke network. The analysis results obtained for different weighted thresholds and different simplex group cycles under four scenarios are given.
4.1. Sample Flight Schedule Table 1
Mechanical problems are detected for aircraft at 9:50, and it is estimated that a 4h period is required for the repair work. In this example, the time band is set to half an hour and the schedule recovery is implemented from 10:00 to 02:00 the next day. In this example, the simplex group cycle partitions are set to and , in order to show the AOCC controllability. Note that these partitions do not correspond to the partition implemented using the simplex group cycle procedure. The weighted value is determined according to the weighted set procedure. We set the disruption threshold ; then the potential disruption flights set is as shown in Table 2.

Subsequently, the flight arc delay costs for the recovery problem are determined using the timeband transformation network [12]. The timeband transformation contains 26 nodes and 57 flight arcs. The delay cost for each flight arc is calculated according to the arc delay, and a 20permin delay cost function is assumed [12] in this paper. The nonzero flight arc delay costs for and 2.0 are listed in Tables 3 and 4, respectively. Finally, we can compute the WTBAM solution using Lingo 14.0. The solutions for and 2.0 are listed in Tables 5 and 6, respectively. The results obtained for the timeband transformation with contain more nodes and flight arcs than that for ; therefore, the WTBAM solution (with ) may yield a superior proximal point than the standard timeband transformation with . Further, the optimal value yielded by the WTBAM is smaller than that for the timeband transformation with .


 
OST: origin stationtime; DST: destination stationtime; D. cost: delay cost. 

4.2. Stochastic Experiments
Our model was also tested using actual data obtained for a specific day in 2015, for a hubandspoke airline network under Sichuan Airlines having approximately 254 daily flights. The WTBAM was implemented in C++ using CPLEX 12.3 on a Lenovo notebook with a 2.16GHz processor. As the arrival time of the disrupted flight was 02:30 on the next day and recovery time interval should contain all the flight schedule, we set the schedule recovery from 07:00 to 03:00 on the next day and the time band to 5 min in order to reduce the gap bound. Then, 1–4 aircraft having mechanical problems or some issue prohibiting them from use between 07:00 to 03:00 on the next day were randomly generated; these events were labelled scenarios 1–4. Combined with the experience of AOCC, the simplex group cycle should not be very large but should not be very small, 60 aircraft in the actual data, roughly divided into 6–15 group. If the simplex group cycle is very large, it may lead to the fact that the calculation time of optimal solution is too long and the potential disruption scope is too large, while if it is too small it may lead to a large potential disruption depth. So we obtain the 7, 12, and 15 three simplex group cycles by setting in the Algorithm. Subsequently, three simplex group cycle partitions are listed in Tables 7–9, and due to the limit of table length, only apart results of WTBAM under these partitions and different disruption thresholds are listed in Tables 12–15. In addition, the results for the simplex group cycle partitions, the scenarios, and the weighted values are presented in Figure 2.
 
SGC : Simplex group cycle number. 


(a) Weighted value under scenario 1
(b) Weighted value under scenario 2
(c) Weighted value under scenario 3
(d) Weighted value under scenario 4
Note that DT 60, DT 30, DT 15, and DT 0 indicate that the flight delay time is more than 60 min, between 30 min and 60 min, between 15 min and 30 min, and between 0 min and 15 min, respectively. Further, GTSN time (s) represents the time required to generate the timespace network in seconds, while CPU time represents the WTBAM calculation time.
Comparison of the data in Tables 12–15 leads to the conclusion that the potential disruption scope , GTSN time, and CPU time increase with increased . Subsequently, by the results in Tables 12–15, we can obtain that GTSN time, CPU time, and potential disruption scope satisfy the function in Table 11 and Figure 3. With the increase of , GTSN time and CPU time are exponential growth and the function relationship is is the real set. In view of this test, GTSN time and CPU time of second simplex group cycle partitions are generally better than the other simplex group cycle partitions. Note that GTSNij and CPUij indicate the function of GTSN time and potential disruption scope in the case of scenario and th simplex group cycle partitions, and .
However, the optimal integer solution primarily decreases with increased , and the gap between the optimal integer and the optimization solution is zero in the examples. Further, when the optimal integer reaches a certain level, it does not necessarily decrease or decreases only slightly, in most cases. These results show that a larger is not necessarily better. In fact, the real extent of the flight disruption decreases with increased ; however, note that this factor, disruption scope, and depth do increase with increased in some cases. Further, the nodes and flight arcs increase with increased , along with the CPLEX 12.3 computational time.
The AOCC controllers can choose the value of by considering the balance between the real flight disruption scope and depth, the CPU time, and the optimal integer value. Subsequently, we first introduce the feasible method to choose the value of (Choose Method) in the following. The AOCC controllers choose the appropriate value of which should be consistent with the four rules in the following.
Choose Method The first rule is to meet the time limit for flight operations recovery. The second rule is about the flight recovery value in the following formula. Let be , be the corresponding threshold of , and be the corresponding threshold of . The third rule is as small as possible to the impact of disruption flights scope and depth. The fourth rule is to ensure the uniqueness flight recovery plan. If there are two flight recovery plans to meet the previous three rules, then select a smaller flight recovery value. And if the flight recovery value is also the same, then select the CPU time which is less than that of the flight recovery plan.
In this paper, we choose which means that the two flight recovery plans can be considered as equivalent when the difference ratio of the two flight recovery values is less than or equal to 0.25, because the AOCC controllers are not just considering about flight recovery values. The airlines can adjust to meet their management objectives. Subsequently, we provide a feasible method to measure the impact of disruption flights scope and depth according to the actual situation in China. Firstly, the weighted value of common delay time band is determined by the AOCC controllers, and then we calculate the total impact value. The delay time band and the weighted value are given out in Table 10. In particular, the weighted value can be adjusted according to the requirements of the airlines.

