• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Airline Disruption Management:A Review of Models and Solution Methods

    2021-08-30 07:05:08YiSuKexinXieHongjinWngZheLingWnprhArtChovlitwongsePnosPrlos
    Engineering 2021年4期

    Yi Su ,Kexin Xie ,Hongjin Wng ,Zhe Ling ,*,Wnprh Art Chovlitwongse ,Pnos M.Prlos

    a School of Economics and Management,Tongji University,Shanghai 200092,China

    b Xiamen Airlines,Xiamen 361006,China

    c Department of Industrial Engineering and Institute for Advanced Data Analytics,University of Arkansas,Fayetteville,AR 72701,USA

    d Center for Applied Optimization,Department of Industrial and Systems Engineering,University of Florida,Gainesville,FL 32611,USA

    Keywords:

    Disruption management Aircraft recovery problem Crew recovery problem Integrated recovery problem

    ABSTRACT The normal operation of aircraft and flights can be affected by various unpredictable factors,such as severe weather,airport closure,and corrective maintenance,leading to disruption of the planned schedule.When a disruption occurs,the airline operation control center performs various operations to reassign resources(e.g.,flights,aircraft,and crews)and redistribute passengers to restore the schedule while minimizing costs.We introduce different sources of disruption and corresponding operations.Then,basic models and recently proposed extensions for aircraft recovery,crew recovery,and integrated recovery are reviewed,with the aim of providing models and methods for different disruption scenarios in the practical implementation of airlines.In addition,we provide suggestions for future research directions in these topics.

    1.Introduction

    The aviation industry has become a major player in the global economy.As people become increasingly dependent on air travel,the number of scheduled flights worldwide grows every year.According to the International Air Transport Association(IATA)[1],civil aviation passenger demand increased by 4.2%while capacity increased by 3.4%worldwide in 2019.Hence,demand appears to grow faster than capacity,suggesting great development potential for the aviation market.With the development of the aviation industry,airline planning and scheduling problems have attracted much attention,and most airlines benefit from advanced optimization methods.Sophisticated models and effective solutions have been developed for each stage of planning,as reported in detailed overviews by Eltoukhy et al.[2]and Zhou et al.[3].

    Under various circumstances,such as aircraft mechanical problems or severe weather conditions,flights cannot be operated as planned.Data provided by the Bureau of Transportation Statistics(BTS)show that approximately 21%of flights in the United States during 2019 experienced more than 15 min of arrival delay.Similarly,the average flight on-time rate in China was 81.43%according to the Civil Aviation Administration of China(CAAC)’s Statistical bulletin of civil aviation industry development in 2019[4].Likewise,the Punctuality League report by the OAG[5]shows that only three airlines achieve a greater than 90%on-time rate worldwide.Flight delays and cancellations have become important factors that affect passengers’airline preferences.Consequently,if irregular operation cannot be properly recovered,the economic and social benefits for airlines may be undermined.Therefore,disruption management has become a major problem in airline operation management.

    There are many differences between airline planning and airline recovery.Planning focuses on optimization,whereas recovery targets a feasible yet possibly suboptimal solution that can be obtained in real time.Furthermore,recovery may be more uncertain than planning,depending on the degree and type of disruption.Unlike flight planning,which can be obtained several months before operation,recovery solutions should be obtained and implemented as quickly as possible after a disruption.In this paper,we focus on the characteristics of recovery and review common models and methods.Details of the models provide the characteristics and applications of recovery options that help airlines to choose appropriate methods to meet their specific requirements.

    The remainder of the paper is organized as follows.Section 2 presents an introduction of disruption management,including sources of disruption and possible recovery operations.Basic models and solution methods for aircraft recovery are provided in Section 3.Section 4 presents the general model and some extensions of crew recovery.Section 5 presents integrated recovery considering multiple resources,including passengers.Section 6 summarizes our work and provides directions for future work.

    2.Disruption management

    Many studies have explored comprehensive recovery methods to handle disruptions efficiently and effectively.Clarke[6]reviews the practices in airline operations control centers during irregular operations.Filar et al.[7]examine disturbance handling in airports.Kohl et al.[8]provide a detailed overview of many aspects of airline disruption management and report on large-scale research on airline disruption management.Clausen et al.[9]review techniques of aviation recovery and introduce models of schedule planning.In Chapter 9 of their book,Belobaba et al.[10]present schedule recovery and robust planning to mitigate the impact of irregular operation.Barnhart and Smith[11]describe basic models of different resource recoveries and airline disruption management tools(Chapter 6.3).Changes to the mathematical models for aviation recovery problems and other related methods are also mentioned by Floudas and Pardalos[12].

    Below,we briefly explain sources of disruption and disruption propagation to the airline flight network.We then describe common airline recovery operations.

    2.1.Sources of disruption

    We classify the sources of flight disruption into two categories:

    (1)Airline resource disruption(e.g.,aircraft,crew).This type of disruption is caused by factors such as additional maintenance due to aircraft mechanical failures and fuel shortage,or when crew members are absent due to illness or personal emergency.

    (2)External environmental disruption(e.g.,weather,air traffic control).Air travel is weather-sensitive.Even minor weather conditions may reduce the airport departure and arrival rates,causing flight delays.Under severe weather conditions,coercive measures such as airport closures and air traffic control are adopted to ensure the safety of passengers and airline assets.

    Fig.1.An illustration of a down-line impact on a network of four flights scheduled with two aircraft and two crews.Crew 1 is assigned to flights 1 and 4,while crew 2 is assigned to flights 2 and 3.

    As various consecutive flights are arranged for each aircraft and crew,and certain connection rules should be satisfied between flights,disruptions are likely to affect subsequent flights,producing a down-line impact.Fig.1 illustrates an example with two aircraft and four flights.It takes at least 30 min for aircraft to transit,and 45 min for crew members to do so.

    Suppose that airport A is closed from 07:00 to 08:30 due to a thunderstorm.Therefore,flight 1 cannot depart until 8:30,rendering aircraft 2 and crew 1 unable to complete the required transference(Fig.2).Flights 2 and 4 are both affected.

    In addition,the capacity of airports and airspaces should be considered during disruptions,because airports and airspaces accept aircraft of various airlines,but the allocated resources(e.g.,boarding gates,runways)are limited at a given time.Overall,even a very minor disruption can cause significant losses due to disagreements in resource allocation among airlines.Therefore,it is important to adopt recovery methods to prevent or mitigate the down-line impact.

    2.2.Recovery operations

    Flight recovery mainly comprises the following operations,which we explain considering the example in Section 2.1.

    (1)Delaying flights.The departure time of affected flights and related flights may be delayed.In Fig.3,the departure of flights 2 and 4 is delayed by at least 15 min to satisfy the minimum transit times of the aircraft and crew members.

    (2)Cancelling flights.During recovery,if the allocated resources to carry out a flight are not feasible,or if the flight can take place but the delays would exceed a limit,the flight is cancelled.As flight cancellation incurs high costs,this operation is usually the last recovery option for airlines.

    (3)Swapping resources(rerouting).When aircraft or crew members are not prepared for the next flight,other aircraft or crew available in the same airport can substitute for the original ones to carry out the flight.The recovered aircraft or crew is then reallocated to other flights when available.For example,when a disruption occurs,aircraft 1 and 2 can be swapped in airport B for flight 2 to be operated using aircraft 1 without delay,as illustrated in Fig.4.

    (4)Using reserved resources(aircraft and crew).Reserved resources are available in airports and do not perform any flight tasks.

    (5)Deadheading and ferrying.Deadheading means that the crew is transported to another airport as passengers,whereas ferrying means that an aircraft is assigned to an unscheduled flight without passengers.Given the high costs incurred by these operations,they are rarely adopted.

    (6)Speed controlling.Various studies have recently addressed speed controlling as a recovery operation that modifies the flight time to reduce the impact of a disruption and its corresponding delay.

    Fig.2.An illustration of a down-line impact on a network when flight 1 cannot depart on time.The delayed flight 1 is represented by flight 1'.Flights 4 and 2 are affected due to the lateness of crew 1 and aircraft 2 which are assigned to flight 1.

    Fig.3.Example of delaying flight 4 and flight 2 to mitigate the total delay.

    Fig.4.Example of swapping aircrafts 1 and 2 to mitigate the total delay.

    (7)Passenger reallocation.If itineraries are disrupted,passengers can be reallocated to itineraries with the same origin and destination.

    Multiple recovery operations can be adopted simultaneously,according to airline preferences and capabilities.Given the complexity of adopting operations,the recovery problem is commonly separated into a sequence of subproblems that are solved in order.Usually,aircraft recovery is solved first,followed by crew recovery and passenger recovery.In the next three sections,we detail specific models and solutions for aircraft recovery,crew recovery,and integrated recovery considering passengers.

    3.Aircraft recovery

    As both crew rescheduling and passenger resettlement depend on the aircraft arrangement,effective reassignment of aircraft is essential for disruption recovery.Compared with that of planning,the time horizon of recovery is relatively short,varying from hours to days.Aircraft recovery aims to reschedule aircraft routes affected by disruptions at minimum cost while ensuring that flights after the recovery period will not be affected by disruptions.In addition,aircraft should be located at specific stations at the end of the recovery period to carry out the subsequent planned flights.

    Aircraft recovery is typically modeled as a network problem.Like many network routing problems,the adopted models are usually arc-based or path-based,as detailed along with various extensions in the following subsections.

    3.1.Arc-based model for aircraft recovery

    Arc-based models allow handling of disruptions such as airport closure and air traffic control caused by severe weather,security issues,military operations,and other factors.The model is usually built in a time-space network,as proposed by Hane et al.[13],and is widely used for fleet assignment problems.The network has three types of nodes and three types of arcs,as detailed below.An example of the time-space network is shown in Fig.5.

    3.1.1.Nodes

    The supply node is the first node that indicates the beginning of the recovery period at each station when a disruption occurs.The demand node is the last node that indicates the end of the recovery period at each station.It indicates,for example,aircraft positioned in the designated airports to proceed with the flight schedule.An intermediate node is a node with time-station information representing the departure or arrival of a specific flight.

    3.1.2.Arcs

    A flight arc is an arc representing a flight with its scheduled departure/arrival times and stations.A ground arc is an arc representing an aircraft staying on the ground.A copied flight arc is an arc representing the delay options(e.g.,10,30,or 60 min)for original flights.Every delay option refers to a flight arc,in which the departure(arrival)time is the scheduled departure(arrival)time plus the delay for the same departure(arrival)station.

    A general arc-based model can be summarized from Refs.[14-21].The aircraft recovery problem can be formulated on a time-space network,like the one illustrated in Fig.5;this method was first introduced by Thengvall et al.[18].

    Assuming a single-fleet problem,the general mathematical model is described below.

    Model 1:Basic arc-based model for aircraft recovery

    Fig.5.Example of slot capacity for aircraft recovery.

    where F is the set of flight,among which Fn+and Fn-denotes the set of flight inbound to and outbound from node n,respectively;N is the set of intermediate nodes;Fb+is the set of flight inbound to node b;B is the set of supply nodes comprising the first node at each station and representing the starting point of the recovery period;Fe-is the set of flight outbound from node e;E is the set of demand nodes comprising the last nodes at each station and representing the end point of the recovery period;Fsis the set of flights involved in slot s(Fscontains the departure flights if slot s is a departure slot and all the flights if slot s restricts both departure and arrival);S is the set of disrupted slots,which may be caused by weather conditions,security issues,military operations,and other factors;I(f)is the set including the original flight arc f and corresponding copied flight arcs,and I(F)=∪fI(f)denotes the set of all the possible flight arcs;ctfis the assignment cost of the t th copy of flight f(the delay cost is included if the t th copy is not an original flight,i.e.,t>0);cfis the cancellation cost of flight f;Asis the restricted number of aircraft allowed in slot s;Numbis the supplied aircraft at supply node b;Numeis the required aircraft for demand node e;xtfis the decision variable with a value of 1 if the t th copy of flight f is chosen and a value of 0 otherwise;yn+is the number of aircraft on the ground after node n;yn-is the number of aircraft on the ground before node n;yb+is the number of aircraft on the ground after node b;ye-is the number of aircraft on the ground before node e;and zfis the decision variable with a value of 1 if flight f is cancelled and a value of 0 otherwise.

    The objective function in Eq.(1)aims to minimize the total assignment cost,delay cost,and cancellation cost.In addition,aircraft flow balance is maintained by the constraints in Eqs.(2)-(4).Specifically,the constraints in Eqs.(3)and(4)indicate that the aircraft are supplied to operate a sequence of flights and reach the demand nodes at the end of the recovery period.The flight coverage constraints in Eq.(5)ensure that each flight is either cancelled or operated according to its scheduled/delay option.The constraints in Eq.(6)require the number of departure and arrival flights involved in a slot to remain below the slot capacity.

    Fig.5 illustrates the reduction of slot capacity,where air traffic control takes place at airport B from 08:30 to 10:30,in which the departure slot capacity decreases from three to two,and flight arcs 1-4 are involved.Copied flight arcs 5 and 6 are set after the slot to represent the delayed flights due to the disruption.

    Under severe weather conditions,such as hurricanes or typhoons,aircraft should be moved into hangars or be fixed by a ground lock.Otherwise,the aircraft should be relocated to other airports to ensure their safety.Therefore,the number of aircraft on the ground with respect to the hangar capacity or ground locks should be guaranteed by the corresponding constraints.The constraints in Eq.(10)limit the number of ground arcs involved in the corresponding slots.where G is the set of ground arcs along with capacity;Agis the restricted number of ground arcs;Ngdenotes the set of some special intermediate nodes;and the number of aircraft on the ground after these nodes is limited.

    Air traffic control can be used for a specific airport and for any set of flights of interest.Therefore,the set of flights involved in Eq.(10)can be extended to any candidate set designated by a controller.

    Thengvall et al.[19]extend the single-fleet model for handling multiple fleets.The main differences between fleets are the configuration,which is mainly reflected in the number of seats and maximum fly distance,and the crew requirements.Swapping between fleets requires more stringent conditions than swapping within one fleet,as a deviation in the preassigned capacity may lead to substantial profit loss.

    The multi-fleet model adopts a set of time-space networks,one for per sub-fleet[19].As in the single-fleet model,a flight can be operated by all aircraft in the scheduled sub-fleet.Furthermore,some flights can be operated by multiple sub-fleets during recovery.As shown in Fig.6,sub-fleet 2 can operate the flights originally scheduled to sub-fleet 1;hence,flights 2 and 3 belong to the network of sub-fleet 2.This model is solved using the IBMCPLEX Optimizer.For the instance with 12 fleets containing 1-6 sub-fleets,1434 flights during the 24 h recovery period,and 10 h of airport closure,the model provides a near-optimal solution after 1838.4 s of computation.In this near-optimal solution,81 flights are operated by sub-fleets different from the original assignments.This indicates that swapping between fleets provides more flexibility and improves the solution quality.

    Aircraft routes cannot be directly obtained from the basic arcbased model,because the flow in the time-space network does not distinguish a specific aircraft.Therefore,the arc-based model cannot handle disruptions such as the unexpected maintenance of a specific aircraft.To solve this problem,Vink et al.[21]extend the arc-based model to incorporate specific aircraft.A set of parallel networks is established,where each network represents a specific aircraft.This set of networks is bundled by the flight coverage constraints.The model is extended by including extra aircraft index for each variable in Model 1.

    Constructing the network for each aircraft and solving the corresponding integer linear program can be very time-consuming.To achieve real-time performance,Vink et al.[21]develop a selection algorithm from the proposal by Vos et al.[20]that comprises three stages.The number of selected aircraft involved in each stage is limited in order to speed up the solving process.If no solution is found by the selected aircraft,the set of candidate aircraft is expanded and moved to the next stage.It takes 22 s on average to solve an instance with 100 aircraft,600 daily flights,and a 16 h recovery period by the selection algorithm,while it takes 10 min for the integer linear program considering the entire set.The selection algorithm finds the global optimal solution in seven of ten scenarios.On average,the gap between the selection algorithm solution and the global optimal solution is 6%.

    Fig.6.Example of a time-space network for a multi-fleet model.

    3.2.General path-based model for aircraft recovery

    The disruption concerning individual aircraft can also be handled using a path-based model.Compared with the arc-based model,the path-based model assigns aircraft to a route that includes detailed information such as flight delay and aircraft swap.Argüello et al.[22],Rosenberger et al.[23],Eggenberg et al.[24],Wu et al.[25],and Liang et al.[26]have developed models based on paths(routes).The path-based model for aircraft recovery can also handle disruptions such as airport closure and air traffic control by adding side constraints to the model.In addition,situations such as unplanned maintenance of an aircraft can be handled by this model.

    The path-based model enables a feasibility check(i.e.,guaranteed minimum turnaround times)when generating routes.Constraints and regulations can be implicitly included during route generation and eliminated from the model.The path-based model is described below.

    Model 2:Basic path-based model for aircraft recovery

    The objective function in Eq.(11)aims to minimize the assignment cost and cancellation cost.The resource utilization constraints in Eq.(12)imply that each aircraft can operate one route.The origin and destination of routes r correspond to the stations where recovery starts and ends.The flight coverage constraints in Eq.(13)indicate that each flight is covered by only one aircraft,and the flight should be cancelled if no aircraft is available.

    Table 1 summarizes the characteristics of the arc-and pathbased models.Compared with the arc-based model,the pathbased model has fewer constraints but a greater number of possible routes.As cancellation and delay are the most common recovery operations,we should analyze their formulations.Cancellation is represented by cancellation variables in both the arc-and pathbased models.However,delay has a variety of representations,as detailed in Section 3.3.

    Table 1 Characteristics of the arc-and path-based models.

    3.3.Recovery with delay

    Flight delay is an important recovery operation.A short delay for recovery can sometimes be absorbed by the buffer time between flights,allowing flight connections to be maintained.Thus,selecting a delay for a minor disruption is effective.

    Using copied flights to represent delay,as described by Model 1,is widely adopted in aircraft recovery.The delay options are given by a discrete set of predetermined delay times.Thengvall et al.[18,19],Eggenberg et al.[24],Vos et al.[20],and Vink et al.[21]adopt copied flights as delay options in their formulations.However,discrete delay presents some disadvantages.①Delay may be overestimated.For example,the actual delay of a flight may be 12 min,but the only available option is 20 min,so an unnecessary delay cost is incurred.②To prevent overestimation,a smaller time interval can be set,but doing so increases the problem size.

    The departure and arrival times of flights can also be decision variables in models.Aktürk et al.[27]propose aircraft rescheduling considering swapping,delaying,and adjusting cruise speed,in which the departure delay and cruise time are decision variables,as detailed below.

    Fig.7 illustrates this model.Case I shows the original schedule of flight f and subsequent flight N(f).Then,a delay combining cruise time decision is added to the original schedule in case II.The relationship of flight f and flight N(f)is easily obtained,as follows:

    Fig.7.Illustration of a relationship of swapping flights.I.Original schedule,in which flight N(f)follows flight f.II.Schedule considering delay and cruise time reduction of flight f.III.Schedule considering delay,cruise time reduction,and swapping,in which flight f is swapped with flight j.df:scheduled departure time of flight f;rf:scheduled arrival time of flight f;TA f:turnaround time of flight f;lf:departure delay time for flight f;lj:departure delay time for flight j;t0 f:original cruise time of aircraft operating flight f;tf:cruise time of aircraft operating flight f;tj:cruise time of aircraft operating flight j;N(f):next flight after operating flight f in the original schedule;dN(f):the scheduled departure time of flight N(f);lN(f):the departure delay time for flight N(f).

    where S(f)is the set of flights that can be swapped with flight f and xf,jis the decision variable with a value of 1 if flights f and j are swapped and a value of 0 otherwise.

    Liang et al.[26]propose an adaptive delay method for airport traffic control,in which the number of departures or arrivals of specific airports are limited in some slots.This can be achieved by introducing the constraints in Eq.(18)to the path-based model.

    This problem is solved using column generation with a multilabel shortest-path algorithm for the corresponding pricing subproblem.Two labels,referring to cost and delay,are considered in the shortest-path problem,and the dominated paths are removed during searching.To ensure subproblem optimality,the affected flights are duplicated at the beginning and end of each slot.The largest evaluated instance under disruptions of air traffic control over 2 h and unexpected maintenance over 8 h with 638 flights and 44 aircraft is optimized within 356 s.

    3.4.Recovery with maintenance

    Maintenance is commonly treated as a fixed activity during aircraft recovery.Eggenberg et al.[24]treat maintenance as a resource-consuming and renewal process.Column generation with the resource-constrained elementary shortest path is used to find the routes with maintenance feasibility in each subproblem.A real-life instance of 16 aircraft and 242 flights over 7 d with a flight-to-aircraft ratio of 18.4 is solved within 3603 s.

    Liang et al.[26]determine maintenance requirements using a path-based model on a connection network.Besides the consuming resource concept used by Eggenberg et al.[24],scheduled maintenance operations can be swapped.Swapping can only take place if aircraft are of the same fleet type and the resource limitation is satisfied.Maintenance swapping improves the flexibility of recovery,increasing the probability of finding better recovery solutions.Again,column generation with the resource-constrained shortest path is used.Three additional labels—namely,total fly time,number of takeoffs and landings,and elapsed time—are considered in the multilabel shortest-path algorithm.Compared with fixed maintenance,the cost of recovery with flexible maintenance decreases by 61.47%.

    An arc-based model considering maintenance is proposed by Vink et al.[21].Two types of maintenance operations are considered.One is planned maintenance,which is difficult to change to another time or station(e.g.,when the maximum elapsed time between two maintenance operations is nearly reached).This maintenance must be performed as scheduled on aircraft.The other type of maintenance provides various time and space options,and at least one option is selected to meet the maintenance requirements in a route.

    Table 2[14-31]summarizes the features and solving methods of different studies on aircraft recovery approaches.

    4.Crew recovery

    Before presenting crew recovery,we briefly introduce crew scheduling to provide a background on some necessary concepts.Crew scheduling consists of arranging a set of flight tasks within a given period(e.g.,15 or 30 d)for the pilot,so that each flight in the schedule can be operated by one or more crews.A crew can only operate the fleet if properly qualified.Therefore,flights are grouped by fleet,and crew scheduling is solved separately for each fleet.In addition,the crew schedule should obey regulations from governing agencies and labor unions to ensure security and operability.

    Compared with crew scheduling,crew recovery involves fewer airports and shorter schedule times,notably reducing the scale of the problem.To fully leverage the small scale of the problem,the flights not involved in the recovery period are usually treated as fixed activities.Hence,crew recovery aims to find a solution witha minimum cost for reassigning the available crews to affected flights,while crews start and end at given fixed activities[32].

    Table 2 Features and solving methods of aircraft recovery approaches.

    Similar to crew scheduling,a crew recovery problem is usually constructed on a single-fleet network.Labor regulations applicable to planning should be complied during recovery.Moreover,preassigned activities,such as vacation and training,should be reflected in the roster of the crew after recovery.

    Teodorovic′and Stojkovic′[33]are the first to develop a mechanism to solve the crew recovery problem during daily operation.They adopt two methods to construct a crew duty—namely,a first-in-first-out scheme and dynamic programming—to minimize the ground time of a crew.However,the optimality of the problem is neglected.Subsequent studies have applied optimization to solve the crew recovery problem.Below,we present a general model for crew recovery as well as an extension and the corresponding method.

    4.1.General model for crew recovery

    Wei et al.[32],Stojkovic′et al.[34],Lettovsky′et al.[35],Yu et al.[36],and Guo et al.[37]use a set partitioning model to solve the crew recovery problem.Below,we detail the formulation proposed by Lettovsky′et al.[35].

    Model 3:Basic model for crew recovery

    This model considers coverage constraints and assignment constraints.The objective function in Eq.(19)aims to minimize the pairing cost,cancellation cost,deadhead cost,and cost of idle crews.Eq.(20)indicates that a flight can be covered more than once,and the surplus coverage identifies the number of flight f used for deadheading.If flight f cannot be covered,the cancellation variable zfhas a value of 1 and incurs the corresponding cancellation cost.Eq.(21)indicates that one crew can operate at most one pairing.

    Wei et al.[32]solve the above formulation using a heuristic algorithm,in which a set of pairings is generated using a shortest-path algorithm on a time-space network.This model aims to find a feasible schedule with minimum deviation from the original plan by setting different costs for flight arcs.An instance of 6 airports,51 flights,and 18 pairings over two days is evaluated.The computation time varies from 1 to 6 s.

    In addition,the structure of crews should be considered for operational solutions.Crew members are usually divided into three ranks:captain,first officer,and second officer.To obtain a realistic solution,Medard and Sawhney[38]consider a multirank extension by modifying Eq.(20)into Eq.(25).

    4.2.Extensions of crew recovery

    Crew recovery can also affect flight delay.Stojkovic′and Soumis[39]consider flight delays when constructing feasible pairings.Moreover,a set of constraints to protect the passenger connection is added to the master problem.

    If arc(w,j)is in the aircraft route,then dw,j=blw+gw,indicating the block time(blw)plus the turnaround time(gw).If there is an important passenger connection on arc(w,j),then dw,j=blw+cw,j,indicating the block time plus the passenger connection time(cw,j).The flight precedence constraints in Eq.(26)protect the aircraft connections and passenger connections by limiting the period between consecutive flights in a pairing.

    Stojkovic′and Soumis[40]develop a multi-rank model by adding synchronization constraints.A set of copies from original flights,defined as tasks,represent the different rank requirements.As these tasks belong to the same flight,Eq.(27)is added to the model of Stojkovic′and Soumis[39].

    Eq.(27)ensures that every task associated with flight f has the same departure time.Column generation is used to solve this multi-rank crew pairing problem and add the synchronization constraints to the master problem.Eq.(27)is rewritten as follows:

    The model is applied to data from domestic flights in the United States.An instance of 190 flights,consisting of 46 and 20 flights with flexible and fixed departure times,respectively,is evaluated.Each flight requires five crew members,and the number of involved crews is 97.This instance is solved in 1237 s.

    Abdelghany et al.[41]use a greedy heuristic to iteratively solve the crew recovery problem.The disrupted flights are grouped into different stages,in which the flights are resource independent,as shown in Fig.8.Resource independence is defined as those flights not sharing the same crew resources in a stage.Hence,in each iteration,the resource-independent flights are assigned to the available crews using the following formulation.

    Fig.8.A rolling framework for crew recovery.

    where dfis the scheduled departure time of flight f;blfis the scheduled block time of flight f;akis the ready time for crew k;vkis the specific time in which maximum on-duty time is met for crew k;xk,f,qis the decision variable with a value of 1 if crew k assigned to flight f for position q;dtfis the actual departure time of flight f;and atfis the actual arrival time of flight f.

    Eqs.(29)and(30)ensure the coverage of flights and that each crew has at most one rank.The constraints in Eqs.(31)and(32)indicate that the actual departure time should not be earlier than the scheduled departure time and the crew-ready time.The constraints in Eq.(33)ensure that the arrival time does not exceed the duty limit.The constraints in Eq.(34)indicate the relationship between departure and arrival times.Data from major airlines in the United States with an 8 h recovery period and 121 crews is solved in 1 min and 51 s.Although this rolling framework cannot retrieve a globally optimal solution,the solution is practical and can be obtained in real time.

    To limit the model size and speed up recovery,Lettovsky′et al.[35]divide the partial pairings in the recovery period into segments.A partial pairing is illustrated in Fig.9.The flights in duty 2 are split into segments that consist of one or more flight legs.Segments are covered instead of flights,thus reducing the number of rows in the model.Using segments may reduce both the number of generated pairings and the computation time.

    4.3.Solving methods for the crew recovery problem

    The general model for crew recovery includes few constraints but a huge number of variables;hence,most recovery approaches are solved using column generation.Column generation is an efficient method to avoid explicitly enumerating all the variables while maintaining the optimality of linear programming relaxation.Under column generation,a path for crew k is generated in the subproblem,and the optimal set of paths for recovery is obtained in the master problem.Although the problems in Refs.[34,35,38,40-42]are solved using column generation,the construction of the networks differs by adopting methods based on elements such as flights,segments,and duties.Table 3[32-44]summarizes the features and solving methods of different studies on crew recovery.

    5.Integrated recovery

    In this section,we present characteristics and considerations for the integration of different resources to perform recovery.

    5.1.Integrated aircraft and crew recovery

    Both aircraft and crews play important roles during recovery,with the former being the scarcest resource for airlines.Complex restrictions and the large scale of the problem hinder the integrated modeling of these two resources.

    5.1.1.Basic link

    The basic link between aircraft and crew recovery problems is provided by flight cancellation and delay decisions[45-51].As detailed in Sections 3 and 4,the basic models of aircraft and crew recovery have the common variable zf,which represents the cancellation of a flight.

    Fig.9.Illustration of segments.Segment 1 contains one flight(i.e.,flight 2),while segment 2 consists of two flights(i.e.,flight 3 and flight 4).

    Table 3 Features and solving methods of crew recovery approaches.

    Maher[50]uses the connection network for recovery by focusing on aircraft and crews.The network uses flight copies to represent different delays.He integrates a path-based model for aircraft recovery and a general model for crew recovery by adding delay consistency constraints.

    The constraints in Eq.(35)ensure that the delay on each flight is consistent for the crews and aircraft.Moreover,the number of delay consistency constraints changes depending on the input paths and duties.

    5.1.2.Aircraft and crew compatibility

    A crew member cannot operate all types of fleets.As a pilot requires the ability to fly a specific aircraft when carrying out a flight,it is important to ensure resource compatibility.Abdelghany et al.[46]and Ar?kan et al.[51]build arc-based models considering these compatibility constraints,which are described as follows in Ref[46]:

    The constraints in Eq.(36)ensure that only compatible aircraft and crew can be assigned to a flight.In addition to these constraints,the model by Abdelghany et al.[46]includes coverage constraints for aircraft and crews,and other constraints to ensure delay time feasibility.A greedy iterative heuristic is used to solve the integrated recovery problem.It takes about 36 s to obtain a recovery plan for an instance with 522 aircraft,1360 pilots,2040 flight attendants,and 1100 daily flights after ten disrupted flights.

    Ar?kan et al.[51]adopt a special flight network that contains four types of nodes—namely,scheduled flight nodes,source nodes,sink nodes,and must-visit nodes—which represent maintenance requirements or scheduled crew rest periods.Then,an arc-based model is established on this network and the resource compatibility constraints are formulated as Eq.(37).

    If a crew is assigned to a flight,the crew should be qualified to operate the corresponding aircraft.By using the problem size control algorithm and passenger aggregation,the total running time for a network containing 1254 flights with 402 aircraft is reduced to less than 12 min when the cruise speed-controlling method is considered.

    It is notable that,as the arc-based model is adopted,it is very difficult—or even impossible—to formulate the complex legality constraints for pairings or duties of crews.

    5.1.3.Connection feasibility

    Different resources,such as aircraft,crews,and passengers,have different transit time requirements.Therefore,the control of the flight interval by delay decisions partly determines whether reasonable and high-quality routings and duties can be generated.Various integrated recovery studies[47,48,50]have incorporated connection considerations into the construction of routings and duties to omit additional constraints that guarantee the feasibility in their models.

    The connection feasibility constraints are included in the model proposed by Ar?kan et al.[51].These constraints can be formulated as Eq.(38).If a resource is assigned to an arc,the period between the arrival of the previous flight and the departure of the next flight should comply with the minimum transit time required for the resources.

    As solving an integrated model is computationally intractable,Zhang et al.[49]use a two-stage heuristic algorithm for the integrated recovery problem.The model mainly decomposes the integrated recovery problem into two models to solve the aircraft and crew recovery problems separately.The two stages are linked by the connection feasibility between consecutive flights for aircraft and crews.

    When crew members fly the same aircraft consecutively,there is no need to obey the crew transit time requirements.This can be considered as a short connection,and considering this situation increases the solution flexibility.Maher[48,50]takes this into account.The corresponding variables are described in Section 5.1.1,and the situation can be formulated as Eq.(39).When an aircraft is assigned to a route containing a short connection,the pairing that contains the same short connection is covered by a crew.

    Nevertheless,it remains difficult to obtain a complete recovery schedule for aircraft and crews.

    5.2.Integrated aircraft recovery considering passengers

    Passenger recovery is important for airlines,as irregular schedules can adversely affect the itineraries of passengers.In serious cases,passengers may fail to reach their destination,and airlines may have to refund the tickets or arrange flights with other airlines for them.Hidden costs include the airline’s loss of credibility,which cannot be easily estimated and which invisibly affects passengers’future travel choices.Therefore,recovering the itineraries of affected passengers quickly and reasonably improves the market competitiveness of airlines.

    5.2.1.Integration of itinerary-based passenger recovery

    Most studies on integrated recovery considering aircraft and passengers are based on itineraries[45,47,48,52-55].In these studies,a binary variable is usually adopted to determine whether an itinerary is affected.An itinerary may be disrupted by either flight cancellation or violation of passenger transit times.This variable and some related constraints in Eqs.(40)and(41)are detailed below.

    The constraints in Eq.(40)ensure that if a flight is cancelled,its itinerary is disrupted.The constraints in Eq.(41)ensure that if the passenger transit time is not satisfied in an itinerary,the itinerary is disrupted.In the objective function of the model proposed by Ar?kan et al.[52],the items related to passengers include delay cost and spill cost.In addition to the disruption of itineraries,some passengers may‘‘spill”(i.e.,be left over)due to capacity shortage when aircraft is swapped.The fuel cost is expressed as a nonlinear function related to cruise time and is included in the objective function.The problem is reformulated into a model with a linear objective function and conic quadratic constraints.Limited by model complexity,passenger reallocation is not part of the formulation in Ref.[52].

    Regarding integrated recovery considering passengers,many studies[47,48,51,54-56]have addressed the rearrangement of passengers.The construction of the corresponding constraints uses a judgment variable in Eqs.(42)-(45).

    The constraints in Eq.(42)guarantee that all passengers on any itinerary are either successfully transported to their destination or refunded.The constraints in Eq.(43)ensure that the number of passengers cannot exceed the capacity of the aircraft assigned to the flight.The constraints in Eq.(44)guarantee that no passenger is allocated to a disrupted itinerary.The constraints in Eq.(45)ensure that passengers follow the original itinerary without any change under no disruption.

    Limited by the scale of the problem and the allowable computation time,although various models have considered passenger reallocation,approximations are usually adopted without considering reallocation constraints but by penalizing the number of disrupted itineraries in the objective function.Marla et al.[53]use a mechanism called flight planning that enables flight speed changes in long-haul flights.Unlike the method in Ref.[52],Marla et al.[53]discretize the dynamic selection of flying time instead of considering the cruise time as a variable.An accurate model and an approximate model have been proposed.As obtaining the solution to the accurate model is very time-consuming,the approximate model is used for case analysis.Similarly,Bratu and Barnhart[45]combine aircraft,crew,and passenger recovery by proposing two models called the disrupted passenger metric model(DPM)and the passenger delay metric model(PDM).The latter can more accurately calculate the cost of passenger delay by allowing passenger rescheduling.In an instance involving the reallocation of more than 80 000 passengers,the average calculation time of the PDM is approximately 25 times that of the DPM.

    5.2.2.Integration of flight-based passenger recovery

    Besides itinerary-based passenger recovery,other methods are available for passenger recovery.The method proposed by Ar?kan et al.[51]treats all resources,including aircraft,crews,and passengers,as the same entities and maintains the connection feasibility and flow balance of each entity in a flight network.However,this method increases the network scale,and entity aggregation should be used to reduce the number of entities without compromising optimality.Hu et al.[56]construct a model based on a time-band network and propose a flight-based passengertransiting mechanism,in which flights with the same departure and arrival airports are arranged for passengers whose itinerary is disrupted.Maher[48]builds a flight-based passenger recovery model based on the point-to-point network,which is considered to be non-multi-itinerary.The formulation of the model extends the airline recovery formulation in Ref.[50]to include passenger recovery.

    5.2.3.Heuristic-based methods

    Various heuristic algorithms exhibit high performance for simultaneously reassigning aircraft and passengers in case of disruptions.

    Bisaillon et al.[57]propose a large neighborhood search heuristic algorithm that is divided into construction,repair,and improvement phases.The first two phases generate a feasible solution.The third phase implements large changes to improve the solution based on the solution from the previous phases,while preserving feasibility.The construction phases are then executed again iteratively.This algorithm provides high-quality solutions and can handle large-scale problems in real time.Sinclair et al.[58]include additional steps for each phase of the algorithm in Ref.[13]to improve performance.Experimental results show that the destroy-and-create step,which is added to the third phase,has the greatest improvement on the quality of the solution.

    Sinclair et al.[59]also propose a column generation postoptimization heuristic algorithm.An integrated linear programming model based on a time-space network is established.The algorithm in Ref.[58]is used to obtain the initial solution and to construct constraints of the restricted master problem(RMP).Then,in order to obtain suitable passenger rearrangement options,the corresponding column generation subproblem is solved.

    Jozefowiez et al.[60]propose the new connection flights heuristic method,which is a three-stage heuristic algorithm based on shortest-path calculation.The first stage integrates various types of disruptions into existing schedules.In the second stage,passengers with a disrupted itinerary are assigned to candidate feasible itineraries by solving the corresponding shortest-path problem.In the third stage,a new sub-rotation is inserted into the existing aircraft rotation to allocate additional passengers.

    Hu et al.[61]use a greedy randomized adaptive search procedure(GRASP)algorithm for the integrated recovery of aircraft and passengers after airline operation disruption.Suitable passenger reassignment can be obtained based on new aircraft routing in each iteration by applying the GRASP algorithm.The iterations proceed until one of the stopping criteria is met.

    5.3.Integrated airline recovery problem

    Thus far,few studies have fully addressed integrated recovery[45,47,48,51],as detailed throughout this section.The model proposed by Bratu and Barnhart[45]is one of the first attempts to solve fully integrated problems,but only considers reserved crews.To obtain real-time solutions for integrated recovery,Ar?kan et al.[51]use two preprocessing methods to considerably reduce the number of constraints and variables,and an algorithm to limit the scope of recovery.

    Petersen et al.[47]integrate five subproblems for recovery of the schedule,aircraft,crew,itinerary,and passengers.They propose a Benders decomposition framework to solve the integration problem.An example with 800 daily flights and two fleets is solved by means of the integrated method and sequential solution,respectively.Although the integrated method shows higher performance,it requires further improvements to solve larger-scale problems within 30 min.

    Maher[48]uses the column-and-row generation method.For integrated recovery,the column generation subproblems are solved for the duty and aircraft routing variables using a shortest-path algorithm.A passenger rearrangement scheme is formulated as a knapsack problem in column generation.The rows,which are added to the model iteratively,represent passenger reallocation constraints and delay consistency on each flight between aircraft and crew.A point-to-point schedule with 262 flights and a hub-and-spoke schedule comprising 441 flights are solved in 427 and 400 s,respectively.

    Some characteristics of existing integrated recovery approaches are summarized in Table 4[45-61].

    6.Conclusions and future work

    In this review,we presented basic models and extensions to disruption management for aircraft and crew recovery,as well as integrated recovery considering passengers.For aircraft recovery,we reviewed the basic arc-based and path-based models and discussed the situations that these models adapt to.We also summarized the formulations related to delay and maintenance.In terms of crew recovery,we reviewed the basic model and the extensions related to crew ranks and delay.For integrated recovery,the key constraints linking different resources and related algorithms wereintroduced.After reviewing the development of disruption management,we provided the following suggestions for future research.In recovery using one stage,more realistic factors should be considered in order for the solution to be more practical and useful.For example,aircraft swapping between different fleets,as a flexible yet more complicated recovery method,can be modeled.In addition,solving the fully integrated problem still remains challenging due to its complexity.Future research should focus on integrating several stages and devising optimization methods to solve the integrated models.In integrated recovery,preprocessing to reduce the size of input data can be adopted in order to reduce the complexity.Decomposition may be used as a preprocessing approach.The loss due to disruption mainly comes from passengers.Therefore,passenger-centered recovery is a promising research topic.For example,passenger preferences should be included for recovery in order to improve service.Likewise,evaluating passenger information through data mining may lead to more accurate optimization.

    Table 4 Characteristics of integrated recovery approaches.

    Acknowledgement

    This study is supported by the National Natural Science Foundation of China(71825001 and 71890973).

    Compliance with ethics guidelines

    Yi Su,Kexin Xie,Hongjian Wang,Zhe Liang,Wanpracha Art Chaovalitwongse,and Panos M.Pardalos declare that they have no conflict of interest or financial conflicts to disclose.

    亚洲精品乱久久久久久| av欧美777| 美女国产高潮福利片在线看| 成人影院久久| 国产成人啪精品午夜网站| 亚洲avbb在线观看| 日韩制服骚丝袜av| 精品久久蜜臀av无| 日韩一卡2卡3卡4卡2021年| 国产欧美日韩综合在线一区二区| 老司机影院毛片| 日韩熟女老妇一区二区性免费视频| 在线天堂中文资源库| 亚洲av美国av| 高清视频免费观看一区二区| 99精国产麻豆久久婷婷| 亚洲精华国产精华精| 视频区图区小说| 欧美国产精品一级二级三级| 亚洲国产欧美在线一区| 亚洲精品久久久久久婷婷小说| 麻豆av在线久日| 男女免费视频国产| 色精品久久人妻99蜜桃| 久久久国产成人免费| 国产精品久久久久久人妻精品电影 | 免费高清在线观看日韩| 日韩视频一区二区在线观看| a级毛片在线看网站| 女警被强在线播放| 男女国产视频网站| bbb黄色大片| 最黄视频免费看| 在线观看免费日韩欧美大片| 久久久久视频综合| 成人手机av| 亚洲黑人精品在线| 久久99一区二区三区| 日韩视频一区二区在线观看| 国产伦人伦偷精品视频| 日韩大码丰满熟妇| 国产91精品成人一区二区三区 | 精品人妻熟女毛片av久久网站| 天天添夜夜摸| cao死你这个sao货| 久久久久久久久久久久大奶| 亚洲av日韩在线播放| 三级毛片av免费| 青春草亚洲视频在线观看| 免费一级毛片在线播放高清视频 | 男人操女人黄网站| 女警被强在线播放| 日韩大码丰满熟妇| 精品国产一区二区久久| 一本—道久久a久久精品蜜桃钙片| 欧美另类亚洲清纯唯美| 亚洲欧美日韩高清在线视频 | av不卡在线播放| 亚洲精品中文字幕在线视频| 无遮挡黄片免费观看| 久久狼人影院| av在线老鸭窝| 国产野战对白在线观看| 亚洲,欧美精品.| www.999成人在线观看| 国产精品一二三区在线看| 精品一区二区三区av网在线观看 | 亚洲黑人精品在线| 国产欧美日韩综合在线一区二区| 人人妻,人人澡人人爽秒播| 国产深夜福利视频在线观看| 人妻一区二区av| 天天添夜夜摸| cao死你这个sao货| av福利片在线| 欧美另类亚洲清纯唯美| 国产精品免费大片| 性少妇av在线| 国产亚洲一区二区精品| 12—13女人毛片做爰片一| 国产高清视频在线播放一区 | 久久国产精品男人的天堂亚洲| 久久精品aⅴ一区二区三区四区| 人妻一区二区av| 国产一区二区激情短视频 | 国产精品一区二区免费欧美 | 欧美黄色片欧美黄色片| 日本黄色日本黄色录像| 亚洲精品美女久久av网站| 久久女婷五月综合色啪小说| 久久天躁狠狠躁夜夜2o2o| 午夜老司机福利片| 自拍欧美九色日韩亚洲蝌蚪91| 各种免费的搞黄视频| 黄色毛片三级朝国网站| 丰满饥渴人妻一区二区三| 天天躁狠狠躁夜夜躁狠狠躁| 他把我摸到了高潮在线观看 | 免费观看人在逋| tocl精华| 亚洲欧洲精品一区二区精品久久久| 人人澡人人妻人| 日韩三级视频一区二区三区| 男女边摸边吃奶| a级片在线免费高清观看视频| 美女脱内裤让男人舔精品视频| 亚洲久久久国产精品| 国产亚洲精品久久久久5区| 国产一区二区激情短视频 | 91精品伊人久久大香线蕉| 99热全是精品| 国产在线一区二区三区精| 99国产极品粉嫩在线观看| a级片在线免费高清观看视频| 大陆偷拍与自拍| 欧美日韩视频精品一区| 夜夜夜夜夜久久久久| 人人澡人人妻人| 亚洲欧美日韩另类电影网站| 午夜免费成人在线视频| 亚洲国产精品999| 69av精品久久久久久 | 可以免费在线观看a视频的电影网站| 窝窝影院91人妻| 婷婷成人精品国产| 宅男免费午夜| 天天操日日干夜夜撸| 欧美亚洲 丝袜 人妻 在线| 国产成人欧美在线观看 | 亚洲国产精品一区二区三区在线| 免费看十八禁软件| 久久久久国产一级毛片高清牌| 久久 成人 亚洲| 日韩中文字幕欧美一区二区| videos熟女内射| 亚洲精品av麻豆狂野| 王馨瑶露胸无遮挡在线观看| 国产视频一区二区在线看| 国产xxxxx性猛交| 国产一区二区激情短视频 | 亚洲精品国产区一区二| 午夜福利乱码中文字幕| 成在线人永久免费视频| 久久香蕉激情| 国产精品 欧美亚洲| 精品国内亚洲2022精品成人 | 国产男女超爽视频在线观看| 久久精品成人免费网站| 日韩大片免费观看网站| 久久国产精品人妻蜜桃| 国产精品久久久久成人av| 婷婷色av中文字幕| 午夜精品国产一区二区电影| 国产xxxxx性猛交| 性高湖久久久久久久久免费观看| 亚洲精品日韩在线中文字幕| 国产精品亚洲av一区麻豆| 久久青草综合色| 亚洲欧美激情在线| 日韩欧美国产一区二区入口| 男女床上黄色一级片免费看| 欧美在线一区亚洲| 首页视频小说图片口味搜索| 久久热在线av| 国产1区2区3区精品| 高清欧美精品videossex| 18禁观看日本| 午夜福利,免费看| 亚洲色图 男人天堂 中文字幕| 久久国产精品人妻蜜桃| a级毛片黄视频| 国产在线一区二区三区精| 日本精品一区二区三区蜜桃| 丝袜美腿诱惑在线| 1024视频免费在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 99国产精品一区二区蜜桃av | 成年人午夜在线观看视频| 飞空精品影院首页| 一个人免费看片子| 日本黄色日本黄色录像| 亚洲五月婷婷丁香| 亚洲avbb在线观看| 黑人操中国人逼视频| 少妇粗大呻吟视频| 国产精品久久久人人做人人爽| 成人免费观看视频高清| 如日韩欧美国产精品一区二区三区| 777米奇影视久久| 91九色精品人成在线观看| 飞空精品影院首页| 高清视频免费观看一区二区| 色精品久久人妻99蜜桃| 久久久久久亚洲精品国产蜜桃av| 精品一区二区三区av网在线观看 | 搡老乐熟女国产| 岛国在线观看网站| 日本欧美视频一区| 大型av网站在线播放| 国产激情久久老熟女| 精品卡一卡二卡四卡免费| 老汉色∧v一级毛片| 久久亚洲精品不卡| 叶爱在线成人免费视频播放| 日本av免费视频播放| 制服人妻中文乱码| 亚洲欧美激情在线| 国产免费av片在线观看野外av| 最近最新中文字幕大全免费视频| 两人在一起打扑克的视频| 国产成人a∨麻豆精品| 美女视频免费永久观看网站| 久久中文看片网| 汤姆久久久久久久影院中文字幕| 欧美久久黑人一区二区| 一区二区三区激情视频| av视频免费观看在线观看| 亚洲av欧美aⅴ国产| 91麻豆av在线| av电影中文网址| 国产在线观看jvid| 成在线人永久免费视频| 日本wwww免费看| 国产免费一区二区三区四区乱码| 亚洲 国产 在线| 国产亚洲欧美在线一区二区| 午夜老司机福利片| 国产成人a∨麻豆精品| 男女午夜视频在线观看| 人成视频在线观看免费观看| 久久久精品94久久精品| 天堂俺去俺来也www色官网| kizo精华| 超色免费av| 午夜免费鲁丝| 每晚都被弄得嗷嗷叫到高潮| 性高湖久久久久久久久免费观看| 免费在线观看黄色视频的| 久久 成人 亚洲| 欧美+亚洲+日韩+国产| www.av在线官网国产| 69精品国产乱码久久久| 91精品伊人久久大香线蕉| 日本91视频免费播放| 久久精品国产亚洲av高清一级| 男女高潮啪啪啪动态图| 久久久久久免费高清国产稀缺| 欧美激情极品国产一区二区三区| 99久久国产精品久久久| 国产成人a∨麻豆精品| 欧美+亚洲+日韩+国产| 男女免费视频国产| 国产成人啪精品午夜网站| 亚洲精品第二区| www.av在线官网国产| 亚洲国产毛片av蜜桃av| 午夜福利一区二区在线看| 国产成人精品久久二区二区91| 欧美国产精品一级二级三级| 国产精品免费视频内射| 蜜桃在线观看..| 亚洲国产欧美一区二区综合| 男女国产视频网站| 午夜两性在线视频| 欧美精品高潮呻吟av久久| 国内毛片毛片毛片毛片毛片| 超色免费av| 国产一区二区 视频在线| 婷婷色av中文字幕| 日韩三级视频一区二区三区| 国产一区二区在线观看av| 高清欧美精品videossex| 国产精品香港三级国产av潘金莲| 1024视频免费在线观看| 午夜免费鲁丝| 日韩制服丝袜自拍偷拍| 香蕉丝袜av| 91麻豆av在线| 91九色精品人成在线观看| 国产主播在线观看一区二区| 欧美在线一区亚洲| 国产日韩一区二区三区精品不卡| 日韩中文字幕欧美一区二区| 亚洲国产日韩一区二区| 久久精品亚洲av国产电影网| 看免费av毛片| 国产精品影院久久| 丝袜美腿诱惑在线| 少妇的丰满在线观看| 无限看片的www在线观看| 午夜两性在线视频| 老汉色av国产亚洲站长工具| 亚洲人成电影免费在线| 啦啦啦在线免费观看视频4| 999精品在线视频| 色94色欧美一区二区| av又黄又爽大尺度在线免费看| 欧美日韩av久久| 久久亚洲精品不卡| 黑丝袜美女国产一区| 欧美精品一区二区免费开放| 国产成人av激情在线播放| 久久久久久久国产电影| 国产av精品麻豆| 美女扒开内裤让男人捅视频| 欧美激情久久久久久爽电影 | 脱女人内裤的视频| 人人妻,人人澡人人爽秒播| 亚洲国产精品一区二区三区在线| 国产精品一区二区在线不卡| 久久久久久人人人人人| 国产精品av久久久久免费| 久久久精品94久久精品| 亚洲精品国产区一区二| 成人免费观看视频高清| 狂野欧美激情性xxxx| 亚洲va日本ⅴa欧美va伊人久久 | 国产福利在线免费观看视频| 老司机午夜十八禁免费视频| 美女扒开内裤让男人捅视频| av电影中文网址| 精品人妻一区二区三区麻豆| 国产成人精品在线电影| 18禁黄网站禁片午夜丰满| 淫妇啪啪啪对白视频 | 国产精品香港三级国产av潘金莲| 国产有黄有色有爽视频| 日本vs欧美在线观看视频| av在线播放精品| 欧美亚洲日本最大视频资源| 久久av网站| 麻豆av在线久日| 午夜久久久在线观看| 欧美人与性动交α欧美软件| 自线自在国产av| 纵有疾风起免费观看全集完整版| av国产精品久久久久影院| 国产日韩欧美亚洲二区| 99国产综合亚洲精品| 国产一区二区三区av在线| netflix在线观看网站| 飞空精品影院首页| 亚洲精品美女久久av网站| 日韩大片免费观看网站| 麻豆乱淫一区二区| 亚洲性夜色夜夜综合| 9色porny在线观看| 涩涩av久久男人的天堂| 美女大奶头黄色视频| 国精品久久久久久国模美| 午夜影院在线不卡| 欧美亚洲 丝袜 人妻 在线| 男男h啪啪无遮挡| 少妇精品久久久久久久| 在线观看免费日韩欧美大片| 亚洲国产欧美网| 人妻 亚洲 视频| 一区二区av电影网| 一本久久精品| 1024香蕉在线观看| 桃花免费在线播放| 国产成+人综合+亚洲专区| 亚洲成人国产一区在线观看| 性少妇av在线| 狂野欧美激情性bbbbbb| 成年人免费黄色播放视频| 亚洲成国产人片在线观看| 欧美日韩黄片免| 久久精品国产亚洲av高清一级| 久久久久久久精品精品| 亚洲第一av免费看| tocl精华| 欧美人与性动交α欧美软件| 欧美激情 高清一区二区三区| 热99re8久久精品国产| 老熟妇乱子伦视频在线观看 | 一进一出抽搐动态| 午夜日韩欧美国产| 在线观看免费视频网站a站| 最黄视频免费看| 水蜜桃什么品种好| 成人亚洲精品一区在线观看| 一级毛片电影观看| 十八禁高潮呻吟视频| 亚洲专区中文字幕在线| 永久免费av网站大全| 亚洲 欧美一区二区三区| 久久精品国产亚洲av香蕉五月 | av超薄肉色丝袜交足视频| 日本撒尿小便嘘嘘汇集6| 黑人巨大精品欧美一区二区mp4| 一二三四社区在线视频社区8| 欧美黄色淫秽网站| 黄网站色视频无遮挡免费观看| 一个人免费在线观看的高清视频 | 一级a爱视频在线免费观看| 久久女婷五月综合色啪小说| 91大片在线观看| 国产在线一区二区三区精| 丝袜人妻中文字幕| av电影中文网址| 美女高潮到喷水免费观看| 多毛熟女@视频| 久久精品国产a三级三级三级| 青草久久国产| 天天躁日日躁夜夜躁夜夜| 精品少妇黑人巨大在线播放| 亚洲综合色网址| 在线亚洲精品国产二区图片欧美| 亚洲国产精品一区二区三区在线| 精品高清国产在线一区| 久久av网站| 在线av久久热| 亚洲国产成人一精品久久久| 亚洲精品日韩在线中文字幕| 国产在视频线精品| 亚洲av美国av| 国产精品免费大片| 十八禁高潮呻吟视频| 亚洲熟女精品中文字幕| 欧美日韩精品网址| 精品国产一区二区三区久久久樱花| 久久久久精品国产欧美久久久 | 两个人免费观看高清视频| 久久香蕉激情| 日韩 亚洲 欧美在线| 免费高清在线观看日韩| 777米奇影视久久| 飞空精品影院首页| 亚洲精品日韩在线中文字幕| 夜夜骑夜夜射夜夜干| 精品少妇黑人巨大在线播放| 人妻久久中文字幕网| 国产在线一区二区三区精| 后天国语完整版免费观看| 黄色视频在线播放观看不卡| 亚洲性夜色夜夜综合| 久久综合国产亚洲精品| 日韩精品免费视频一区二区三区| 久久久精品国产亚洲av高清涩受| 丝袜脚勾引网站| 美女脱内裤让男人舔精品视频| 久久青草综合色| 久久久欧美国产精品| 亚洲国产欧美在线一区| 99久久人妻综合| 制服诱惑二区| 亚洲精品第二区| 国产精品久久久久久人妻精品电影 | 成人免费观看视频高清| 亚洲欧洲精品一区二区精品久久久| 亚洲国产毛片av蜜桃av| 免费观看a级毛片全部| 国产精品一区二区精品视频观看| 丝袜美腿诱惑在线| 看免费av毛片| 亚洲男人天堂网一区| 欧美大码av| 午夜福利乱码中文字幕| 最新在线观看一区二区三区| 免费在线观看影片大全网站| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲人成电影观看| 无限看片的www在线观看| 久久久国产精品麻豆| 黄色a级毛片大全视频| 欧美变态另类bdsm刘玥| 久久久久久亚洲精品国产蜜桃av| 中文字幕制服av| 国产成人影院久久av| a 毛片基地| 亚洲欧美一区二区三区黑人| 精品高清国产在线一区| 精品人妻熟女毛片av久久网站| 国产高清国产精品国产三级| 久久狼人影院| 免费不卡黄色视频| 成人av一区二区三区在线看 | 老司机亚洲免费影院| 国产精品一区二区在线观看99| 在线观看免费高清a一片| 亚洲国产日韩一区二区| 亚洲av电影在线观看一区二区三区| 亚洲第一av免费看| 在线av久久热| 精品亚洲成a人片在线观看| 99久久精品国产亚洲精品| 老熟妇仑乱视频hdxx| 国产av又大| 国产视频一区二区在线看| 久久免费观看电影| 男人爽女人下面视频在线观看| 丝袜美腿诱惑在线| 日日爽夜夜爽网站| 日本精品一区二区三区蜜桃| 搡老乐熟女国产| 日韩 欧美 亚洲 中文字幕| 中文字幕色久视频| videos熟女内射| 少妇人妻久久综合中文| 天堂俺去俺来也www色官网| 国产黄频视频在线观看| 亚洲熟女毛片儿| 91国产中文字幕| 动漫黄色视频在线观看| 黄片大片在线免费观看| 欧美成人午夜精品| 亚洲精品自拍成人| 国产精品免费视频内射| 色精品久久人妻99蜜桃| 亚洲专区国产一区二区| 一级,二级,三级黄色视频| 成年人午夜在线观看视频| 天堂中文最新版在线下载| 老司机影院成人| av网站免费在线观看视频| 色老头精品视频在线观看| 69av精品久久久久久 | 美女扒开内裤让男人捅视频| 国产又爽黄色视频| 国产区一区二久久| av天堂在线播放| 菩萨蛮人人尽说江南好唐韦庄| 精品人妻1区二区| 人人妻人人澡人人看| 精品人妻1区二区| 99久久精品国产亚洲精品| 精品人妻一区二区三区麻豆| 桃花免费在线播放| 欧美 日韩 精品 国产| 波多野结衣一区麻豆| 悠悠久久av| bbb黄色大片| av天堂在线播放| 一本—道久久a久久精品蜜桃钙片| 国产野战对白在线观看| 18禁裸乳无遮挡动漫免费视频| 亚洲成人免费av在线播放| 亚洲国产欧美在线一区| 亚洲精品久久成人aⅴ小说| 欧美少妇被猛烈插入视频| 国产精品久久久人人做人人爽| 中文字幕最新亚洲高清| 午夜福利,免费看| 黄色怎么调成土黄色| 一二三四社区在线视频社区8| 亚洲精品av麻豆狂野| 人人妻人人澡人人爽人人夜夜| 91九色精品人成在线观看| 国产有黄有色有爽视频| 一级毛片女人18水好多| 老司机深夜福利视频在线观看 | 亚洲色图 男人天堂 中文字幕| 欧美日韩亚洲高清精品| 久久精品亚洲av国产电影网| 老司机在亚洲福利影院| 91国产中文字幕| 亚洲全国av大片| 国产男人的电影天堂91| 悠悠久久av| 一本—道久久a久久精品蜜桃钙片| 亚洲国产看品久久| 欧美国产精品一级二级三级| 一二三四社区在线视频社区8| 啦啦啦视频在线资源免费观看| 99精国产麻豆久久婷婷| www日本在线高清视频| 一级,二级,三级黄色视频| 丰满饥渴人妻一区二区三| 国产av又大| 亚洲av成人不卡在线观看播放网 | videos熟女内射| 美女扒开内裤让男人捅视频| 亚洲精品一二三| 久久 成人 亚洲| 精品免费久久久久久久清纯 | 欧美一级毛片孕妇| 亚洲视频免费观看视频| 精品少妇久久久久久888优播| 美女高潮到喷水免费观看| 啦啦啦在线免费观看视频4| 午夜福利免费观看在线| 又黄又粗又硬又大视频| 亚洲国产欧美一区二区综合| 另类精品久久| 99精品久久久久人妻精品| 国产在线一区二区三区精| 夜夜夜夜夜久久久久| 欧美精品高潮呻吟av久久| 亚洲精品国产色婷婷电影| 亚洲精品粉嫩美女一区| 午夜福利影视在线免费观看| 69av精品久久久久久 | 成年av动漫网址| 老司机午夜福利在线观看视频 | 成人国产一区最新在线观看| 欧美黄色淫秽网站| 少妇被粗大的猛进出69影院| 脱女人内裤的视频| 国产有黄有色有爽视频| 久久久久国产精品人妻一区二区| 久久久国产成人免费| 久久久国产欧美日韩av| 多毛熟女@视频| 成年人免费黄色播放视频| 欧美 亚洲 国产 日韩一| 欧美精品人与动牲交sv欧美| 久久人人爽人人片av| 99精品久久久久人妻精品| 国产精品麻豆人妻色哟哟久久| 久久午夜综合久久蜜桃| 国产男女超爽视频在线观看| 水蜜桃什么品种好| 国产熟女午夜一区二区三区| 国产精品久久久av美女十八|