5.1. Test Cases
The benchmark data sets of Cordeau [
43] are well-known in MDVRPTW. In this study, we choose the narrow time window problem pro01–pro10 from the typical database. Then, they are applied to the MDOVRPTW and used to test the two-phase algorithm. The distance in the test instances is the Euclidean distance and supposing that between two nodes the travel time of the vehicle is equal to the Euclidean distance. This article uses the literature [
44] classification method to divide the test cases into three categories according to the number of customers: small-scale study 1–100 (pro01, pro02 and pro07), medium-scale study 101–200 (pro03, pro04 and pro08), large-scale study 201–300(pro05, pro06, pro09 and pro10).
Table 4 shows part of information about the test instances, which includes the number of customers
C, depots
D, vehicles
H, maximum vehicle load
Q0 and longest duration of each route
L. And we set the parameters of this model according to the former studies [
16,
45,
46], which are shown in
Table 5. According the
Table 2,
can be calculated as 2.623. We set the carbon trading price
as 25 CNY/t, which refers to the literature [
37].
5.3. Effectiveness of Two-Phase Algorithm
As we all know, the carbon quota is generally set by the government. Here, we initially assume that the carbon quota is 0 and we set the carbon trading price as 0.025 CNY/kg. To investigate the effectiveness of the algorithm, we propose the traditional particle swarm optimization to compare with the proposed two-phase algorithm. For each of the following experiments, we perform them 20 times and record the best value as the optimal result.
Table 7 shows the detailed computational results of PSO and PSO-TS, which include total costs, length of route, carbon emission and optimization rate of total costs.
We can easily see from
Table 7, compared with PSO algorithm, total costs, length of route, carbon emissions of the PSO-TS algorithm all have a great improvement in the quality of solution. And the average optimization rate of total costs reaches to 35.21%.
In this paper, we propose three scale cases: small, medium and large and we calculate the ten cases (pro01–pro10) with the PSO-TS algorithm.
Table 8 shows the detailed calculation results, which include vehicles’ number, drivers’ salary, penalty costs, fuel costs (fuel consumption costs), carbon costs (carbon emission trading costs), total costs, carbon emissions and length of route.
From the
Table 8, we find the drivers’ salary of these cases (pro01, pro02 and pro07) are all zero and those cases all belong to small scale. We propose an assumption that the two-phase algorithm is more suitable for small-scale cases.
To verify the applicability of the PSO-TS algorithm, we calculate the average optimization rate of objective function and the average running time of the three scales. The results are shown in the
Figure 5. The average running time demonstrates the computational complexity of the algorithm. The average optimization rate of objective function illustrates the performance of the algorithm.
It can be seen in
Figure 5, with the increase in the size of cases, the PSO-TS algorithm needs to spend more time to get the optimal result and the best average optimization rate of objective function 43.73% is of the small-scale case. Considering the speed of convergence and the quality of solution, the proposed two-phase algorithm is more suitable for small-scale problems.
5.4. Experimental Results
In the carbon trading environment, carbon quotas and carbon trading prices have a direct bearing on environmental protection costs (carbon emission trading costs) and they can indirectly change the vehicle arrangements and route planning, which will further affect economic costs (penalty costs and drivers’ salary costs) and energy costs (fuel consumption costs). In the following, we conduct a detailed study on carbon trading prices and carbon quotas respectively.
In the following research, we only study the three cases (pro01, pro02 and pro07), because the proposed algorithm is more applicable to small-scale cases. What is more, we perform 20 times for each of the following experiments and record the best value as the optimal result. Then decoding the optimal solution can obtain total costs, each sub-cost, carbon emissions and so on.
First, in order to study the impact of carbon trading price on the objective function and carbon emission, we set up a comparative experiment for small-scale cases. Here, we just consider the carbon trading price, thus the carbon quotas are 0. We set the carbon trading price to 0, calculate the values associated with the objective function and compare it with the previous experiment results shown in
Table 8 when the carbon trading price is 0.025.
Table 9 shows the detailed values of carbon costs, fuel costs, total costs, carbon emissions when the carbon trading price is 0 and 0.025 respectively.
According to the results in
Table 9, the following findings can be observed. When carbon trading prices are 0, the carbon costs are also 0. When carbon trading prices is 0.025, the carbon costs are bigger than 0. That the carbon trading price is 0 means there is not a carbon prices constraint. Compared with the case in which the carbon trading price is 0, the fuel costs and total costs of the three cases (pro01, pro02 and pro07) are all lower when carbon trading price is 0.025. In addition, carbon emissions of the three cases (pro01, pro02 and pro07) have also been effectively improved when carbon trading prices are 0.025. In conclusion, it shows that under the constraints of carbon trading prices, even if the carbon costs increase, the total costs as the objective function and fuel costs can reduce and the carbon emissions can also reduce.
Next, we study the effect of carbon quotas on the objective function and carbon emissions for fixed carbon trading prices. Here, the fixed carbon trading prices are 0.025. The carbon emissions quotas are difficult to estimate the exact value. In this study, for each case, we set 3 new carbon quotas (
) around the carbon emissions (pro01:615.50, pro02:1411.18, pro07:1208.70), which can be obtained from
Table 8 when the carbon trading price is 0.025 and the carbon quota is 0. Aimed at the new different carbon quotas, the values of total costs, carbon costs, carbon emissions and the difference between
and carbon emissions (CE) can be calculated respectively, which are shown in
Table 10.
According to the results in
Table 10, the following findings can be observed. With the increase of carbon quotas, for each case, the total costs and carbon costs also increase. Since the carbon quotas are a small increase, the increase of total costs and carbon costs is slight. However, as the carbon quotas increase, the carbon emissions are a fixed value for each case (pro01:615.50, pro02:1411.18, pro07:1208.70). We know that the calculation of carbon costs is related to the difference between carbon quotas and carbon emissions. The increase in carbon costs is due to the increase in the difference between CE and
. For each case, the change of this difference (
−CE) can be seen from
Table 10.
On the basis of the findings, when carbon quotas change, carbon emissions are the fixed value for each case. We know that carbon emissions are related to the distribution paths. Therefore, we propose an assumption that when the carbon quota changes, the distribution paths are unchanged. Then, we use the case of pro01 to verify this hypothesis. According to
Table 10, we set the carbon quotas (
) as 600, 650 and 700. Then, the optimal distribution paths for pro01, which are obtained by solving the model, are shown in
Figure 6,
Figure 7 and
Figure 8.
According to
Figure 6,
Figure 7 and
Figure 8, we can easily see that the optimal distribution paths for pro01 are same when the carbon quotas are 600, 650 and 700. Therefore, the changed carbon quotas have no effect to the optimal distribution paths, which supports our hypothesis.
At last, in order to further study the impact of carbon trading on the objective function and carbon emission, we add two new carbon prices 0.015 and 0.035 based on the study of fixed carbon prices 0.025. Thus, a further study of changing carbon prices and carbon quotas is formed. The computational results of the changing carbon prices and carbon quotas are shown in
Table 11.
Figure 9 shows the changing trend of total costs under different carbon trading price and carbon quotas and
Figure 10 shows the changing trend of carbon costs under different carbon trading price and carbon quotas.
- (1)
When the carbon trading prices are 0.015, 0.025 and 0.035, the corresponding carbon emissions are also fixed at 661.27, 615.50 and 595.42, respectively. This further validates that for a fixed carbon trading price, carbon emissions and optimal distribution paths will not change when the carbon quota changes. As the carbon trading prices increase, the carbon emission decrease. It also further validates that carbon trading prices directly affect carbon emissions.
- (2)
When the carbon trading prices and carbon quotas increase within a certain range, the total costs and carbon costs have the downward trend. When the carbon trading price is 0.035 and the carbon quotas are 700, the total costs have the lowest value 1681.53 and the carbon costs have the lowest value −3.66.
5.5. Discussion and Analysis
In this study, the proposed PSO-TS algorithm is proven to be more applicable to small-scale problems (pro01, pro02 and pro07), which can get the better solution quickly by this method. For the small-scale MDOVRPTW, we study the impact of carbon trading on total costs, carbon costs and carbon emissions from both carbon trading prices and carbon quotas. The main summings-up are listed as follows:
- (1)
When the carbon trading prices exist, the carbon costs, fuel costs and total costs are all lower and the carbon emissions are also effectively improved compared with the trading prices as zero.
- (2)
When the carbon quota is increasing and the carbon trading price is fixed, the total cost and carbon costs also increase but the carbon emissions fixed on a certain value. We propose and prove that the changed carbon quotas don’t affect the optimal distribution paths, which also validates that carbon trading prices directly affect carbon emissions.
- (3)
When the carbon trading prices and carbon quotas both increase within a certain range, the total costs and carbon costs are both decreasing. For the fixed carbon quotas, when the carbon trading prices increase, the carbon emissions will decrease. For the fixed carbon trading prices, when the carbon quotas increase, the carbon emissions will be a fixed value.
Based on these results, some advice to third-party logistics enterprises and the government is forthcoming.
For logistics enterprises, they use the technological means and path optimization methods to reduce carbon emissions to create low-carbon logistics. Technological means can develop new energy or use energy-saving cars. However, these require a huge investment compared to path optimization. Therefore, this research is very necessary for logistics companies. It can quickly optimize the distribution route for government-specific carbon trading policies. According to this study, the logistics enterprises can also take some optimization measures to balance economic costs, energy costs and environmental costs. First, they should have the awareness of environmental protection and introduce the carbon emissions into the distribution path optimization. Secondly, from the operation and management level, it can reduce business operating costs and establish a good corporate image when the carbon emissions are taken into consideration. Thirdly, logistics enterprises should respond positively to the carbon trading policies, which are proposed by the government. They can choose low-carbon strategies to reduce total costs.
For the government, in order to promote low-carbon development, they can use carbon trading policies to reach the goal of fewer carbon emissions. Firstly, the government should encourage the enterprises to save energy and reduce carbon emissions. Secondly, from the management insight, they should form a relatively complete carbon trading policies and strengthen the regulation of carbon emissions. Research shows that the government can increase carbon prices in a certain range to reduce carbon emissions. The government can give logistics companies corresponding carbon quotas as much as possible to save the company’s costs without increasing carbon emissions. It can encourage the companies to respond to the carbon trading policies.