Optimization of VRR for Cold Chain with Minimum Loss Based on Actual Traffic Conditions

Recently, fresh agricultural cold-chain logistics have been greatly developed with the increasing needs of people’s life. Reducing costs of cold-chain distribution has become the main object of loss control in logistics enterprises. The objective of this research is to find a set of optimal routes that minimize the total loss, including fuel cost, refrigeration cost, soft time window penalty cost, and cargo damage cost over transit time. In this paper, the definition and model construction of vehicle routing problem (VRP) with multiobjective minimum lost are introduced first. Then, an ant colony optimization (ACO) algorithm combined with Pareto local search (PLS) is put forward to solve the minimum loss model. In order to avoid the influence of complex road conditions during distribution, the distance matrix and the transit time matrix are both derived from the recommended navigation road based on E-map API. At last, a compare experiment between the traditional method and our proposed method is performed. The results indicate that our method has strong applicability and potential advantages in coldchain logistic and has important practical significance and application value.


Introduction
According to relevant media reports, rotting fruit, vegetables, and other foods transported by truck alone are worth about 70 billion yuan each year, causing a huge economic waste. Agricultural products especially for fruit, vegetables, fish, meat, etc. require strictly limited temperature, humidity, and time in the process of transportation and storage. As a branch of the logistics industry, cold-chain logistics provides a guarantee for the safe transportation of fresh agricultural products. In the process of transportation from cold storage to customers, a complete cold-chain logistics realizes the temperature control of the whole process of refrigerated and frozen food, as well as the closed environment, storage, and transportation during loading and unloading of goods. By selecting the optimal path in the process of cold-chain logistics transmission, the circulating rate of fruits and vegetables, meat, and aquatic products is reduced; the waste of resources and the cost of logistics costs are reduced too.
The vehicle routing problem (VRP) introduced by Dantzig and Ramser [1] plays a central role in the optimization of distribution networks.
Recently, the loss problem of cold-chain distribution has attracted the attention of many scholars and experts. For distribution problem of fresh agricultural products, fuel cost, refrigeration cost, cargo damage cost, and soft time window penalty cost have become the important components of loss control in distribution companies [2]. Li et al. [3] proposed a low-carbon model for fresh food and a genetic simulated annealing algorithm to solve the model in the cold-chain distribution. Yao et al. [4] proposed a minimizing fuel consumption solution to timedependent VRP with time window. Kim et al. [5] proposed a Markov decision process method to solve a dynamic VRP (DVRP) model with nonstationary transit times under actual traffic congestion. Chen et al. [6] developed a hybrid heuristic algorithm including harmony search and neighborhood descent to solve DVRP with time window. Abidi et al. [7] proposed a GA with a simple heuristic to solve a variant of RVPR (Rich VPR) with time widows and dynamically changing orders. Dongdong and Yinzhen [8] proposed a green multitype VRP with time windows to reduce the wastes of fuel consumption and carbon emission and use an improved tabu search algorithm to solve the G-MVRPTW. Fan et al. [9] pointed out that in fresh agricultural product cold-chain logistics, the total costs are composed of five kinds: fixed, transportation, damage, penalty, and energy consumption. Fang and Ai [10] proposed a mathematical model with soft time window penalty cost, refrigeration cost, cargo damage cost, and a hybrid ant colony algorithm to minimize the total costs. The DVRP mainly deal with the time-variation information of customer demands and road conditions. To share the traffic information and classify traffic conditions, Big Data and classification techniques are used in logistics and transportation. If the drivers receive the information of the traffic congestion or poor weather, they can change their way to reduce the costs and save time. Dimensionality reduction must be performed first in Big Data transmitted and stored process. Thippa et al. [11] proposed a machine learning (ML) algorithm with PCA to reduce the dimension of Big Data when the data sets are high. Gadekallu et al. [12] have studied the hybrid PCA-whale optimization algorithm to extract features and used a deep neural network to classify the diseases of tomato. Guha et al. [13] proposed an ANN-based content classification in combination with n -grame TF-IDF feature descriptor to classify the documents with accurate, sensitive information. The ant colony optimization (ACO) proposed by Colorni et al. [14] has been widely used in solving NP-hard vehicle routing problems. Many scholars used ACO algorithms to solve multiobjective combinatorial optimization problems [15][16][17][18][19]. These papers considered the costs of fuel, refrigeration, cargo damage, and delay simultaneously. However, their studies are not comprehensive enough in factors affecting the loss, such as the cost of energy consumption and rotting consumption when the door of compartment is opening.
We construct a total cost model, including soft time window penalty cost, fuel cost, refrigeration cost in transit and during unloading, and cargo damage cost in transit and during unloading. In our model, to take into account the actual traffic condition, the distance matrix and the transit time matrix are both archived from the navigation functions based on E-Map API. The traditional heuristic algorithm ACO has some limitations, such as easy stagnation in the initial stage and slow search speed. We perform an ACO algorithm with Pareto local search (PLS) on ants to obtain the uniform Pareto-optimal frontier and keep the diversity of Pareto solution set.
The paper is organized as follows. Section 1 introduces the definition and structure of VRP for cold-chain logistics. Section 2 constructs a minimum cost model based on actual traffic conditions. Then, an improved ACO algorithm with PLS is presented in Section 3. In Section 4, a compare experiment is performed to indicate our proposed model and the improved method. At last, we make a conclusion of our contributions in Section 5.

Problem Description
A cold-chain supplier has a warehouse with a certain number of transportation vehicles that delivers a variety of fresh products to a certain number of customers. The capacity of the transportation vehicles is limited, and the vehicles are all the same type with a total capacity at the starting point. The object of vehicle routing problem is to distribute goods to each customer with correct goods, limited time, minimum cost, and so on. VRP with time window means that a vehicle has to visit a customer within a certain time window. So, the service time, the transit time, and the total time must be calculated when dealing with the VRP. The mathematical model of VRPTW in cold-chain distribution is defined as follows [20].
There are K transportation trucks in this distribution center. The maximum load and the maximum transit time of per vehicle are Q and T. The kth distribution vehicle is responsible for path k. N is the total customers. The needs and the service time of customer point i are q i and u i ; t ij means the delivery time from customer points i to j. x k i is a 0-1 variable, and x k i = 1 means the point i in the path k. x k ij is a 0-1 variable, and x k ij = 1 means that the kth distribution vehicle travels from i to j.

Model
Minimizing the loss cost is the target of this paper. Therefore, the following paper introduces the structure of the loss cost of cold chain.   In our literature, the navigation function of E-map is used to avoid the mentioned shortcomings. The distance matrix and the transit time matrix under actual traffic conditions among customer points (including the distribution center) are derived from the E-map public platform [21].
As Figure 1 shows, there are two points with their latitudes and longitudes on the E-map.   Figure 2 shows, the distance and the transit time between the two points are shown by the navigation function of the E-map platform.
In our study, t ij represents transit time between customer i and customer j, and d ij is the distance between customer i and customer j. As abovementioned, t ij and d ij are both achieved from the navigation result based on E-map API.  11 is the penalty coefficient when vehicle k i arrives at customer i earlier than Et i , and C 12 is the penalty coefficient when vehicle k i leaves from customer i later than Lt i . The arrival time of vehicle k at customer i is T ki . According to the analysis above, the penalty cost of cold chain can be defined as the following:

Loss Cost Target
3.2.2. Fuel Cost. The fuel consumption of the distribution vehicle is inevitable for completing the distribution task. We assume that the fuel cost is proportional to the load of the vehicle pert unit distance. Let Q k be the load of the k vehicle and P k be the load rate of the kth vehicle.
Sets C 21 and C 22 represent no load and full load fuel consumption costs per unit distance of the distribution vehicle. C k 2 means the fuel consumption cost per unit distance of the kth vehicle.
The total cost of fuel consumption C 2 can be expressed as Initialize the coefficient variables and pheromone;  for iteration 1,…,M{  Initialize a taboo table and a start position for each ant. for ant, 1,…,N{ Select the next node according to rules The selected node is stored in the taboo table If all nodes are stored in taboo tables, the iteration is completed, break; The total cost is calculated The local pheromone is updated } The global pheromone is updated All ants' total costs are compared The current optimal solution is stored } The process is stopped and the current path with the shortest cost is output.  The temperature of the cold-chain logistic must be maintained at a certain low level to keep the freshness of goods. The refrigeration function of vehicles must be operated to achieve the required temperature, which causes the refrigeration cost immediately.

Mathematical Model of Minimum Loss Analysis.
According to the above analysis, the minimum loss model of the cold-chain distribution problem is as follows: Equation (12) is the objective optimization function which is aimed at finding the minimum total cost during the whole distribution process. Equation (13) Figure 5: The optimal path plan based on the traditional method.  (14) is a 0-1 variable too; the value 1 means the given route is travelled by the given vehicle. Equation (15) means each customer point is serviced by one and only one vehicle for distribution. Equation (16) means that the total demand served by the vehicle cannot exceed its maximum load. Equation (17) means that the vehicle can only start or end once.

Ant Colony Optimization Algorithm
4.1. Ant Colony Optimization Algorithm. The ant colony optimization (ACO) algorithm with parallelism, positive feedback, and strong robustness is used to solve the NPhard and highly constrained problem. In the 1990s, Italian scholars Dorigo and Maniezzo found that ants will release a certain amount of substance called pheromone in the path when they are searching for food. When they need to choose a path, they prefer to choose the path with high concentration of pheromone. At last, the route with the highest concentration will be selected as an optimist route between their home and the food source.
The ACO algorithm is defined as follows: where p k ij is probability of ant k transferring from point i to point j at time t; η ij ðtÞ is the heuristic function, and β is the heuristic function important factor; τ ij ðtÞ means the pheromone concentration on the path at time t, and α is the pheromone important factor. allowed k represents points that ant k is allowed to select in the next step; τ ij ðt + nÞ is the pheromone update over time; ρ is the pheromone factor; and Δτ k ij represents the pheromone enhancer. The rough process of the algorithm is as follows: 4.2. Combining ACO Algorithm with PLS. We applied Pareto local search (PLS) in our paper. After finding a path, PLS is used to further optimize the path. PLS proposed by Paquete et al. is a heuristic algorithm for tackling NP-hard multiobjective combinatorial optimization problems in the Pareto sense [22]. Pareto dominance defines a partial order on the set of feasible solutions. The goal when tackling the multiobjective problems in the Pareto sense is to find the set of Paretooptimal solutions. The weak component-wise ordering is used as a mutually nondominated criterion in PLS. Let a solution s ∈ A and its neighborhood s ′ ∈ NðsÞ. If each s ′ in NðsÞ is Table 5: Distribution paths of 4 vehicles based on the traditional method.

Vehicles
The route of each vehicle The realistic delivery on E-map (minutes) The      [23]. An improved ACO with PLS is proposed in our research to deal with the mode. In our approach, ants are used to optimize solutions and generate new Pareto solutions, and the new pheromone updating strategy is used to control the path selection. The procedure is described as follows: Step 1. Compute p k ij of n customers and update the trail level τ ij new for each k.
Step 2. Run PLS, and update the current solution according to calculating results.
Step 4. For each move (i, j) in BestSolution to update the trail level τ ij new.
Step 5. Repeat step 1 until the maximum computation time reaches.
The proportion of the pheromone update model of the ACO algorithm with PLS is as follows: where SðtÞ means information entropy value and p ij ðtÞ is transition probability. Equation (19) means that the proportion of the pheromone update will decreases with the increase of the pheromone value. k 1 is the positive constant, and its value range is [0,1]. The value 1 means the proportion of the pheromone update is high influenced by the information entropy value, and 0.75 is usually used as the value of k 1 . This definition combined a sequence of arithmetic with information entropy to adjust the pheromone adaptively.

Data Set Description.
In our experiments, the distribution situation of JingKeLong in Beijing City with one distribution center and 20 customer points is selected to verify our proposed mathematical model and algorithm. The distribution center has 4 vehicles of the same type with an average speed of 60 km/h, a maximum driving time of 1200 minutes, and a maximum load of 2000 kg. Table 1 shows the settings of parameters in our given minimum loss model and ACO+PLS algorithm. The data of the distribution center and customers are shown in Table 2. The distribution center is defined as 0, and customer points are defined as 1 to 20. Figure 3 shows the latitude and longitude of each point in the XY system, and Figure 4 shows the actual locations of them in E-map.
Our proposed method needs an actual transit time matrix and distance matrix both derived from E-map API. We supposed that the coordinates of points are stored into vertex [N] [2]; the core codes with E-map API to achieve the transit time matrix and distance matrix are described as follows. Table 3 shows the transit time matrix among customer points and distribution center obtained by E-map API. Table 4 shows the distance matrix among customer points and distribution center obtained by E-map API.

Experimental Result and Discussion.
In the traditional method, the distance matrix is achieved by the Euclidean distance formula. The Euclidean distance formula of two points on earth is defined as the following: where A and B are two points on the surface of earth and R is the radius of the earth. LatA and lonA represent the longitude and latitude of point A.
First, the traditional method ACO with the Euclidean distance matrix is performed to solve the VRP with minimum cost. Figures 5 and 6 are the optimal distribution routes and the trend of cost during iteration processed. The optimal distribution routing scheme, the total transit time, and the total cost based on the traditional method are shown in Table 5. Table 6: Distribution paths of 3 vehicles based on the proposed method.

Vehicles
The route of each vehicle The realistic delivery on E-map (minutes) The total time 593 The total cost 3673 9 Wireless Communications and Mobile Computing Second, our proposed method ACO+PLS with the actual transit time matrix and travel distance matrix is performed to solve the VRP with minimum cost. Figures 7 and 8 are the calculated distribution routes on map and the trend of cost during iteration processed. The optimal distribution routing scheme, the total transit time, and the total cost based on the proposed method are shown in Table 6.
From the above compare experimental results, we can see that ACO+PLS with actual traffic conditions can get smaller total cost and save more transit time. Our proposed method can be used to get the optimal solution effectively for VRP of cold-chain logistics.

Conclusion
This paper studies the model of VPR with minimum cost of cold-chain distribution. Minimizing the total cost of distribution can maximize the economic benefit of distribution enterprises. To consider the dynamic change of transit speed under actual traffic conditions, the transit time matrix and distance matrix are both derived from the navigation function based on E-map API. We proposed a heuristic approach ACO+PLS to solve the minimum loss model of cold-chain logistics. The experiments show that our proposed method has strong applicability and potential advantages in coldchain distribution. In future studies, a combination of multiple spatial information technologies such as geographic information technology and remote sensing technology can be realized to make the problem more practical.

Data Availability
The data used to support the findings of this study are included within the article.

Conflicts of Interest
The authors declare that they have no competing interest.