CN111291888A - Scheduling optimization method for airport special vehicles - Google Patents
Scheduling optimization method for airport special vehicles Download PDFInfo
- Publication number
- CN111291888A CN111291888A CN202010068912.8A CN202010068912A CN111291888A CN 111291888 A CN111291888 A CN 111291888A CN 202010068912 A CN202010068912 A CN 202010068912A CN 111291888 A CN111291888 A CN 111291888A
- Authority
- CN
- China
- Prior art keywords
- flight
- service
- vehicle
- time
- ferry
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000005457 optimization Methods 0.000 title claims abstract description 18
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 15
- 210000000349 chromosome Anatomy 0.000 claims description 30
- 230000010006 flight Effects 0.000 claims description 18
- 230000035772 mutation Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000000295 fuel oil Substances 0.000 claims description 4
- 238000010187 selection method Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000003921 oil Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Genetics & Genomics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- Development Economics (AREA)
- Evolutionary Computation (AREA)
- Game Theory and Decision Science (AREA)
- General Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Entrepreneurship & Innovation (AREA)
- Biomedical Technology (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Artificial Intelligence (AREA)
- Physiology (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Traffic Control Systems (AREA)
Abstract
The invention relates to the field of airport vehicle scheduling, in particular to an airport special vehicle scheduling optimization method, which is characterized by comprising the following steps: and constructing a combined scheduling model of the flight refueling service and the customer service based on the time constraint relation of the refueling service and the customer service, and solving the model by using an NSGA-II algorithm. Compared with manual scheduling, the invention has the following advantages: the oil consumption of the special vehicle and the vehicle circulation time are greatly saved; the use cost of the vehicle is reduced; compared with the single scheduling, the joint scheduling of the refueling truck and the ferry vehicle is as follows: the using quantity of the vehicles is greatly reduced, the running distance is relatively reduced, and the difference of the arrival time difference is not large; the vehicle scheduling optimization method is not only limited to vehicle scheduling in airports, but also can be popularized to vehicle scheduling in other fields. Therefore, the dispatching scheme designed by the combined dispatching of the refueling truck and the ferry vehicle is superior to manual dispatching and single dispatching in comprehensive consideration.
Description
Technical Field
The invention relates to the field of airport vehicle scheduling, in particular to an airport special vehicle scheduling optimization method.
Background
With the increasing expansion of the size and the traffic volume of airports, the bottleneck problems of low operation efficiency and insufficient collaborative decision capability of busy airports are increasingly highlighted. How to quickly respond to basic business and provide decision support for related departments becomes a key problem for realizing efficient operation of airports. The airplane punctuality rate is one of important index indexes for evaluating the aviation service quality. The air transportation punctuality rate can be influenced by a plurality of factors, wherein airport scheduling and other factors account for nearly 30%, and the punctuality rate of the flight is greatly influenced.
At present, the dispatching mode of special vehicles at airports by most airports in China is mainly manual dispatching and is a dispatching mode of single-vehicle single-flight service. The scheduling method has extremely low efficiency, and the possibility of flight delay caused by the fact that path optimization is not considered is very high. Because special vehicles are generally high in cost, an airport only provides limited vehicles to complete service, and the method of one-time service of one vehicle also causes waste of resource cost. For large hub airports, or for medium and small airports at peak transportation, the take-off and landing of flights can reach higher densities within a short period of time. Therefore, an effective airport special vehicle dispatching mode needs to be found in colleges and universities, and the method has great economic and social significance for improving the utilization rate of airport resources and improving the service quality of airports. The existing research is developed aiming at single type guarantee vehicle scheduling, Anna Norin, Mao X, Xing Z and the like sequentially research the deicing vehicle scheduling, and Anna Norin and Di Yuan adopt a greedy random self-adaptive search algorithm to search the optimal running path for the deicing vehicle to work, so that flight delay caused by overlong deicing time is reduced; mao X and the like propose a traditional centralized deicing scheduling method which cannot solve the problem of autonomy of different interest-related parties, and a new scheduling method based on Agent is proposed for the problem, so that the method can better solve the problem of airport deicing scheduling; and Xing Z and the like research the scheduling problem of the deicing vehicle by using related knowledge of the game theory, and verify that the research of the deicing scheduling problem based on the game theory method is feasible through experiments. JiaYan Du et al studied the problem of trailer dispatch in flight ground support services. Establishing a MIP model related to trailer dispatching; angus Cheung et al studied the individual dispatch of three different airport special vehicles, namely, clean water vehicles, tractors and cleaning vehicles; the problem of dynamic refueling truck scheduling and the problem of individual baggage truck scheduling were studied in succession by the paddle army, yangxedong et al; the yellowbird poems are subjected to simulation experiments by means of SIMIO simulation software, and the dispatching of the refueling truck and the ferry vehicle is respectively researched.
Most of the existing research aims at single type of guarantee vehicles, and service constraint relations among different types of guarantee vehicles are not considered. The rigor of the airport special vehicle service process cannot be better reflected.
Disclosure of Invention
Aiming at the defects, the invention provides a method for jointly scheduling a refueling truck and a ferry vehicle.
The technical scheme adopted by the invention for solving the technical problems is as follows: a scheduling optimization method for airport special vehicles is characterized in that a combined scheduling model of flight refueling service and passenger service is constructed based on a time constraint relation of refueling service and passenger service, and the model is solved by using an NSGA-II algorithm.
The combined scheduling model of flight refueling service and customer service is specifically described as follows: has m1Refuelling vehicle and m2Ferry vehicles having n flights, p, parked at different positions to be servediIndicating the fueling service time, q, for flight iiRepresents the service time of flight i, [ a ]i,bi]Time window representing flight i fueling service start time, [ c ]i,di]A time window representing the start time of the customer service on flight i, based on the time constraint relationship between the fueling service and the customer service, [ a ]i,bi]Should be earlier than [ c ]i,di]。
The time constraint relation between the refueling service and the customer service is specifically as follows:
when it is satisfied withAndwhen the flight number is larger than the preset value, each flight number is represented to have service only by one refueller and one ferry;
when s is satisfiedi∈[ai,bi]Namely ai≤si≤biAnd si+pi≤ti≤diRepresenting the time relation between each flight refuelling vehicle and the ferry vehicle; when simultaneously satisfying
si+pi+hij+M(xik+xjk-2)≤sj+M(1-μij)
ti+qi+h′ij+M(yil+yjl-2)≤tj+M(1-νij)
When, represents the constraint relationship between fueling service and service time of flight i and flight j (i < j), where M is a sufficiently large constant, if and only if xik=xjk=μijWhen 1, there is a constraint si+pi+hij≤sjSimilarly, if and only if yil=yjl=νijWhen 1, there is a constraint ti+qi+h′ij≤tjOtherwise this is not true;
when x is satisfiedik+xjk-1≤μij+μji≤1、yil+yjl-1≤νij+νjiAnd when the value is less than or equal to 1, representing the relation between decision variables.
Solving the model by using an NSGA-II algorithm, comprising the following specific steps:
step one, performing chromosome coding by using a real number code and a vehicle serial number, namely, a service vehicle corresponding to all flights is arranged on a chromosome, the length of the chromosome is K ═ M + D, M is the number of objective functions, D is the number of total flights, and the coding form is as follows:
s11s22s33…snk…sflightnumNj+l11l22…lnk…llightnumNb+D1…D3wherein s isnkThe tank service van for the nth flight is denoted as k cars, lnkThe ferry vehicle serving the nth flight is represented as a k vehicle, positions D +1 to k are objective function values, the objective function reads one chromosome each time, only a decision variable part is calculated, and the objective function values are returned, and the objective function values are stored at the tail ends of the chromosomes, so that the calculation and the data processing are convenient;
calculating the fitness value of each individual in the population according to the airport special vehicle scheduling model, and then performing non-inferior layering according to the fitness value of the individual to calculate the crowdedness of the individual at the same level;
step three, performing the following operations in each generation: selecting parents suitable for reproduction; performing cross mutation operations between selected parents; performing a selection operation between a parent and a descendant; replacing the non-adapted individuals with the adapted individuals to ensure that the population number is a constant;
step four, generating a parent population: according to the principle that the parent population and the offspring population are combined in the third step, the parent population and the offspring population are combined to obtain a new parent population;
judging whether a termination condition is met, wherein the termination condition of the algorithm is the maximum iteration algebra, and if the termination condition is met, terminating the iteration; otherwise, returning to the step two.
The chromosome coding in the step one comprises the following specific steps: the total number of the refuelers and the ferry vehicles is initially set, the vehicles are respectively numbered, the vehicle numbers from 0 to the total number of the special vehicles are randomly generated for each flight match until all N flights are matched, chromosomes of the refuelers and the ferry vehicles are spliced into a one-dimensional vector, and the objective function value is stored at the tail end of the chromosome to finally form a chromosome.
The fitness value of each individual is respectively as follows:
the total number of the refueling truck and the ferry vehicle is minimum
Shortest vehicle form path
The time difference between the vehicle arrival time and the earliest permitted service starting time is minimum
min(∑f1+∑f2)
The specific process of the cross mutation operation in the third step is as follows: adopting a single-point cross variation mode, matching every two individuals in a group into pairs randomly, wherein the exchange mode of each pair of individuals is as follows: randomly selecting two positions, exchanging elements between k1 and k2 in the two parents, calculating the fitness of the two new filial generations after exchange, and finishing the exchange to generate two filial generations if the fitness of the two filial generations is greater than that of the parents; if the fitness of the child is larger than that of the parent and the fitness of the child is smaller than that of the parent, the large child is reserved, and the small child is restored to be the parent; if the fitness of the offspring is smaller than that of the parent, the exchange is cancelled, namely the fitness of the two new offspring after the exchange is calculated, and the exchange mode is determined according to the value of the fitness smaller than that of the parent.
The specific process of selecting operation in the third step is as follows: and (3) adopting a binary tournament selection method, randomly selecting two individuals and comparing the fitness of the two individuals during selection, wherein the individual with the best fitness is selected as a parent until the number of the pairing pools is reached.
The invention has the beneficial effects that: according to the invention, through the combined scheduling of the refueling truck and the ferry vehicle, based on the time constraint relation of the refueling service and the boarding service, a combined scheduling model of flight refueling service and boarding service is provided, and the NSGA-II algorithm is used for solving the model, and experimental results show that the provided model can better solve the problem of cooperative scheduling of the refueling truck and the ferry vehicle.
(1) Compared with manual scheduling, the invention has the following advantages: the oil consumption of the special vehicle and the vehicle circulation time are greatly saved; the use cost of the vehicle is reduced.
(2) Compared with single scheduling, the combined scheduling of the refueling vehicle and the ferry vehicle has the advantages that the number of vehicles used is greatly reduced, the running distance is relatively reduced, and the difference of arrival time is small.
(3) The vehicle scheduling optimization method is not only limited to vehicle scheduling in airports, but also can be popularized to vehicle scheduling in other fields.
And comprehensively considering, the scheduling scheme designed by the joint scheduling of the refueling truck and the ferry vehicle is superior to manual scheduling and independent scheduling.
Drawings
FIG. 1 is a flow chart of the present invention for solving a model using the NSGA-II algorithm;
FIG. 2 is a diagrammatic illustration of a selected Tianjin airport stand layout of the present invention;
FIG. 3 is a value diagram of two objective functions corresponding to the Pareto optimal solution of the present invention;
fig. 4 is a schematic diagram of objective function values corresponding to the Pareto optimal solution of the present invention, and the abscissa and ordinate axes respectively correspond to three objective function values.
Detailed Description
To further explain the technical means and effects of the present invention adopted to achieve the predetermined object, the following detailed description of the embodiments, structures, features and effects according to the present invention will be given with reference to the accompanying drawings and preferred embodiments.
The technical scheme adopted by the invention for solving the technical problems is as follows: a scheduling optimization method for airport special vehicles is characterized in that a combined scheduling model of flight refueling service and passenger service is constructed based on a time constraint relation of refueling service and passenger service, and the model is solved by using an NSGA-II algorithm.
The combined scheduling model of flight refueling service and customer service is specifically described as follows: has m1Refuelling vehicle and m2Ferry vehicles having n flights, p, parked at different positions to be servediIndicating the fueling service time, q, for flight iiRepresents the service time of flight i, [ a ]i,bi]Time window representing flight i fueling service start time, [ c ]i,di]Indicating the start time of a customer service on flight iTime window, based on the time constraint relationship between the fueling service and the customer service, [ a ]i,bi]Should be earlier than [ c ]i,di]。
The time constraint relation between the refueling service and the customer service is specifically as follows:
when it is satisfied withAndwhen the flight number is larger than the preset value, each flight number is represented to have service only by one refueller and one ferry;
when s is satisfiedi∈[ai,bi]Namely ai≤si≤biAnd si+pi≤ti≤diRepresenting the time relation between each flight refuelling vehicle and the ferry vehicle; when simultaneously satisfying
si+pi+hij+M(xik+xjk-2)≤sj+M(1-μij)
ti+qi+h′ij+M(yil+yjl-2)≤tj+M(1-νij)
When, represents the constraint relationship between fueling service and service time of flight i and flight j (i < j), where M is a sufficiently large constant, if and only if xik=xjk=μijWhen 1, there is a constraint si+pi+hij≤sjSimilarly, if and only if yil=yjl=νijWhen 1, there is a constraint ti+qi+h′ij≤tjOtherwise this is not true;
when x is satisfiedik+xjk-1≤μij+μji≤1、yil+yjl-1≤νij+νjiAnd when the value is less than or equal to 1, representing the relation between decision variables.
Wherein m is1: total number of refuelers; m is2: the total number of ferry vehicles;n: the total number of flights to be serviced at different positions; a isi: flight i receives the earliest service starting time allowed by the service of the refuelling truck; bi: flight i receives the latest service starting time allowed by the service of the refuelling car; [ a ] Ai,bi]: a tank service time window; c. Ci: flight i receives the earliest permitted service starting time of the last customer service; di: flight i receives the latest starting service time allowed by the last customer service; [ c ] isi,di]: a ferry service time window; p is a radical ofi: the fuel oil filling time is 30 min; q. q.si: the service time of the ferry vehicle is 15 min; h isij: the time required for the refueller to reach the station of the flight j from the station of the flight i; dij: the distance from the refueller to the station of the flight j from the station of the flight i; h'ij: the time required for the ferry vehicle to reach the location of the flight j from the location of the flight i; d'ij: the distance from the place where the flight i is located to the place where the flight j is located by the ferry vehicle; d0i: the distance from the tank truck to the flight i stand; d'0i: the distance of the ferry car from the yard to flight i; x is the number ofik: whether a refueller k is assigned to flight i; y isil: whether ferry i is assigned to flight i; z is a radical ofk: whether the fuelling vehicle k is in use; z'l: whether the ferry vehicle is used or not; si: flight i's fueling service start time; t is ti: the service start time of the last passenger of flight i; s'i: the time when the refueller arrives at the stop of the flight i; t'i: and the time when the ferry vehicle arrives at the stand where the flight i is located.
As shown in fig. 1, the solution of the model by using the NSGA-II algorithm includes the following specific steps:
step one, performing chromosome coding by using a real number code and a vehicle serial number, namely, a service vehicle corresponding to all flights is arranged on a chromosome, the length of the chromosome is K ═ M + D, M is the number of objective functions, D is the number of total flights, and the coding form is as follows:
s11s22s33…snk...sflightnumNj+l11l22...lnk...llightnumNb+D1...D3wherein s isnkThe tank service van for the nth flight is denoted as k cars, lnkThe ferry vehicle serving the nth flight is represented as a k vehicle, positions D +1 to k are objective function values, the objective function reads one chromosome each time, only a decision variable part is calculated, and the objective function values are returned, and the objective function values are stored at the tail ends of the chromosomes, so that the calculation and the data processing are convenient;
calculating the fitness value of each individual in the population according to the airport special vehicle scheduling model, and then performing non-inferior layering according to the fitness value of the individual to calculate the crowdedness of the individual at the same level;
step three, performing the following operations in each generation: selecting parents suitable for reproduction; performing cross mutation operations between selected parents; performing a selection operation between a parent and a descendant; replacing the non-adapted individuals with the adapted individuals to ensure that the population number is a constant;
step four, generating a parent population: according to the principle that the parent population and the offspring population are combined in the third step, the parent population and the offspring population are combined to obtain a new parent population;
judging whether a termination condition is met, wherein the termination condition of the algorithm is the maximum iteration algebra, and if the termination condition is met, terminating the iteration; otherwise, returning to the step two.
The chromosome coding in the step one comprises the following specific steps: the total number of the refuelers and the ferry vehicles is initially set, the vehicles are respectively numbered, the vehicle numbers from 0 to the total number of the special vehicles are randomly generated for each flight match until all N flights are matched, chromosomes of the refuelers and the ferry vehicles are spliced into a one-dimensional vector, and the objective function value is stored at the tail end of the chromosome to finally form a chromosome.
Flight acceptance special car service number as represented by chromosome
The chromosome coded by the method can be regarded as spliced by two child gene strings, and the code has the advantages that all flights have the same probability of receiving the services of each refueller and ferry vehicle, the generated initial population is rich, and the value of the objective function is stored at the tail end of the chromosome, so that the calculation and the data processing are convenient.
The fitness value of each individual is respectively as follows:
The time difference between the vehicle arrival time and the earliest permitted service starting time is minimum
min(∑f1+∑f2)
The specific process of the cross mutation operation in the third step is as follows: adopting a single-point cross variation mode, matching every two individuals in a group into pairs randomly, wherein the exchange mode of each pair of individuals is as follows: randomly selecting two positions, exchanging elements between k1 and k2 in the two parents, calculating the fitness of the two new filial generations after exchange, and finishing the exchange to generate two filial generations if the fitness of the two filial generations is greater than that of the parents; if the fitness of the child is larger than that of the parent and the fitness of the child is smaller than that of the parent, the large child is reserved, and the small child is restored to be the parent; if the fitness of the offspring is smaller than that of the parent, the exchange is cancelled, namely the fitness of the two new offspring after the exchange is calculated, and the exchange mode is determined according to the value of the fitness smaller than that of the parent.
The specific process of selecting operation in the third step is as follows: and (3) adopting a binary tournament selection method, randomly selecting two individuals and comparing the fitness of the two individuals during selection, wherein the individual with the best fitness is selected as a parent until the number of the pairing pools is reached.
Further, the model is constructed based on the following assumptions: the number of airport special vehicles is limited; each flight receives the refueling service and the boarding service once respectively, and only needs one refueling vehicle and one ferry vehicle; the total service time of the refueling truck and the ferry vehicle is not limited; service is started once and is not interrupted until completion.
Experimental results and analysis:
1. experimental data and pretreatment thereof
The invention selects 126 outbound flights and departure flight data of the Tianjin airport in 1 day of 9 months to carry out an experiment on the algorithm. This section is tested refueling truck and ferry vehicle coordinated scheduling.
(1) Airport stand layout used in the experiments herein, which now has 63 passenger aircraft stands, as shown in fig. 2, is generally distributed in three areas: t1 terminal, T2 terminal, and remote gate. At civil aviation airports, special vehicles at airports must travel on a prescribed route, i.e. the connecting lines between the stations in fig. 1 cannot freely enter other areas to travel. All 63 stand positions are numbered in sequence according to the adjacency relation: 409. 410, 411, 412, 413, 414, 415, 416, 417, 418, 419, 101, 102, 103, 104, 105, 106, 107, 108, 109, 501, 502, 503, 504, 110, 111, 112, 113, 114, 115, 116, 117, 118, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, the distance between adjacent stands being about 40 meters. The parking lot of the special vehicle is positioned between the No. 501 parking stand and the No. 109 parking stand, the number is D, and the distance matrix (unit: meter) between the parking lot and partial parking stands is shown in the table 1:
TABLE 1
(2) Time window for station-passing flight and departure flight
According to the regulations of civil aviation bureaus: for fuel oil filling service, an airport refueling truck needs to finish the fuel oil filling service of an airplane five minutes before the start of boarding. The passenger refueling, or special case refueling exceptions, should also be completed at least five minutes before the expected departure time. Table 2 shows the allowable starting service time windows for some of the departure flights requiring fueling and receiving passenger service from the ferry.
TABLE 2
2. Theory of experimental results
The invention mainly utilizes MATLAB to carry out experiments. Through multiple experimental debugging, the main parameters of the genetic algorithm are set as follows: the population size is 100, the maximum iteration algebra is 800, the cross probability is 0.8, and the variation probability is 0.3.
Many leading edge solutions can be generated in the solving process, and the Pareto optimal solution can be finally obtained on the premise of meeting the algorithm termination condition. Fig. 2 and fig. 3 show objective function values corresponding to the Pareto optimal solution, where fig. 3 shows values of two objective functions corresponding to the Pareto optimal solution, fig. 4 shows objective function values corresponding to the Pareto optimal solution, and the horizontal and vertical axes respectively correspond to three objective function values.
The above experimental results for the airport refueller and ferry vehicle dispatch show that: compared with manual scheduling: the total distance of the special vehicle in the airport is reduced from 112.96km to 60.4km, the amplitude reduction reaches 46.5 percent, the special vehicle runs for less 52.56km, and the oil consumption of the special vehicle and the vehicle circulation time are greatly saved; the total number of the special vehicles is only 24 at least, so that the use cost of the vehicles is reduced. Compared with single scheduling, the combined scheduling of the refueling truck and the ferry vehicle has the advantages that the number of vehicles is greatly reduced, the number of the used combined scheduling vehicles is only 20 at least, the running distance is relatively reduced, the difference of the arrival time difference is small, comprehensive consideration is given, and the scheduling scheme designed by the combined scheduling of the refueling truck and the ferry vehicle is superior to manual scheduling and single scheduling.
The vehicle scheduling optimization method is not only limited to vehicle scheduling in airports, but also can be popularized to vehicle scheduling in other fields, and experiments prove that the scheduling scheme of the invention is superior to manual scheduling and independent scheduling, so that the vehicle use cost can be greatly reduced, and the time can be saved.
Claims (8)
1. The scheduling optimization method for the special type vehicles in the airport is characterized by comprising the following steps: and constructing a combined scheduling model of the flight refueling service and the customer service based on the time constraint relation of the refueling service and the customer service, and solving the model by using an NSGA-II algorithm.
2. The airport special vehicle scheduling optimization method according to claim 1, wherein: the combined scheduling model of flight refueling service and customer service is specifically described as follows: has m1Refuelling vehicle and m2Ferry vehicles having n flights, p, parked at different positions to be servediIndicating the fueling service time, q, for flight iiRepresents the service time of flight i, [ a ]i,bi]Time window representing flight i fueling service start time, [ c ]i,di]Time window indicating the start time of the customer service on flight i, [ a ]i,bi]Before [ c ]i,di]。
3. The airport special vehicle scheduling optimization method according to claim 1, wherein: the time constraint relation between the refueling service and the customer service is specifically as follows:
when it is satisfied withAndin time, it means that each flight has and only consists of one refueller and one ferryVehicle service;
when s is satisfiedi∈[ai,bi]Namely ai≤si≤biAnd si+pi≤ti≤diRepresenting the time constraint relation between each flight refuelling vehicle and the ferry vehicle; when simultaneously satisfying
si+pi+hij+M(xik+xjk-2)≤sj+M(1-μij)
ti+qi+h′ij+M(yil+yjl-2)≤tj+M(1-νij)
When, represents the constraint relationship between fueling service and service time of flight i and flight j (i < j), where M is a sufficiently large constant, if and only if xik=xjk=μijWhen 1, there is a constraint si+pi+hij≤sjSimilarly, if and only if yil=yjl=νijWhen 1, there is a constraint ti+qi+h′ij≤tjOtherwise this is not true;
when x is satisfiedik+xjk-1≤μij+μji≤1、yil+yjl-1≤νij+νjiWhen the decision variables are less than or equal to 1, representing the relation between the decision variables;
wherein m is1: total number of refuelers; m is2: the total number of ferry vehicles; n: the total number of flights to be serviced at different positions; a isi: flight i receives the earliest service starting time allowed by the service of the refuelling truck; bi: flight i receives the latest service starting time allowed by the service of the refuelling car; [ a ] Ai,bi]: a tank service time window; c. Ci: flight i receives the earliest permitted service starting time of the last customer service; di: flight i receives the latest starting service time allowed by the last customer service; [ c ] isi,di]: a ferry service time window; p is a radical ofi: the fuel oil filling time is 30 min; q. q.si: the service time of the ferry vehicle is 15 min; h isij: the refueller is from the position of the flight iThe time required to reach the flight j at the flight level; dij: the distance from the refueller to the station of the flight j from the station of the flight i; h'ij: the time required for the ferry vehicle to reach the location of the flight j from the location of the flight i; d'ij: the distance from the place where the flight i is located to the place where the flight j is located by the ferry vehicle; d0i: the distance from the tank truck to the flight i stand; d'0i: the distance of the ferry car from the yard to flight i; x is the number ofik: whether a refueller k is assigned to flight i; y isil: whether ferry i is assigned to flight i; z is a radical ofk: whether the fuelling vehicle k is in use; z'l: whether the ferry vehicle is used or not; si: flight i's fueling service start time; t is ti: the service start time of the last passenger of flight i; s'i: the time when the refueller arrives at the stop of the flight i; t'i: and the time when the ferry vehicle arrives at the stand where the flight i is located.
4. The airport special vehicle scheduling optimization method according to claim 1, wherein: solving the model by using an NSGA-II algorithm, comprising the following specific steps:
step one, performing chromosome coding by using a real number code and a vehicle serial number, namely, a service vehicle corresponding to all flights is arranged on a chromosome, the length of the chromosome is K ═ M + D, M is the number of objective functions, D is the number of total flights, and the coding form is as follows: s11s22s33...snk...sflightnumNj+l11l22...lnk...llightnumNb+D1...D3Wherein s isnkThe tank service van for the nth flight is denoted as k cars, lnkThe ferry vehicle serving the nth flight is represented as a k vehicle, positions D +1 to k are objective function values, the objective function reads one chromosome each time, only a decision variable part is calculated, and the objective function values are returned, and the objective function values are stored at the tail ends of the chromosomes, so that the calculation and the data processing are convenient;
calculating the fitness value of each individual in the population according to the airport special vehicle scheduling model, and then performing non-inferior layering according to the fitness value of the individual to calculate the crowdedness of the individual at the same level;
step three, performing the following operations in each generation: selecting parents suitable for reproduction; performing cross mutation operations between selected parents; performing a selection operation between a parent and a descendant; replacing the non-adapted individuals with the adapted individuals to ensure that the population number is a constant;
step four, generating a parent population: according to the principle that the parent population and the offspring population are combined in the third step, the parent population and the offspring population are combined to obtain a new parent population;
judging whether a termination condition is met, wherein the termination condition of the algorithm is the maximum iteration algebra, and if the termination condition is met, terminating the iteration; otherwise, returning to the step two.
5. The airport special vehicle scheduling optimization method according to claim 4, wherein: the chromosome coding in the step one comprises the following specific steps: the total number of the refuelers and the ferry vehicles is initially set, the vehicles are respectively numbered, the vehicle numbers from 0 to the total number of the special vehicles are randomly generated for each flight match until all N flights are matched, chromosomes of the refuelers and the ferry vehicles are spliced into a one-dimensional vector, and the objective function value is stored at the tail end of the chromosome to finally form a chromosome.
6. The airport special vehicle scheduling optimization method according to claim 4, wherein: the fitness value of each individual is respectively as follows:
The time difference between the vehicle arrival time and the earliest permitted service starting time is minimum
min(∑f1+∑f2)。
7. The airport special vehicle scheduling optimization method according to claim 4, wherein: the specific process of the cross mutation operation in the third step is as follows: adopting a single-point cross variation mode, matching every two individuals in a group into pairs randomly, wherein the exchange mode of each pair of individuals is as follows: randomly selecting two positions, exchanging elements between k1 and k2 in the two parents, calculating the fitness of the two new filial generations after exchange, and finishing the exchange to generate two filial generations if the fitness of the two filial generations is greater than that of the parents; if the fitness of the child is larger than that of the parent and the fitness of the child is smaller than that of the parent, the large child is reserved, and the small child is restored to be the parent; if the fitness of the offspring is smaller than that of the parent, the exchange is cancelled, namely the fitness of the two new offspring after the exchange is calculated, and the exchange mode is determined according to the value of the fitness smaller than that of the parent.
8. The airport special vehicle scheduling optimization method according to claim 4, wherein: the specific process of selecting operation in the third step is as follows: and (3) adopting a binary tournament selection method, randomly selecting two individuals and comparing the fitness of the two individuals during selection, wherein the individual with the best fitness is selected as a parent until the number of the pairing pools is reached.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010068912.8A CN111291888A (en) | 2020-01-21 | 2020-01-21 | Scheduling optimization method for airport special vehicles |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010068912.8A CN111291888A (en) | 2020-01-21 | 2020-01-21 | Scheduling optimization method for airport special vehicles |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111291888A true CN111291888A (en) | 2020-06-16 |
Family
ID=71023489
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010068912.8A Pending CN111291888A (en) | 2020-01-21 | 2020-01-21 | Scheduling optimization method for airport special vehicles |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111291888A (en) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112070355A (en) * | 2020-08-05 | 2020-12-11 | 北京交通大学 | Distribution scheduling method for airport ferry vehicle |
CN112163790A (en) * | 2020-10-29 | 2021-01-01 | 广东机场白云信息科技有限公司 | Airport ground resource scheduling method, electronic device and computer readable storage medium |
CN112734302A (en) * | 2021-02-04 | 2021-04-30 | 成都国翼电子技术有限公司 | Airport ferry vehicle scheduling method, device and equipment based on man-vehicle separation and storage medium |
CN112862260A (en) * | 2021-01-18 | 2021-05-28 | 华中科技大学 | Airport ferry service scheduling method and system |
CN113065692A (en) * | 2021-03-22 | 2021-07-02 | 南京航空航天大学 | Prediction method for service time of airport water truck |
CN113205216A (en) * | 2021-05-07 | 2021-08-03 | 哈尔滨工业大学 | Dynamic dispatching method and system for ferry vehicle at hub airport |
CN113344368A (en) * | 2021-05-31 | 2021-09-03 | 中国航空油料集团有限公司 | Refueling operation scheduling method and system |
CN114491317A (en) * | 2022-04-18 | 2022-05-13 | 中国民航大学 | Centralized deicing operation method and system for airplane, storage medium and computer equipment |
CN115130787A (en) * | 2022-08-29 | 2022-09-30 | 深圳市城市公共安全技术研究院有限公司 | Configuration method, system, terminal equipment and medium of emergency resource scheduling scheme |
CN115796554A (en) * | 2023-01-31 | 2023-03-14 | 中国航空油料有限责任公司 | Method and device for creating offline refueling order |
CN116485091A (en) * | 2023-02-14 | 2023-07-25 | 航安云创数据科技(上海)有限公司 | Airport vehicle scheduling method, airport vehicle scheduling device, electronic equipment and storage medium |
CN116777145A (en) * | 2023-05-29 | 2023-09-19 | 威海广泰空港设备股份有限公司 | Method for optimizing airport vehicle operating rate based on Internet of vehicles |
CN117391401A (en) * | 2023-12-07 | 2024-01-12 | 中国民用航空飞行学院 | Dispatching method of airport electric ground service vehicle |
CN118134207A (en) * | 2024-05-06 | 2024-06-04 | 青岛民航凯亚系统集成有限公司 | Production scheduling method and system in automatic distribution scene of civil airport guided vehicle |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016165392A1 (en) * | 2015-04-17 | 2016-10-20 | 华南理工大学 | Genetic algorithm-based cloud computing resource scheduling method |
CN107301510A (en) * | 2017-06-26 | 2017-10-27 | 北京首都国际机场股份有限公司 | A kind of tank service truck and ferry bus coordinated dispatching method based on genetic algorithm |
-
2020
- 2020-01-21 CN CN202010068912.8A patent/CN111291888A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016165392A1 (en) * | 2015-04-17 | 2016-10-20 | 华南理工大学 | Genetic algorithm-based cloud computing resource scheduling method |
CN107301510A (en) * | 2017-06-26 | 2017-10-27 | 北京首都国际机场股份有限公司 | A kind of tank service truck and ferry bus coordinated dispatching method based on genetic algorithm |
Non-Patent Citations (4)
Title |
---|
冯霞等: "基于遗传算法的加油车和摆渡车协同调度研究", 《交通运输系统工程与信息》, no. 02, pages 2 - 5 * |
李秦渝等: "随机配送时间车辆路径优化模型及求解算法", 《自动化与仪器仪表》, no. 03 * |
罗佳: "客户资源共享及需求可拆分的车辆路径问题", 广西民族大学学报, pages 2 * |
高伟;王俊义: "机场特种服务保障车辆优化调度研究", 《计算机仿真》, pages 3 - 4 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112070355A (en) * | 2020-08-05 | 2020-12-11 | 北京交通大学 | Distribution scheduling method for airport ferry vehicle |
CN112070355B (en) * | 2020-08-05 | 2024-04-05 | 北京交通大学 | Distribution scheduling method of airport ferry vehicle |
CN112163790A (en) * | 2020-10-29 | 2021-01-01 | 广东机场白云信息科技有限公司 | Airport ground resource scheduling method, electronic device and computer readable storage medium |
CN112862260A (en) * | 2021-01-18 | 2021-05-28 | 华中科技大学 | Airport ferry service scheduling method and system |
CN112862260B (en) * | 2021-01-18 | 2023-10-03 | 华中科技大学 | Airport ferry service scheduling method and system |
CN112734302A (en) * | 2021-02-04 | 2021-04-30 | 成都国翼电子技术有限公司 | Airport ferry vehicle scheduling method, device and equipment based on man-vehicle separation and storage medium |
CN112734302B (en) * | 2021-02-04 | 2023-08-29 | 成都国翼电子技术有限公司 | Airport ferry vehicle scheduling method, device and equipment based on man-vehicle separation and storage medium |
CN113065692A (en) * | 2021-03-22 | 2021-07-02 | 南京航空航天大学 | Prediction method for service time of airport water truck |
CN113205216A (en) * | 2021-05-07 | 2021-08-03 | 哈尔滨工业大学 | Dynamic dispatching method and system for ferry vehicle at hub airport |
CN113344368A (en) * | 2021-05-31 | 2021-09-03 | 中国航空油料集团有限公司 | Refueling operation scheduling method and system |
CN113344368B (en) * | 2021-05-31 | 2024-07-23 | 中国航空油料集团有限公司 | Oiling operation scheduling method and system |
CN114491317B (en) * | 2022-04-18 | 2022-06-21 | 中国民航大学 | Centralized deicing operation method and system for airplane, storage medium and computer equipment |
CN114491317A (en) * | 2022-04-18 | 2022-05-13 | 中国民航大学 | Centralized deicing operation method and system for airplane, storage medium and computer equipment |
CN115130787A (en) * | 2022-08-29 | 2022-09-30 | 深圳市城市公共安全技术研究院有限公司 | Configuration method, system, terminal equipment and medium of emergency resource scheduling scheme |
CN115796554A (en) * | 2023-01-31 | 2023-03-14 | 中国航空油料有限责任公司 | Method and device for creating offline refueling order |
CN116485091A (en) * | 2023-02-14 | 2023-07-25 | 航安云创数据科技(上海)有限公司 | Airport vehicle scheduling method, airport vehicle scheduling device, electronic equipment and storage medium |
CN116777145B (en) * | 2023-05-29 | 2024-02-20 | 威海广泰空港设备股份有限公司 | Method for optimizing airport vehicle operating rate based on Internet of vehicles |
CN116777145A (en) * | 2023-05-29 | 2023-09-19 | 威海广泰空港设备股份有限公司 | Method for optimizing airport vehicle operating rate based on Internet of vehicles |
CN117391401A (en) * | 2023-12-07 | 2024-01-12 | 中国民用航空飞行学院 | Dispatching method of airport electric ground service vehicle |
CN117391401B (en) * | 2023-12-07 | 2024-03-05 | 中国民用航空飞行学院 | Dispatching method of airport electric ground service vehicle |
CN118134207A (en) * | 2024-05-06 | 2024-06-04 | 青岛民航凯亚系统集成有限公司 | Production scheduling method and system in automatic distribution scene of civil airport guided vehicle |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111291888A (en) | Scheduling optimization method for airport special vehicles | |
Wang et al. | Solving multi-depot electric vehicle scheduling problem by column generation and genetic algorithm | |
CN109508751B (en) | Deep neural network model modeling method for high-speed railway train late time prediction | |
CN108805344A (en) | A kind of high-speed railway network train running scheme optimization method considering time-dependent demand | |
CN111582691A (en) | Double-layer planning-based transportation capacity matching method for multiple transportation modes of passenger transport hub | |
CN106530180A (en) | High-cold region charging service network planning method | |
CN113393073B (en) | Dispatching system and method for sharing transfer and connection between automobile and rail transit | |
CN105046379A (en) | City bus departure timetable optimization method | |
CN107301510A (en) | A kind of tank service truck and ferry bus coordinated dispatching method based on genetic algorithm | |
CN114925995A (en) | Road trunk vehicle and cargo matching and path planning method based on single parent genetic algorithm | |
CN110942193A (en) | Vehicle scheduling method and storage medium | |
CN114444965B (en) | Single-yard multi-line electric bus collaborative scheduling method | |
CN116070422A (en) | Flexible marshalling train turnover plan decomposition optimization method based on space-time network | |
CN109522584B (en) | Tramcar departure timetable design method | |
Alisoltani et al. | Optimal fleet management for real-time ride-sharing service considering network congestion | |
CN117057158A (en) | Electric automobile charging load simulation method considering cellular automata and space-time travel | |
Zhang et al. | Optimizing the strategic decisions for one-way station-based carsharing systems: A mean-CVaR approach | |
CN116739213A (en) | Subway connection bus optimization method based on agent model auxiliary algorithm | |
Hao et al. | Max-min ant system for bus transit multi-depot vehicle scheduling problem with route time constraints | |
Bi et al. | [Retracted] Multiobjective Optimization of Airport Ferry Vehicle Scheduling during Peak Hours Based on NSGA‐II | |
Zhaohua et al. | Research on airport special vehicle scheduling problem | |
CN115759567A (en) | Flight refueling scheduling scheme determining method and device | |
CN110627010B (en) | User classification guiding system and method for gas station | |
Zhao et al. | Research on optimal scheduling method of Airport Rapid Transit system | |
CN115034522B (en) | Dynamic dispatching method for commuting regular bus based on employee off-duty time and off-duty station |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |