| Home | E-Submission | Sitemap | Contact Us |  
Environ Eng Res > Volume 28(1); 2023 > Article
Fan: Routing optimization method of waste transportation vehicle using biological evolutionary algorithm under the perspective of low carbon and environmental protection


Reasonably and effectively formulating the best route for urban waste transportation vehicle is particularly important for realizing low carbon and environmental protection of Green China construction concept. However, the current path planning method has shortcomings such as local optimization. In order to solve this problem, this paper aims at low carbon and environmental protection construction needs and proposes a routing optimization method of waste transportation vehicles based on improved ant colony algorithm. Firstly, the mathematical model of Vehicle Routing Problem (VRP) is constructed by considering transportation distance and carbon emissions cost. Then, network parameters in traditional ant colony algorithm are combined and optimized to realize self-adaptive update in the whole cycle. Furthermore, the neighborhood algorithm is used to iteratively optimize improved algorithm to avoid the defects of local optimization when solving VRP problem. Finally, the simulation results based on an actual dataset in North China show that the proposed method can achieve efficient and accurate optimal routing optimization for complex samples, and its solution stability index Pbest is 0.87 and the average deviation Var is 0.011, the lowest distribution cost after optimization is about 3,080 yuan, which are better than the comparison methods.

1. Introduction

With the rapid development of social economy and acceleration of urbanization, the amount of various types for waste has increased sharply [13]. According to incomplete statistics, the total amount of domestic waste in China’s cities and towns ranks first in the world. In recent years, with the acceleration of my country’s urbanization process and the improvement of people’s living standards, urban domestic waste is still increasing at an annual rate of 5 to 8% [4].
At present, on the one hand, the amount of waste is huge, on the other hand, there are phenomena such as backward transportation management methods and high equipment operation costs in the process of waste collection and transportation [5, 6]. This will lead to excessive waste of resources and low work efficiency problems in the current waste transportation links.
According to statistics, 80–85% of waste disposal system cost is attributed to the transportation system [7]. At the same time, the carbon dioxide emissions of my country’s transportation industry are showing an increasing trend, and transportation energy consumption accounts for about 9% of the total energy consumption of whole society [8]. With the transfer of waste treatment facilities to the outer suburbs and expansion of main urban area, the coverage of entire waste collection and transportation system is getting larger and larger. Therefore, scientific and reasonable planning of waste collection and transportation routes is an important guarantee for the realization for concept of “low energy consumption, low pollution and low emissions”, which is particularly urgent.
In recent years, the academic research on waste collection path is mainly based on the corresponding analysis of Traveling Salesman Problem (TSP) and the Chinese Postman Problem (CPP) [9, 10]. But in fact, we can divide waste collection and transportation routes into two aspects: waste collection routes and waste transfer routes. If focusing on the two aspects of modeling ideas and solving goals, the waste collection routing problem should belong to the category of Vehicle Routing Problem (VRP) [11, 12].
At present, some researchers have discussed this based on intelligent algorithms. Reference [13] established cost-oriented urban distribution optimization model, adopted an improved genetic algorithm and designed successive correction operators on both sides to improve the solution quality and efficiency of algorithm. Reference [14] established mixed integer programming model and designed a joint optimization genetic algorithm to solve the model. It used the proposed joint optimization genetic algorithm and two-stage genetic algorithm for simulation calculation and analysis. Reference [15] proposed a hybrid method based on sine cosine algorithm and paper swarm algorithm to solve the problem of heterogeneous vehicle routing. Reference [16] was based on the gray wolf spatial integer coding and solution of routing before grouping. It proposed an adaptive genetic gray wolf optimization algorithm to solve VRP with capacity constraints. However, it should be noted that the current routing optimization methods based on intelligent algorithms generally have the problems of weak adaptive ability and weak global optimization convergence ability [17]. This make the problem of searching for the best path in complex space often lose the optimal solution, and it will cause network to diverge in severe cases.
Aiming at the problems existing in current VRP planning method, this paper proposes a routing optimization method for waste transportation vehicles that meets the needs of low carbon environmental protection based on improved ant colony algorithm. The innovative points of proposed method are:
In order to further improve the optimization characteristics of the best route for waste vehicle transportation, this paper optimizes the network parameters based on traditional ant colony algorithm, which improves the global convergence and solution efficiency of the ant colony algorithm.
The pheromone setting method in the network is improved to realize the full-period self-adaptive update of network parameters and avoids losing the global optimal solution.
The remaining content of this paper is arranged as follows: Section 2 describes and models VRP problem accordingly. Section 3 describes the corresponding improvement content of ant colony algorithm. Section 4 is the simulation verification of proposed method based on actual sample dataset. Section 5 is the conclusion of this paper.

2. Problem Description and Model Construction

2.1. Problem Description

Waste management requires not only government investment and management, but also residents’ awareness of environmental protection. The government needs to increase waste collection facilities in each community under its jurisdiction and arrange waste vehicles to recycle waste in each community in time.
In an area, the process of waste collection can be described as follows: multiple waste vehicles depart from parking lot at the same time and begin to collect waste in various communities. Each community must be served and can only be served once (assuming that the amount of waste generated in each community is less than or equal to the maximum load capacity of waste vehicles). When the waste vehicle reaches or is close to the maximum load, it will first drive to transfer stations, unload the waste and then return to the depot. The goal of solving the problem is to find a shortest path that minimizes the total distance traveled by all waste vehicles. Fig. 1 is a simple schematic diagram of waste collection path.
As can be seen in Fig. 1, there are a total of 7 communities, and 2 waste vehicles can be used to complete the waste collection work in the community. Each waste vehicle departs from parking lot and starts to serve the community. After the waste collection volume reaches or approaches the maximum load capacity, waste vehicles go to the transfer station (treatment plant) to unload waste. Finally it returns to the parking lot.
In order to better describe the mathematical model of problem, the model elements are elaborated accordingly. The components of waste collection path problem include: parking lot, waste vehicle, waste, community, transfer station (processing plant), constraint conditions and objective function.

2.1.1. Depot

The parking lot is similar to logistics distribution center. It is the starting point and ending point of the waste vehicle’s driving process. According to the number of parking lots, it can be divided into single parking lot and multiple parking lot. However, considering the actual situation, usually a parking lot in an area can park required waste vehicles. So this paper studies the situation of bicycle yard.

2.1.2. Waste vehicle

Waste vehicle is a tool for transporting waste, which can be divided into single vehicle type and multiple vehicle types. The difference of model often determines the difference in maximum load. The maximum load capacity refers to the upper limit of waste vehicle loading capacity. The type of waste loaded is divided into solid and non-solid. The load condition refers to full load or partial load. The number of vehicles refers to the number of waste vehicles, which are divided into single vehicles and multiple vehicles according to their differences. The maximum mileage limit refers to whether to consider the restriction of waste vehicle mileage [18]. This paper mainly studies the problem of recycling path of community domestic waste. Most domestic waste is solid waste, so a waste vehicle loaded with solid waste is selected for research. Considering comprehensively, multiple waste vehicles of the same model (the same maximum load) are used for waste collection. The waste vehicles are fully loaded and mileage limit is not considered.

2.1.3. Waste

Waste is similar to goods in logistics transportation, with attributes such as weight and collection location. These two attributes are the basis for waste vehicles to choose the driving route.

2.1.4. Community

Each community is a service object corresponding to a waste collection location, and must be served once.

2.1.5. Transfer station (processing plant)

Taking into account the actual situation of community waste treatment by county, this element may be a transfer station or a treatment plant. According to their different numbers, it can be divided into single transfer station (treatment plant) and multiple transfer station (treatment plant). In order to reduce the complexity of solving the problem, this paper conducts research on the premise of a single transfer station (processing plant).

2.1.6. Constraint conditions

Waste vehicles must meet the following constraints when reclaiming waste: (1) Each waste vehicle must start from the yard; (2) The loading capacity is close to but not exceeding its own maximum load; (3) It must be transported before returning to the yard. The station (treatment plant) unloads the waste; (4) Every community must be served and only served once.
The time for waste collection by waste vehicles is arranged by the government every day, and it will not be shifted by the will of a certain community. There is no strong internal relationship between service and time, and the time window constraint has no practical significance in solving this problem. Thus, this paper does not consider the time window limitation.

2.1.7. Objective function

The commonly used objective functions for solving VRP problems are: (1) The total mileage of vehicles is the shortest; (2) The time-consuming is the least; (3) The cost is the least. This paper uses objective function (1) when solving the waste collection path problem, that is, the sum of all waste vehicles has the shortest mileage.

2.2. Model Construction

The routing optimization problem of waste transportation vehicles can be represented by an undirected graph G = (V, E) [19]. For ease of expression, the following symbols are introduced: point set V = {1, 2, ……, n} of parking lot, community, transfer station (processing plant), where point i = 1 represents parking lot, point i = n represents transfer station (processing plant), the remaining points represent communities, and the number of communities is n − 2. The element in edge set E is dij, which represents the distance between point i and point j. The point may be a parking lot, a community or a transfer station (processing plant). The waste discharge of each community is wi(i = 2, 3, ……, n − 1). There are r(k = 1, 2, ……, r)waste vehicles for waste collection in these n − 2 communities. Each vehicle has the same model, and the maximum load capacity is W and Wwi(i = 2, 3, ……, n − 1). Define the following variables:
The goal of solving the problem is to find an optimal path so that the total mileage Z of all waste vehicles is taken to the minimum.
The constraints that the model should meet are as follows:
In the above constraints, Eq. (3) indicates that waste vehicles must and can only return to the depot after arriving at transfer station (processing plant). Eq. (4) means that waste vehicles must arrive at the transfer station (processing plant) to unload the waste before returning to the parking lot. Eq. (5) means that every waste vehicle must depart from the parking lot and drive to transfer station (processing plant) after completing tasks. From the first two constraints, the driving route of each waste vehicle can be obtained as: parking lot-each community node-transfer station (processing plant)-parking lot. Eq. (6) and (7) indicate that every community can be served and can only be served once. Eq. (8) indicates that the total amount of waste collected by each waste vehicle in the community must not exceed its own maximum load.
However, a large amount of carbon emissions will be generated in the process of waste transportation, which will not only pollute the environment, but also increase the transportation cost of responsible department. Therefore, we also introduce carbon emissions into the routing optimization analysis.
The cost of carbon emissions mainly refers to the cost of carbon emissions generated in the transportation of waste vehicles. Carbon emissions are the product of fuel consumption and CO2 emission coefficient. The fuel consumption of vehicles is related to the load capacity and driving distance of vehicles. The relationship between the fuel consumption per unit distance of a transport vehicle and the load capacity can be expressed as:
Among them, ρ0 represents the fuel consumption per unit distance when the transport vehicle is empty; ρ* represents the fuel consumption per unit distance when the transport vehicle is fully loaded; W represents the full load capacity of vehicles.
Then, the target formula of routing optimization problem of waste transportation vehicles from the perspective of low carbon environmental protection is rewritten as:
where ξ represents the carbon tax; ɛ represents the carbon emission coefficient.
The community waste collection routing problem belongs to the category of VRP problem, which is similar to Capacitated VRP (CVRP) [20, 21]. Compared with the CVRP model, the biggest difference in the mathematical model of this paper is that each waste vehicle must drive to the transfer station (treatment plant) to unload waste before returning to the depot, and then return from the transfer station (treatment plant) to depot.

3. Routing Optimization of Waste Transportation Vehicles Based on Improved Ant Colony Algorithm

Ant colony algorithm has a certain memory function, which can continuously transfer from an old state to a new state when solving combinatorial optimization problems. But there are defects of local optimization and low search ability [22, 23]. For this reason, this paper improves traditional ant colony algorithm to realize the model solution of waste transportation vehicle routing optimization.

3.1. Mathematical Model of Ant Colony Algorithm

Before ant colony algorithm completes an iteration, the ants can only visit each node once, and return to the initial node after visiting all nodes. Ants secrete pheromone along the path they traveled. The amount and time of pheromone secretion are determined by algorithm.
Assuming that the number of communities to be visited is n, the size of ant colony is m, and the distance between nodes i and j is dij, then dij=(xi-xj)2+(yi-yj)2, where xi, xj, yi and yi represent the horizontal and vertical coordinates of nodes i and j respectively or latitude and longitude coordinates. At the beginning of each iteration, m ants will be randomly placed in these n nodes, making the starting point of ants random. In order to prevent the same ant from accessing the same node multiple times in one iteration, we set up a taboo table to record nodes that each ant has visited. tabuk records the names and order of the nodes that ant k have visited. τij(t) represents the pheromone concentration on the path <i, j> at time t. At the initial time (at time t = 0), the pheromone concentration on all paths is set to the same constant. The probability that ant k moves from current node i to node j at time t is represented by pijk(t). We call this state transition rule a random proportional rule, and its formula is as follows:
where the heuristic factor ㅋ ηij represents the visibility of path <i, j>, which reflects the heuristic degree of ants moving from node i to j, usually its value is:
Parameter α and parameter β respectively reflect the relative importance of pheromone concentration and heuristic information on the path when ants choose the path. The value of α affects the ant’s choice of repeated paths. The larger the value, the greater the ant’s ability to walk, which will make the algorithm prone to stagnation. The node numbers that ants have not yet visited are stored in allowedk, and the calculation formula is as follows.
When ant k selects an accessible node c according to the random ratio strategy and arrives, we need to update the taboo table, that is, add the number of node c to tabuk. If |tabuk| < n, then ant k continues to select the next node and visit according to the random ratio strategy. After repeating the above process, until tabuk contains all the communities. At this time, ant k returns to the starting point and adds starting community node to the end of tabuk. When m ants have completed their visit tasks, calculate the total distance traveled by each ant according to the visit order of each ant in the taboo table, and update pheromone on the path according to the following formula.
where δ represents the pheromone volatilization coefficient, which simulates the volatilization of pheromone secreted by ants in nature over time. The value range of δ is (0, 1). If the value is too large, the pheromone on the path will volatilize too fast, which will increase the possibility of re-selecting the previously searched path, thereby greatly reducing the understanding space and weakening the ability of algorithm to search globally. If the value is too small, although the possibility of searching for global optimal solution is increased, it also weakens the convergence of algorithm. Δτijk represents the number of pheromone left by ant k on the path <i, j> in this iteration.

3.2. Algorithm Design

3.2.1. Parameter selection optimization

In the realization of ant colony algorithm, there are many parameters that need to be initialized. Pheromone and heuristic function, information quantity-heuristic function product, and cooperative behavior between ants will seriously affect the convergence of algorithm [24]. At the same time, the parameters of ant colony algorithm are also key factors affecting the efficiency of its solution performance. Pheromone residual factor, information heuristic factor, expected heuristic factor, pheromone intensity, number of ants, etc. are all very important parameters. The selection method and selection principle directly affect the global convergence and solution efficiency of ant colony algorithm. The choice of different parameters can have a vital impact on the performance of ant colony algorithm. The influence of parameter α on the performance of ant colony algorithm

The information heuristic factor α reflects the relative importance of the amount of information accumulated by ants in the process of movement in guiding the search of ant colony. The larger the value, the more likely the ants will choose the path they have walked before, and the randomness of the search will be weakened. When the value of the information heuristic factor α is too small, the search of the ant colony is likely to fall into the local optimum prematurely.
It can be seen from Fig. 2 that when the value of α is 1.5–4.5, the comprehensive solution performance of ant colony algorithm is better. The influence of parameter β on the performance of ant colony algorithm

The expected heuristic factor β reflects the relative importance of heuristic information in guiding ant colony search process, and its size reflects the strength of a priori and deterministic factors in the ant colony optimization process. The larger the value is, the more likely the ant will choose the local shortest path at a certain local point. Although the convergence speed of algorithm is accelerated at this time, the randomness of ant colony’s search for the optimal path is weakened, and it is easy to fall into the local optimal solution.
It can be seen from Fig. 3 that when the parameter β is 2–6, the comprehensive solution performance of ant colony algorithm is better. The influence of parameter δ on the performance of ant colony algorithm

The size of the pheromone volatilization factor δ is directly related to the global search ability of ant colony algorithm and its convergence speed. The pheromone residual factor (1 − δ) reflects the strength of individual interaction between ants. Due to the existence of the parameter δ, when the scale of the problem to be dealt with is relatively large, the amount of information on the path that has never been searched will be reduced to close to 0, thus affecting the global search ability of algorithm.
It can be seen from Fig. 4 that under the same conditions of other initialization parameters, when the pheromone residual factor is between 2 and 4, the convergence performance and convergence speed of ant colony algorithm are better, and the calculation performance is also relatively stable. The influence of parameter τ on the performance of ant colony algorithm

The pheromone intensity τ is the total amount of pheromone released on the path that the ant circulates for a week. Its role is to make full use of the global information feedback on the directed graph, so that the algorithm can search for the global optimal solution of problem at a reasonable evolution speed under the action of positive feedback mechanism. When the Q value is large, although the convergence speed of ant colony algorithm is very fast, the global search ability of algorithm will deteriorate at this time, it is easy to fall into the local optimal solution, and calculation performance becomes very unstable.
It can be known from experience that the model used in this paper is ant week model, and the total amount of information has no obvious effect on the performance of ant colony algorithm of ant week model. Generally take τ = 100. When the size of city becomes larger, the value of τ should be relatively increased, otherwise the total amount of information remains unchanged. When the scale of the city is larger, the amount of information volatilized to each path becomes smaller, and it is not easy to find the global optimal solution.
The model used in this paper is ant week model, and the total amount of information has no obvious effect on the performance of ant colony algorithm of ant week model. Based on the experimental analysis of parameter selection law for ant colony algorithm, the specific steps of effective method for selecting the optimal combination parameters of ant colony algorithm are as follows:
  1. Determine the number of ants, that is, the proposed selection strategy of city size/number of ants≈2 to determine the total number of ants;

  2. Coarse parameter adjustment, that is, to adjust the information heuristic factor α, expected heuristic factor β, and pheromone intensity τ with a large value range to obtain a more ideal solution;

  3. Parameter fine-tuning, that is, adjusting the pheromone volatilization factor δ with a smaller value range.

3.2.2. Pheromone update strategy

When m ants have completed a cycle, that is, when all ants have constructed a feasible path R, compare it with the global optimal path R*, where R* = min(Rk | k = 1, 2, ……, m).
In order to make each ant’s iterative transportation scheme produce positive and negative feedback on the increase and decrease of pheromone, according to the rule of Fig. 5, the pheromone on all ants R in each iteration is updated globally:
  1. The pheromone update formula is:

    τijnew=θ2τijold         (i,j)R
  2. The pheromone update formula is:

    τijnew=θτiiold         (i,j)R
  3. The pheromone update formula is:

    Δτij*=QL*         (i,j)R
Q is a constant, which represents the total amount of pheromone released by the ant to complete a search; L* is the optimal path length explored by the ant. θ is the degree of pheromone retention. In the early stage of algorithm, the diversity of understanding can take a relatively large value. In the middle and late stages, for the rapid convergence of the optimal solution, the pheromone is more concentrated, and it can be smaller.
In order to avoid falling into a local optimal solution in the early stage, MMAS can be used for reference to limit the pheromone within a certain range to avoid losing the possibility of exploring new paths. which is

3.2.3. Optimize the optimal solution using neighborhood algorithm

In order to improve the convergence speed of algorithm, neighborhood algorithm is used to optimize the feasible solution produced by each iteration of ant colony algorithm [25]. There are three ways to find the neighborhood solution: insertion, reversal and exchange. In this paper, the exchange method is used to find the neighborhood solution, that is, the solution produced by each iteration of ant colony algorithm is first obtained. Then randomly select different delivery points to exchange the positions of the two in the path where the solution is obtained, perform random perturbation and recalculate to obtain a new solution. If the new solution is better than the original solution, the solution is updated, otherwise it remains unchanged. Repeat the above operations until the optimal solution is obtained when the maximum number of iterations is reached.

3.3. Algorithm Flow

  • Step 1: Define the initial value of each algorithm parameter. The pheromone value on the road section between any two points is initially τmax, and NC _ max is set as the maximum number of iterations of algorithm.

  • Step 2: Put m ants in the parking lot, and put the point in the current solution, using the method of constructing the solution in order.

  • Step 3: Calculate the transition probability, and require all constraints to be met at the same time, select another demand point in addition to the current solution set, and then place it in the current solution. If the next node that satisfies the constraint conditions is not found, return to the parking lot.

  • Step 4: Cycle Step 3 until m ants have visited all points, and get the corresponding loop that starts from the parking lot and meets the constraints. The several loops of each ant are equivalent to the distribution path formed by several transportation vehicles issued by vehicle site, and the shortest path R* is calculated and saved. Step 5: Perform global update and adjust θ according to the method described in this paper. The adjustment rules are:

    θi={0.95θt-1if         0.95θt-1θminθminotherwise
  • Step 6: Judge whether the number of iterations reaches the maximum value, if yes, terminate the transportation, otherwise clear the taboo table and go to Step 2.

4. Case Study and Verification

The verification and simulation of the proposed method are all carried out on a high-performance computer, so that the comparison method also runs in a unified simulation environment. The experimental hardware and software configurations are shown in Table 1. The experimental codes of the routing optimization method based on improved biological evolution algorithm are all realized by Matlab language.

4.1. Network Parameters and Application Cases

In the improved ant colony routing optimization proposed method, the chromosome coding uniformly adopts the natural number coding method. The parameter settings involved in ant colony algorithm: m = 30, α = 2.5, β = 3 and δ = 3; the number of iterations is 300 as the termination condition.
The experimental example uses an actual dataset in a certain area of North China, which governs 31 villages and uses a waste vehicle with a maximum load of 7,000 kg. Table 2 shows the serial number, name, coordinates and daily waste output of 1 depot, 31 villages, and a waste treatment plant.

4.2. Case Analysis

In order to verify the feasibility and engineering applicability of the proposed method, we choose reference [13], reference [14] and reference [16] as the comparison method, and our proposed method is analyzed separately under a unified experimental scenario.
Fig. S1 is a schematic diagram of the optimal path proposed in this paper based on improved ant colony algorithm. When the simulation experiment is about the 40th iteration, the best path of waste transportation vehicles is:
  • Waste vehicle No. 1: 1-19-20-21-25-24-22-23-28-29-33-1

  • Waste vehicle No. 2: 1-2-3-9-10-12-17-16-15-18-31-33-1

  • Waste vehicle No. 3: 1-4-5-6-7-8-11-14-13-33-1

  • Waste vehicle No. 4: 1-26-27-30-32-33-1

The best path of garbage transportation vehicles and the lowest cost is shown in Fig. S2. It can be seen from Fig. S2 that when the number of iterations reaches 40, the improved ant colony algorithm can converge quickly, the proposed model tends to be stable, and the minimum cost of distribution after optimization is about 3,080 yuan.
At the same time, we also analyze the cost of carbon emissions. The carbon emissions under different methods are shown in Fig. S3.
Fig. S3 shows the carbon emissions under different methods. It can be seen from Fig. S3 that the proposed method in this paper achieves effective convergence in about 40 iterations, and the cost is roughly controlled at 3,080 yuan. Compared with reference [13] and reference [14], carbon emission costs are reduced by 100 and 70 yuan, respectively, and the convergence speed of comparison algorithm is significantly slower than our proposed method. However, within the specified number of iterations-300 times, the problem is still effectively solved, and it is still in a divergent state in the reference [16]. The reason is that this paper improves the pheromone update strategy, which can realize the self-adaptive update of pheromone in the whole cycle of solution. This limit avoids the routing optimization problem being vulnerable to the risk of local optimization and realizes the optimization analysis of waste vehicle transportation path.
In addition, in order to more intuitively quantify the analysis performance of various methods, this paper presents the optimal path problem analysis capabilities under different methods in the form of Table S1. In Table S1, Gbest is the optimal value found, Gavg is the average value, Pbest is the probability of finding the optimal solution, and Var is the average deviation.
As shown in Table S1, when the waste vehicle transportation path is optimally solved, as the solution performance stability index, the value of our proposed method is 0.87, which is an improvement of 0.08 compared to the reference [14]. The average deviation value of the indicators analyzed in this paper is 0.011, which is 0.03 lower than that in the reference [14]. This shows that the proposed method can effectively analyze complex path analysis problems under the same experimental conditions. The comparison of indicators GB and G illustrates this phenomenon. The indicator values obtained by the proposed method are significantly lower than the comparison method.
In summary, this paper proposes a path analysis method for waste vehicle transportation based on improved biological evolutionary algorithm. This method is oriented to low carbon and environmentally-friendly transportation requirements, which can realize efficient and accurate optimal routing optimization analysis.

5. Conclusions

Aiming at the low efficiency of current VRP problem analysis, this paper considers the concept of low carbon environmental protection and proposes a new method for optimizing waste transportation vehicles based on improved ant colony algorithm. This method improves pheromone by combining and optimizing network parameters, and realizes the full-cycle adaptive update of algorithm. Moreover, it introduces neighborhood algorithm to iteratively optimize the biological evolution algorithm, which further improves the solution performance of VRP algorithm. Experimental results show that the proposed method can achieve optimal solution of complex paths, and it has better execution efficiency than other methods. The future research direction will make VRP problem more concrete and study the routing optimization method of waste transportation vehicles under fuzzy demands. In addition, due to the high timeliness requirements of waste transportation, the driving time in the distribution process will generally be affected by factors such as weather conditions, road conditions, emergency time and vehicle status. Therefore, the uncertainty of these factors can be more considered in the future research on waste transportation.

Supplementary Information



The author declares that he has no conflict of interest.

Author Contributions

The main idea of this paper is proposed by F.L.S. (engineer). All the work of the article was completed by F.L.S. himself.


1. Elgarej M, Mansouri K, Youssfi M, Nezha B, Fazazi EL. Distributed Swarm Optimization Modeling for Waste Collection Vehicle Routing Problem. Int J Adv Comput SC. 2017;8(9)306–312.

2. Vu HL, Bolingbroke D, Ng KTW. Assessment of waste characteristics and their impact on GIS vehicle collection route optimization using ANN waste forecasts. Waste Manage. 2019;88(1)118–130.

3. Rybnytska O, Burstein F, Rybin AV. Decision support for optimizing waste management. J Decis Syst. 2018;27(1)68–78.

4. Zhang SQ, Mu D, Wang C. A Solution for the Full-load Collection Vehicle Routing Problem with Multiple Trips and Demands: An Application in Beijing. IEEE Access. 2020;8(1)89381–89394.

5. Codu MK, Yilmaz M, Codur MY. Arc Routing Problem Approach for Reducing Exhaust Gas Emission in Road Transportation: A Case Study of Erzurum. Fresen Environ Bull. 2019;28(10)7196–7205.

6. Benrahou F, Tair A. Capacitated Vehicle Routing Problem for Collection Waste Lube Oil in Algiers. Fresen Environ Bull. 2019;28(6)4500–4505.

7. Banyai T, Tamas P, Illes B, Živilė S, Ágota B. Optimization of Municipal Waste Collection Routing: Impact of Industry 4.0 Technologies on Environmental Awareness and Sustainability. Int J Environ Res Pub HE. 2019;16(4)1–26.

8. Zhao HX, Liu GS, Li Y. Optimal path planning for classified waste recovery based on random walk. J Transp Eng Inf. 2018;16(3)103–108.

9. Vu HL, Ng K, Fallah B, Richter A, Kabir G. Interactions of residential waste composition and collection truck compartment design on GIS route optimization. Waste Manage. 2019;102(1)613–623.

10. Wohlk S, Laporte G. A districting-based heuristic for the coordinated capacitated arc routing problem. Comput Oper Res. 2019;111(1)271–284.

11. Yang J, Qiu K. Routing Optimization of White Pollution Garbage Cleaning and Transportation Vehicles using MMKCA Algorithm under the Perspective of Low Carbon and Environmental Protection. Fresen Environ Bull. 2021;30(4)4402–4410.

12. Amal L, Le HS, Chabchoub H, Lahiani H. Analysis of municipal solid waste collection using GIS and multi-criteria decision aid. Appl Geomat. 2019;12:193–208.

13. Ge XL, Zhang H. Research on vehicle routing optimization of urban real-time traffic network. Ind Eng Manage. 2018;23:140–149.

14. Li ZP, Zhang YW. Multi demand vehicle routing problem with time windows and service sequence constraints. Control Decis. 2019;34:1565–1570.

15. Bansal S, Wadhawan S. A Hybrid of Sine Cosine and Particle Swarm Optimization (HSPS) for Solving Heterogeneous Fixed Fleet Vehicle Routing Problem. Int J Appl Metaheur. 2021;12:41–65.

16. Huang GW, Cai YG. Adaptive genetic gray wolf optimization algorithm for vehicle routing problem with capacity constraint. Acta Electron Sinica. 2019;47:2602–2610.

17. Rossit DG, Nesmachnow S, Toutouh J. A bi-objective integer programming model for locating garbage accumulation points a case study. Rev Fac Ing. 2019;93:70–81.

18. Shen Y, Liu M, Yang J, Shi Y, Middendorf M. A Hybrid Swarm Intelligence Algorithm for Vehicle Routing Problem with Time Windows. IEEE Access. 2020;8:93882–93893.

19. Guo YN, Cheng J, Luo S, Gong DW. Robust Dynamic Multi-objective Vehicle Routing Optimization Method. IEEE ACM T Comput BI. 2018;15:1891–1903.

20. Vega-Mejía CA, Montoya-Torres JR, Islam S. A nonlinear optimization model for the balanced vehicle routing problem with loading constraints. Int T Oper Res. 2019;26:794–835.

21. Saeheaw T, Charoenchai N. A comparative study among different parallel hybrid artificial intelligent approaches to solve the capacitated vehicle routing problem. Int J Bio-Inspir Com. 2018;11:171–191.

22. Zhang S, Gajpal Y, Appadoo SS. A meta-heuristic for capacitated green vehicle routing problem. Ann Oper Res. 2018;269:753–771.

23. Wu Y, Pan F, Li S, Chen Z, Dong M. Peer-induced fairness capacitated vehicle routing scheduling using a hybrid optimization ACO–VNS algorithm. Soft Comput. 2019;24:2201–2213.

24. Huang YH, Blazquez CA, Huang SH, Paredes-Belmar G, Latorre-Nuñez G. Solving the Feeder Vehicle Routing Problem using ant colony optimization. Comput Ind Eng. 2019;127:520–535.

25. Wang JH, Weng TY, Zhang QF. A Two-Stage Multiobjective Evolutionary Algorithm for Multiobjective Multidepot Vehicle Routing Problem With Time Windows. IEEE T Cybernetics. 2019;49:2467–2478.

Fig. 1
Schematic diagram of waste vehicle path.
Fig. 2
Sensitivity analysis of parameter α.
Fig. 3
Sensitivity analysis of parameter β.
Fig. 4
Sensitivity analysis of parameter δ.
Fig. 5
Global pheromone update process.
Table 1
Experimental Simulation Scenarios
Software environment Operating system Windows10
Program editor Matlab2016a
Hardware environment CPU 3.2GHz Intel Core i5-6500 CPU
GPU GeForce GTX 1050Max-Q
Running memory 8 GB
Table 2
Actual Sample Dataset of a Certain Area
No. Name Waste output (kg)
1 Parking lot 0
2 Linglong goushang Village 828
3 Shanqian Village 517
4 Hengzhang Jiangjia Village 925
5 Lugezhuang Village 967
6 Mu Lin Zhuang Village 840
7 Shengjia Village 879
8 Wangjia Village 872
9 Qianhuayuan Village 696
10 Houhuayuan Village 585
11 Panjiaji Village 828
12 Oujiakuang Village 853
13 Lugezhuang Village 515
14 Lujia village 638
15 Linglong Taishang Village 523
16 Dajiangjia Village 912
17 Xiaojiangjia Village 548
18 Linglong Gaojia Village 848
19 Fengjia Village 658
20 Liujia Village 776
21 Zhaizi Village 517
22 Dongzhuangtou Village 719
23 Xizhuangtou Village 669
24 Dongtuanbu Village 883
25 Xituanbu Village 792
26 Huwangzhuang Village 593
27 Longquanzhuang Village 769
28 Gaojiatuan Village 824
29 Guanjiahe Village 775
30 Kuang Lujia village 591
31 Lijia village 628
32 Yuquanzhuang Village 857
33 Transfer station 0
PDF Links  PDF Links
PubReader  PubReader
Full text via DOI  Full text via DOI
Download Citation  Download Citation
Supplement  Supplement
Web of Science
Editorial Office
464 Cheongpa-ro, #726, Jung-gu, Seoul 04510, Republic of Korea
FAX : +82-2-383-9654   E-mail : eer@kosenv.or.kr

Copyright© Korean Society of Environmental Engineers.        Developed in M2PI
About |  Browse Articles |  Current Issue |  For Authors and Reviewers