• <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香蕉五月| 黄片wwwwww| 国产精品不卡视频一区二区| 成熟少妇高潮喷水视频| 久久精品影院6| 搡老岳熟女国产| 国产三级中文精品| 91午夜精品亚洲一区二区三区 | 联通29元200g的流量卡| 99久久无色码亚洲精品果冻| 亚洲美女黄片视频| 自拍偷自拍亚洲精品老妇| 亚洲国产高清在线一区二区三| 欧美人与善性xxx| 在线看三级毛片| 在线免费十八禁| 免费看av在线观看网站| 欧美性猛交╳xxx乱大交人| 国产免费av片在线观看野外av| 美女免费视频网站| 久久热精品热| 变态另类成人亚洲欧美熟女| 久99久视频精品免费| 少妇的逼好多水| 最新中文字幕久久久久| 国内少妇人妻偷人精品xxx网站| 国产 一区 欧美 日韩| a在线观看视频网站| 深夜精品福利| 无遮挡黄片免费观看| h日本视频在线播放| 国产色爽女视频免费观看| 国产精品久久久久久av不卡| 国产精品,欧美在线| 国产真实伦视频高清在线观看 | 亚洲精品成人久久久久久| 欧美日韩国产亚洲二区| 干丝袜人妻中文字幕| 日韩av在线大香蕉| 国产精品,欧美在线| 国产伦精品一区二区三区视频9| 国产高潮美女av| 成人国产综合亚洲| 欧美中文日本在线观看视频| 亚洲av熟女| 欧美性猛交╳xxx乱大交人| 国内毛片毛片毛片毛片毛片| 色哟哟哟哟哟哟| 精品久久久久久久人妻蜜臀av| 99热这里只有是精品在线观看| 一区二区三区免费毛片| 日本 欧美在线| 国产成人福利小说| 午夜福利欧美成人| 少妇丰满av| 中文亚洲av片在线观看爽| 国产亚洲精品av在线| 亚洲国产精品久久男人天堂| 亚洲av成人av| 91午夜精品亚洲一区二区三区 | 村上凉子中文字幕在线| 国产精品亚洲一级av第二区| 免费观看在线日韩| 日韩国内少妇激情av| 欧美绝顶高潮抽搐喷水| 午夜视频国产福利| 观看免费一级毛片| av黄色大香蕉| 美女 人体艺术 gogo| 18禁黄网站禁片午夜丰满| 久久午夜亚洲精品久久| 亚洲欧美日韩高清专用| 久久天躁狠狠躁夜夜2o2o| 亚洲av成人av| av在线蜜桃| av天堂中文字幕网| 成年版毛片免费区| 亚洲av二区三区四区| 免费无遮挡裸体视频| 国产久久久一区二区三区| 久久久久性生活片| 18禁裸乳无遮挡免费网站照片| 一个人免费在线观看电影| 亚洲内射少妇av| 丰满的人妻完整版| 亚洲自拍偷在线| 波多野结衣高清无吗| 熟妇人妻久久中文字幕3abv| 日日摸夜夜添夜夜添av毛片 | 成人国产综合亚洲| 少妇的逼水好多| 亚洲不卡免费看| 51国产日韩欧美| 欧美成人a在线观看| 欧美日韩中文字幕国产精品一区二区三区| avwww免费| 久久欧美精品欧美久久欧美| 全区人妻精品视频| 欧美日韩国产亚洲二区| 欧美潮喷喷水| 中文字幕人妻熟人妻熟丝袜美| 国产成人aa在线观看| 国产精品爽爽va在线观看网站| www.色视频.com| 久久人人精品亚洲av| 中文字幕av成人在线电影| 久久久久国产精品人妻aⅴ院| 在线国产一区二区在线| 国产精品不卡视频一区二区| 国产精品亚洲美女久久久| 高清在线国产一区| av女优亚洲男人天堂| 性插视频无遮挡在线免费观看| 99久久无色码亚洲精品果冻| 精品久久久久久久末码| 欧美xxxx黑人xx丫x性爽| 国内揄拍国产精品人妻在线| 高清日韩中文字幕在线| 极品教师在线免费播放| 99riav亚洲国产免费| 最近最新免费中文字幕在线| 日日干狠狠操夜夜爽| 黄片wwwwww| 日日啪夜夜撸| av天堂中文字幕网| 亚洲中文字幕日韩| 夜夜看夜夜爽夜夜摸| 成人av一区二区三区在线看| 欧美激情国产日韩精品一区| 亚洲国产日韩欧美精品在线观看| 精品99又大又爽又粗少妇毛片 | 可以在线观看毛片的网站| 久久精品国产亚洲av天美| 精品99又大又爽又粗少妇毛片 | 3wmmmm亚洲av在线观看| 色综合亚洲欧美另类图片| 欧美日韩精品成人综合77777| 日本黄色视频三级网站网址| 动漫黄色视频在线观看| 午夜福利在线观看免费完整高清在 | 午夜久久久久精精品| 亚洲不卡免费看| 99热这里只有精品一区| 九九在线视频观看精品| 婷婷精品国产亚洲av| 少妇高潮的动态图| 黄色配什么色好看| 一区二区三区激情视频| 国产黄片美女视频| 国产伦一二天堂av在线观看| 窝窝影院91人妻| 日本一本二区三区精品| 国产综合懂色| 中文在线观看免费www的网站| 成人亚洲精品av一区二区| 久久精品国产99精品国产亚洲性色| 欧美xxxx性猛交bbbb| 免费观看的影片在线观看| 毛片一级片免费看久久久久 | 性色avwww在线观看| 成人无遮挡网站| 91久久精品国产一区二区成人| 午夜亚洲福利在线播放| 一区二区三区免费毛片| 春色校园在线视频观看| 简卡轻食公司| 亚洲精华国产精华精| 男女边吃奶边做爰视频| 久久精品夜夜夜夜夜久久蜜豆| 日本 av在线| 久久久久久久午夜电影| 美女cb高潮喷水在线观看| 国产精品美女特级片免费视频播放器| 日本欧美国产在线视频| 国产麻豆成人av免费视频| 欧美日韩亚洲国产一区二区在线观看| 欧美潮喷喷水| 国产视频一区二区在线看| 永久网站在线| 搡老熟女国产l中国老女人| 亚洲三级黄色毛片| 十八禁国产超污无遮挡网站| 国产亚洲精品综合一区在线观看| 小说图片视频综合网站| 91久久精品电影网| 五月伊人婷婷丁香| 欧美日本亚洲视频在线播放| 草草在线视频免费看| 欧美一区二区国产精品久久精品| 十八禁网站免费在线| 亚洲精品国产成人久久av| 午夜老司机福利剧场| 蜜桃亚洲精品一区二区三区| 国内精品久久久久精免费| 国产精品不卡视频一区二区| 成人永久免费在线观看视频| 天堂影院成人在线观看| 久久久精品大字幕| 亚洲无线观看免费| 91午夜精品亚洲一区二区三区 | 日韩精品有码人妻一区| 男人舔女人下体高潮全视频| 丰满的人妻完整版| 国产精品人妻久久久影院| or卡值多少钱| 日本a在线网址| 亚洲欧美清纯卡通| bbb黄色大片| 成人二区视频| 两人在一起打扑克的视频| 精品午夜福利视频在线观看一区| 亚洲人成网站在线播| 亚洲av不卡在线观看| 欧美日韩亚洲国产一区二区在线观看| 日本三级黄在线观看| 精品久久久久久久久久久久久| 中文字幕人妻熟人妻熟丝袜美| 国产成人a区在线观看| 日韩一区二区视频免费看| 亚洲国产日韩欧美精品在线观看| 特级一级黄色大片| 日韩精品中文字幕看吧| 国产精品98久久久久久宅男小说| 日本黄色视频三级网站网址| 成人毛片a级毛片在线播放| 日本五十路高清| 男人舔女人下体高潮全视频| 国产精品一及| 长腿黑丝高跟| 精品乱码久久久久久99久播| 亚洲精品国产成人久久av| 国产伦在线观看视频一区| 大型黄色视频在线免费观看| 丰满的人妻完整版| 一区二区三区免费毛片| 精品午夜福利视频在线观看一区| 伦精品一区二区三区| 精品不卡国产一区二区三区| 国产精品一区二区性色av| 一进一出好大好爽视频| 男人的好看免费观看在线视频| 一本精品99久久精品77| 久久久久免费精品人妻一区二区| 日本黄大片高清| 久久久久久伊人网av| 欧美xxxx黑人xx丫x性爽| 日韩欧美在线二视频| 哪里可以看免费的av片| 一级毛片久久久久久久久女| 99精品久久久久人妻精品| 成年人黄色毛片网站| 久久久久性生活片| 18禁黄网站禁片午夜丰满| 黄色欧美视频在线观看| 婷婷精品国产亚洲av| 男女做爰动态图高潮gif福利片| 国产精品精品国产色婷婷| 99精品久久久久人妻精品| 日本撒尿小便嘘嘘汇集6| 国产久久久一区二区三区| 久久久精品大字幕| 一个人看的www免费观看视频| 亚洲美女黄片视频| 俄罗斯特黄特色一大片| 国产日本99.免费观看| 精品国内亚洲2022精品成人| 日本成人三级电影网站| 日韩一本色道免费dvd| 亚洲av.av天堂| а√天堂www在线а√下载| 啦啦啦啦在线视频资源| 亚洲人成网站高清观看| h日本视频在线播放| 深夜a级毛片| 搡老妇女老女人老熟妇| 51国产日韩欧美| 国产色婷婷99| 成人永久免费在线观看视频| 少妇熟女aⅴ在线视频| 内射极品少妇av片p| 中文字幕av在线有码专区| 老司机深夜福利视频在线观看| 日韩一本色道免费dvd| 国产亚洲91精品色在线| 日韩欧美免费精品| 麻豆国产av国片精品| 免费看日本二区| 国产伦精品一区二区三区四那| 婷婷精品国产亚洲av| 自拍偷自拍亚洲精品老妇| 国产亚洲91精品色在线| 天天躁日日操中文字幕| 夜夜夜夜夜久久久久| 精品久久久噜噜| 久久精品国产亚洲av涩爱 | 人人妻人人澡欧美一区二区| www.色视频.com| 少妇的逼好多水| 少妇猛男粗大的猛烈进出视频 | 午夜精品在线福利| 51国产日韩欧美| 老女人水多毛片| 精品福利观看| 午夜影院日韩av| 久久久久久伊人网av| 啪啪无遮挡十八禁网站| 蜜桃久久精品国产亚洲av| 成人永久免费在线观看视频| 少妇丰满av| 日韩av在线大香蕉| 在线免费观看不下载黄p国产 | 国产免费一级a男人的天堂| 国产一区二区三区av在线 | 日日摸夜夜添夜夜添av毛片 | 成年女人永久免费观看视频| 女人十人毛片免费观看3o分钟| 精品无人区乱码1区二区| 欧美一级a爱片免费观看看| 白带黄色成豆腐渣| 日日啪夜夜撸| 97超视频在线观看视频| 丝袜美腿在线中文| 成人鲁丝片一二三区免费| 亚洲精华国产精华液的使用体验 | 美女高潮喷水抽搐中文字幕| 99在线人妻在线中文字幕| 免费观看人在逋| 最新在线观看一区二区三区| 成人高潮视频无遮挡免费网站| 美女免费视频网站| 女人十人毛片免费观看3o分钟| 国产又黄又爽又无遮挡在线| 国产高清不卡午夜福利| eeuss影院久久| 免费看美女性在线毛片视频| 直男gayav资源| 国产白丝娇喘喷水9色精品| 美女大奶头视频| 搡老妇女老女人老熟妇| 午夜福利成人在线免费观看| 黄色日韩在线| 男女做爰动态图高潮gif福利片| 国产高清有码在线观看视频| 两人在一起打扑克的视频| 国产精品99久久久久久久久| 亚洲精品在线观看二区| 午夜亚洲福利在线播放| 国产精品久久电影中文字幕| 国产成人一区二区在线| 三级男女做爰猛烈吃奶摸视频| 久久精品人妻少妇| 精品乱码久久久久久99久播| 别揉我奶头~嗯~啊~动态视频| 真人做人爱边吃奶动态| 亚洲av中文av极速乱 | 日本撒尿小便嘘嘘汇集6| 久久亚洲精品不卡| 免费搜索国产男女视频| 日韩大尺度精品在线看网址| 国产三级中文精品| 亚洲熟妇熟女久久| 麻豆av噜噜一区二区三区| 日韩在线高清观看一区二区三区 | 国产精品国产三级国产av玫瑰| 日韩高清综合在线| 亚洲av不卡在线观看| 少妇熟女aⅴ在线视频| 日本 av在线| 69av精品久久久久久| 51国产日韩欧美| 两性午夜刺激爽爽歪歪视频在线观看| 丰满的人妻完整版| 久久久久九九精品影院| 国产大屁股一区二区在线视频| 三级男女做爰猛烈吃奶摸视频| av国产免费在线观看| 偷拍熟女少妇极品色| 国产精品伦人一区二区| 熟女电影av网| 精品久久久久久久久av| 亚洲精品456在线播放app | 一区二区三区四区激情视频 | 少妇丰满av| 久久国内精品自在自线图片| 赤兔流量卡办理| 久久热精品热| 国国产精品蜜臀av免费| 美女高潮的动态| 亚洲成人免费电影在线观看| 国产一区二区亚洲精品在线观看| 亚洲成a人片在线一区二区| 久久精品国产亚洲av涩爱 | 人妻少妇偷人精品九色| 国产色婷婷99| 午夜福利高清视频| 欧美日韩乱码在线| 国产在视频线在精品| 大型黄色视频在线免费观看| 91在线精品国自产拍蜜月| 亚洲一区高清亚洲精品| 国产精品不卡视频一区二区| 日韩欧美精品v在线| 欧美zozozo另类| 尤物成人国产欧美一区二区三区| 国产高清不卡午夜福利| 亚洲av免费在线观看| 亚洲va在线va天堂va国产| 亚洲精华国产精华液的使用体验 | 国产亚洲av嫩草精品影院| 久久久久久久久久久丰满 | 欧美又色又爽又黄视频| 嫩草影院新地址| 欧美日韩亚洲国产一区二区在线观看| 欧美区成人在线视频| a级毛片免费高清观看在线播放| 久久精品国产亚洲网站| 欧美+日韩+精品| 3wmmmm亚洲av在线观看| 一本一本综合久久| 欧美激情国产日韩精品一区| 一a级毛片在线观看| 亚洲欧美清纯卡通| 99热精品在线国产| 草草在线视频免费看| 两人在一起打扑克的视频| 国产精品98久久久久久宅男小说| 欧美精品国产亚洲| 成人亚洲精品av一区二区| 国产一区二区三区视频了| 超碰av人人做人人爽久久| 级片在线观看| 国产乱人视频| 精品午夜福利在线看| av女优亚洲男人天堂| 日本熟妇午夜| 亚洲熟妇中文字幕五十中出| 久久草成人影院| 不卡一级毛片| 免费观看人在逋| 亚洲天堂国产精品一区在线| 成人毛片a级毛片在线播放| 久久久久久大精品| 精品日产1卡2卡| 免费观看在线日韩| 露出奶头的视频| 亚洲人成网站在线播| 永久网站在线| 欧美一区二区亚洲| 日韩欧美在线二视频| 色哟哟哟哟哟哟| 在线观看美女被高潮喷水网站| 成年女人毛片免费观看观看9| 成年女人看的毛片在线观看| 18禁在线播放成人免费| 蜜桃久久精品国产亚洲av| 永久网站在线| 成人特级av手机在线观看| 一级av片app| 国产白丝娇喘喷水9色精品| 国产精品人妻久久久久久| 婷婷亚洲欧美| 五月玫瑰六月丁香| 精品午夜福利视频在线观看一区| 亚洲国产高清在线一区二区三| 乱人视频在线观看| 色尼玛亚洲综合影院| 国产极品精品免费视频能看的| av在线天堂中文字幕| 高清在线国产一区| 欧美又色又爽又黄视频| 亚洲最大成人av| 精品日产1卡2卡| 日本黄大片高清| 亚洲综合色惰| 国产人妻一区二区三区在| 人妻制服诱惑在线中文字幕| 99久久中文字幕三级久久日本| 自拍偷自拍亚洲精品老妇| 欧美xxxx黑人xx丫x性爽| 成人精品一区二区免费| 男女做爰动态图高潮gif福利片| 可以在线观看的亚洲视频| 亚洲最大成人中文| 99在线人妻在线中文字幕| 午夜福利成人在线免费观看| 欧美一区二区精品小视频在线| 最近中文字幕高清免费大全6 | 欧美不卡视频在线免费观看| 国产一区二区在线av高清观看| 九九爱精品视频在线观看| 午夜福利在线在线| 国产美女午夜福利| 国产亚洲精品久久久久久毛片| 我要看日韩黄色一级片| 亚洲图色成人| 国产爱豆传媒在线观看| 国产精品精品国产色婷婷| 亚洲三级黄色毛片| 久久久久久久久大av| 蜜桃亚洲精品一区二区三区| 婷婷精品国产亚洲av在线| 久久天躁狠狠躁夜夜2o2o| 亚洲va日本ⅴa欧美va伊人久久| 综合色av麻豆| av专区在线播放| АⅤ资源中文在线天堂| 精品99又大又爽又粗少妇毛片 | 国模一区二区三区四区视频| 少妇熟女aⅴ在线视频| 亚洲狠狠婷婷综合久久图片| 成人一区二区视频在线观看| 亚洲精品日韩av片在线观看| 99久久精品国产国产毛片| 免费看光身美女| 日韩强制内射视频| 亚洲三级黄色毛片| 丰满的人妻完整版| 成年免费大片在线观看| 身体一侧抽搐| 国产精品久久久久久av不卡| 久久精品国产亚洲av涩爱 | 国产午夜精品论理片| 欧美性猛交黑人性爽| 99热精品在线国产| 九色成人免费人妻av| 如何舔出高潮| 欧美日本亚洲视频在线播放| 国产真实乱freesex| 精品一区二区三区人妻视频| 最后的刺客免费高清国语| h日本视频在线播放| 久久6这里有精品| 搡老熟女国产l中国老女人| 国产一区二区三区在线臀色熟女| 美女被艹到高潮喷水动态| 色吧在线观看| 精品国内亚洲2022精品成人| 久久久久久久午夜电影| 欧美日韩国产亚洲二区| 舔av片在线| 欧美最黄视频在线播放免费| 亚洲一级一片aⅴ在线观看| 欧美另类亚洲清纯唯美| 深夜a级毛片| 在线看三级毛片| 一区二区三区高清视频在线| 嫁个100分男人电影在线观看| 免费无遮挡裸体视频| 日本a在线网址| 国产伦在线观看视频一区| 成人特级av手机在线观看| 精品久久久噜噜| 99久久中文字幕三级久久日本| 亚洲图色成人| 久久久成人免费电影| 亚洲男人的天堂狠狠| 亚洲成人精品中文字幕电影| 午夜激情福利司机影院| 国产一区二区三区视频了| 男插女下体视频免费在线播放| 国产精品1区2区在线观看.| 天美传媒精品一区二区| 欧美日本视频| 久久精品影院6| 国内精品久久久久精免费| 亚洲精品粉嫩美女一区| 欧美潮喷喷水| 久久精品影院6| 九色国产91popny在线| 亚洲综合色惰| 天堂影院成人在线观看| 久久99热这里只有精品18| 国产精品久久久久久久久免| 婷婷精品国产亚洲av在线| 国产麻豆成人av免费视频| 欧美日韩黄片免| 99在线视频只有这里精品首页| 1000部很黄的大片| 国产淫片久久久久久久久| 国产极品精品免费视频能看的| 久久天躁狠狠躁夜夜2o2o| 国产白丝娇喘喷水9色精品| 亚洲国产精品成人综合色| 国产精品国产高清国产av| 欧美一区二区亚洲| 中国美白少妇内射xxxbb| 九九热线精品视视频播放| 天堂av国产一区二区熟女人妻| 亚洲精品国产成人久久av| 欧美日本视频| 最近中文字幕高清免费大全6 | 精品一区二区三区人妻视频| 亚洲最大成人手机在线| 不卡一级毛片| 亚洲国产色片| 日韩人妻高清精品专区| 午夜精品久久久久久毛片777| 欧美三级亚洲精品| 嫁个100分男人电影在线观看| 亚洲国产高清在线一区二区三| 亚洲最大成人手机在线| 免费在线观看影片大全网站| 亚洲欧美日韩高清专用| 亚洲中文字幕日韩| 国产亚洲精品久久久com| 日本黄大片高清| 午夜爱爱视频在线播放| 亚洲精品亚洲一区二区| 国内精品久久久久久久电影| 日本五十路高清| 国产精品久久久久久久久免| 国产免费av片在线观看野外av| 国产在线精品亚洲第一网站| 久久精品综合一区二区三区| 在线观看舔阴道视频|