Dynamic allocation method for airport apron removal
Technical Field
The invention belongs to the field of scheduling distribution, and particularly relates to a dynamic distribution method for airport deicing plateaus.
Background
In winter, the ice on the surface of the airplane can be frozen in cold ice and snow weather, the air transportation safety is damaged, and the deicing of the ground of the airplane is an important guarantee for the safe operation of flights in winter.
The distribution problem of the airport deicing plateaus is a highly real-time dynamic distribution problem, the domestic airports basically still adopt a manual distribution method, however, with the continuous expansion of the airport scale, flights needing deicing in cold ice and snow weather in winter grow exponentially, which brings huge workload to workers, and the manual distribution is very easy to make mistakes during the departure peak period on duty, so that the intellectualization and the convenience of the distribution of the airport deicing plateaus are urgently realized. The problem of airport deicing apron distribution is solved, airport resources can be fully utilized, and airport operation benefits and service quality are improved.
In the current stage of research on the distribution problem of the airport ice lawn removal, most of the methods such as the traditional genetic algorithm, the queuing theory, the game theory and the like are adopted, and the method stays in the static distribution stage. However, in actual work, the distributor is highly dependent on information transmitted by the air traffic control in real time and information obtained from the radar of the field monitor in real time when the ice terrace is distributed. When a flight applies for deicing, an allocator allocates one deicing apron in real time according to the current occupancy of each deicing apron and the slide way, and continuously adjusts the deicing aprons according to real-time information, which is a highly real-time dynamic allocation process.
Disclosure of Invention
Aiming at the defects of the prior art, the invention aims to provide a dynamic allocation method of airport deicing aprons, and aims to solve the problem of low static allocation efficiency of airport deicing aprons.
In order to achieve the above object, the present invention provides a dynamic allocation method for airport ice apron removal, comprising:
s1: screening the ice terrace removing set Q selectable by the flight to be distributed according to the hard service constraint0;
S2: obtaining a set Q for removing ice aprons which do not meet the time constraint by calculating the predicted time for entering the aprons and the predicted time for leaving the aprons of the flights to be distributed1;
S3: judging ice terrace removing set Q1If it is empty, if the ice terrace is removed, set Q1If the set is not an empty set, the set Q is in the ice terrace removing set1Screening the deicing plateaus with the earliest estimated plateau time to obtain a distribution result, and turning to the step S4; otherwise, after delaying the assignment of the flight's time to enter the plateau, divideIceground set Q0Screening the deicing plateaus with the earliest estimated plateau time to obtain a distribution result, and turning to the step S4;
s4: before the flight leaves the airport, judging whether the flight state changes in real time, if so, turning to the step S2, otherwise, maintaining the distribution result;
preferably, the hard traffic constraint comprises: the number of the ice aprons corresponds to the model constraint, the unit constraint of the flights to which the units of the ice aprons belong corresponds to the unit constraint of the flights to which the units of the ice aprons belong, and the position constraint of the ice aprons to which the positions of the flight berths correspond to the positions of the ice aprons.
The step S2 obtains the ice terrace detachment set Q1The method comprises the following steps:
s2.1 calculating set Q of ice terrace removing0The flight sets which are distributed in all the ice aprons and are not empowered are determined;
s2.2, judging whether the flight set which is distributed in each deicing station and has no station is empty, and if the flight set is not empty, turning to the step S2.3;
s2.3, respectively calculating the estimated apron entering time and the estimated apron exiting time of the flights to be distributed and the distributed and undistributed flight sets;
s2.4, according to the calculation in the step S2.3, judging whether the flight to be distributed and the flight without taking the flight have overlapping of the occupied time, if so, removing the ice pad corresponding to the flight without taking the flight to obtain an ice pad set Q1。
Preferably, the method for calculating the estimated entering time and the estimated exiting time of the flight to be allocated comprises the following steps:
(1) judging the current position of the flight, and if the flight sends an deicing application but does not slide out of the parking space, turning to the step (2); if the flight has already slipped out of the parking space but has not entered the floor, go to step (3); if the flight has entered the level, go to step (4);
(2) according to the current time T and the application time T for sending deicing to the flightaCalculating the estimated time T to enter the plateau of the flight to be allocatedeAnd the estimated plateau time Tl;
Wherein, Tl=Te+T3;
If Ta+T1T, then Te=Ta+T1+T2(ii) a If Ta+T1If < T, then Te=t;
(3) According to the current time T and the time T of the flight sliding out of the standsCalculating the estimated time T to enter the plateau of the flight to be allocatedeAnd the estimated plateau time Tl;
Wherein, Tl=Te+T3;
If Ts+T2T, then Te=Ts+T2(ii) a If Ts+T2If < T, then Te=t;
(4) Calculating the predicted entering time T of the flight to be distributed according to the current time TeAnd the estimated plateau time Tl;
Wherein if Te+T3T, then Tl=Te+T3(ii) a If Te+T3If < T, then Tl=t;
Wherein, T1Time from applying for de-icing to docking for the flight; t is2Time for flights from berthing to each de-icepad; t is3The flight takes time to remove the ice pad.
The method for delaying the predicted entering time of the distributed flights comprises the following steps:
when removing plateau set Q1When the flight is empty, the estimated entering time of the flight to be allocated for removing each deicing plateau is the estimated exiting time of the previous flight allocated to each deicing plateau +1 min.
The change in flight status includes: the flight status changes from normal to urgent and the flight is not settled, or the urgent changes to delayed and the flight is not settled.
The method for judging the flight state comprises the following steps:
if the estimated time of taking out the flight to be distributed is later than the estimated taking-off time and the time difference is less than the threshold time, the flight is a normal flight;
if the current time is later than the expected takeoff time and the time difference is greater than the threshold time, the flight is a delayed flight;
otherwise, the flight is an emergency flight.
Preferably, the threshold time is 25 minutes.
Through the technical scheme, compared with the prior art, the invention can obtain the following advantages
Has the advantages that:
(1) on one hand, the invention obtains the ice terrace removing set Q by considering hardware business constraint and time constraint1Then, in order to improve distribution efficiency, a ice terrace removing set Q is further screened1The deicing plateau with the earliest plateau time is predicted; on the other hand, the estimated time for entering and leaving the flight of the flight to be distributed is updated by monitoring the change of the actual state of the flight to be distributed, so that the assignment scheme is timely adjusted while the flight to be distributed is monitored in real time, and the working efficiency of an assignor is greatly improved.
(2) Firstly, comprehensively considering hard service constraint and time constraint, selecting the optimal ice terrace corresponding to the flight to be distributed, continuously screening and determining an ice terrace set for three times, and completing a pre-distribution process; and then, the estimated time for entering and leaving the flight corresponding to the flight to be distributed is fed back and updated in real time in combination with the actual flight state, and the ice terrace corresponding to the flight to be distributed is redistributed.
Drawings
FIG. 1 is a flow chart of a method for airport de-icing apron assignment provided by the present invention;
FIG. 2 shows a set Q of screened deicing plateaus provided by the present invention1Is a schematic flow diagram.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
As shown in fig. 1, the present invention provides a dynamic allocation method for airport ice apron removal, comprising:
s1: screening the ice terrace removing set Q selectable by the flight to be distributed according to the hard service constraint0;
S2: obtaining a set Q for removing ice aprons which do not meet the time constraint by calculating the predicted time for entering the aprons and the predicted time for leaving the aprons of the flights to be distributed1;
S3: judging ice terrace removing set Q1If it is empty, if the ice terrace is removed, set Q1If the set is not an empty set, the set Q is in the ice terrace removing set1Screening the deicing plateaus with the earliest estimated plateau time to obtain a distribution result, and turning to the step S4; otherwise, delaying the time of allocating flights to enter the flight level, and collecting Q in the ice level removing set0Screening the deicing plateaus with the earliest estimated plateau time to obtain a distribution result, and turning to the step S4;
generally, if the plateau is removed, Q1When the flight is empty, adjusting the estimated entering time of the flight to be distributed to be one minute after the departure time of the flight on each ice removal plateau;
s4: and before the flight leaves the airport, judging whether the flight state changes in real time, if so, turning to the step S2, otherwise, maintaining the distribution result. Preferably, the hard traffic constraint comprises: the number of the ice aprons corresponds to the model constraint, the unit constraint of the flights to which the units of the ice aprons belong corresponds to the unit constraint of the flights to which the units of the ice aprons belong, and the position constraint of the ice aprons to which the positions of the flight berths correspond to the positions of the ice aprons.
As shown in fig. 2, the above step S2 obtains the ice pad detachment set Q1The method of (1), comprising:
s2.1 calculating set Q of ice terrace removing0The flight sets which are distributed in all the ice aprons and are not empowered are determined;
s2.2, judging whether the flight set which is distributed in each deicing station and has no station is empty, and if the flight set is not empty, turning to the step S2.3;
s2.3, respectively calculating the estimated apron entering time and the estimated apron exiting time of the flights to be distributed and the distributed and undistributed flight sets;
s2.4, according to the calculation in the step S2.3, judging whether the occupied space time of the flight to be distributed and the flight without leaving the flight are overlapped, if so, removing the ice terrace corresponding to the distributed flight without leaving the flight, and acquiring an ice terrace set Q1。
Preferably, the method for calculating the estimated entering time and the estimated exiting time of the flight to be allocated comprises the following steps:
(1) judging the current position of the flight, and if the flight has a deicing application but does not slide out of the parking space, turning to the step (2); if the flight has already slipped out of the parking space but has not entered the floor, go to step (3); if the flight has entered the level, go to step (4);
(2) according to the current time T and the application time T for sending deicing to the flightaCalculating the estimated time T to enter the plateau of the flight to be allocatedeAnd the estimated plateau time Tl;
Wherein, Tl=Te+T3;
If Ta+T1T, then Te=Ta+T1+T2(ii) a If Ta+T1If < T, then Te=t;
(3) According to the current time T and the time T of the flight sliding out of the standsCalculating the estimated time T to enter the plateau of the flight to be allocatedeAnd the estimated plateau time Tl;
Wherein, Tl=Te+T3;
If Ts+T2T, then Te=Ts+T2(ii) a If Ts+T2If < T, then Te=t;
(4) Calculating the predicted entering time T of the flight to be distributed according to the current time TeAnd the estimated plateau time Tl;
Wherein if Te+T3T, then Tl=Te+T3(ii) a If Te+T3If < T, then Tl=t;
Wherein, T1Time from applying for de-icing to sliding out of the stand for the flight; t is2The time for the flight from the slide-out flight level to each ice pad; t is3The flight takes time to remove the ice pad.
T1The method comprises the steps of calculating the average time from applying for deicing to stopping of a flight to be distributed in one month to obtain the average time; t is3The initial value is the time of taking the flight of the first flight in each ice terrace on the day, and each ice terrace in one month is calculatedAverage take-over time for flights in the ice terrace, starting with the second flight, T3The mean of the time of flight in the plateau is the previous flight.
Preferably, the method for delaying the estimated approach time of the flight to be allocated comprises: when removing plateau set Q1When the flight is empty, the estimated entering time of the flight to be allocated for removing each deicing plateau is the estimated exiting time of the previous flight allocated to each deicing plateau +1 min.
The change in flight status includes: the flight status changes from normal to urgent and the flight is not settled, or the urgent changes to delayed and the flight is not settled.
The method for judging the flight state comprises the following steps:
if the estimated time of taking the flight to be distributed is later than the estimated taking-off time and the time difference is less than the threshold time, the flight is a normal flight;
if the current time is later than the expected takeoff time and the time difference is greater than the threshold time, the flight is a delayed flight;
otherwise, the flight is an emergency flight.
Preferably, the threshold time is 25 minutes.
It will be understood by those skilled in the art that the foregoing is only a preferred embodiment of the present invention, and is not intended to limit the invention, and that any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the scope of the present invention.