Next Article in Journal
Financial Aspects of Sustainable Rainwater Management in Small-Scale Urban Housing Communities
Previous Article in Journal
A New Land Use Dataset Fusion Algorithm for the Runoff Simulation Accuracy Improvement: A Case Study of the Yangtze River Basin, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Route Optimization for Hazardous Chemicals Transportation under Time-Varying Conditions

School of Management, Shanghai University, Shanghai 200444, China
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(2), 779; https://doi.org/10.3390/su16020779
Submission received: 13 December 2023 / Revised: 7 January 2024 / Accepted: 9 January 2024 / Published: 16 January 2024
(This article belongs to the Section Hazards and Sustainability)

Abstract

:
Since accidents of hazardous chemicals transportation will cause serious loss to the surrounding environment and lives and properties, this paper studies the transportation route optimization problem of hazardous chemicals under dynamic time-varying conditions. Combined with the goal of green sustainable development, a multiobjective nonlinear optimization model is constructed to minimize the transportation risk, transportation cost, and carbon emissions generated in the transportation. The model is solved by the improved Fast Non-Dominated Sorting Genetic Algorithm with Elite Strategy (NSGA-II) algorithm. The effectiveness of the model and the algorithm are tested on the Sioux Falls network. The experimental results show that under time-varying conditions, a vehicle’s departure at different times will generate different transportation costs and risks. Therefore, enterprises need to rationally arrange the departure time of vehicles according to the time windows of customer nodes and road conditions. In additio, from the relationship between the optimization objectives, in order to achieve green, sustainable and low-risk transportation, enterprises should first reduce their transportation costs.

1. Introduction

Hazardous chemicals refer to highly toxic chemicals and other chemicals with the properties of poison, corrosion, explosion, combustion, and other hazards to human beings, facilities, and the environment. China is the world’s largest producer of hazardous chemicals. As of the end of December 2021, there were as many as 26,947 chemical enterprises in China, with an annual total profit of CNY 1.16 trillion [1]. The chemicals industry is characterized by “high transportation demand, wide distribution of production capacity, and rapid growth in production”. From the perspective of China’s economic development and the demand of the petrochemicals industry, most of the chemicals processing enterprises are concentrated in the eastern coastal areas, while the production areas of chemical raw materials are relatively concentrated in the western areas. This phenomenon of uneven distribution of production and marketing increases the demand for the transportation of hazardous chemicals, with road transportation accounting for 70%.
Compared with developed countries which have established complete hazardous chemicals facilities and regulatory systems, China’s hazardous chemicals transportation industry is in a stage of rapid development. The supporting facilities and infrastructure of the hazardous chemicals transportation industry are relatively weak, the professional skills of the practitioners need to be improved, and the relevant laws and regulations also need to be improved. These reasons are intertwined and result in accidents of hazardous chemicals road transportation. According to relevant statistics, a total of 344 safety accidents in the transportation of hazardous chemicals occurred from January to October in 2021 [2,3]. Once an accident occurs during the transportation of hazardous chemicals, it will cause serious damage to the surrounding population, environment, and economic property [4]. On 13 June 2020, an explosion of a liquefied petroleum gas tanker occurred on the G15 Shenhai Expressway near Liangshan Village, Daxi Town, Wenling City, Zhejiang Province. The explosion caused the extensive collapse of surrounding residential and factory buildings, resulting in a total of 20 deaths and 175 injuries, with a direct economic loss of CNY 94.77815 million [5].
The vehicle routing problem (VRP) was first proposed by Dantzig and Ramser in 1959 [6]. Afterwards, the vehicle routing problem with time windows (VRPTW), multidepot vehicle routing problem (MDVRP), vehicle scheduling problem with full load (VRPFL), and other derivatives of VRP have been proposed one after another. Desrochers et al. [7] solved the VRPTW with an unconstrained number of vehicles by utilizing the branch-and-bound and dynamic programming methods. Thangiah et al. [8] were the first to use genetic algorithms for solving the VRPTW. Wren and Holliday [9] proposed a scan-based algorithm and applied it to the multivehicle route problem. Golden et al. [10] proposed two heuristic algorithms for solving the MDVRP. The hazardous chemicals transportation route optimization problem belongs to a special branch of the VRP. The transportation of hazardous chemicals involves multiple subjects, including hazardous chemicals productions and transportation enterprises, governments, and the public. The objectives pursued by each subject in the process of hazardous chemicals transportation are different. Among them, the hazardous chemicals production and transportation enterprises usually prioritize minimizing the transportation costs from the perspective of economic interests; however, governments want to minimize the risk of hazardous chemicals transportation. So, the hazardous materials transportation route optimization problem is usually a multiobjective vehicle route optimization problem (MOVRP).
Transportation route optimization is an effective way to reduce the occurrence of accidents. Therefore, in order to reduce risk in the road transportation of hazardous chemicals and realize sustainable development, this paper optimizes the road transportation routes of hazardous chemicals by constructing a mathematical optimization model under time-varying conditions, solved by using the improved NSGA-II algorithm, which can expand the research scope of the optimization theory of hazardous chemicals road transportation. In addition, while minimizing transportation costs and transportation risks, the reduction of carbon emissions is incorporated into the optimization objective, which is important for the reduction of carbon emissions and the decision-making and operational management of hazardous chemicals transportation enterprises. The structure of this paper is arranged as follows: Section 2 gives a brief introduction to the MOVRP of hazardous materials transportation, and the optimization problem of hazardous chemicals transportation under time-varying conditions. Section 3 describes the multiobjective time-varying route optimization model for hazardous chemicals transportation. Section 4 proposes an improved NSGA-II algorithm and introduces the implementation process of the algorithm. Section 5 validates the effectiveness of the algorithm on the Sioux Falls network and analyzes the experimental results. Finally, Section 6 makes a conclusion and discusses the prospects for future research directions.

2. Literature Review

2.1. MOVRP for Hazardous Chemicals

Due to the high risk from hazardous chemicals during road transportation, the hazardous chemicals transportation route optimization problem usually needs to minimize the transportation risks and transportation costs. Konstantinos and Androutsopoulos [11] proposed a biobjective hazardous chemicals road transportation optimization model with a time window to minimize the transportation risk and time, and the biobjective optimization problem was transformed into a single-objective optimization problem, solved using the weighted summation method. Zou and Zhang [12] proposed five major hazardous chemicals transportation route selection indexes from the perspective of the transportation subject, and constructed a multiobjective hazardous chemicals transportation route optimization model with a hybrid time window, which was solved by using an improved genetic algorithm. Yuan et al. [13] established a biobjective optimization model considering both transportation cost and transportation risk, introduced a population hybridization strategy to improve the particle swarm algorithm, and tested it on a benchmark algorithm. Bula et al. [14] proposed a biobjective optimization model to minimize the transportation cost and the transportation risk from the point of view of different interested parties and solved it by using a neighborhood search algorithm. Jiang et al. [15] investigated the heterogeneous vehicle problem of hazardous materials transportation with time windows to minimize the transportation cost, transportation risk, and the average number of vehicle redundancies, and used a hybrid multiobjective algorithm to solve the problem. Chai [16] proposed a new driving risk during hazardous chemicals transportation taking into account the risk evaluation model, and constructed a biobjective optimization model for minimizing the transportation risk and transportation cost, which was solved by using a modified nondominated sequential genetic algorithm. Esmaeilidoukia et al. [17] proposed a fuzzy planning model for minimizing the time and risk, which was transformed into a deterministic model and then solved using an invasive weed algorithm.
The growth of carbon emissions and global climate change has become an important factor affecting the survival and development of mankind, and taking the path of green, low-carbon and sustainable development has become a global consensus. Along with the development of a global low-carbon economy, the green vehicle route transportation optimization problem (GVRP) has also attracted a lot of attention. Rahbari et al. [18] investigated the negative impacts of greenhouse gas emissions on the environment caused by hazardous materials and waste materials from transportation. They expressed the greenhouse gas emissions as a function related to the transportation distance, and established a location–inventory–routing problem, which minimizes the transportation cost, risk, and carbon emissions, and was solved using a multiobjective black widow optimization algorithm. Zhao and Cao [19] expressed carbon emissions as a function related to fuel consumption during transportation in the VRP problem, incorporated it into the transportation cost as the minimization objective, and constructed a biobjective optimization model to minimize the cost and risk. Wang et al. [20,21] investigated the effect of traffic flow uncertainty on transportation in the vicinity of an intra-city transportation route, expressed the carbon emissions as a function of vehicle travel speed, distance, and vehicle loading, and constructed a multiobjective optimization model minimizing the transportation risk, transportation cost, and carbon emissions, which was solved by using the improved NSGA-II algorithm. Lyu and He in [22] studied a multijourney heterogeneous vehicle route problem for hazardous chemicals, constructed a single-objective optimization model for minimizing the cost in two stages by adding the cost of carbon emissions to the transportation cost, and expressed it as a function of fuel consumption and vehicle loading. The model was solved using a two-stage mixing and summing metaheuristic algorithm. Sun et al. [23] investigated the multimodal transportation problem of hazardous materials by rail and public transport, considering the uncertainty of population exposure, constructed an uncertain optimization model to minimize the transportation risk and transportation cost, and used the carbon emissions related to the transportation distance and loading as the constraints of the optimization model to limit the maximum carbon emissions generated during transportation.

2.2. The Application of Spatiotemporal Data and Information Technology in Transportation

Vehicles will be affected by the combination of various factors such as time, weather and crowd flow in transportation. In the era of big data, with the development of information technology and the easier acquisition of dynamic time-varying data, many scholars use modern information technology to study dynamic spatiotemporal data in transportation. Meng et al. [24] studied the use of web-based real-time traffic data to estimate vehicle carbon emissions from the perspective of reducing urban carbon emissions. Experiments were performed in Chengdu, China, and results showed roadway-based emissions matched well spatio-temporally with traffic conditions. Xiao et al. [25] analyzed the relevant literature on spatiotemporal crowd flow prediction in Web of Science, from the perspectives of traffic control and public safety, using scientometric methods, social network analysis, and stochastic actor-oriented models (SAOMs), and proposed a knowledge evolution path construction process based on SAO structures. Zhang et al. [26] studied the application of dynamic data analysis in pedestrian movement risk and constructed a complete pedestrian movement risk identification process using video recognition technology and a support vector model (SVM). Li et al. [27] proposed a data-driven traffic shock wave detection technique, which combines machine learning techniques with traffic features, and conducted experiments on collecting vehicle trajectory data for highways in the United States, demonstrating the efficiency and accuracy of the proposed four-step method. Pourmoradnasseri et al. [28] studied the application of real-time traffic data flow on improving transportation systems to address the problems caused by traffic congestion in cities, and used road Internet of Things (IoT) data as the data source for a two-layer optimization model. The time dependence of traffic status was ensured by splitting the calculation process into a short-term framework during solving, and the efficiency of this method in processing real-time data flow was verified in a road network experiment in Tartu city. Chen et al. [29] proposed a framework that integrates a feature-enhanced scale aware descriptor, Kalman filter, Hungarian algorithm, and perspective projection theory to extract vehicle speeds from port surveillance videos for analyzing automated guided vehicle trajectories.

2.3. Hazardous Chemicals Dynamic and Time-Varying Route Optimization Problem

The road transportation of hazardous chemicals is a complex process related to time, and changes in the environment under different time conditions will have different effects on the transportation risk. Liu and Zhu [30] considered the effects of uncertain population density and terminal demand, expressed the accident rate during transportation as a time-dependent function, constructed a stochastic optimization model, and solved it using a hybrid heuristic algorithm combining particle swarm algorithms and genetic algorithms. Ke [31] explored the impacts caused by system disruptions in the hazardous material emergency logistics system, and proposed a risk metric model under discrete time-varying conditions. Ouertani et al. [32] investigated the dynamic route optimization problem of hazardous materials transportation with time windows, and considered that the customer’s demand changes dynamically over time, so constructed a biobjective optimization model that minimizes both cost and risk, and solved it using a hybrid meta-heuristic algorithm. Xu et al. [33,34] studied the real-time route planning problem for hazardous chemicals transportation under the condition of changing road traffic service level. They first used the Dijkstra algorithm for initial route global planning, and then updated the local route in case of changes in traffic service level. He et al. [35] studied the hazardous chemicals transportation route optimization problem under discrete time-varying conditions and established an optimization model to minimize the conditional value-at-risk, and the results showed that the optimal route for transportation under different departure times would be different.
In summary, most of the current research on road transportation of hazardous chemicals focuses on the minimization of transportation risk and transportation cost under static conditions, and the carbon emissions during transportation are usually expressed as a function of transportation distance or transportation fuel consumption. During the road transportation of hazardous chemicals, spatiotemporal data has not been effectively utilized. The time factors of road transportation of hazardous chemicals are less frequently studied and are usually expressed as the time windows of the customers. However, the transportation cost, the incidence of transportation accidents, and the consequences of transportation accidents during the transportation of hazardous chemicals are all related to time. Therefore, this paper proposes a multiobjective optimization model for the transportation of hazardous chemicals under discrete time-varying conditions. The model considers the different accident rates and consequences of hazardous chemicals transportation under different time conditions. In addition to minimizing transportation risk and cost, the model also incorporates minimizing carbon emissions as an optimization objective.

3. Problem Formulation

In the multiobjective route optimization model under time-varying conditions for hazardous chemicals transportation, G = (N, A) is an undirected road network, in which N represents the collection of nodes in the network, including the starting node O where the transportation enterprise is located, the collection of demand nodes NC and the collection of the general nodes NO. A represents the collection of road segments. Now, the transportation enterprise at O needs to deliver a certain kind of hazardous chemical to each customer node through the transportation network. The enterprise should be in the control of transportation cost, while minimizing the transportation risk. In addition, in order to achieve the goal of sustainable development, it is required to realize the minimization of carbon emissions. Under discrete time-varying conditions, the maximum speed, population density and accident rate of each road segment in the road network change dynamically over time, and the optimal transportation route of hazardous chemicals is planned according to the minimization objective.

3.1. Assumptions

(1)
In the transportation network, there is only one type of hazardous chemicals transportation vehicle. The maximum loading of each vehicle is the same, and the interaction between vehicles is not considered.
(2)
The vehicle leaves from the starting node and must return to it after providing services to different customer nodes.
(3)
Each customer node is limited by a soft time window, and when a vehicle arrives outside the time window specified for that customer node it can be unloaded normally, but there will be a penalty cost.
(4)
The demand of each customer node is indivisible and does not exceed the maximum loading of vehicle.
(5)
In order to reduce the amount of traffic in the network, the vehicle serves each customer node in order from smallest to largest demand, and the actual loading of the vehicle is the sum of the customer demand it needs to serve.
(6)
Each customer node has a fixed service time and can be converted to an ordinary node in the road network when the demand of customer is satisfied.
(7)
The maximum speed limit for transportation vehicles is different at each time period. Considering the transportation cost of the enterprise, transport vehicles travel at a constant speed, at the maximum speed limit, during each time period.
(8)
The population density around the road segment is divided into the population density on the road and the population density around the road, where the population density on the road changes dynamically over time.
(9)
The accident rate in transportation of hazardous chemicals changes dynamically over time.
(10)
The carbon emissions during transportation are related to the distance traveled and vehicle loading.

3.2. Symbols Definition

G: road network, G = (N, A)
N: node set, N =   { 0 ,   1 ,   2 ,   ,   n } , node 0 is the starting node.
N C : customer node set, N C N .
N O : general node set, N O N .
A: road segment set among nodes, A  = ( i ,   j ) :   i ,   j N .
K: vehicle set, K  = 0 ,   1 ,   2 ,   ,   k .
T: time period set, T  = t 1 ,   t 2 ,   t 2 ,   t 3 ,   t 3 ,   t 4 ,   ,   t n 1 ,   t n .
V: speed set, V  = v t 1 ,   t 2 , v t 2 ,   t 3 ,   v t 3 ,   t 4 ,   ,   v t n 1 ,   t n .
dij: the length of road segment   ( i ,   j ) , ( i ,   j ) A .
d t n 1 ,   t n : the distance that can be traveled in the time period t n 1 ,   t n .
fc: fixed cost per transport vehicle.
tc: cost per unit of distance traveled.
Q: maximum loading of the vehicle.
e: carbon emission factor per unit of distance during transportation.
q c : demand at customer node c, c     N C .
η c 1 : penalty factor when the vehicle at customer node c arrives earlier than the time window, c     N C .
η c 2 : penalty factor when the vehicle at customer node c arrives later than the time window, c     N C .
b t i k : departure time of vehicle k  ( k   K ) from node j  ( j   N ) .
a t j k : time for vehicle k  ( k   K ) to arrive at node i  ( i N ) .
e t c ,   l t c : time window at customer node c  ( c     N C ) .
s t c : service time at customer node c   ( c     N C ) .
l i j k : the loading of vehicle k  ( k   K ) when traveling through the road segment   ( i ,   j ) , ( i ,   j ) A .
v i j t n 1 ,     t n k : the speed of vehicle k  ( k   K ) on the segment   ( i ,   j ) in the time period t n 1 , t n , ( i ,   j ) A , t n 1 ,   t n   T .
ρ i j t n 1 ,   t n l : population density on the road segment ( i ,   j ) in the time period t n 1 , t n , ( i ,   j ) A , t n 1 ,   t n   T .
ρ i j o : population density around the road segment   ( i ,   j ) , ( i ,   j ) A .
p i j k : accident rate of hazardous chemicals transportation when vehicle k  ( k   K ) is traveling through the road segment   ( i ,   j ) , ( i ,   j ) A .
A R i j : the rate of traffic accidents on the road segment   ( i ,   j ) , ( i ,   j ) A .
P r i j : conditional release probability of hazardous materials after a traffic accident on the road segment   ( i ,   j ) , ( i ,   j ) A .
α : fixed parameters related to transport vehicle.
β : fixed parameters related to transport vehicle.
ε m : fuel consumption per unit of distance when the transport vehicle is fully loaded.
ε 0 : fuel consumption per unit of distance when the transport vehicle is empty.
r i j k : radius of impact after an accident when the vehicle k  ( k   K ) is traveling through the road segment   ( i ,   j ) , ( i ,   j ) A .
C c k : penalty cost of vehicle k  k   K at customer node c  ( c     N C )
P i j c k : consequence of an accident after the vehicle k  k   K passes through the road segment ( i ,   j ) at risk while traveling to the customer node c   ( c     N C ) , ( i ,   j ) A .
C i j c k : transportation cost of vehicle k  ( k   K ) passing through the road segment ( i ,   j ) on its way to customer node c  ( c     N C ) , ( i ,   j ) A .
R i j c k : transportation risk of vehicle k  ( k   K ) passing through the road segment ( i ,   j ) on its way to customer node c   ( c     N C ) , ( i ,   j ) A .
E i j c k : carbon emissions of vehicle k  ( k   K ) passing through the road segment i , j on its way to customer node c  ( c     N C ) , ( i ,   j ) A .
x i j c k : decision variable, x i j c k = 1 when vehicle k  ( k   K ) passes through the road segment ( i ,   j ) on its journey to customer c  ( c     N C ) ; else x i j c k = 0 , ( i ,   j ) A .
y c k : decision variable, y c k = 1 when customer node c  ( c     N C ) is served by vehicle k  ( k   K ) ; else y c k = 0 .

3.3. Mathematical Optimization Model

3.3.1. Transportation Cost

The cost involved in the transportation of hazardous chemicals is the combination of the fixed cost per vehicle dispatched, the cost of the vehicle during transportation, and the penalty cost that may be incurred at each customer node. In a time-varying road network, the time taken by a vehicle to arrive at node j is related to the departure time from the previous node i as well as the vehicle speed. When transport vehicle k departs from node i at b t i k through the road segment   ( i ,   j ) to arrive at node j, the arrival time a t j k is shown as below, where t n b t i k v i j t n 1 ,   t n k represents the distance the vehicle traveled from b t i k to t n , and x = 1 i 1 d t n +   x ,   t n + x + 1 represents the accumulation of the distance traveled by vehicles during the time period t n ,   t n + i , and so the numerator of the formula represents the total distance traveled by vehicles during the time period t n + i ,   t n + i   + 1 :
a t j k = t n + i + d i j t n b t i k v i j t n 1 ,   t n k x = 0 i 1 d t n +   x ,   t n + x + 1 v i j t n + i ,   t n + i + 1 k ,   b t i k t n 1 ,   t n ,   a t j k t n + i ,   t n + i   + 1
then the time for vehicle k to arrive at any demand node c from node m at b t m k through the road segment ( m ,   c ) is
a t c k = t n + i + d m c t n b t m k v m c t n 1 , t n k x = 0 i 1 d t n + x , t n + x + 1 v m c t n + i , t n + i + 1 k ,   b t m k   t n 1 , t n ,   a t c k   t n + i , t n + i + 1
Each customer node c has a soft time window, and the vehicle can arrive at the customer node at any time and complete normal unloading. When the time of arrival at the customer node is not within the time window of the customer node, the penalty cost is incurred for arriving earlier or later than the time window. The penalty cost incurred by vehicle k at customer node c is
C c k = η c 1 max e t c a t c k ,   0   +   η c 2 max a t c k l t c ,   0
and the total cost of vehicle k on the segment ( i ,   j ) can be composed of the fixed dispatch cost of the vehicle, the driving cost of the vehicle, and the penalty cost incurred at the customer nodes:
C i j c k = y c k f c +   y c k x i j c k d i j t c + η c 1 max e t c a t c k ,   0 + η c 2 max a t c k l t c ,   0

3.3.2. Transportation Risk

The risk of hazardous chemicals transportation is usually expressed as the product of the probability of the accident and the consequences of the accident [36]. The probability of hazardous chemicals road transportation accidents varies dynamically with the loading of the transport vehicle as well as with different road segments and times [15,37,38]. In Equation (5), α and β are constant values that depend on the type of hazardous chemicals [39]:
p i j k = A R i j P r i j d i j α l i j k β
The population exposure model is used to measure the consequences caused by the occurrence of transportation accidents [40]; the population density is divided into the population density on the road where the accident occurs and the population density around the road. Figure 1 shows the population exposure area after a traffic accident occurs. The circular area in the middle represents the road area directly affected by the accident, and the remaining area is the surrounding area. The population density of the road segment where the accident occurred varies dynamically with different time periods, and the consequences of road transportation accidents involving hazardous chemicals can be expressed as the number of people in the exposed area:
P i j c k = ρ i j o π r i j k 2 + 2 π r i j k d i j + ρ i j t n 1 , t n l π r i j k 2
The transportation risk incurred by vehicle k  ( k   K ) on road segment ( i ,   j ) can be expressed as the product of p i j k and P i j c   k :
R i j c k = y c k x i j c k p i j k P i j c   k = y c k x i j c k A R i j P r i j d i j α l i j k β ρ i j o π r i j k 2 + 2 π r i j k d i j + ρ i j t n 1 , t n l π r i j k 2 .

3.3.3. Carbon Emissions

The carbon emissions of vehicle k  ( k   K ) on road segment ( i ,   j ) are related to the transportation distance d i j and the vehicle loading l i j k . According to [22,41], the calculation formula of carbon emissions is
E i j c k =   x i j c k d i j e ε m ε 0 Q l i j k + ε 0

3.3.4. Mathematical Model

In summary, the multiobjective road transportation route optimization model for hazardous chemicals under time-varying conditions is formulated as follows:
min   f 1 = k   K c   N c i ,   j   A { y c k f c +   y c k x i j c k d i j t c + η c 1 max e t c a t c k ,   0 + η c 2 max a t c k l t c ,   0 }
min   f 2 = k     K c     N c i ,   j   A y c k x i j c k A R i j P r i j d i j α l i j k β ρ i j o π r i j k 2 + 2 π r i j k d i j + ρ i j t n 1 , t n l π r i j k 2
min   f 3 = k     K c     N c i ,   j   A x i j c k d i j e ε m ε 0 Q l i j k + ε 0
s.t.
k K c   N c i N x i j c k = 1 ,   j   N O
k     K c     N c j     N x i j c k = 1 ,   i     N O
j     N x 0 j c k = i     N x i 0 c k = 1 ,     k     K ,   c     N C
i     N j   N x i j c k = i     N j     N x j i c k ,   k     K ,   c   N C
j   N l j c k j N l c j k = q c ,     k     K ,   c     N C
i     N x i j c k c   N c q c   l i c k ,     k     K ,   c   N C
c   N c y c k q c     Q ,     k     K
x i j c k 0 ,   1 ,     k     K ,   c   N C ,   ( i ,   j )   A
y c k 0 ,   1 ,     k     K ,   c   N C
y c k = 0 , i     N j     N x i j c k = 0 1 , i     N j     N x i j c k = 1 , k     K , c N C
In this multiobjective nonlinear integer optimization model, Equations (9)–(11) are objective functions that represent the minimization of transportation costs of hazardous chemicals, the minimization of transportation risk, and the minimization of carbon emissions, respectively. Equations (12)–(20) are constraints of the model. Equations (12) and (13) indicate that each customer node must and can be served only once. Equation (14) indicates that each transport vehicle must depart from the starting node and return to the starting node after serving each customer node. Equation (15) provides the vehicle flow balance constraint for each node in the transportation network. Equation (17) indicates that the loading of the vehicle must meet the demand of the customer node to be served. Equation (18) indicates that the loading of each vehicle must not exceed the maximum loading of the transport vehicle. Equations (19) and (20) indicate that the decision variables are 0–1 variables. Equation (21) indicates the interrelation between the decision variables.

4. Solution Methodology

4.1. NSGA-II Algorithm

The multiobjective hazardous chemicals road transportation route optimization model under time-varying conditions belongs to the NP-Hard problem, which is not applicable to large-scale arithmetic cases using the exact algorithm. There may be conflicts between the optimization objectives, which make it impossible to find a deterministic optimal solution. The NSGA-II, as a heuristic algorithm, has been shown to have better results in solving multiobjective optimization problems [42,43], which utilizes the computation of nondominated ordering and congestion distance for population optimization and uses an elite strategy to select the optimal solution generated during the iteration process. In this paper, NSGA-II is improved to solve the multiobjective nonlinear optimization model.

4.2. Encoding

In the optimization model proposed in this paper, firstly the number of transport vehicles needs to be determined according to the demand of each customer node, and then the transportation route planning is carried out according to the customer nodes to be served by each transport vehicle. A chromosome needs to form multiple transportation routes based on the number of required transport vehicles, so the chromosome adopts priority integer coding. Each integer gene represents the priority of the corresponding node, and the length of the chromosome is determined by the number of nodes in the road network. The sequence of integers starting from 0 is randomly generated. Take the transportation road network shown in Figure 2 as an example; there are eight nodes in this road network, which contains the starting node A. Then the chromosome is an integer list from 0 to 7.

4.3. Decoding

Decoding is the process of mapping the generated integer chromosomes into multiple hazardous chemicals road transportation routes. In this paper, the decoding process is divided into two parts. The first part is the allocation of transport vehicles, which arranges transport vehicles for each customer node according to its demand and defines key nodes for each vehicle. The key nodes are some nodes that this transport vehicle must pass through during traveling, including the customer nodes it serves and the starting node to which it will eventually return. The second part is the generation of a transportation route, which consists of decoding the chromosome to generate a transportation route based on the key nodes. The specific decoding method is as follows:
Part 1—Transport vehicle allocation
Step 1.1: Obtain the demand of each customer node in the road network.
Step 1.2: Accumulate the demand of each customer node sequentially from the first customer node until it meets the demand of customer node that does not exceed the maximum loading of the transport vehicle. Allocate the same transport vehicle to these customer nodes and generate a list of key nodes for this vehicle.
Step 1.3: Repeat the above process until all the demand points can be served.
Part 2—Transportation route generation
Step 2.1: Mark the key nodes of the transport vehicle in the road network.
Step 2.2: Extract the adjacency node matrix of the node where the transport vehicle is currently located. Determine whether the next key node is in the adjacency matrix of the current node; if so, select the next key node to join the transportation route, update the node where it is currently located, and go to Step 2.4; otherwise, go to the Step 2.3.
Step 2.3: Find out the priority of each node in the adjacency node matrix of the current node according to the chromosome encoding and select the node with the highest priority. Determine whether it already exists in the current transportation route, and the node with the highest priority that is not in the current transportation route joins the transportation route; then update the node where it is currently located and go to the Step 2.4. If there is no node that meets the conditions, then determine that this route is an unfeasible route and the current chromosome is invalid.
Step 2.4: Judge whether the current node is the last key node. If so, end this route planning process and turn to Step 2.5; otherwise turn to Step 2.3.
Step 2.5: Judge whether all transportation vehicles have completed the route planning. If so, end the chromosome decoding process and output all transportation routes corresponding to this chromosome; otherwise turn to Step 2.1.
Using the road network in Figure 2 as an example, a chromosome is randomly generated: 8–3–4–7–6–1–2–5. If node G is a customer node, the list of key nodes during vehicle transportation is [A, G, A]. The adjacency node matrix of node A is [B, C], where node C has higher priority, so the next route node is chosen to be point C. The adjacency node matrix of node C is [A, B, D, E], because nodes A and B are already in the departure route, so the next route node selects point E. The adjacency node matrix of node E is [C, F, H], where node H has higher priority and is not in the current route, so node H is selected as the next node. The adjacency node matrix of node H is [F, E, G], which contains the customer node, so it goes directly to node G. After serving the customer node at node G, it needs to return to the starting node, and the adjacency node matrix of node G is [D, H], of which node D is of higher priority, so the next node of the route is D. Among the adjacent nodes of node D, node C has a higher priority, and the starting node A is in the adjacent nodes of node C. Therefore, the complete transportation route is A → C → E → H → G → D → C → A.

4.4. Crossover Operator

In this paper, the sequential crossover operator is used. As shown in Figure 3, taking the generation of subindividual C1 as an example. On two randomly selected individuals, two crossover points are randomly selected to form the gene fragments, and the gene fragments in P1 are copied to the corresponding positions in the subindividual C1. Scanning each gene fragment of the individual P2 from left to right, the gene fragments that are not in the subindividual C1 are filled into the empty positions in the individual C1 in order to generate the subindividual C1. The same steps are used to generate another subindividual C2.

4.5. Mutation Operator

The two kinds of mutation operators used in this paper are two-point mutation and sequential mutation. As shown in Figure 4, two crossover points are randomly selected on individual P, and the gene fragments at the corresponding positions are exchanged to form a new individual P’.
As shown in Figure 5, a gene fragment on individual P is intercepted, and then the gene fragment is arranged in reverse order to form a new gene fragment, which is then placed in the corresponding position on the original chromosome to form a new individual P”.

4.6. Repair Operator

During the process of decoding, unfeasible routes are generated, resulting in invalid chromosomes. In order to reduce the number of invalid chromosomes, this paper introduces the repair operator. When a chromosome is judged as an invalid chromosome, firstly, a repair probability is randomly generated for each gene fragment on the chromosome and a repair threshold is set. Compare the repair probability of each gene fragment with the repair threshold from left to right, and if the repair probability is greater than the repair threshold, extract the node of the road network corresponding to the gene fragment as well as the adjacency matrix. Then, generate a new node priority by arranging the priorities of nodes in the node’s adjacency matrix in inverse order. Replace the gene fragments in the existing chromosome with the new priorities, until all the gene fragments in the chromosome are traversed. A repair threshold of 0.5 is selected as shown in Figure 6.

4.7. Termination Condition

In this paper, the iteration of the algorithm is terminated by setting the maximum number of iterations. When the number of iterations of the algorithm reaches the preset iteration threshold, the algorithm is aborted and the Pareto optimal solution is outputted; otherwise, the algorithm is continued.

4.8. The Improved NSGA-II Algorithm Flow

The specific flow of the improved NSGA-II algorithm is as follows, and the flowchart is shown in Figure 7.
Step 1: Initialize the population and set the maximum number of iterations.
Step 2: Determine if it is an invalid chromosome. If it is, perform chromosome repair; otherwise, proceed to the next step.
Step 3: Determine whether the first subpopulation has been generated, and if so, proceed to the next step; otherwise, perform fast nondominated sorting on the initial population and generate offspring population through elite selection, crossover, and mutation operations.
Step 4: Update the number of iterations.
Step 5: Combine the parent population and the offspring population.
Step 6: Determine if a new parent population has been generated, and if so, proceed to the next step; otherwise, calculate the fitness function value of individuals in the population, perform fast nondominated sorting, calculate crowding distance, and obtain a new parent population through elite selection, crossover, and mutation operations.
Step 7: Determine if it is an invalid chromosome. If it is, perform chromosome repair; otherwise, proceed to the next step.
Step 8: Determine whether the preset number of termination iterations is reached, and end the algorithm operation; otherwise, return to Step 4.

5. Case Analysis

5.1. Road Network

In this section, Sioux Falls network [44] is used to verify the feasibility and effectiveness of the time-varying model and the improved NSGA-II algorithm. Sioux Falls network was derived from Sioux Falls, the largest city in South Dakota, USA, and this network has been widely used in route optimization problems. As shown in Figure 8, the network has 24 nodes and 38 road segments. Based on the optimization model of the hazardous chemicals transportation network under time-varying conditions proposed, the length of each road segment, the population density, hazardous chemicals leakage rate, and accident rate at different time periods are randomly generated, as shown in Table 1. Node 1 is the starting node of the hazardous chemicals transport vehicle, and the customer nodes, as well as the demands, time windows, and penalty cost, are randomly generated, as shown in Table 2. The day is divided into seven time periods based on light visibility as well as road traffic density, as shown in Table 3 [45]. The other parameters involved in this calculation example are set as shown in Table 4.
The experimental running environment consists of an Intel i5 processor, 3.20 GHz main frequency, 16 GB RAM, and the Windows 11 operating system, and the improved NSGA-II algorithm was tested using the Python 3.10.0 program on the VS Code platform. The parameter settings for the algorithm operation, according to the problem size and the analysis, are as follows: population of 200, maximum number of iterations of 100, crossover probability of 0.6 and variance probability of 0.8. For each experiment, run ten times and use the optimal results.

5.2. Experiment 1: Transportation Route Optimisation under Different Departure Times

In the time-varying model of this paper, the traveling speed of the vehicle and the population density around the road segment are affected by time and each customer node in the road network has a soft time window, so we firstly explore the time-varying characteristics of hazardous chemicals road transportation. On the basis of dividing a day into different discrete time periods, different vehicle departure times can change the time period of the entire transportation process. The experimental design randomly generates seven vehicle departure times in seven time periods and performs optimal route planning at each departure time. The Pareto solutions under each department time are shown in Figure 9, and the optimal routes responding to the Pareto solutions, as well as the objective function values, are shown in Table 5.
As can be seen from Figure 9, running the improved NSGA-II algorithm proposed in this paper at different departure times can obtain the corresponding multiple Pareto solutions. And there may be overlap between multiple different Pareto optimal solutions, which suggests that the proposed algorithm can obtain the Pareto optimal solutions in many iterations, and the verification demonstrates the feasibility and validity of the improved NSGA-II algorithm in this paper.
Solving the multiobjective optimization model under each departure time obtains multiple sets of Pareto optimal solutions and randomly intercepts three optimal solutions in each group, as shown in Table 5. Three different customer nodes (14, 17, and 18) are set in the experiment and the sum of the overall demand of customer node 14 and customer node 17 does not exceed the maximum loading of the transport vehicle, so customer node 14 and customer node 17 are transported by one transport vehicle. Each Pareto optimal solution obtained by the optimization experiment when the customer node is the same at different departure times contains the optimization route of two different transport vehicles.
As can be seen from the experimental results in Table 5, vehicles also have different objective function values when transporting hazardous chemicals along the same transportation routes at different departure times. Taking the routes r 1 1 and r 2 2 as examples, at the different departure times of 4:20 and 7:20, the transportation routes are both [1, 3, 4, 11, 14, 15, 10, 17, 10, 11, 4, 3, 1] and [1, 3, 4, 11, 10, 16, 18, 16, 8, 6, 2, 1], but in terms of transportation cost, the transportation cost of r 1 1 is CNY 2646.93 higher than r 2 2 and the transportation cost of r 2 2 is CNY 2592.35. In terms of transportation risk, the transportation risk of r 1 1 is 20.47, which is higher than the transportation risk of r 2 2 . The carbon emissions when departing at 7:20 are 226.63 kg, while the carbon emissions when departing at 7:20 were slightly lower at 240.55 kg. When the departure time of the vehicle at node 1 is 4:20, the time interval from the time window of each customer node is longer than that of the departure time of 7:20, resulting in an increase in the penalty cost at the earlier departure time and thus an increase in the total transportation cost. When the driving speed during transportation is low, the overall transportation time of the vehicle increases, resulting in an increase in carbon emissions. In addition, at different departure times, speed limits on different road segments and the surrounding population density will affect the transportation risk. When the speed is low, the travel time increases and the time it takes for the transporter to reach the customer node will vary, resulting in different penalty costs at each customer node.
In Table 5, when the departure time is 7:20, there are two optimal transportation solutions r 1 2 and r 2 2 . In these two transportation route plans, transport vehicle 2 has the same driving path, while transport vehicle 1 has different driving paths in the road network, resulting in differences in transportation cost, transportation risk, and carbon emissions between the two route schemes. Specifically, the vehicles in transportation path r 2 2 have a shorter travel distance, and the difference in travel distance between the two route schemes is mainly reflected in the transportation distance when returning to the starting node after serving the customer node. Therefore, the longer travel distance results in higher transportation cost and carbon emissions for route r 1 2 . In terms of transportation risk, by comparing the probability of conditional leakage, accident occurrence rate, and population density of each road segment, it can be seen that the various indicators of transportation path r 2 2 are relatively high, so the transportation risk of transportation path r 2 2 is relatively high.
In conclusion, in the road transportation of hazardous chemicals, time factors affect the optimization of transportation route by affecting the driving speed and accident rate of transport vehicles. If transportation enterprises want to minimize cost during transportation, they need to choose a departure time closer to the customer’s time window to minimize the penalty cost incurred at the customer node. If transportation enterprises want to minimize transportation risk, they should choose roads with lower accident rates among multiple transportation route schemes.

5.3. Experiment 2: Transportation Route Optimization with Customer Nodes Changing

In order to verify the scalability of the time-varying model constructed and explore the impact of changes in customer nodes on the optimization of hazardous chemical road transportation routes, on the basis of the existing customer nodes (14, 17, and 18), the customer node 10 and customer node 22 shown in the Table 6 were added, and the departure time was fixed at 10:30 to carry out the route optimization experiment.
Figure 10 shows the Pareto optimal solution obtained by optimization experiments under the condition that the number of customer nodes in the road network increases at the same departure time, and Figure 10a–c are the Pareto optimal solutions at three customer nodes, four customer nodes and five customer nodes, respectively. As can be seen from the results in Figure 10, the number of Pareto optimal solutions gradually increases as the number of customer nodes increases. This indicates that the improved algorithm proposed in this paper has good scalability, which can obtain good solutions under the condition of changing numbers of customer nodes, and meet the requirements of the optimization experiment.
Table 7 shows each optimized route and the corresponding cost, risk and carbon emissions obtained by performing the optimization experiments at different customer nodes at the same departure time. As can be seen from the results in Table 7, transportation cost, transportation risk and carbon emissions rise as the number of customer nodes in the road network gradually increases. When the customer node is (14, 17, or 18) the transportation cost, transportation risk, and carbon emissions are between [2000, 3000], [10, 20] and [200, 300], respectively. When the customer node is (14, 10, 17, or 18), the transportation cost, transportation risk and carbon emissions are between [3000, 4000], [15, 30] and [300, 500], respectively. After adding node 22 into the network as the customer node, the objective function is increased to between [5000, 6000], [20, 35] and [400, 600], respectively. This is because the increase in customer nodes will first lead to an increase in the overall carrying capacity of transporters and an increase in the number of transport vehicles, which in turn leads to an increase in the transportation cost. Since the risk during transportation is directly related to the vehicle loading, when the number of customer nodes increases and the distance traveled by the vehicle increases, so the transportation risk also increases.

5.4. Green Transportation

When the customer node in the network is (14, 10, 17, 18, or 22), excluding carbon emissions, the transportation cost of transportation route r 2 10 is the lowest at CNY 4882.80, the transportation cost of transportation route r 8 10   is CNY 5823.06, while the transportation risk of transportation route r 1 10 is 23.58, and the transportation risk of transportation route r 2 10 is 30.69. Figure 11 shows the relationship between transportation cost and carbon emissions. It can be seen that in the road transportation route optimization problem of hazardous chemicals, there is often a contradictory relationship between transportation cost and transportation risk, and the reduction of transportation cost often means the increase of transportation risk. Therefore, the hazardous chemical transportation enterprises need to weigh the relationship between the two in the actual path planning. In the case that the transportation cost and transportation risk cannot be optimal at the same time, the optimal road transportation route should be selected combined with enterprises’ risk tolerance and economic benefits.
Figure 12 shows that carbon emissions will increase when transportation cost gradually increases, and there is a positive correlation between transportation cost and carbon emissions, indicating that when the transportation cost of vehicles gradually increase, it is often accompanied by an increase in carbon emissions. For example, the transportation cost of route r 2 8 in Table 7 is higher than r 4 8 , but the carbon emissions are lower, indicating that the enterprise can indirectly reduce carbon dioxide emissions by reducing vehicle travel distance, setting reasonable departure times, and other means to reduce transportation cost, and promote the green development of hazardous chemicals transportation.

6. Conclusions and Future Work

Road transportation of hazardous chemicals is characterized by high risk, so the risk should be minimized during the transportation process. This paper innovatively proposes that from the perspective of transportation route schemes, enterprises should optimize transportation routes based on the time windows of customer nodes and dynamic road conditions by studying the time-varying factors in hazardous chemical road transportation. At the same time, by studying the interrelationships between optimization objectives, it is believed that the reduction of transportation cost should be the primary consideration during transportation. The specific research content of this paper is as follows:
(1)
A multiobjective route optimization model is constructed to minimize transportation cost, transportation risk, and carbon emissions based on the risk of hazardous chemicals road transportation, the economy of transportation enterprises, and the “double carbon” goal. At the same time, the time factor is added into the model, and a soft time window is set for each customer node. The day is divided into different time periods. The speed limit of the road as well as the population density and accident rate are different during each time period.
(2)
The NSGA-II algorithm is used for solving the problem, and the NSGA-II algorithm is improved to reduce the occurrence of a large number of invalid solutions during the optimization problem solving process. Based on the model proposed in this paper, road network optimization experiments in classical Sioux Falls networks are performed to verify the feasibility and effectiveness of the algorithm.
(3)
The experimental results of route optimization show that the transportation routes obtained at different departure times are very different. When transport vehicles depart at different times, the same transportation route will generate different transportation costs, transportation risks, and carbon emissions. In addition, as the number of customer nodes increases, transportation cost, transportation risk, and carbon emissions during transportation will correspondingly increase.
(4)
When optimizing the transportation route, enterprises need to reasonably plan the department time and transportation route according to the time-varying factors If enterprises want to minimize transportation cost, it needs to choose a departure time that is closer to the customer’s time window and a transportation route that is shorter. If enterprises pay more attention to safety, they need to choose transportation routes with lower accident rates and population density based on road conditions.
(5)
In the process of optimizing the road transportation route of hazardous chemicals, from the perspective of transportation enterprises, transportation cost is the primary factor to consider. Based on the negative correlation between cost and risk, as well as the positive correlation between cost and carbon emissions, multiple measures should be taken to reduce transportation cost, and then machine learning, artificial intelligence, and other technologies should be applied to construct a more reasonable and effective optimization model for minimizing transportation risk. Firstly, it is necessary to improve the construction of information infrastructure for road transportation of hazardous chemicals, construct a big data system for hazardous chemical transportation, achieve the circulation of transportation big data, and reduce the cost of data acquisition for transportation enterprises. Secondly, through the application of in-vehicle sensors, the IoT, and video recognition detection technology, intelligent supervision of transportation can be achieved, reducing the cost of supervision during transportation. Finally, the government can allocate financial subsidies to small and medium-sized transportation enterprises to reduce their cost burden.
Road transportation of hazardous chemicals is a process that is affected by time, weather, road conditions, and other factors, and the goals pursued by each participant in the transportation process are also in conflict. Firstly, this paper conducts route optimization experiments based on discrete time-varying conditions, and preliminarily concludes that the departure time of vehicles and the time window of customer nodes will have an impact on the cost and risk of the transportation. However, there are still shortcomings in specific time path planning, and more comprehensive experiments are needed to explore the impact mechanism of changes in time factors on transportation. Secondly, the time factor is included in the transportation route optimization model of hazardous chemicals, but the road information is known with certainty. In the future, the uncertain road transportation optimization problem based on time-varying factors can be considered to study the influence of uncertain external environmental conditions on the transportation route. Thirdly, the interaction between vehicles and the overall riskiness of the road network when multiple transport vehicles are traveling in the road network at the same time are not considered, and further research is needed on the interaction between vehicles and the overall risk and fairness of the transportation network. Meanwhile, the model and algorithm in the paper also require a larger network for testing.

Author Contributions

Conceptualization, Z.Z. and S.K.; methodology, S.K.; software, S.K.; validation, S.K. and Z.Z.; formal analysis, S.K.; investigation, S.K.; resources, Z.Z.; data curation, S.K.; writing—original draft preparation, S.K.; writing—review and editing, S.K.; visualization, S.K.; supervision, Z.Z.; project administration, S.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this published article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zhao, G.; Li, H. Economic operation report of China’s petroleum and chemical industry in 2021. Mod. Ind. 2022, 42, 249–251. [Google Scholar]
  2. Zhao, J.; Yuan, J.; Zhao, R.; Huang, H.; Su, B. Multi-Agent emergency cooperative dispatching system for dangerous chemicals leakage on highways. China Saf. Sci. J. 2020, 30, 116–121. [Google Scholar]
  3. Yao, J.; Xie, B.; Wu, X.; Zhang, C. Two-level programming model based on cooperative operation study of stakeholders in hazardous chemical storage. Sustainability 2023, 15, 1221. [Google Scholar] [CrossRef]
  4. Li, Y.; Wang, Y.; Lai, Y.; Shuai, J.; Zhang, L. Monte Carlo-based quantitative risk assessment of parking areas for vehicles carrying hazardous chemicals. Reliability Eng. Syst. Saf. 2023, 231, 109010. [Google Scholar] [CrossRef]
  5. Hou, S.; Wang, Z.; Luan, X.; Liu, Y.; Gong, Y.; Guo, C.; Wang, S.; Zhang, B. Analysis of tank truck explosion accident in Wenling. Zhejiang Prov. J. Nanjing Univ. Technol. 2021, 43, 144–149. [Google Scholar]
  6. Dantzig, G.B.; John, H.R. The truck dispatching problem. Manag. Sci. 1959, 6, 80–91. [Google Scholar] [CrossRef]
  7. Desrochers, M.; Jacques, D.; Marius, S. A new optimization algorithm for the vehicle routing problem with time windows. Oper. Res. 1992, 40, 342–354. [Google Scholar] [CrossRef]
  8. Thangiah, S.R.; Nygard, K.E.; Juell, P.L. GIDEON: A genetic algorithm system for vehicle routing with time windows. 1991. Proceedings. In Proceedings of the 7th IEEE Conference on Artificial Intelligence Application, Miami Beach, FL, USA, 24–28 February 1991; pp. 322–328. [Google Scholar]
  9. Wren, A.; Holliday, A. Computer scheduling of vehicles from one or more depots to a number of delivery points. J. Oper. Res. Soc. 1972, 23, 333–344. [Google Scholar] [CrossRef]
  10. Golden, B.L.; Magnanti, T.L.; Nguyen, H.Q. Implementing vehicle routing algorithms. Networks 1977, 7, 113–148. [Google Scholar] [CrossRef]
  11. Androutsopoulos, K.N.; Zografos, K.G. A bi-objective time-dependent vehicle routing and scheduling problem for hazardous materials distribution. EURO J. Transp. Logist. 2012, 1, 157–183. [Google Scholar] [CrossRef]
  12. Zou, Z.; Zhang, B. Route optimization of hazardous chemicals transportation with mixed time windows. China Saf. Sci. J. 2012, 22, 83–89. [Google Scholar]
  13. Yuan, W.; Wang, J.; Wu, J.; Li, J. A new model of hazardous material vehicle routing problem and its algorithm. J. Syst. Sci. Math. Sci. 2017, 37, 393–406. [Google Scholar]
  14. Bula, G.A.; Afsar, H.M.; Gonzalez, F.A.; Prodhon, C.; Velasco, N. Bi-objective vehicle routing problem for hazardous materials transportation. J. Clean. Prod. 2019, 206, 976–986. [Google Scholar] [CrossRef]
  15. Jiang, P.; Men, J.; Xu, H.; Zheng, S.; Kong, Y.; Zhang, L. A variable neighborhood search-based hybrid multiobjective evolutionary algorithm for HazMat heterogeneous vehicle routing problem with time windows. IEEE Syst. J. 2020, 14, 4344–4355. [Google Scholar] [CrossRef]
  16. Chai, H.; He, R.; Kang, R.; Jia, X.; Dai, C. Solving bi-objective vehicle routing problems with driving risk consideration for hazardous materials transportation. Sustainability 2023, 15, 7619. [Google Scholar] [CrossRef]
  17. Esmaeili-Douki, A.; Mahzouni-Sani, M.; Jahromi, A.N.; Jolai, F. A novel fuzzy bi-objective vehicle routing and scheduling problem with time window constraint for a distribution system: A case study. Sci. Iran. 2021, 28, 2868–2889. [Google Scholar] [CrossRef]
  18. Rahbari, M.; Khamseh, A.A.; Sadati-Keneti, Y.; Jafari, M.J. A risk-based green location-inventory-routing problem for hazardous materials: NSGA II, MOSA, and multi-objective black widow optimization. Environ. Dev. Sustain. 2022, 24, 2804–2840. [Google Scholar] [CrossRef]
  19. Zhao, L.; Cao, N. Fuzzy random chance-constrained programming model for the vehicle routing problem of hazardous materials transportation. Symmetry 2020, 12, 1208. [Google Scholar] [CrossRef]
  20. Wang, X.; Dai, C.; Li, J.; Wu, R. Optimization of green distribution routes for hazardous materials under uncertain transportation risk. J. Comput. Eng. Appl. 2023, 59, 323–332. [Google Scholar]
  21. Dai, C.; Li, Y.; Ma, C.; Chai, H.; Mu, H. Multi-criteria optimization for hazardous materials distribution routes under uncertain conditions. J. Jilin Univ. 2018, 48, 1694–1702. [Google Scholar]
  22. Lyu, J.; He, Y. A two-stage hybrid metaheuristic for a low-carbon vehicle routing problem in hazardous chemicals road transportation. Appl. Sci. 2021, 11, 4864. [Google Scholar] [CrossRef]
  23. Sun, Y.; Li, X.; Liang, X.; Zhang, C. A bi-objective fuzzy credibilistic chance-constrained programming approach for the hazardous materials road-rail multimodal routing problem under uncertainty and sustainability. Sustainability 2019, 11, 2577. [Google Scholar] [CrossRef]
  24. Meng, X.; Zhang, K.; Pang, K.; Xiang, X. Characterization of spatio-temporal distribution of vehicle emissions using web-based real-time traffic data. Sci. Total Environ. 2020, 709, 136227. [Google Scholar] [CrossRef] [PubMed]
  25. Xiao, G.; Chen, L.; Chen, X.; Jiang, C.; Ni, A.; Zhang, C.; Zong, F. A hybrid visualization model for knowledge mapping: Scientometrics, SAOM, and SAO. IEEE Trans. Intell. Transp. Syst. 2023, 1–14. [Google Scholar] [CrossRef]
  26. Zhang, Y.; Lu, L.; Liu, Q.; Hu, M. Modeling of low-risk behavior of pedestrian movement based on dynamic data analysis. Transp. Res. Part A 2023, 168, 103576. [Google Scholar] [CrossRef]
  27. Li, J.; Guo, F.; Zhou, Y.; Yang, W.; Ni, D. Predicting the severity of traffic accidents on mountain freeways with dynamic traffic and weather data. Transp. Saf. Environ. 2023, 5, tdad001. [Google Scholar] [CrossRef]
  28. Pourmoradnasseri, M.; Khoshkhah, K.; Hadachi, A. Leveraging IoT data stream for near-real-time calibration of city-scale microscopic traffic simulation. IET Smart Cities 2023, 5, 269–290. [Google Scholar] [CrossRef]
  29. Chen, X.; Wang, Z.; Hua, Q.; Shang, W.; Luo, Q.; Yu, K. AI-empowered speed extraction via port-like videos for vehicular trajectory analysis. IEEE Trans. Intell. Transp. Syst. 2023, 24, 4541–4552. [Google Scholar] [CrossRef]
  30. Liu, Y.; Zhu, X. Multi-objective optimization of hazardous materials transportation under double uncertainty conditions. Comput. Integr. Manuf. Syst. 2020, 26, 1130–1141. [Google Scholar]
  31. Ke, G.Y. Managing reliable emergency logistics for hazardous materials: A two-stage robust optimization approach. Comput. Oper. Res. 2022, 138, 105557. [Google Scholar] [CrossRef]
  32. Ouertani, N.; Ben-Romdhane, H.; Krichen, S. A decision support system for the dynamic hazardous materials vehicle routing problem. Oper. Res. 2022, 22, 551–576. [Google Scholar] [CrossRef]
  33. Xu, W.; Liang, J.; Bian, W.; Dai, B.; Tao, G.; Liu, C. R Spontaneous path selection for hazardous chemical transportation based on quasi-decision tree pruning. CIESC J. 2018, 69, 324–329. [Google Scholar]
  34. Xu, W.; Bian, W.; Wang, W.; Liu, C.; Zhuang, J. Spontaneous path selection for hazardous chemical transportation based on real-time traffic condition. CIESC J. 2018, 69, 1136–1140. [Google Scholar]
  35. He, Z.; Zhang, C.; Fang, J. Variable path-ways selection for the transportation of hazardous goods based on the conditional value. J. Saf. Environ. 2018, 18, 428–433. [Google Scholar]
  36. Russo, F.; Vitetta, A. Risk evaluation in a transportation system. Int. J. Sustain. Dev. Plan. 2006, 1, 170–191. [Google Scholar] [CrossRef]
  37. Men, J.; Chen, G.; Zhou, L.; Chen, P. A pareto-based multi-objective network design approach for mitigating the risk of hazardous materials transportation. Process Saf. Environ. Prot. 2022, 161, 860–875. [Google Scholar] [CrossRef]
  38. Guo, J.; Luo, C. Risk assessment of hazardous materials transportation: A review of research progress in the last thirty years. J. Traffic Transp. Eng. 2022, 9, 571–590. [Google Scholar] [CrossRef]
  39. Ronza, A.; Vílchez, J.A.; Casal, J. Using transportation accident databases to investigate ignition and explosion probabilities of flammable spills. J. Hazard. Mater. 2007, 146, 106–123. [Google Scholar] [CrossRef]
  40. Erkut, E.; Ingolfsson, A. Transport risk models for hazardous materials: Revisited. Oper. Res. Lett. 2005, 33, 81–89. [Google Scholar] [CrossRef]
  41. Zhang, L.-Y.; Tseng, M.-L.; Wang, C.-H.; Xiao, C.; Fei, T. Low-carbon cold chain logistics using ribonucleic acid-ant colony optimization algorithm. J. Clean. Prod. 2019, 233, 169–180. [Google Scholar] [CrossRef]
  42. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef]
  43. Luo, Q.; Fan, Q.; Deng, Q.; Guo, X.; Gong, G.; Liu, X. Solving bi-objective integrated scheduling problem of production, inventory and distribution using a modified NSGA-II. Expert Syst. Appl. 2023, 225, 120074. [Google Scholar] [CrossRef]
  44. Zahiri, B.; Suresh, N.C. Hub network design for hazardous-materials transportation under uncertainty. Transp. Res. Part E-Logist. Transp. Rev. 2021, 152, 102424. [Google Scholar] [CrossRef]
  45. Lu, Y.; Wu, J.; Shao, S.; Shi, S.; Zhou, R.; Wang, W. Prediction model for road transport accidents of hazardous chemicals based on Bayesian network. China Saf. Sci. J. 2022, 32, 174–182. [Google Scholar]
Figure 1. Population exposure area after the accident.
Figure 1. Population exposure area after the accident.
Sustainability 16 00779 g001
Figure 2. Example of a transportation network.
Figure 2. Example of a transportation network.
Sustainability 16 00779 g002
Figure 3. Example of crossover operator.
Figure 3. Example of crossover operator.
Sustainability 16 00779 g003
Figure 4. Example of two-point mutation.
Figure 4. Example of two-point mutation.
Sustainability 16 00779 g004
Figure 5. Example of sequential mutation.
Figure 5. Example of sequential mutation.
Sustainability 16 00779 g005
Figure 6. Example of repair operator.
Figure 6. Example of repair operator.
Sustainability 16 00779 g006
Figure 7. The flowchart of the improved NSGA-II algorithm.
Figure 7. The flowchart of the improved NSGA-II algorithm.
Sustainability 16 00779 g007
Figure 8. Sioux Falls network.
Figure 8. Sioux Falls network.
Sustainability 16 00779 g008
Figure 9. The Pareto fronts of Experiment 1.
Figure 9. The Pareto fronts of Experiment 1.
Sustainability 16 00779 g009
Figure 10. The Pareto optimal solution: (a) optimization results for three customer nodes; (b) optimization results for four customer nodes; and (c) optimization results for five customer nodes.
Figure 10. The Pareto optimal solution: (a) optimization results for three customer nodes; (b) optimization results for four customer nodes; and (c) optimization results for five customer nodes.
Sustainability 16 00779 g010
Figure 11. The relationship between transportation cost and transportation risk. The green circles represent the experimental results when the customer nodes are 14, 17, and 10, while the blue triangles represent the experimental results when the customer nodes are 14, 17, 10, and 18, and the yellow square represents the experimental results when the customer nodes are 14, 17, 10, 18, and 22.
Figure 11. The relationship between transportation cost and transportation risk. The green circles represent the experimental results when the customer nodes are 14, 17, and 10, while the blue triangles represent the experimental results when the customer nodes are 14, 17, 10, and 18, and the yellow square represents the experimental results when the customer nodes are 14, 17, 10, 18, and 22.
Sustainability 16 00779 g011
Figure 12. The relationship between transportation cost and carbon emissions. The green circles represent the experimental results when the customer nodes are 14, 17, and 10, while the blue triangles represent the experimental results when the customer nodes are 14, 17, 10, and 18, and the yellow square represents the experimental results when the customer nodes are 14, 17, 10, 18, and 22.
Figure 12. The relationship between transportation cost and carbon emissions. The green circles represent the experimental results when the customer nodes are 14, 17, and 10, while the blue triangles represent the experimental results when the customer nodes are 14, 17, 10, and 18, and the yellow square represents the experimental results when the customer nodes are 14, 17, 10, 18, and 22.
Sustainability 16 00779 g012
Table 1. Transportation road segment.
Table 1. Transportation road segment.
No.Road Segment Length / k m Conditional Release ProbabilityAccident Rate Population   Density / Person · k m 2
ρ i j o ρ i j 1 l ρ i j 2 l ρ i j 3 l ρ i j 4 l ρ i j 5 l ρ i j 6 l ρ i j 7 l
1[1, 2]18.510.00570.0014228180141316106132133141
2[1, 3]15.190.00660.004612931128951445254438387
3[2, 6]23.510.00370.008767181415496435353448385
4[3, 4]23.910.00620.004539620390198342107278474
5[4, 5]18.900.00760.0066381316473325282122264456
6[5, 6]14.020.00570.004136117238730312224227661
7[3, 12]27.950.00460.001017917823716627237990389
8[4, 11]13.290.00750.003948443631067206161131113
9[5, 9]26.550.00690.0022359440136271141489149161
10[6, 8]27.190.00660.006627011344984248271121192
11[9, 8]16.920.00740.0026379361149472272165271351
12[8, 7]26.420.00680.0022169111469130497329197442
13[9, 10]14.940.00380.0096227266166437442210445358
14[8, 16]16.250.00620.003544279232371483411158412
15[7, 18]18.650.00510.004640433641532041044165412
16[12, 11]12.570.00470.009029621919249622716629356
17[11, 10]17.620.00500.0052196430204125274137271315
18[10, 16]17.260.00500.00257526075297493370339177
19[16, 18]18.580.00600.009724510585110468437364465
20[12, 13]16.310.00420.0065209378107446383439253302
21[11, 14]16.160.00440.00544321391046856121101127
22[10, 15]25.680.00640.008235729443425246940346961
23[10, 17]20.270.00490.0034482198421139468245269490
24[16, 17]15.810.00630.00736038816335726214150106
25[17, 19]10.110.00740.00284567523219821035924493
26[18, 20]28.350.00510.008919449529613832934724198
27[14, 15]17.530.00550.0047163106300393437211187407
28[15, 19]29.070.00580.0013315450314156495156394442
29[14, 23]27.110.00430.001324239714024044397123495
30[15, 22]16.370.00660.0056218350179295176251236472
31[19, 20]18.920.00620.00464082395349461271426161
32[23, 22]22.060.00750.00822461635614086375425386
33[23, 24]22.540.00660.006123011829818389321170431
34[22, 21]15.670.00790.0023386343313226299372462155
35[22, 20]21.530.00420.0027434247263162127244435447
36[13, 24]12.870.00400.0045387235335442104330469264
37[24, 21]11.520.00750.006320314585347125167122425
38[21, 20]14.630.00600.0092403417125477245477371345
In ρ i j n l , n represents the time periods.
Table 2. Customer nodes.
Table 2. Customer nodes.
Customer NodeDemand/tTime WindowService Time/h η c 1 /CNY · h 1 η c 2 /CNY · h 1
142[8:30, 10:00]0.6911
175[9:00, 10:00]0.8610
189[13:00, 16:00]0.9514
Table 3. Time periods.
Table 3. Time periods.
No.StartEnd Speed / k m · h 1
10:006:0060
26:007:0070
38:0011:0080
411:0013:0070
513:0018:0080
618:0020:0070
720:000:0060
Table 4. Other parameters.
Table 4. Other parameters.
ParameterValue
Fix cost180 CNY
Transportation cost5 CNY/km
α0.1
β0.2
r0.5 km
E2.61 kg/L
ε m 0.255 L/km
ε 0 0.165 L/km
Maximum loading10 t
Table 5. The results of Experiment 1.
Table 5. The results of Experiment 1.
Departure TimeRoad No.Vehicle No.RouteCost/CNYRiskCarbon Emission/kg
4:20 r 1 1 11, 3, 4, 11, 14, 15, 10, 17, 10, 11, 4, 3, 12646.9320.47240.55
21, 3, 4, 11, 10, 16, 18, 16, 8, 6, 2, 1
r 2 1 11, 3, 12, 11, 14, 15, 19, 17, 16, 8, 6, 2, 12700.9110.15242.34
21, 3, 12, 11, 10, 16, 18, 16, 8, 6, 2, 1
r 3 1 11, 3, 4, 11, 14, 15, 10, 17, 10, 11, 14, 15, 22, 21, 20, 18, 7, 8, 6, 2, 12668.6913.69239.26
21, 3, 4, 11, 10, 15, 22, 21, 20, 18, 7, 8, 6, 2, 1
7:20 r 1 2 11, 3, 4, 11, 14, 15, 19, 17, 16, 8, 6, 2, 12614.058.72235.87
21, 3, 4, 11, 10, 16, 18, 16, 8, 6, 2, 1
r 2 2 11, 3, 4, 11, 14, 15, 10, 17, 10, 11, 4, 3, 12592.3512.37226.63
21, 3, 4, 11, 10, 16, 18, 16, 8, 6, 2, 1
9:20 r 1 3 11, 3, 4, 11, 14, 15, 10, 17, 10, 11, 4, 3, 12564.2317.56228.83
21, 3, 4, 11, 10, 16, 18, 16, 10, 11, 4, 3, 1
r 2 3 11, 3, 12, 11, 14, 23, 22, 20, 19, 17, 10, 11, 12, 3, 12808.389.45248.79
21, 3, 12, 11, 10, 16, 18, 16, 10, 11, 12, 3, 1
r 3 3 11, 3, 4, 11, 14, 23, 22, 20, 19, 17, 10, 11, 4, 3, 12742.289.18253.70
21, 3, 4, 11, 10, 16, 18, 16, 10, 11, 4, 3, 1
12:20 r 1 4 11, 3, 12, 11, 14, 15, 19, 17, 19, 20, 21, 24, 13, 12, 3, 13048.2312.51273.31
21, 3, 12, 11, 10, 9, 8, 7, 18, 20, 21, 24, 13, 12, 3, 1
r 2 4 11, 3, 4, 11, 14, 15, 10, 17, 10, 11, 4, 3, 12545.8123.12218.83
21, 3, 4, 11, 10, 16, 18, 16, 8, 6, 2, 1
r 3 4 11, 3, 12, 11, 14, 15, 19, 17, 16, 8, 6, 2, 12600.5512.72222.32
21, 3, 12, 11, 10, 16, 18, 16, 8, 6, 2, 1
15:20 r 1 5 11, 3, 4, 11, 14, 15, 10, 17, 10, 11, 4, 3, 12601.4117.87219.71
21, 3, 4, 11, 10, 16, 18, 16, 10, 11, 4, 3, 1
r 2 5 11, 3, 12, 11, 14, 15, 19, 17, 19, 20, 21, 24, 13, 12, 3, 13214.3711.48291.88
21, 3, 12, 11, 10, 9, 8, 7, 18, 7, 8, 9, 10, 11, 12, 3, 1
r 3 5 11, 3, 4, 11, 14, 15, 19, 17, 16, 8, 6, 2, 12612.7511.77223.60
21, 3, 4, 11, 10, 16, 18, 16, 8, 6, 2, 1
19:20 r 1 6 11, 3, 4, 11, 14, 15, 19, 17, 19, 20, 21, 24, 13, 12, 3, 13163.5511.57270.43
21, 3, 4, 11, 10, 9, 8, 7, 18, 7, 8, 9, 10, 11, 4, 3, 1
r 2 6 11, 3, 4, 11, 14, 15, 10, 17, 10, 11, 4, 3, 12748.1718.94218.83
21, 3, 4, 11, 10, 16, 18, 16, 8, 6, 2, 1
r 3 6 11, 3, 12, 11, 14, 15, 19, 17, 16, 8, 6, 2, 12802.2210.33222.32
21, 3, 12, 11, 10, 16, 18, 16, 8, 6, 2, 1
21:20 r 1 7 11, 3, 4, 11, 14, 15, 10, 17, 10, 11, 4, 3, 12756.2820.73227.78
21, 3, 4, 11, 10, 16, 18, 16, 8, 6, 2, 1
r 2 7 11, 3, 12, 11, 14, 15, 19, 17, 19, 20, 22, 15, 10, 11, 12, 3, 13354.5513.27274.43
21, 3, 12, 11, 10, 17, 19, 20, 18, 16, 8, 6, 2, 1
r 3 7 11, 3, 12, 11, 14, 15, 10, 17, 10, 11, 12, 3, 12807.2914.34221.91
21, 3, 12, 11, 10, 16, 18, 16, 8, 6, 2, 1
Table 6. New customer nodes.
Table 6. New customer nodes.
Customer NodeDemand/tTime WindowService Time/h η c 1 /CNY · h 1 η c 2 /CNY · h 1
104[10:00, 13:00]0.649
228[11:00, 14:00]1312
Table 7. The results of Experiment 2.
Table 7. The results of Experiment 2.
Customer NodeRoad No.Vehicle No.RouteCost
/CNY
RiskCarbon Emission/kg
14, 17, 18 r 1 8 11, 3, 4, 11, 14, 15, 10, 17, 10, 11, 4, 3, 12642.9823.67230.23
21, 3, 4, 11, 10, 17, 16, 18, 16, 10, 11, 4, 3, 1
r 2 8 11, 3, 4, 11, 14, 15, 19, 17, 16, 8, 6, 2, 12655.7516.70228.29
21, 3, 4, 11, 10, 17, 16, 18, 16, 8, 6, 2, 1
r 3 8 11, 3, 12, 11, 14, 15, 19, 17, 16, 8, 6, 2, 12595.8713.26247.44
21, 3, 12, 11, 10, 16, 18, 16, 8, 6, 2, 1
r 4 8 11, 3, 4, 11, 14, 15, 10, 17, 10, 11, 4, 3, 12550.3321.46232.74
21, 3, 4, 11, 10, 16, 18, 16, 10, 11, 4, 3, 1
r 5 8 11, 3, 4, 11, 14, 15, 19, 17, 16, 8, 6, 2, 12563.0914.49230.81
21, 3, 4, 11, 10, 16, 18, 16, 8, 6, 2, 1
14, 17, 10, 18 r 1 9 11, 3, 4, 11, 14, 15, 10, 11, 4, 3, 13521.0626.28314.08
21, 3, 4, 11, 10, 17, 10, 11, 4, 3, 1
31, 3, 4, 11, 10, 17, 16, 18, 16, 8, 6, 2, 1
r 2 9 11, 3, 12, 11, 14, 23, 22, 20, 19, 17, 10, 11, 12, 3, 13939.2718.89399.89
21, 3, 12, 11, 10, 17, 16, 8, 6, 2, 1
31, 3, 12, 11, 10, 16, 18, 16, 8, 6, 2, 1
r 3 9 11, 3, 12, 11, 14, 23, 22, 20, 19, 17, 10, 11, 12, 3, 13926.5718.89403.10
21, 3, 12, 11, 10, 17, 10, 11, 12, 3, 1
31, 3, 12, 11, 10, 16, 18, 16, 10, 11, 12, 3, 1
r 4 9 11, 3, 4, 11, 14, 15, 10, 11, 4, 3, 13428.4024.07316.60
21, 3, 4, 11, 10, 17, 10, 11, 4, 3, 1
31, 3, 4, 11, 10, 16, 18, 16, 8, 6, 2, 1
r 5 9 11, 3, 4, 11, 14, 23, 22, 20, 19, 17, 10, 11, 4, 3, 13826.8220.75371.16
21, 3, 4, 11, 10, 17, 10, 11, 4, 3, 1
31, 3, 4, 11, 10, 16, 18, 16, 10, 11, 4, 3, 1
r 6 9 11, 3, 12, 11, 14, 15, 10, 11, 12, 3, 13511.1322.29346.58
21, 3, 12, 11, 10, 17, 10, 11, 12, 3, 1
31, 3, 12, 11, 10, 16, 18, 16, 8, 6, 2, 1
14, 17, 10, 18, 22 r 1 10 11, 3, 12, 13, 24, 21, 22, 20, 19, 17, 10, 11, 14, 15, 10, 11, 12, 3, 15956.8723.58574.91
21, 3, 12, 13, 24, 21, 22, 20, 19, 17, 10, 11, 12, 3, 1
31, 3, 12, 13, 24, 21, 22, 21, 24, 13, 12, 3, 1
41, 3, 12, 13, 24, 21, 22, 20, 18, 20, 21, 24, 13, 12, 3, 1
r 2 10 11, 3, 4, 11, 14, 15, 10, 11, 4, 3, 14882.8030.69446.20
21, 3, 4, 11, 10, 17, 10, 11, 4, 3, 1
31, 3, 4, 11, 10, 16, 17, 19, 20, 22, 21, 24, 13, 12, 3, 1
41, 3, 4, 11, 10, 16, 18, 16, 10, 11, 4, 3, 1
r 3 10 11, 3, 12, 11, 14, 15, 10, 11, 12, 3, 14998.8328.29486.45
21, 3, 12, 11, 10, 17, 10, 11, 12, 3, 1
31, 3, 12, 11, 10, 16, 17, 19, 20, 22, 21, 24, 13, 12, 3, 1
41, 3, 12, 11, 10, 16, 18, 16, 10, 11, 12, 3, 1
r 4 10 11, 3, 12, 13, 24, 21, 20, 19, 15, 14, 11, 10, 11, 12, 3, 15481.9926.28523.81
21, 3, 12, 13, 24, 21, 20, 19, 17, 10, 11, 12, 3, 1
31, 3, 12, 13, 24, 21, 22, 21, 24, 13, 12, 3, 1
41, 3, 12, 13, 24, 21, 20, 18, 20, 21, 24, 13, 12, 3, 1
r 5 10 11, 3, 4, 11, 14, 23, 22, 20, 19, 17, 10, 11, 4, 3, 15363.2727.37508.45
21, 3, 4, 11, 10, 17, 10, 11, 4, 3, 1
31, 3, 4, 11, 10, 16, 17, 19, 20, 22, 23, 14, 11, 4, 3, 1
41, 3, 4, 11, 10, 16, 18, 16, 10, 11, 4, 3, 1
r 6 10 11, 3, 12, 11, 14, 23, 22, 20, 19, 17, 10, 11, 12, 3, 15475.6224.90548.23
21, 3, 12, 11, 10, 17, 10, 11, 12, 3, 1
31, 3, 12, 11, 10, 16, 17, 19, 20, 22, 23, 24, 13, 12, 3, 1
41, 3, 12, 11, 10, 16, 18, 16, 10, 11, 12, 3, 1
r 7 10 11, 3, 12, 13, 24, 21, 20, 19, 17, 10, 11, 14, 15, 10, 11, 12, 3, 15618.4724.89542.24
21, 3, 12, 13, 24, 21, 20, 19, 17, 10, 11, 12, 3, 1
31, 3, 12, 13, 24, 21, 22, 21, 24, 13, 12, 3, 1
41, 3, 12, 13, 24, 21, 20, 18, 20, 21, 24, 13, 12, 3, 1
r 8 10 11, 3, 12, 13, 24, 21, 22, 20, 19, 15, 14, 11, 10, 11, 12, 3, 15823.0624.81565.13
21, 3, 12, 13, 24, 21, 22, 20, 19, 17, 10, 11, 12, 3, 1
31, 3, 12, 13, 24, 21, 22, 21, 24, 13, 12, 3, 1
41, 3, 12, 13, 24, 21, 22, 20, 18, 20, 21, 24, 13, 12, 3, 1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zou, Z.; Kang, S. Route Optimization for Hazardous Chemicals Transportation under Time-Varying Conditions. Sustainability 2024, 16, 779. https://doi.org/10.3390/su16020779

AMA Style

Zou Z, Kang S. Route Optimization for Hazardous Chemicals Transportation under Time-Varying Conditions. Sustainability. 2024; 16(2):779. https://doi.org/10.3390/su16020779

Chicago/Turabian Style

Zou, Zongfeng, and Shuangping Kang. 2024. "Route Optimization for Hazardous Chemicals Transportation under Time-Varying Conditions" Sustainability 16, no. 2: 779. https://doi.org/10.3390/su16020779

APA Style

Zou, Z., & Kang, S. (2024). Route Optimization for Hazardous Chemicals Transportation under Time-Varying Conditions. Sustainability, 16(2), 779. https://doi.org/10.3390/su16020779

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop