Research Article  Open Access
A Hybrid Tabu Search Algorithm for a RealWorld Open Vehicle Routing Problem Involving Fuel Consumption Constraints
Abstract
Outsourcing logistics operation to thirdparty logistics has attracted more attention in the past several years. However, very few papers analyzed fuel consumption model in the context of outsourcing logistics. This problem involves more complexity than traditional open vehicle routing problem (OVRP), because the calculation of fuel emissions depends on many factors, such as the speed of vehicles, the road angle, the total load, the engine friction, and the engine displacement. Our paper proposed a green open vehicle routing problem (GOVRP) model with fuel consumption constraints for outsourcing logistics operations. Moreover, a hybrid tabu search algorithm was presented to deal with this problem. Experiments were conducted on instances based on realistic road data of Beijing, China, considering that outsourcing logistics plays an increasingly important role in China’s freight transportation. Open routes were compared with closed routes through statistical analysis of the cost components. Compared with closed routes, open routes reduce the total cost by 18.5% with the fuel emissions cost down by nearly 29.1% and the diver cost down by 13.8%. The effect of different vehicle types was also studied. Over all the 60 and 120node instances, the mean total cost by using the lightduty vehicles is the lowest.
1. Introduction
Many developing countries are confronted with two problems: how to reduce economic costs and how to develop an environmentfriendly society. As a result, researchers as well as companies with integrated logistics expertise tend to adopt new transportation models to counter these problems. Outsourcing logistics operations to the thirdparty logistics would reduce costs through better resource utilization and operations efficiency in freight transportation. It plays an increasingly important role in freight transportation. In that case, a company can hire vehicles from other companies to deliver its goods. Vehicles do not need return to the company as usual. They are usually described as the open vehicle routing problem (OVRP) [1]. Compared with the vehicle routing problem (VRP), routes in the OVRP model are open; see Figure 1.
The OVRP is NPhard [2], so heuristics or metaheuristics methods are usually used to deal with it [3], such as the tabu search [4], the neighborhoodbased search [5], particle swarm optimization [6], ant colony optimization [7], or evolutionary computing [8]. Hybrid metaheuristic algorithms were also designed to solve the OVRP [9]. Moreover, several variants of the OVRP were studied to model specific practical problems. OVRP with time windows (OVRPTW) was presented by [10] to model the delivery of multiproduct newspapers. The multidepot OVRP (MDOVRP) was first proposed by [11] to model the distribution of fresh meat. The open vehicle routing problem with decoupling points (OVRPDP) was introduced in [12] to describe open routes performed by more than one carrier. The OVRP with uncertain demands was investigated to deal with nondeterministic customer demands and avoid unsatisfied demands or more extra operation cost [13]. However, among models studied in the past there are rare formulations considering the impact of fuel and carbon emission [14]. This paper is a contribution to this line of research.
Carbon emissions along with freight transportation have hazardous impacts on the environment. How to minimize the fuel consumption becomes a hot topic. In [15], six fuel consumption models were compared with each other and several factors affecting carbon emissions in road transportation were analyzed. In [16], the pollutionrouting problem (PRP) was introduced to evaluate the greenhouse gas emissions for the VRP. The total energy consumed on each road can be directly translated into greenhouse gas emissions with the PRP model. Reference [17] considered a heterogeneous vehicle fleet and extended the PRP model. Reference [18] presented a recent review of studies on green freight transportation problems.
This study has been motivated by a real problem in Beijing, China. As outsourcing logistics is playing a more and more important role in China’s freight transportation, many companies hire vehicles from the thirdparty logistics. Vehicles do not need return to the depot after delivering products to customers. Each route departs from the depot and ends up with one of the customers. These companies are facing the problem of minimizing the total cost, including both fuel consumption cost and driver salaries.
The main contributions of our paper are summarized as follows. (1) Fuel consumption was analyzed in the context of the OVRP; and the green open vehicle routing problem (GOVRP) model was constructed. With the purpose of optimizing the fuel emissions cost in the model, we introduced the comprehensive modal emission model (CMEM) into the OVRP model and extend the OVRP with an objective that accounts for the amount of greenhouse emissions. (2) A hybrid tabu search algorithm involving several neighborhood search strategies was proposed to solve the GOVRP. A modified nearest neighborhood heuristic (mNNH) was proposed to get the initial solution. It considered two factors load and distance when searching for a suitable initial solution. And then four neighborhood search operators were designed to produce the neighborhood of the current solution. (3) Computational experiments were conducted on instances derived from real geographical data of customers in Beijing, China. Both the effect of open routes and the effect of vehicle types were analyzed in respects of reducing total cost.
2. Mathematical Model of the GOVRP
The CMEM was developed by [19, 20] to calculate fuel emissions. In this section, we introduced the CMEM to the OVRP and formulated the mathematical model of the GOVRP. The fuel consumption (in liters) of vehicle type is given by [21] where , , , and are constraints. is the total vehicle weight, is the distance, and is the vehicle speed. Three terms of are referred to as the engine module, the weight module, and the speed module, respectively. Notations and their default values are listed in the following two tables [17, 22]. Common parameters of vehicles are listed in Table 1. Specific parameters for different vehicle types are listed in Table 2.


Let be the set of customers and the depot, and let be the set of arcs. Node 0 denotes the depot, so the customer set is . represents the capacity of vehicles. Variables means the demand of customer . Variables is the distance from node to node . Variables are the total amount of flow on arc . The mathematical model of GOVRP is defined as follows:
Formula (2)–(5) describes the objective of the GOVRP, where terms (2)–(4) calculate the total fuel consumption cost and term (5) measures the total driver wages. More specifically, is the fuel emissions cost per liter, and terms (2)–(4) calculate the fuel consumption costs induced by the engine module, the weight module, and the speed module, respectively. denotes the total time spent on a route where is the last served customer.
The binary variable is equal to 1 if there is a vehicle traveling on arc ; otherwise, is equal to 0. Constraints (6)–(9) ensure that each customer is served by one vehicle and it is served only once. Vehicles do not return to the deport. If is the last served customer, constraint (8) can be written as . Otherwise, it can be written as . Constraints (10) and (11) define flows. The binary variable is equal to 1 when a vehicle travels on arc at speed .
3. Hybrid Tabu Search Algorithm for the GOVRP
In this study, a hybrid tabu search algorithm including several neighborhood search strategies was designed to solve the GOVRP. The flowchart of this algorithm is shown in Figure 2. The detailed steps are listed as follows.
Step 1. The initial setup includes initializing an empty tabu list and an empty candidate list, creating an initial solution, setting that initial solution as the best solution to date. A modified nearest neighborhood heuristic (mNNH) is used to gain the initial solution. We will elaborate on it specifically in Section 3.1.
Step 2. If the stopping condition is satisfied, the search process stops, and the best solution is returned. Otherwise, it turns to Step 3.
Step 3. The solution space is explored each step from a current solution to the accepted solutions in the neighborhood by using four efficient neighborhood search algorithms. For each neighbor of the current solution , the optimal speed on each arc of their routes is set by using the speed optimization algorithm according to [23].
Step 4. The neighboring solutions are checked for tabu elements. We search for the best solution and update the candidate list.
Step 5. If the best local candidate is better than the current best solution, update the current best solution. The local optimal solution is added to the tabu list. Then, turn to Step 2.
In this paper, the stopping condition is the maximum number of iterations . When our algorithm is iterated times, the algorithm is terminated.
3.1. The Initial Solution
We use a modified nearest neighborhood heuristic (mNNH) to obtain an initial solution. The mNNH considers two factors load and distance when searching for a suitable initial solution. In a delivery system, a vehicle can reduce load after servicing customer . We define where is defined as the unrouted customer set.
The mNNH builds routes one by one according to the following steps. The first route begins with an unrouted customer , which can be calculated by the following equation:Then, we calculate between the current node and the other unrouted nodes. We choose the node with the greatest current value as the next node if the current route does not violate the vehicle capacity constraint. We update the current node and search for the next node in the same way. When no customer can be assigned to the route, a new route is started. When all the customers are already routed, the process stops. The algorithm is described as Algorithm 1.

3.2. The Neighborhood Search
The neighborhood of the current solution is first obtained by four neighborhood search operators. They cannot violate the capacity constraint during this process as before. Then, the optimal speed on each route is set according to [23].
(1) Random Operator (RO). It randomly selects one node from the solution and then randomly finds a possible position for it.
(2) HighCostNode Improvement Operator (HCNIO). As illustrated in Figure 3, the operator tries to reassign the high cost node , because the total distance between its preceding customer and its following customer is the longest. Customer can be calculated as follows: where is the preceding customer and is the following customer.
We compute the best possible position for node based on where is the cost for inserting node between customer and customer and is the set of currently routed customers.
(3) ShortRoute Improvement Operator (SRIO). The operator tries to combine or delete short routes by reassigning nodes belonging to those routes. A threshold is set to evaluate the length of a route. Only if the length of a route is shorter than the threshold , the nodes in that route can be reassigned to other positions. The process can be illustrated in Figure 4 (). First, arcs in the shorter routes are deleted and then isolated nodes are assigned to other longer routes.
The best possible position for each unrouted node is also calculated by (19). By using (20), the most suitable customer is selected to be inserted in the route repeatedly. Unrouted nodes will be inserted to the current routes one by one.
(4) RandomRouteImprovement Operator (RRIO). Each route in the current solution has a chance of being destroyed. In that case, customers in that route should be reassigned. The best possible position for each node is calculated by (19). Then best customer to be inserted in the route is selected by using (20) repeatedly. All other unrouted customers will be inserted into the sequence one by one. As shown in Figure 5, two routes included in the dotted oval are destroyed; nodes in those routes are reassigned in possible positions.
For each neighbor of the current solution , the optimal speed on each route is set. According to [23], the optimal speed is
4. Computational Analysis
Experiments are run with data derived from real geographical distances of customers in Beijing. Three smaller classes with 10, 20, and 30 customers and four larger classes with 60, 80, 100, and 120 customers are generated. Each class includes 10 instances. All the customers involved in this study are shown in Figure 6. The algorithmic parameters values are given in Table 3.

4.1. Effect of Objectives Min TD, Min FEC and Min DC
Three different objectives, Min TD (total distance), Min FEC (fuel emissions cost), and Min DC (driver cost), are used to minimize the total distance, fuel emissions cost, and driver cost, respectively. We obtained solutions with different cost components on each performance measure and analyzed the effect of different objectives. Experiments were conducted on 10, 20, and 30node instances by using lightduty vehicles. The mean result of each instance collected over ten runs is reported; see Tables 4, 5, and 6. The columns display total distance (TD), execution time (ET), fuel emissions cost (FEC), driver cost (DC), total cost (TC), and CO_{2} emissions (CE). Table 7 presents the average results of each class.




On average of 10node class, the fuel emissions cost represents about 30% of the total cost. For 20 and 30node classes, the fuel emissions cost accounts for approximately 26% of the total cost. The carbon dioxide emissions increase with the number of customers. Only considering the driver cost in the objective leads to the most carbon dioxide emissions and the poorest total cost performance. Only considering the fuel emissions gets solutions with the lowest carbon dioxide emissions, but a higher total cost. Min TC as an objective yields the lowest total cost and the shortest total distance.
4.2. Effect of Open Routes
In this subsection, we compared the total cost of open routes with that of closed routes. Experiments were also conducted on smaller instance classes by using lightduty vehicles. The mean results of closed routes collected over ten runs are listed in Table 8. The average results of each class are reported in Table 9 to illustrate the effect of open routes. As shown in Table 9, open routes reduce the total cost by 18.5% with the fuel emissions cost down by nearly 29.1% and the diver cost down by 13.8% over all the instances.


4.3. Results for Larger Scale Instances
Experiments were also conducted on the four largest sets with the objective Min TC. The results collected over ten runs are reported in Table 10 by using lightduty vehicles. Notations BS, MS, WS, SD, and ET in columns represent the best solution, mean solution, the worst solution, standard deviation, and execution time, respectively. As shown in Table 10, the total cost increases with the number of customers.

4.4. Results with Different Vehicle Types
Experiments were conducted on the 60 and 120node instances by using different types of vehicles. In this subsection, the analysis focuses on which type of vehicle is the most suitable to minimize the total cost. The specific vehicle parameters are listed in Table 2. Results with different vehicle types are listed in Tables 11 and 12. Over all the instances in the two classes, the mean total cost by using the lightduty vehicles is the lowest. For 120node instances, the mean total cost by using mediumduty vehicles increases by 4.8% and the result by using heavyduty vehicles rises 23.6%. For 60node instances, the mean total cost by using mediumduty vehicles increases by 7.7% and the value by using heavyduty vehicles rises 28%.


5. Conclusions
In this study, the GOVRP was introduced as a formulation considering the impact of fuel emissions for the thirdparty logistics. The problem was to construct open routes for vehicles to visit all customers with the vehicle capacity constraints. The objective was to minimize the total cost composing of the fuel emissions cost and the driver cost.
A hybrid tabu search algorithm was designed to deal with the GOVRP instances. Experiments were conducted on 60 instances derived from real geographical data of customers in Beijing. We analyzed the cost components with different objectives and compared open routes with closed routes in respects of reducing the total cost. A homogenous fleet was considered in this paper. Computational results showed that vehicle type influenced the total cost and it should be changed according to the instance’s size. In the future, a heterogeneous fleet of vehicles can be used to minimize the total cost. Other types of the OVRP, such as OVRPTW and MDOVRP, can be solved by our algorithm with some changes.
Nowadays, several novel computing techniques are used to deal with complex problems. Some of them are bioinspired models, such as membraneinspired evolutionary algorithms [24–26] and probe machine [27]. Their nondeterministic distributed parallel frameworks have been proved to improve the performance of optimization algorithms [28]. Most of them can be used to solve reallife problems [29]. We hope that more competitive results for our GOVRP instances can be obtained by using those algorithms in the future.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grant nos. 61502012, 61772290, and 61373066); the Beijing Natural Science Foundation (Grant no. 4164096); the Humanity and Social Science Youth Foundation of Ministry of Education of China (Grant no. 13YJC630010); the Science and Technology Development Strategy Research Program of Tianjin (Grant no. 16ZLZXZF00030); and the Collaborative Innovation Center for China Economy.
References
 L. Schrage, “Formulation and structure of more complex/realistic routing and scheduling problems,” Networks, vol. 11, no. 2, pp. 229–232, 1981. View at: Publisher Site  Google Scholar
 S. Raff, “Routing and scheduling of vehicles and crews: the state of the art,” Computers and Operations Research, vol. 10, no. 2, pp. 63–211, 1983. View at: Publisher Site  Google Scholar
 J. Brandao, “A tabu search algorithm for the open vehicle routing problem,” European Journal of Operational Research, vol. 157, no. 3, pp. 552–564, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 U. Derigs and K. Reuter, “A simple and efficient tabu search heuristic for solving the open vehicle routing problem,” Journal of the Operational Research Society, vol. 60, no. 12, pp. 1658–1669, 2009. View at: Publisher Site  Google Scholar
 E. E. Zachariadis and C. T. Kiranoudis, “An open vehicle routing problem metaheuristic for examining wide solution neighborhoods,” Computers and Operations Research, vol. 37, no. 4, pp. 712–723, 2010. View at: Publisher Site  Google Scholar
 S. A. Mirhassani and N. Abolghasemi, “A particle swarm optimization algorithm for open vehicle routing problem,” Expert Systems with Applications, vol. 38, no. 9, pp. 11547–11551, 2011. View at: Publisher Site  Google Scholar
 X.Y. Li, P. Tian, and S. C. H. Leung, “An ant colony optimization metaheuristic hybridized with tabu search for open vehicle routing problems,” Journal of the Operational Research Society, vol. 60, no. 7, pp. 1012–1025, 2009. View at: Publisher Site  Google Scholar
 S. Yu, C. Ding, and K. Zhu, “A hybrid GATS algorithm for open vehicle routing optimization of coal mines material,” Expert Systems with Applications, vol. 38, no. 8, pp. 10568–10573, 2011. View at: Publisher Site  Google Scholar
 P. P. Repoussis, C. D. Tarantilis, O. Braysy, and G. Ioannou, “A hybrid evolution strategy for the open vehicle routing problem,” Computers and Operations Research, vol. 37, no. 3, pp. 443–455, 2010. View at: Publisher Site  Google Scholar
 R. Russell, W.C. Chiang, and D. Zepeda, “Integrating multiproduct production and distribution in newspaper logistics,” Computers & Operations Research, vol. 35, no. 5, pp. 1576–1588, 2008. View at: Publisher Site  Google Scholar
 C. D. Tarantilis and C. T. Kiranoudis, “Distribution of fresh meat,” Journal of Food Engineering, vol. 51, no. 1, pp. 85–91, 2002. View at: Publisher Site  Google Scholar
 R. Atefi, M. Salari, L. C. Coelho, and J. Renaud, “The open vehicle routing problem with decoupling points,” European Journal of Operational Research, vol. 265, no. 1, pp. 316–327, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 C. Erbao, L. Mingyong, and Y. Hongming, “Open vehicle routing problem with demand uncertainty and its robust strategies,” Expert Systems with Applications, vol. 41, no. 7, pp. 3569–3575, 2014. View at: Publisher Site  Google Scholar
 K. Fleszar, I. H. Osman, and K. S. Hindi, “A variable neighbourhood search algorithm for the open vehicle routing problem,” European Journal of Operational Research, vol. 195, no. 3, pp. 803–809, 2009. View at: Publisher Site  Google Scholar
 E. Demir, T. Bektaş, and G. Laporte, “A comparative analysis of several vehicle emission models for road freight transportation,” Transportation Research Part D: Transport and Environment, vol. 16, no. 5, pp. 347–357, 2011. View at: Publisher Site  Google Scholar
 T. Bektaş and G. Laporte, “The pollutionrouting problem,” Transportation Research Part B: Methodological, vol. 45, no. 8, pp. 1232–1250, 2011. View at: Publisher Site  Google Scholar
 Ç. Koç, T. Bektaş, O. Jabali, and G. Laporte, “The fleet size and mix pollutionrouting problem,” Transportation Research Part B: Methodological, vol. 70, pp. 239–254, 2014. View at: Publisher Site  Google Scholar
 E. Demir, T. Bektaş, and G. Laporte, “A review of recent research on green road freight transportation,” European Journal of Operational Research, vol. 237, no. 3, pp. 775–793, 2014. View at: Publisher Site  Google Scholar
 M. Barth, T. Younglove, and G. Scora, “Development of a Heavyduty Diesel Modal Emissions and Fuel Consumption Model,” Tech. Rep. UCBITSPRR20051, California PATH Program, Institute of transportation Studies, University of California at Berkeley, California, Calif, USA, 2005. View at: Google Scholar
 M. Scora and G. Barth, “Comprehensive Modal Emission Model (CMEM), Version 3.01, User Guide,” Tech. Rep., 2006, http://www.cert.ucr.edu/cmem/docs/CMEM_User_Guide_v3.01d.pdf. View at: Google Scholar
 M. Barth and K. Boriboonsomsin, “Realworld carbon dioxide impacts of traffic congestion,” Transportation Research Record, no. 2058, pp. 163–171, 2008. View at: Publisher Site  Google Scholar
 Man, Trucks in Distribution Transport, http://www.mantruckandbus.co.uk/en/trucks/start_trucks.html.
 E. Demir, T. Bektas, and G. Laporte, “An adaptive large neighborhood search heuristic for the PollutionRouting PROblem,” European Journal of Operational Research, vol. 223, no. 2, pp. 346–359, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 G. Zhang, H. Rong, F. Neri, and M. J. PérezJiménez, “An optimization spiking neural P system for approximately solving combinatorial optimization problems,” International Journal of Neural Systems, vol. 24, no. 5, pp. 1–16, 2014. View at: Google Scholar
 G. Zhang, H. Rong, J. Cheng, and Y. Qin, “A populationmembranesysteminspired evolutionary algorithm for distribution network reconfiguration,” Chinese Journal of Electronics, vol. 23, no. 3, pp. 437–441, 2014. View at: Google Scholar
 X. Zhang, J. Li, and L. Zhang, “A multiobjective membrane algorithm guided by the skin membrane,” Natural Computing, vol. 15, no. 4, pp. 597–610, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 J. Xu, “Probe machine,” IEEE Transactions on Neural Networks and Learning Systems, vol. 27, no. 7, pp. 1405–1416, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 G. Zhang, J. Cheng, M. Gheorghe, and Q. Meng, “A hybrid approach based on differential evolution and tissue membrane systems for solving constrained manufacturing parameter optimization problems,” Applied Soft Computing, vol. 13, no. 3, pp. 1528–1542, 2013. View at: Publisher Site  Google Scholar
 G. Zhang, M. J. PrezJimnez, and M. Gheorghe, RealLife Applications with Membrane Computing, Springer, Berlin, Germany, 2017.
Copyright
Copyright © 2018 Yunyun Niu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.