Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2016 / Article

Research Article | Open Access

Volume 2016 |Article ID 3201490 | https://doi.org/10.1155/2016/3201490

Haiwen Xu, Songchen Han, "Weighted Time-Band 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 Time-Band Approximation Model for Flight Operations Recovery considering Simplex Group Cycle Approaches in China

Academic Editor: Leonid Shaikhet
Received15 Apr 2016
Revised30 Jul 2016
Accepted18 Sep 2016
Published27 Oct 2016

Abstract

The time-band 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 hub-and-spoke 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 hub-and-spoke 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 time-band 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 time-band 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 management-level 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 minimum-cost 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 time-band 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 time-based two-commodity 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 time-space network representation, Yan and Lin [7] have developed network flow models to handle station closures, while Yan and Tu [8] have proposed a multiple-fleet network flow model. Finally, Cao and Kanafani [9, 10] have developed a minimum-cost network flow aircraft recovery model, along with network flow algorithms to determine solutions.

However, the mathematical solution derived from a time-based 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 time-band 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 time-based 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 single-commodity integral minimum-cost 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 trade-off 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 time-band 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 [1619]. 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 [2125].) 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 hub-and-spoke 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 time-band approximation model on flight operations recovery problems is constructed based on the time-band network and approximated delay costs. Further, a series of numerical experiments show that the time-band approximation model on flight operations recovery with the random flying time may be more efficient than the time-band 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 fleet-assignment 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 hub-and-spoke network structure. A hub-and-spoke 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, multi-hub-and-spoke airlines are associated with multibase airlines in China. Therefore, the hub-and-spoke 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 hub-and-spoke airlines. Hence, we describe the weighted time-band approximation model (WTBAM) for flight operations recovery using the simplex group cycle approach, referring to the concept of the time-band 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 time-band approximation model and the weighted time-band 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 time-band 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 lowest-cost 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 multi-hub-and-spoke 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 multi-hub-and-spoke 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 multi-hub-and-spoke 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 multi-hub-and-spoke airline is comprised of two or more hub-and-spoke structures. Thus, multi-hub-and-spoke airlines have more maintenance bases and a greater number of feasible recovery plans. In fact, the group of flights between one or more hub-and-spoke 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 hub-and-spoke 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 hub-and-spoke 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 Tail From To Dep. Arr. Ass. A C. cost

CA419311 CTU PEK 7:00 9:30 11380
CA410212 PEK CTU 11:30 14:00 15150
CA121621 SHA PEK 8:00 10:00 12985
CA121522 PEK SHA 12:30 14:30 16991
CA194531 SHA CTU 8:05 10:35 11260
CA194632 CTU SHA 12:20 14:50 15785
CA42533 SHA CTU 17:35 20:05 13031
CA195034 CTU SHA 21:00 23:30 13820
CA120341 PEK XIY 07:05 9:00 11650
CA120442 XIY PEK 10:00 11:55 10256
CA410751 CTU PEK 09:00 11:30 11000
CA410852 PEK CTU 15:00 17:30 15236
CA420153 CTU XIA 20:45 22:00 9325
CA420254 XIA CTU 23:00 00:15 9128

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.

Input:     = Original flight schedule
   Aircraft = All the executive aircraft associated with the flights
Output: Airport = An ordered list of all airports according to the number of related flights
      = The aircraft group in simplex group cycle
      = The flight group in simplex group cycle
begin
 for each airport
  Count the related flights
  mark;
  Sort according to ;
end    /According to the node degree at the airport-flight adjacency graph, we sort the airport from large to small./
  ;     /where is the approximate simplex group cycle number/
   = The related aircraft in airport ;
for
   = The related aircraft in airport
  ; /Ensure that each aircraft is in the only one
  ;
  ;
end
if ;
   = The related aircraft in airport
  ;
  ;   ;
else output .
end
  output ,
   = The flight group comprised of all flights in all simplex group cycles;
end
2.3. Weighted Set and Potential Disruption Set

When disruptions occur in a hub-and-spoke 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 hub-and-spoke 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 multi-hub-and-spoke 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 time-band network incorporating is constructed using the time-band network transformation procedure found in [11, 12]. The time-band network mainly consists of station-time nodes, station-sink nodes, and flight arcs, so we refer to some concepts of the literature [11, 12] in the followingFlight arcs: arcs representing specific scheduled flights.Station-time nodes: nodes representing aggregate activity for a specific time segment at a station.Station-sink nodes: nodes representing the end of the recovery period at stations.

Input:
Output: , potential disrupted aircraft set and
begin
 for each
  compute the average available flight time in the ;
  / is the set of aircraft in the simplex group cycle
 end
 for each simplex group cycle
  for each simplex group cycle
   compute
  end
 end  / is the aircraft partition related to .
    is the aircraft partition related to
for each .
  compute
end
;
;
;
 for
 if is the threshold./
  ;  ;
 end
 end
end

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 time-band 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 time-band approximation model.

3. Time-Band 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 station-time node. The station-sink 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 station-time node. The station-sink 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 time-band approximation model (LPRWTBAM) is expressed as follows:Naturally, and , , , . We can set Further, we assume that and satisfy the following properties.
Well-Ordered 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 well-ordered 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 well-ordered 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 hub-and-spoke 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 4-h 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.


Tail From To Dep. Arr. Ass. A C. Cost

12 PEK CTU 11:30 14:00 15150
22 PEK SHA 12:30 14:30 16991
32 CTU SHA 12:20 14:50 15785
33 SHA CTU 17:35 20:05 13031
34 CTU SHA 21:00 23:30 13820
42 XIY PEK 10:00 11:55 10256
52 PEK CTU 15:00 17:30 15236
53 CTU XIY 20:45 22:00 9325
54 XIY CTU 23:00 00:15 9128

Subsequently, the flight arc delay costs for the recovery problem are determined using the time-band transformation network [12]. The time-band transformation contains 26 nodes and 57 flight arcs. The delay cost for each flight arc is calculated according to the arc delay, and a 20-per-min 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 time-band 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 time-band transformation with . Further, the optimal value yielded by the WTBAM is smaller than that for the time-band transformation with .


Tail Origin node Destination node Arc delay cost

12 3 7 3000
21 3 6 1800
32 4 8 3200
32 7 9 6200
32 10 14 10400
32 11 14 11600
33 8 11 1100
33 9 13 4100
34 11 14 1200


Tail Origin node Destination node arc delay cost

12 4 10 1200
12 5 11 1800
12 6 13 1800
22 5 9 600
22 6 12 3420
32 7 14 3200
32 10 16 4400
32 11 17 5000
32 13 18 6200
32 15 20 7400
32 19 27 10400
32 21 27 11600
32 22 27 12800
32 24 27 14000
33 14 21 1100
33 16 22 2300
33 17 24 2900
33 18 26 4100
33 20 29 5300
34 21 27 1200
34 22 27 2400
34 24 27 3600
42 23 28 29640
42 25 28 31920
53 19 23 570
53 21 25 1800
53 22 30 5130
53 24 30 6270
53 26 30 8550
54 25 29 2280



Tail From To OST node DST node Ass. A D. cost C. cost

12 PEK CTU 2 4 0
21 PEK SHA 1 5 0
32 CTU SHA 15758
33 SHA CTU 5 10 0
34 CTU SHA 4 14 0

OST: origin station-time; DST: destination station-time; D. cost: delay cost.


Tail FromToOST nodeDST nodeAss. AD. costC. cost

12 PEK CTU 2 7 0
22 PEK SHA 4 8 0
32 CTU SHA 3 8 0
33 SHA CTU 8 19 0
34 CTU SHA 19 27 0
42 XIY PEK 1 5 0
52 PEK CTU 6 15 0
53 CTU XIY 15 23 0
54 XIY CTU 23 29 0

4.2. Stochastic Experiments

Our model was also tested using actual data obtained for a specific day in 2015, for a hub-and-spoke 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.16-GHz 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 79, and due to the limit of table length, only apart results of WTBAM under these partitions and different disruption thresholds are listed in Tables 1215. In addition, the results for the simplex group cycle partitions, the scenarios, and the weighted values are presented in Figure 2.


SGC Tail

16 17 18 21 31 32 33 3539 42
1 22 2526 27 30 3848 52 54 55 5760
2 3 8 910 11 34 40 41 43 44
29 36 37 47 4953
5 615 20 5658 59
4 712 19 2324 46
1314 28 45 50 51

SGC : Simplex group cycle number.

SGC Tail

17 18 32 35 39 42
22 25 2638 57
2 89 1011 34 40 41
29 36 37 49 53
5 615 20 5659
12 2324 46
2743 52 58
1 30 33 47 48
13 14 2845 50 51
3 16 44 5560
4 19 21
7 3154


SGC Tail

18 32 35 39 42
22 25 2638
2 89 1011 34 40 41
29 36 37 49
5 615 20 5659
23 24
4358
1 30 47 48
14 28 4551
3 16 44 5560
4 19 21
7 3154
12 46 50
13 52 57
15 1727 53

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 time-space network in seconds, while CPU time represents the WTBAM calculation time.

Comparison of the data in Tables 1215 leads to the conclusion that the potential disruption scope , GTSN time, and CPU time increase with increased . Subsequently, by the results in Tables 1215, 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 MethodThe 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.


Delay time band
DT 0 DT 15 DT 30 DT 60 Cancelled flight

Weighted value 1 2 4 6 8


First SGP Second SGP Third SGP First SGP Second SGP Third SGP

Scenario 1Scenario 2

GTSN
CPU

Scenario 3Scenario 4

GTSN
CPU