Route Selection of Multimodal Transport Based on China Railway Transportation

The advantage of multimodal transport is that it can deliver the goods to their destination in a reasonable combination of transport modes while ensuring security and punctuality. Multimodal transportation can eﬀectively reduce logistics costs, improve logistics eﬃciency, and reduce environmental pollution. In the process of multimodal transportation, due to the interference of natural factors (weather, terrain, etc.) and some special human factors, it may have diﬀerent degrees of impact on the transportation time and transportation safety of diﬀerent transportation modes. Therefore, when choosing a transportation method, it is necessary to consider the transportation time and transportation safety under the interference. However, the current research on multimodal transport has not considered the impact of external interference on transportation time and transportation safety. Compared with other modes of transportation, external interference has a relatively small impact on railway transportation. Railways can safely deliver goods to their destinations on time. Under the background of China’s huge railway network and advanced heavy-duty technology, this paper establishes a multimodal transport route selection model for considering railway as the core, introduces time penalty cost and damage compensation cost, and takes the lowest comprehensive transportation cost as the model objective under the premise of considering transportation reliability and transportation safety. Finally, taking a multimodal transport network in China as an example, an improved ant colony algorithm is used to solve the model and the results verify the rationality of the model.


Introduction
Multimodal transport [1] is a mode of cargo transportation through the connection and transshipment of two or more transportation modes (such as railway transportation, water transportation, road transportation, and air transportation) to complete the transportation process. At present, China is at an important stage of economic structural transformation and upgrading, and the freight transportation market is facing historical changes. e establishment of a new and efficient logistics channel based on multimodal transport can connect China's inland economy and coastal economy, which is the key to the transformation and upgrading of China's economic structure. is shows that promoting the development of multimodal transport is of great significance. According to statistics, in 2017, China's multimodal freight volume was 1.368 billion tons, accounting for 2.9% of the total freight volume of the society. According to the "Notice on Further Encouraging Multimodal Transport Work" by 18 departments including the Ministry of Transport, China's multimodal transport freight volume will reach 3.02 billion tons in 2020. By then, the scale of multimodal transport freight will account for 6% of the total freight volume%, the proportion is still low. e multimodal transport freight volume accounts for a small proportion of the whole society's freight volume, which is caused by the incomplete construction of China's current integrated transportation system. With the increasing market requirements, more efficient, more resource-efficient, and more intensive transportation services are required to achieve a rational division of labor and organic linkage of various modes of transportation, and to give full play to the overall efficiency and effectiveness of the integrated transportation system. However, at present, all kinds of transportation modes in the integrated transportation system are self-contained, which are only physically connected through hubs, and the limited joint command is carry out under the emergency condition.
ere are some problems, such as imperfect mechanism, market segmentation, prominent transport structural contradictions, weak hub service capacity, poor operational coordination, and low safety and service level, which make the great potential of integrated transportation not fully released. It affects the full play of the efficiency of various modes of transport and the public's sense of gain (quoted: the 13th Five-Year Plan of Integrated Transport Services), which cannot meet the economic, reliable, and efficient freight demand.
Compared with other modes of transportation, railway transportation is more economical and environmentally friendly, and the advantages of transportation reliability and transportation security are more obvious. At the same time, railway transportation as a participant of multimodal transportation has good transportation conditions. Data show that in 2018, the national railways completed 3.19 billion tons of cargo shipments, a year-onyear increase of 9.3%, which increased transportation by 272 million tons. With the formation of the "four vertical and four horizontals" networks of high-speed railways and the promotion of the "eight vertical and eight horizontals" high-speed rail network, the operating mileage of the railway reached 131,000 kilometers, of which the operating mileage of the high-speed railway was 31,000 kilometers. It is concluded that China has a huge railway network, and the operating mileage of railways covers the entire country. Although railways have a nationwide road network and huge transportation capacity, they have not been fully utilized. In 2018, China's highway freight volume was 39.568 billion tons, accounting for 76.79% of the country's total freight; China's railway freight volume was 4.022 billion tons, accounting for 7.81% of the country's total freight. So, making full use of railway capacity has become a key goal for China to promote the development of multimodal transport. In response to this situation, relevant Chinese departments have formulated a "road to rail transportation" policy, and railway management is included in the Ministry of Transport. ese policies make it possible for railways to play a backbone role in the integrated transportation system. So, it is urgent to study the cooperative operation theory and key technologies of the railway-based integrated transportation system. is paper will deeply study the optimization of multimodal transport path selection with railway transportation as the main body, and realize the rational division of labor and coordinated development of various transportation modes.
In multimodal transport network design, Christine and Sabine [2] conducted a scenario-based assessment of the future position of intermodal transport within Belgiumbased case studies and discussed service network design models for consolidation-based freight transport systems; Demir et al. [3] studied the design problem of green multimodal transport network with uncertain transportation time [4][5][6] and established a stochastic optimization model that can generate robust transportation plans according to different objectives (cost, time, and carbon emissions); Wang and Meng [7] considered a discrete multimodal transport network design problem and decided whether the network planner needs to establish or extend a path to minimize operating costs. Cheng and Jin [8] constructed the multimodal transportation route selection models with congestion considered under different low-carbon polices; Demir et al. [3] introduced a green intermodal service network design problem with travel time uncertainty (GISND-TTU) for combined offline intermodal routing decisions of multiple commodities. In the path optimization of multimodal transport, many scholars have made in-depth research based on different aspects. Fazayeli et al. [9] studied the problem of multimodal transport path optimization under fuzzy demand, established a mixed-integer fuzzy mathematical model, and solved it by genetic algorithms [10]; Fotuhi studied the optimization problem of multimodal transport path with uncertain network topology [11]; Idri et al. [12] introduced a search method for the shortest path algorithm in parallel distributed architecture to solve the path optimization problem of dynamic multimodal transport network with time dependence [13,14]; Liu et al. [15] established a super transportation network to predict the generalized cost of multimodal transportation, which provides technical support for the comprehensive transportation planning. Li constructed a route selection model for cargo multimodal transportation and used an improved ant colony algorithm to be solved. Tian et al. [16] explored the influence of multiple time windows on the route selection of fresh products and designed genetic algorithms to solve it. Zhu et al. established an evaluation model of container multimodal transport based on BP neural network and solved the optimal path. Russo et al. [17] research object is a variable of the operating cycle of a rail transit terminal, the variable considered is the overall average time of the terminal operation cycle relative to the number of trucks and vehicles entering the terminal, which solved the problem of the efficiency of public rail transport; Abderrahman et al. [18] proposed a robust optimization model for the transportation cost of multimodal freight and the uncertainty of network terminal capacity [19]; Liu et al. [20] in the context of global containerized transportation of soybeans considers the transportation and purchase costs of soybeans, establishes a multimodal transport network model, and optimizes the soybean transportation path [21]; Verma et al. [22] designed a dual-objective mixed-integer linear programming model for the problem of joint transportation path selection for dangerous goods with minimum transportation costs and risks of transportation and designed a modelsolving algorithm based on taboo search; Toumazi and Kwon [23] introduced the conditional value-at-risk (CVaR) theory and built a CVaR-based road transport route optimization model for dangerous goods under a time-varying network [16,24]. e optimal dangerous goods road transport can be selected based on the risk aversion of decision makers' path; Sun and Long [25] took into account the customer's need for timeliness of transportation when planning the multimodal transport path [26] and obtained the optimal solution through Pareto optimality.
At present, the research on multimodal transport does not consider the influence of natural factors such as weather, terrain, and some special human factors on the transport [15]. In actual transportation, these factors may lead to the delay [27] of transport time, and even the occurrence of dangerous accidents [28] and damage of goods. At this time, the cost and time of the original optimal path will change, and the optimal path will change accordingly [29,30].
e multimodal transport research based on railways mostly focuses on transportation policy, and there are few research projects on path selection. erefore, this paper establishes a multimodal transport route selection model based on railway transportation, introduces time penalty cost and damage compensation cost, considering the transportation reliability and transportation safety factors in the transportation process, aiming at the minimum transportation cost, and finds the optimal path that is in line with reality.

Problem Description.
e route selection problem of multimodal transportation mainly based on railway transportation can be described as follows: In a multimodal transport network G � (V, A), Vis the node set of networks G and Ais the path set of networks G. e shipper requires that a batch of goods be transported from the origination o ∈ A to the destination d ∈ A. e carrier needs to consider the reliability and safety of the transportation routes between the origination and destination, so as to select the transportation route and mode that can reach the destination within the time specified by the shipper and ensure that the goods are intact.

Variable Definitions.
e notation of parameters and variables is shown in Table 1.

Objective Functions
e goal of multimodal transport is to maximize the transport efficiency and spend as little cost and time as possible to make the goods arrive at the destination safely. Multimodal transport is the combined transport of multiple modes of transport, and transportation cost and transportation time for different modes of transport are different. erefore, when calculating the model target, it is necessary to choose a transportation mode with lower transportation cost and shorter transportation time. However, in the actual transportation process, the transportation cost and transportation time of various transportation modes are not fixed. ey will change with the interference of natural factors and some special human factors. Different transportation modes have different antijamming capabilities, and transportation costs and transportation time have different degrees of variation. For example, in foggy weather conditions, various modes of transportation will be affected, but rail transportation will be less affected, and the transportation time will be slightly extended; the speed of road transportation will be greatly reduced, and safety accidents may occur due to operational errors, not only the safety of life is involved but also the transported goods will be damaged. e carrier will have to compensate the shipper for losses, and the transportation cost and transportation time will increase significantly.
When the transportation cost and transportation time are dynamically changing, it will dynamically affect the choice of transportation mode. erefore, the model built in this paper focuses on the change of transportation cost and transportation time during transportation. In this paper, the degree of change in transportation time is expressed as transportation reliability, and the time penalty cost is established. e degree of change in transportation cost is expressed as transportation safety, and the damage compensation cost is established. At the same time, considering the global low-carbon environmental protection at the current stage, low-carbon transportation is the top priority for the country. Multimodal operation is an emerging efficient transportation mode. It is also necessary to reduce carbon emissions as much as possible in response to national calls. erefore, this paper will start from the cost, energy conservation, time, transport reliability, and transport safety to build a path optimization model with the minimum comprehensive cost: (1) In the formula, Z is the comprehensive transport cost; C 1 , C 2 , C 3 , C 4 , and C 5 are transportation cost, transit cost, carbon emission cost, time penalty cost, and damage compensation cost, respectively; and ω 1 , ω 2 , and ω 3 are the weights of each cost in the objective function.

Transportation Cost and Transit Cost.
e cost of freight transportation between adjacent nodes is the transportation cost, and the cost of replacing the transportation mode inside the transferable node is the transit cost:

Carbon Emission Cost.
In order to develop a low-carbon economy, countries all over the world have studied energy conservation and emission reduction and proposed a carbon tax scheme. Carbon tax is a tax imposed on carbon emissions Journal of Advanced Transportation 3 and is an important means of limiting carbon emissions. e calculation formula of carbon tax is as follows: In the formula, D t is the damage value of thet-year caused by the unit carbon emissions and λ is the discount rate. In fact, the carbon tax is not levied as high as possible. If the carbon tax is higher, it will restrict the economic development of the country. If the carbon tax is too low, it will not play a good role in energy conservation and emission reduction. According to the research of experts and scholars in China, the proposal of the Ministry of environmental protection for carbon tax is 20 $/ton. erefore, in this paper, the carbon tax is assumed to be 20 $/ton and 0.02 $/ kg. e carbon emission cost is the product of carbon emissions and carbon tax:

Transportation Reliability and Time Penalty Cost.
Transport reliability herein refers to the on-time arrival rate from the starting point to the end point within a specified time using either mode of transport. In the actual transportation process, various modes of transportation will be affected by natural factors and some special human factors, resulting in delays in transportation, such as traffic jams caused by road conditions; trains will stop and wait during the meeting, and long ramps will have a great impact on the running speed and so on. Transportation delay will prolong the transportation time and may exceed the delivery time limit, so if the model does not take into account the transportation reliability factors, the obtained optimal path will lose its significance in the actual transportation context. In order to quantify the impact of transportation reliability on transportation time, this paper defines the time delay coefficient: taking into account the complex influencing factors between adjacent nodes and taking values within 0-1 so that the transportation time of the section is extended linearly on the basis of the original transportation time, as shown in Figure 1.
In the diagram, t represents the transportation time between two nodes without considering the delay, and μ is the time delay coefficient. When μ � 0, there is no time delay in the transportation process of goods, and the actual transportation time of goods is the initial transportation time t. When μ � 1, there is time delay in the transportation process of goods, and the actual transportation time is the initial time t plus the delay time quantified as t. When μ is between 0 and 1, the delay time of goods is related to the delay coefficient, defining that the delay time has a linear relationship with the delay coefficient.
A certain penalty is imposed for transportation schemes that exceed the time limit. In order to quantify the impact of the time limit of delivery on the transportation plan selection, the time penalty function is introduced in this paper. Taking into account the changes in the psychological state of the consignee while waiting for the delivery of the goods, this paper adopts the fuzzy convenience function to express and defines 2 hours as a state change gradient in case of exceeding the arrival time limit and assigns different degrees of punishment to different state gradients, as shown in Figure 2. Total quantity of goods to be transported in the multimodal transport model π Unit value of goods that need to be transported ψ s ij Accident damage coefficient: from node city i to node city j, transport modesis used for transportation d s ij Transportation mileage: from node city i to node city j, transport modesis used for transportation μ s ij Time delay coefficient: from node city i to node city j, transport modesis used for transportation c s ij Unit transportation cost: from node city i to node cityj, transport modesis used for transportation c ss' i Unit transfer cost: at the transferable node i, the mode of transport is changed fromsto s′ T ss ′ Journal of Advanced Transportation Formula (5) is the transportation time calculation formula for any transportation scheme: the first part is the sum of transportation time and delay time, and the second part is the transit time of the transportable node. Formula (6) is the time penalty cost calculation formula, when the transportation time is from T to T + 2, the penalty function is a linear increasing function; when the time is from T + 2 to T + 4, the penalty function is a constant (y 1 ); when the transportation time exceeds T + 4, the penalty function is a constant (y 2 ).

Transportation Safety and Damage Compensation Costs.
Transport safety in this paper refers to the value loss rate of goods transported during transportation from the starting point to the end point. e influencing factors mainly include road conditions and traffic management status, vehicle technical performance, and warranty quality, as well as the driver's operating skill level and responsibility. If a safety accident occurs during transportation, the transportation goods will be damaged. e carrier shall compensate the consignor for the corresponding loss according to the degree of damage to the goods because the liability of the accident is that the carrier has not planned the transportation route. If the safety factors in the transportation process are not considered when optimizing the route, the final result is likely to lose more. In order to facilitate the calculation, this paper assumes that the degree of damage to the goods is determined by the severity of the accident when a safety accident occurs during transportation. at is to say, if the severity of the accident is relatively low, the goods may only be partially damaged rather than completely damaged. In order to quantify the impact of the severity of the accident on the damage of the goods, this paper proposes and defines the accident damage coefficient: the value is within 0-1, and the change from 0 to 1 indicates that the degree of damage of the goods increases gradually and the amount of compensation required by the carrier increases gradually, and the amount of compensation is based on the value of the goods:

Transportation Process Constraint
s,s′∈S Constraint (8) makes sure that only one mode of transport can be selected between any adjacent nodes. Constraint (9) ensures that only one transit can occur in any transit node. Constraint (10) guarantees the continuity of transportation mode during transit. Constraint (11) is the flow equilibrium constraint of nodes in the network.

Carbon Emission Constraint
Cd ≤ E.
Constraint (12) is the formula of total carbon emissions, and for this paper, which is a known constant. Constraint (13) ensures that the transportation plan meets the total carbon emission limit.  Journal of Advanced Transportation multiobjective optimization problem in discrete space; however, the basic ant colony algorithm has the problems of slow running efficiency, premature convergence, and easy to fall into a local optimal solution. is paper combines the solution requirements of the multimodal transport path selection model. In order to improve the effectiveness of the algorithm, the basic algorithm is improved accordingly.

Node Transition Probability.
Although the objective function of the model in this paper is the total cost value in the transportation process, the time penalty cost is determined by the length of the transportation time. e model target expectation cost and time are as small as possible, so in the improved ant colony algorithm, a cost heuristic function and a time heuristic function need to be set at the same time. It will mainly rely on the influence of path pheromone concentration, cost heuristic experience, and time heuristic experience when choosing to reach the node: Formula (14) is the node transition probability formula in the multiobjective ant colony algorithm, where α is the information heuristic factor, β is the cost expectation heuristic factor, λ is the time expectation heuristic factor, and τ ij is the pheromone concentration on path (i, j). Formula (15) is a cost heuristic function, whose value is inversely proportional to the sum of transportation cost, carbon emission cost, and damage compensation cost on path 1 and transit cost of inode. Formula (16) is a time heuristic function, whose value is inversely proportional to the sum of transportation time on path (i, j) and transit time of inode.

Pheromone Update Rule.
e principle of the positive feedback mechanism of ant colony algorithm is based on the continuous update of pheromone, so the pheromone update rule is the key to solving the problem of speed and accuracy of ant colony algorithm. is paper considers the basic idea of ant colony algorithm with elite strategy: at the end of each cycle, additional pheromone enhancement is given to the optimal solution found in this cycle in order to make the optimal solution found so far more attractive to the ants in the next cycle. e pheromone is updated according to the following formula: j) is part of the optimal path in this cycle, 0, otherwise.
Formula (17) is the pheromone concentration on path (i, j) at time t + n, and ρ is the pheromone global volatilization coefficient, Δτ ij (t, t + n) is the amount of pheromone increase on path (i, j) from time t to time t + n. Formula (18) is the amount of pheromone increase on path (i, j) caused by elite ants from time t to time t + n, q is the pheromone intensity, σ is the number of elite ants, and C * is the optimal solution value found in this cycle.
Using the elite strategy can make the ant system find a better solution through fewer cycles, and the solution speed is faster. However, after using the elite strategy, the attractiveness of the current optimal solution increases, so that the search process is concentrated near the optimal solution found so far, thereby preventing further searches for a better solution. e problem of premature convergence is more likely to occur. In order to maintain the selection pressure, this paper extends the concept of ranking in genetic 6 Journal of Advanced Transportation algorithms to the elite mechanism and calculates Δτ ij (t, t + n) by a weighted method. After completing a cycle, the ants are sorted according to the size of the target value found in the cycle, and the ant's contribution to the pheromone update is considered to be inversely proportional to the target value of the path the ant passes. e contribution of the ant to the pheromone update is weighted according to the ranking of the ant (ξ). In order to ensure the speed of the solution while expanding the solution space, it is considered that only the top (m/4) ants can contribute to the pheromone update, and the amount of pheromone obtained by the path it passes is proportional to the ant's ranking. Δτ ij (t, t + n) is calculated as follows: Formula (20) represents the contribution of the ant ξ to the increase in the amount of pheromone on the path (i, j), ξ is the ranking of the ants, C av represents the average value of the target of this cycle, and C ξ represents the target value found by the ant ξ.

Pheromone Limit Interval.
Although the aforementioned improved pheromone update rule can solve the problems of solution speed and precocious convergence, it is easy to fall into the dilemma of local optimization and the algorithm stops searching. erefore, in this paper, the idea of MMAS is introduced in the process of pheromone updating. e upper and lower limits are imposed on the pheromone, and the limiting interval of the pheromone is set to[τ min , τ max ]. When performing pheromone update according to the aforementioned pheromone update rules, the following is obtained: is can effectively avoid that the amount of pheromone on one path is much larger than the rest, so that all ants are concentrated on the same path, resulting in a local optimal situation.

Algorithm Flowchart.
In summary, the multiobjective ant colony algorithm flowchart is as shown in Figure 3.

Case Study
Take China's multimodal transport network with 10 cities as an example, it involves 10 cities, the origination is Chengdu and the destination is Shanghai, as shown in Table 2. Since there are many inland cities in this multimodal transport network and the water transport is poor, we considered three transport modes: road, railway, and air in this case. e node cities in the multimodal transport network are numbered sequentially, and the schematic diagram of the multimodal transport network of this case is made accordingly, as shown in Figure 4.
is case involves three transportation modes. We need to consider how to clearly express the ant's choice behavior for transportation mode when using the ant colony algorithm. erefore, we have simply dealt with the multimodal transport network. Each transport node is expanded into multiple virtual subnodes according to the number of transport modes. For example, transfer node 3 is extended to railway node, highway node, and airport node. In the processed network diagram, there is a transport problem among multiple subnodes of the same node. But it is not the ant to select the next arrival node, that is to say, there is no sequence between the subnodes, which is the same level parallel relationship. e network deformation diagram is shown in Figure 5. e transportation mileage, time delay coefficient, and accident damage coefficient of different transportation modes among node cities are shown in Table 3. e time delay coefficient is given according to the road condition of each node city and the different characteristics of transportation mode. e accident damage coefficient is given by the traffic management status of each node city and the technical performance of each transportation mode. e requirements for freight volume, transportation time limit, and unit value of goods by freight forwarders and the requirements for carbon emission limit by government officials are shown in Table 4. e unit transit cost and unit transit time between different modes of transport are shown in Table 5. e transportation cost, average speed, and carbon emission factors of each transportation mode are shown in Table 6.

Results' Analysis.
e multiobjective ant colony algorithm is used to write the algorithm code on the MATLAB R2014a platform. e algorithm parameters are set as follows: λ � 3, α � 1, β � 3, q � 2, ρ � 0.9, the number of ants is 100, and the number of iterations is 200, y 1 � 10000, y 2 � 20000. Taking the weight of the objective function as 1/3 as an example, the optimal path of multimodal transport is Chengdu to Chongqing to Wuhan to Hefei to Nanjing, and the total cost is 331750. To explore the results of multimodal transport route optimization considering different cost factors, the weight of each objective function is assigned in turn, and the following conclusions are drawn in this paper.

Objective Weight Sensitivity Analysis.
In order to explore the multimodal path optimization results when considering different cost factors, this paper analyzes the sensitivity of the objective function weights as shown in Table 7 and draws the following conclusions: In Figure 6, Experiment (1,5) and (2,4) are compared, it is found that when the weight of damage compensation cost Journal of Advanced Transportation is the same, and the weight of the sum of transportation cost, transshipment cost, and carbon emission cost is high, the optimal path is the railway-road combined transportation. e railway-road combined transport pays more attention to the cost of transport. It pays more attention to the size of the cost incurred in the transportation process. When the weight     of time penalty cost is high, the optimal route is mainly airrailway combined transportation, and more attention is paid to the length of transportation time. As the carbon emission factor of road transportation is relatively high and the transport time of railway transportation is relatively long, the carbon emission cost of Experiment 1 and 2 is relatively high and the time penalty cost is relatively high. Due to the poor safety of air transport compared with other modes of transport, the cost of compensation for damages in Experiment 4 and 5 is relatively high. Generally speaking, in the comprehensive consideration of cost, time, energy saving, and transportation safety, whether it is air-railway combined transport or road-railway combined transport, railway occupies a large proportion, which also proves that multimodal transport based on railway transportation can indeed maximize the transport benefits.    Weight assignment

Analysis of Transportation Reliability and Transportation
Safety. In order to explore the practicability of the model, all cost weights are set to 1/3, and the ant colony algorithm is used to solve the problem without considering the transportation reliability and transportation safety. From Table 8, it can be concluded that the highway has a higher proportion in the optimal path solved without considering the transportation reliability and transportation safety. However, the transportation reliability and transportation safety of road transportation are poor, so the total cost of this optimal path in the model constructed in this paper exceeds the total cost of multimodal transportation optimal path based on railway transportation. is result is also in line with China's policy of transitioning from road transport to railway transport because railway transport is more stable and safer, and it can transport goods safely and in a timely manner. With the continuous development of other modes of transportation, in the future, China or other countries may have a multimodal transport mode based on road transport or other modes of transport. However, transportation reliability and transportation safety are the factors that must be considered in the multimodal path optimization problem; therefore, the model established in this paper is still applicable to the route selection and optimization of other transportation modes.

Conclusions
e paper analyzes and demonstrates the important position of railway transportation in multimodal transport and establishes a multimodal transport route selection model based on railway transportation. In view of some problems that often occur during transportation, this paper focuses on transportation reliability and transportation safety, introduces time penalty cost and damage compensation cost, and defines time delay coefficient and accident damage coefficient, respectively. Taking a multimodal transport network in China as a case, and the network is transformed according to the characteristics of the ant colony algorithm. Sensitivity analysis of cost weights proves that the multimodal transport based on railway transportation can maximize transportation efficiency. e analysis of transportation reliability and transportation safety shows that the model constructed in this paper has strong practicability and wide application range.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.