Next Article in Journal
Vibrational Communication of Scolypopa australis (Walker, 1851) (Hemiptera: Ricaniidae)—Towards a Novel Sustainable Pest Management Tool
Previous Article in Journal
Performance of Concrete Pavement Incorporating Portland Limestone Cement in Cold Weather
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improving the Resilience of Port–Hinterland Container Logistics Transportation Systems: A Bi-Level Programming Approach

School of Management, Zhejiang University, Hangzhou 310058, China
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(1), 180; https://doi.org/10.3390/su14010180
Submission received: 19 November 2021 / Revised: 16 December 2021 / Accepted: 21 December 2021 / Published: 24 December 2021
(This article belongs to the Section Sustainable Transportation)

Abstract

:
Port–hinterland container logistics transportation systems (PHCLTSs) are significant to economic and social development. However, various kinds of unconventional emergency events (UEEs), such as natural or human-caused disasters, threaten PHCLTSs. This study aims to measure and improve the resilience of PHCLTSs. Bi-level programming models with two different lower level models are established to help PHCLTSs recover their capacity efficiently in the face of UEEs. In the upper level model, the government makes immediate recovery decisions about a damaged PHCLTS with the goal of improving the resilience of the PHCLTS. In the lower level models, truck carriers make decisions about transportation routes and freight volume in the recovered PHCLTS. They cooperate fully to pursue the maximization of total profit and are coordinated by a central authority, or they make their own decisions to pursue maximization of their own profit noncooperatively. An algorithm combining particle swarm optimization (PSO) and traditional optimization algorithms is proposed to solve the bi-level programming models. The numerical experimental results show the validity of the proposed models.

1. Introduction

World trade has maintained a growth trend since the 21st century. The World Trade Organization [1] reported that trade in goods expanded by 4.6% on average between 2005 and 2017, which enables the continuous growth of transportation demand. China’s “One Belt One Road” initiative promotes world trade development and creates higher requirements for transportation infrastructure.
Container transportation is an advanced transportation mode in the modern transportation industry, which is an important tool to realize multimodal transportation and to complete commodity transportation and delivery for world trade efficiently and safely. Maritime transportation systems are important infrastructures to enable transportation activities [2]. This study focuses on a crucial subsystem of maritime transportation systems called the “port–hinterland container logistics transportation system” (PHCLTS). It is the landside component of maritime transportation systems [3], which realizes container transportation between a port and its hinterlands to support the national economy and people’s lives.
A PHCLTS is an extremely complex transportation system that contains diverse infrastructures, transportation modes, and participants. Figure 1 shows the infrastructure structure of a PHCLTS, and it is composed of several nodes and arcs. The nodes of the PHCLTS represent a port, several multimodal connection points, and numerous hinterland points. Only one port exists in a PHCLTS, and the port is an important hub connecting maritime process and port–hinterland transportation process. The transportation mode of containers can be changed at the multimodal connection points, including multimodal terminals, dry ports, and inland ports. The hinterland points correspond to the demand and supply areas, including warehouses, distribution centers, and factories. The arcs of the PHCLTS represent roads, railways, and inland rivers, which connect the different nodes of the PHCLTS. A PHCLTS contains diverse transportation modes, including road, rail, inland river, and multimodal transportation. It also contains diverse participants, such as the government, shippers, port companies, truck carriers, rail carriers, inland-river carriers, and logistics service providers.
However, various kinds of unconventional emergency events (UEEs), such as natural or human-caused disasters, threaten PHCLTSs. Consequently, the performance of PHCLTSs will decline and exert negative effects on the economy and society. For example, in August 2005, the United States of America was hit by Hurricane Katrina, which stopped the operation of highway and railway trunk lines between the Port of New Orleans and its hinterland. The supply chains of many industries were interrupted because of the severe damage in infrastructure, people’s lives and the production processes of enterprises were affected, and the economic growth rate of the United States of America slowed down.
Improving the recovery capability of PHCLTSs in the face of UEEs is critical. The concept of “resilience” is applied to describe such a capability; this concept refers to the capacity of a logistics transportation system to recover to an acceptable operational service level after being affected by a disruption, and the level of efforts under resource restrictions, such as immediate recovery activities implemented in the logistics transportation system, should be considered when measuring the resilience of the logistics transportation system [4]. This study aims to provide an optimal decision-making reference for governments to improve the resilience of PHCLTSs in the face of UEEs to recover the operational capacity efficiently and reduce the negative effects of UEEs.
Logistics transportation system participants, including organizations and individuals, will affect the resilience of the system through their actions and interactions [5]. To study the resilience problem of PHCLTSs, the behaviors and interactions of associated participants should be analyzed.
Therefore, a bi-level programming model consisting of two submodels is applied to analyze the resilience problem of PHCLTSs. The upper and lower level models have decision variables, objective functions, and constraints of their own. Furthermore, it contains the hierarchical relationship and interaction between decision makers of the two submodels and can deal well with their two conflicting goals [6]. For the upper level model, the decision maker is the planner of PHCLTSs, i.e., the government, which as the leader makes decision on immediate recovery activities to improve the resilience of PHCLTSs. For the lower level model, the decision makers are the users of PHCLTSs, i.e., shippers and carriers. The shippers make decisions about the demands for containers in hinterlands, and the carriers complete the actual transportation process of containers in the recovered PHCLTS that are from port to hinterland points in accordance with the shippers’ demand. In this study, truck carriers are selected because road transportation is a critical mode of freight transportation at present. Truck carriers, which are the followers, make decisions about transportation routes and freight volume in the recovered PHCLTS to pursue profit maximization [7], affecting the achievement of the goal of the government.
The contributions of this paper are the following:
(1)
Time resilience is adopted to measure and improve the resilience of a PHCLTS with considerations of PHCLTS user behaviors and the hierarchical relationship and interaction of government and truck carriers.
(2)
Bi-level programming models with two different lower level models are established according to different behavior patterns of PHCLTS users, supporting the government in selecting a reasonable budget level and implementing efficient recovery decisions to aid the PHCLTS to recover the operational capacity in the face of UEEs.
(3)
It expands the prearranged planning about UEEs and provides a decision-making reference and scientific theoretical support for an efficient response to UEEs.
The remainder of this paper is structured as follows. An overview of the relevant literature is provided in Section 2. The bi-level programming models with two different lower level models are established in Section 3. The solution algorithm is constructed in Section 4. Numerical experiments are conducted to illustrate the validity of the models in Section 5. Lastly, managerial insights and conclusions are summarized in Section 6.

2. Literature Review

This study analyzes the resilience problem of PHCLTSs by using the bi-level programming approach. It aims to learn about three aspects from the relevant literature: the resilience of transportation systems corresponding to the goal of the government, logistics transportation system assignment and equilibrium corresponding to the behaviors of PHCLTS users, and the bi-level programming model for logistics transportation system problems corresponding to the approach used in this study.

2.1. Resilience of Transportation Systems

Scholars discussed characteristics of resilient systems and sought to optimize system resilience. Bruneau et al. [8] defined resilience as the capacity of a system to alleviate the negative influence of a shock, absorb it, and recover promptly after it happens. Murray-Tuite [9] listed 10 dimensions of transportation resilience and provided several measurement methods for transportation resilience. Ta et al. [5] discussed the properties of the resilience of a freight transportation system and argued that the resilience is related to physical infrastructure, managing organization, and users. Quantitative studies are conducted on the basis of qualitative studies on this issue. Scholars established mathematical models to measure the resilience of transportation systems while accounting for the limitation of time, cost, and other resources. Omer et al. [10] defined resilience as the capacity of a maritime transportation system to recover to the normal level in the face of a shock, and tonnage resiliency, time resiliency, and cost resiliency are proposed as metrics. Faturechi and Miller-Hooks [11] established a bi-level and three-stage stochastic programming model to optimize travel time resilience of passenger transportation systems.
Regarding the resilience of port–hinterland logistics and multimodal transportation systems, Nair et al. [12] employed resilience as a quantitative measure to determine the optimal actions to enhance security at the nodal facilities of a freight transportation network. Chen and Miller-Hooks [4] indicated that the level of effort should be considered when measuring the resilience of an intermodal freight transportation system. They proposed an indicator that the system resilience should be equal to the ratio of the transportation demand satisfied post-disaster to the transportation demand satisfied pre-disaster. Subsequently, Miller-Hooks et al. [13] maximized the resilience of freight transportation networks while considering the implementation of pre-disaster preparation activities and post-disaster recovery activities. Chen et al. [14] introduced a quantitative method to measure the resilience of a port–hinterland container transportation network from the perspective of shippers. Chen et al. [3] used the network game theory approach to analyze the strategic investment decision of each participant in a port–hinterland container transportation network to enhance the network’s resilience.

2.2. Logistics Transportation System Assignment and Equilibrium

The behaviors of the transportation system users will affect the transportation system resilience. Existing research on modeling behaviors of logistics transportation system users mostly follows Wardrop’s two principles. Wardrop’s first principle indicates that each transportation system user pursues cost minimization noncooperatively. A “user equilibrium (UE)” assignment in accordance with Wardrop’s first principle occurs when no one in the network can decrease his cost by unilateral action [15]. Using Wardrop’s first principle, Jones et al. [16] established an equilibrium model with capacity constraints of links and congestion delays to analyze import and export container flows of the U.S. Meng and Wang [17] studied the intermodal transportation system design problem with multiple stakeholders and multitype containers by using a mathematical program with an equilibrium constraint model. Corman et al. [18] proposed an assignment model which can obtain the equilibrium flows and generalized costs to optimize the performance of transportation systems.
Wardrop’s second principle indicates a central authority coordinates all transportation system users to cooperate fully to achieve the total transportation cost minimization of the entire network. A “system optimum (SO)” assignment in accordance with Wardrop’s second principle occurs when the total marginal costs of transportation alternatives are equal [15]. Using Wardrop’s second principle, Zhang et al. [19] adopted an SO model to minimize the total cost of a trucking company in the upper level of a bi-level programming model to analyze the urban freight transportation problem. Jiang et al. [20] established a bi-level programming model to study the design problem of a regional multimodal logistics system with the goal of emission reduction under uncertain demands, and in the lower level model, the carrier arranged the transportation routes for freight demand based on the SO principle to minimize the total generalized transportation cost. Yu and Jiang [21] constructed a bi-level model to study the logistics transportation system design and delivery scheme problem while adopting an SO model to achieve cost minimization in the lower level model.

2.3. Bi-Level Programming Model for Logistics Transportation System Problems

Research on logistics transportation system optimization that applies a bi-level programming model remains limited. In general, a planning authority in the upper level seeks to maximize the total utility of the entire transportation system. Logistics transportation system users in the lower level determine the transportation process according to the authority’s decision. Moreno-Quintero et al. [22] established a bi-level programming model to characterize the interaction between the road planner and carriers and studied the road maintenance and overloading control problem. Lee et al. [23] constructed a bi-level programming model to capture interactions between shippers and different types of carriers in maritime logistics transportation systems. Qiu et al. [24] employed a bi-level programming model to analyze the carbon emission problem in road freight transportation. Li et al. [25] used the goal programming approach to study the problem of logistics infrastructure investment and CO2 emission reduction.

2.4. Summary of Literature Review

“Resilience” is a relatively new concept for logistics transportation systems. Previous research mostly studied the resilience problems of logistics transportation systems from the perspective of demand realization [4,12,13,14], and such research quantitatively characterized the resilience of logistics transportation systems by using the ratio of the freight demand that can be satisfied post-disaster to the freight demand that can be satisfied pre-disaster. Their perspective focused on measuring and improving the resilience of logistics transportation systems under the following single transportation period after UEEs.
In this study, a medium- and long-term perspective, which focuses on the resilience of logistics transportation systems under multiple transportation periods after UEEs, is applied. The recovered state of a PHCLTS will last for a period of time after the implementation of immediate recovery activities, and the logistics transportation process will be repeated for multiple transportation periods during this time. Hence, time resilience is adopted to measure and improve the resilience of a PHCLTS. It is represented as the ratio of the total transportation time of the PHCLTS in each transportation period pre-disaster to the total transportation time of the PHCLTS in each transportation period post-disaster.
Previous studies rarely considered the hierarchical relationship and interaction among different participants when analyzing the resilience problem of transportation systems. For instance, Faturechi and Miller-Hooks [11] optimized the travel time resilience of a passenger transportation system while employing a bi-level programming approach to describe the behaviors and interactions of different participants. User behaviors were formulated as a partial user equilibrium in the lower level model. In the current study, the objectives are to explore PHCLTSs, which are freight transportation systems, and to improve the resilience of PHCLTSs, while the hierarchical relationship and interactions among different participants of a PHCLTS are considered. Specifically, two different behavior modes of PHCLTS users are considered on the basis of the SO and UE principles.
Optimization problems of logistics transportation systems have been successfully solved using bi-level programming methods under conventional situations. Nevertheless, research on logistics transportation systems with different decision makers in the face of UEEs is still lacking. This study fills the research gap in this field to analyze how to improve the resilience of PHCLTSs under UEE situations while considering PHCLTS user behaviors, and the hierarchical relationship and interaction of the PHCLTS planner and users. It provides decision support to governments to make highly effective decisions about immediate recovery activities after UEEs, aiding PHCLTSs to recover their operational capacity efficiently.

3. Model

3.1. Model Description and Assumptions

Only one port exists in the targeted PHCLTS. All containers are transported from the port to hinterland points via single mode transportation. The research target is a single-mode PHCLTS, i.e., road transportation by truck carriers. The reason is that road transportation is the critical mode of freight transportation at present. Some PHCLTSs are relatively small in size in the early stage of their development and are connected only by roads. For example, Jones et al. [16] mentioned that few railway facilities had access to the ports in South Florida at that time, and almost all freight transportation between the ports and hinterlands was realized through road transportation.
The participants of a PHCLTS include the PHCLTS planner and users. The PHCLTS planner is the government, which is the leader and the decision maker in the upper level model. The PHCLTS users include shippers and carriers, which are the followers and decision makers in the lower level model. The decision framework is shown in Figure 2.
UEEs affect the arcs of PHCLTSs. After the occurrence of a UEE, the government makes decisions on immediate recovery activities for a damaged PHCLTS within the given budget constraint, with the goal of improving the resilience of the PHCLTS. The recovery activities refer to the relevant activities implemented immediately to reduce the negative influence of a UEE and restore the capacity of the PHCLTS. In this study, the recovery activities are filling and repairing the damaged part of an arc to recover a proportion of the transportation capacity of the arc. The recovery activities are implemented only once for each attacked arc in the PHCLTS. The recovery state of the PHCLTS will remain for a period. To complete the recovery activities on each damaged arc efficiently, no limitation is imposed on the implementation time of the recovery activities. The government fully understands the container demands of shippers in the PHCLTS and the behavior modes of truck carriers in the recovered PHCLTS. After the government completes the implementation of recovery activities, the transportation capacity and post-disaster, free-flow travel time of each arc of the recovered PHCLTS will be determined. These two factors will have an influence on truck carriers in arranging the transportation process of containers to achieve the goal of profit maximization.
There are diverse decision makers in the freight transportation process [26]. The demands for containers in hinterland points are determined by shippers, which is one type of PHCLTS user. Container transportation activities from the port to hinterland points are repeated every day. The daily transportation demands of hinterland points are determined by shippers, which is one type of PHCLTS user, and the container demands are constants and unaffected by UEEs, i.e., the daily transportation demands of hinterland points are regarded as exogenous constants. This assumption has also been adopted in the existing literatures [14,18]. Before a UEE happens, all the transportation demands in hinterland points can be satisfied. The sizes of the containers are the same in this study, and all containers are forty-foot equivalent unit (FEU) containers. The value of goods loaded in each container is equal, i.e., the urgency for transporting each container is the same. Each truck transports one FEU container at one time.
Another type of PHCLTS user is truck carriers. The truck carriers fully understand the recovery state of the PHCLTS. In the recovered PHCLTS, in accordance with the container demands of the shippers, the truck carriers make decisions about transportation routes and freight volume on each arc to pursue profit maximization. Then, the freight volume and the actual transportation time on each arc are determined. These two factors will affect the operation level of the PHCLTS, thus having an influence on the government in achieving the goal of improving the resilience of the PHCLTS. Considering that the container demands of the shippers in the hinterland points are constants, the actual PHCLTS users to be analyzed are the truck carriers. Nonetheless, the shippers are indispensable participants in the PHCLTS because freight transportation is a derivative demand. First, an original demand must exist. Then, actual the transportation process of the containers by the truck carriers will occur.
The transportation fee of one container transport between the same origin–destination (O–D) pair paid by the shippers to truck carriers is constant and unaffected by UEEs. Given that the demands and the unit transportation fee of containers are constant, the revenue of truck carriers is a fixed value, and the goal of profit maximization can be transformed into the goal of transportation cost minimization. Road impedance exists in PHCLTSs. As the flow on one arc increases, the transportation time on this arc increases accordingly [18,23].
In accordance with the above description, the models of this study are based on the following assumptions:
(1)
The daily transportation demands of containers are constant and unaffected by UEEs.
(2)
All containers are homogeneous, are the same size, and have the same value of goods loaded.
(3)
A recovery activity is implemented only once for each attacked arc in a PHCLTS.
(4)
The government has sufficient time to complete the recovery activities.
(5)
Road impedance exists in the PHCLTS.
The notations employed in the proposed models are indicated in Table 1.

3.2. Time Resilience of a PHCLTS

Time resilience is employed to measure the resilience of a PHCLTS, which represents the influence of UEEs on the total transportation time. Owing to the randomness of UEEs, the resilience level of a PHCLTS is defined as the expected value of the ratio of total transportation time pre-disaster to total transportation time post-disaster after the implementation of recovery activities under different UEE scenarios:
R T = E ( T t o t a l s i n g l e m o d e / T t o t a l s i n g l e m o d e ( ξ ) )  
The total transportation time pre-disaster, T t o t a l s i n g l e m o d e , can be calculated as:
T t o t a l s i n g l e m o d e = a A h a · t a  
Before the PHCLTS is affected by a UEE, the operation state of the transportation system is stable. The transportation arrangement of truck carriers for containers is the same in each transportation cycle. Therefore, h a and t a are constants, then T t o t a l s i n g l e m o d e is a constant. The total transportation time post-disaster after the implementation of recovery activities, T t o t a l s i n g l e m o d e ( ξ ) , will be solved using the bi-level programming models established in this study.

3.3. Model Formulation

3.3.1. Upper Level Model

The government implements recovery activities within the given budget constraint of cost to recover the transportation capacity of each affected arc. The goal of the government is to seek the maximization of RT. In accordance with the PHCLTS resilience measure expression, Equation (1), the total transportation time pre-disaster, T t o t a l s i n g l e m o d e , is a deterministic value. The goal of the government can be replaced with minimizing the total transportation time post-disaster after recovery activities, T t o t a l s i n g l e m o d e ( ξ ) . Therefore, under a realized UEE scenario ξ , the upper level model is:
min   T t o t a l s i n g l e m o d e ( ξ ) = a A h a ( ξ ) · t a ( ξ )  
s . t .
a b a · l a ( ξ ) · ( c a c a u ( ξ ) ) c a B  
0 < l a ( ξ ) 1 ,   if   a A u  
l a ( ξ ) = 0 ,   if   a A n u  
a A  
Objective function (3) seeks to minimize the total transportation time post-disaster T t o t a l s i n g l e m o d e ( ξ ) under each UEE scenario to improve the resilience of the PHCLTS. h a ( ξ ) is determined by the container transportation arrangement of truck carriers in the lower level models. t a ( ξ ) is jointly determined by the decision on recovery activities of the government in the upper level model and the container transportation arrangement of truck carriers in the lower level models. The calculation formula of t a ( ξ ) will be described in the next section.
Constraint (4) is the budget constraint. Constraint (5) indicates that if UEEs affect arc a, then a recovery activity may be taken on it. Constraint (6) indicates that if UEEs do not affect arc a, then no recovery activity will be taken on it.
After the implementation of recovery activities by the government, the capacity of arc a after recovery, c a ( ξ ) , can be calculated as:
c a ( ξ ) = c a u ( ξ ) + l a ( ξ ) · ( c a c a u ( ξ ) )  
The transportation capacity of an arc cannot exceed its original capacity pre-disaster. The change in transportation time of an arc is determined as a function of change in capacity [4]. In this study, if z% decreases in capacity along with an arc, then the free-flow travel time of this arc increases by c h u · z %. For arc a A u , the free-flow travel time on arc a post-disaster and after recovery under a UEE scenario ξ , t a r 0 ( ξ ) , can be calculated as:
t a r 0 ( ξ ) = t a 0 · ( 1 + c h u · ( c a c a ( ξ ) c a ) )  
c a ( ξ ) and t a r 0 ( ξ ) will affect the truck carriers to arrange the container transportation process to achieve the goal of minimizing transportation costs in the lower level models.

3.3.2. Lower Level Model 1: The SO Model

Two different lower level models are established on the basis of an SO model and a UE model. The first lower level model supposes that only one organization manages all trucks in a PHCLTS. That is, only one trucking company exists in the PHCLTS, or two or more trucking companies belonging to one truck association operate in the PHCLTS. For the convenience of description, it is collectively referred to as “trucking company” in the subsequent contents of this paper, which refers to the unique trucking company in the PHCLTS or the truck association that uniformly manages all trucking companies in the PHCLTS. This trucking company is the decision maker of the first lower level model, and it is a central authority coordinating all truck drivers to cooperate fully to minimize total costs. It can be described as an SO model, which follows Wardrop’s second principle.
The lower level SO model is established with respect to the arc flows. The reason is that the lower level SO model with respect to the arc flows is a strict convex optimization problem which has a unique optimum solution. However, the lower level SO model is not convex with respect to path flows, so the corresponding optimum solutions for the path flows may not be unique [27]. The uniqueness of the optimal solution of the lower SO level model is crucial to solving the process of the bi-level programming model. The objective function is to minimize the total costs of the trucking company in accordance with the recovered PHCLTS:
min   G C t o t a l s i n g l e m o d e S O ( ξ ) = a A ( r a + q · t a ( ξ ) ) · h a ( ξ )  
It is a generalized cost function [28], which includes transportation time. The value of time is adopted to translate time into money. In accordance with the model assumption that all containers are homogeneous, the value of time for all containers is equal (parameter q).
Road impedance exists in the PHCLTS. The “Bureau of Public Roads (BPR) function” is applied. The basic mathematical expression of the BPR function is:
T a = T a 0 · ( 1 + α · ( X a C a ) β )  
Ta and Xa are the actual travel time and the flow of arc a, respectively. T a 0 is the free-flow travel time. Ca is the value of flow at which the transportation time on arc a is higher than T a 0 . In this study, such a value is denoted as c a ( ξ ) . The quantities α and β are model parameters [27]. In each transportation period, the value of freight flow on arc a is equal to the freight volume on arc a, i.e., h a ( ξ ) , because each truck transports one container (FEU) at one time. Therefore, the post-disaster transportation time, t a ( ξ ) , in objective function (9) can be calculated as
t a ( ξ ) = t a r 0 ( ξ ) · ( 1 + α · ( h a ( ξ ) c a ( ξ ) ) β )  
The constraints of the lower level SO model include the following:
a A w h a ( ξ ) = D w  
a ( : , k ) A h a ( ξ ) a ( k , : ) A h a ( ξ ) = 0  
h a ( ξ ) c a ( ξ )  
h a ( ξ ) 0  
a A ,   k K  
Constraint (12) means that all freight demands are satisfiable. Constraint (13) requires flow conservation on arcs. Constraint (14) indicates that the freight flow on arc a cannot exceed the transportation capacity of arc a after recovery. Constraint (15) denotes the attributes of the lower level decision variables.
Combined with the upper level model, the bi-level model with the lower level SO model is shown as follows:
Bi-level model with the lower level SO model
Upper   level   model { Objective   function   ( 3 ) s . t . Constraints   ( 4 ) ( 6 ) a A Lower   level   SO   model { Objective   function   ( 9 ) s . t . Constraints   ( 12 ) ( 15 ) a A ,   k K  

3.3.3. Lower Level Model 2: The UE Model

The second lower level model supposes that no central authority exists to coordinate the actions of truck drivers. Under this circumstance, the truck drivers are the decision makers, and each of them decides on his own transportation routes to minimize costs noncooperatively. When all truck drivers cannot reduce their transportation costs by changing personal decisions, the PHCLTS reaches equilibrium. This condition can be described as a UE model, which follows Wardrop’s first principle.
The lower level UE model is established with respect to the arc flows, too, and it is also a strict convex optimization problem which has a unique optimum solution [27]. The generalized cost function is also employed in this model. The objective function of the lower level UE model is for a particular function to reach UE assignment:
min a A 0 h a ( ξ ) ( r a + q · t a ( ξ , v ) ) d v  
The lower level UE model has the same Constraints (12)–(15) as the lower level SO model, which have been illustrated in Section 3.3.2. Combined with the upper level model, the bi-level model with the lower level UE model is shown as follows:
Bi-level model with the lower level UE model
Upper   level   model { Objective   function   ( 3 ) s . t . Constraints   ( 4 ) ( 6 ) a A Lower   level   UE   model { Objective   function   ( 16 ) s . t . Constraints   ( 12 ) ( 15 ) a A ,   k K  

4. Solution Algorithm

Bi-level programming models are NP-hard problems [29], and it is difficult to solve them by using traditional optimization algorithms. Therefore, many scholars have tried to utilize heuristic algorithms to solve them [30,31,32]. No special requirements are imposed for the differentiability and convexity of the objective function and constraints in the models when using heuristic algorithms [32]. In this study, particle swarm optimization (PSO), one of the heuristic algorithms, is adopted to solve the proposed bi-level programming models.

4.1. Description of Standard PSO

PSO was first proposed by Kennedy and Eberhart [33]. It is a swarm intelligence stochastic optimization algorithm inspired by bird foraging behaviors. The PSO algorithm is suitable for solving large-scale, nonlinear, nondifferentiable, and other complex models [32].
The notations employed in the PSO algorithm are presented in Table 2. The standard PSO algorithm is briefly introduced below:
(1)
The velocity vi and position Pi of particle i of the particle swarm are initialized through random generation, i [ 1 ,   n ] . Each particle i corresponds to a solution of the model needed to be solved. A fitness function is adopted to evaluate the particle. The position Pi of each particle i in the initialized population is set to pbest, i [ 1 ,   n ] . The particle with the best fitness function value is set to gbest.
(2)
In each iteration, Equations (17) and (18) are adopted to update the velocity vi and position Pi of particle i to generate a newswarm, i [ 1 ,   n ] .
v i ( t + 1 ) = w v i ( t ) + c 1 r 1 ( p b e s t P i ( t ) ) + c 2 r 2 ( g b e s t P i ( t ) )  
P i ( t + 1 ) = P i ( t ) + v i ( t + 1 )  
where w is the inertia weight. According to Shi [34], when w decreases linearly as the iteration increases, the convergence performance of the PSO algorithm will be improved. The following equation is adopted to calculate the value of w in each iteration:
w = w M a x i t e r · w M a x w M i n i t e r M a x  
The evaluation method is conducted. If particle i is better than pbest, pbest is updated to particle i. If particle i is better than gbest, gbest is updated to particle i.
(3)
The iteration process is repeated until the algorithm termination condition is reached, then the gbest is output to be solution of the model.

4.2. An Algorithm Combining PSO and Traditional Optimization Algorithms

In consideration of the characteristics of the bi-level programming model, a hybrid algorithm, which is a combination of PSO and traditional optimization algorithms, is constructed to solve the model. It can find out an approximate optimal global solution within a reasonable computing time [31,35]. The upper level model is solved using the PSO algorithm. The lower level model is solved using the traditional optimization algorithm, which is realized using MATLAB software (The MathWorks, Inc., Natick, MA, USA) and the TOMLAB/CONOPT solver. The TOMLAB/CONOPT solver uses the generalized reduced gradient (GRG) method and relevant extended methods to solve the nonlinear programming problems [36]. As mentioned in Section 3.3.2 and Section 3.3.3, the two lower level models are strictly convex optimization problems. For such problems, the local optimal solution of the lower level model obtained using the TOMLAB/CONOPT solver is the unique global optimal solution [37].
The solution procedure of the algorithm is illustrated in Figure 3. The steps of the algorithm are described below:
Step 1: The parameters for the PHCTN are set as follows: (1) original transportation capacity of arc a ( c a ), free-flow travel time on arc a pre-disaster ( t a 0 ), transportation cost on arc a (ra), and recovery cost on arc a (ba); (2) daily transportation demand for warehouse w (Dw); (3) maximum budget for recovery activities (B); (4) value of time (q); (5) parameters for the BPR function (α and β).
Step 2: The parameters for a UEE scenario are set as follows: (1) UEE type (u); (2) set of arcs affected by UEEs (Au) and set of arcs unaffected by UEEs (Anu); (3) capacity of arc a post-disaster and before recovery ( c a u ( ξ ) ).
Step 3: The parameters for PSO are set as follows: (1) swarm size (n); (2) maximum number of iterations (iterMax); (3) learning factors c1 and c2; (4) maximum inertia weight (wMax) and minimum inertia weight (wMin); (5) penalty factor in fitness function (M).
Step 4: For the upper level model, velocity vi and position Pi for each particle i, i [ 1 ,   n ] , are initialized.
(1)
The initial velocity vi for each particle i, i [ 1 ,   n ] , is randomly generated with a restriction of −4 to 4. In the subsequent iterations, the restrictions of particle velocities are removed due to the adopted inertia weight w [38].
(2)
For arc a A u , pi,a is a uniform random number in the interval [0, 1] (Figure 4a). For arc a A n u , pi,a = 0.
Step 5: The position Pi may be an infeasible solution when it does not meet the constraints of the upper level model. In this case, referring to Alves et al. [39], a repair method is applied to the infeasible Pi (Figure 4b):
(1)
For Constraint (5), a A u , if pi,a > 1, then it is adjusted to pi,a = 1; if pi,a < 0, then it is adjusted to pi,a = 0.
(2)
For Constraint (6), a A n u , if p i , a 0 , then it is adjusted to pi,a = 0.
(3)
For Constraint (4), if position Pi does not meet this budget constraint, the corresponding total implementation cost B P i of recovery activities for position Pi is calculated as
B P i = a b a · p i , a · ( c a c a u ( ξ ) ) c a
For a A u , p i , a = p i , a · B / B P i . The particle position Pi after adjustment will satisfy all the constraints of the upper level model.
Step 6: The particle position Pi, i.e., the solution of the upper level model, is substituted into the lower level model. The TOMLAB/CONOPT solver is applied to solve the corresponding lower level model (Figure 4a).
(1)
The capacity and free-flow travel time of arc a after recovery, c a ( ξ ) and t a r 0 ( ξ ) , are calculated.
(2)
The optimal solution of the lower level model Hi ( h i ,   a H i , a A ) is obtained by the TOMLAB/CONOPT solver. Vector Hi corresponds to particle position Pi.
Step 7: Particle position Pi and its corresponding vector Hi are combined to form a solution of the bi-level programming model. (Pi, Hi) is substituted into the bi-level programming model. The evaluation method is conducted, and pbest and gbest are updated (Figure 4c).
(1)
The evaluation method for different particles is as follows: First, whether (Pi, Hi) is a feasible solution is determined. Small errors may exist when using the TOMLAB/CONOPT solver, which lead to the possibility that equality constraints of lower level models, Constraints (12) and (13), cannot be accurately satisfied. In this case, with reference to Soares et al. [31], a small positive tolerance ε is set for Constraints (12) and (13) as follows:
| a A w h a ( ξ ) D w | ε
| a ( : , k ) A h a ( ξ ) a ( k , : ) A h a ( ξ ) | ε
If (Pi, Hi) satisfies all constraints of the bi-level programming model, it is a feasible solution. Then, the fitness function value f i t n e s s ( P i ,   H i ) of (Pi, Hi) is calculated. The objective function of the upper level model is adopted as the fitness function:
f i t n e s s ( P i , H i ) = a A h i ,   a · t a ( ξ )
The evaluation criteria for different particles are presented in Table 3.
(2)
For the initial population, the position Pi of each particle i is set to pbest, i [ 1 ,   n ] , and the corresponding vector Hi is also recorded. In accordance with the evaluation criteria in Table 3, the best particle is set to gbest, and the corresponding vector Hi is also recorded.
(3)
For noninitial populations, in accordance with the evaluation criteria in Table 3, if particle i is better than pbest, then pbest is updated to the current location Pi of particle i, and the corresponding vector Hi is also updated. If particle i is better than gbest, then gbest is updated to the current location Pi of particle i, and the corresponding vector Hi is also updated.
Step 8: If iter does not reach iterMax, the velocity vi and position Pi of particle i, i [ 1 ,   n ] , are updated in accordance with Equations (17) and (18), then step 5 is repeated (Figure 4b). If iter reaches iterMax, step 9 is performed.
Step 9: The iteration is terminated. Outputting gbest, the corresponding vector Hi, and the value of the objective function of the upper level model as the ultimate solution of the bi-level programming model.

5. Numerical Studies

5.1. Description of Numerical Studies

A numerical case is presented to test the bi-level programming models and solution procedure. The models are solved on a test PHCLTS consisting of one port and five warehouses in hinterland with 16 arcs between them, as illustrated in Figure 5. Table 4 shows the arc parameters.
The daily transportation demands are 1000, 500, 800, 500, and 1200 for warehouses 1–5, respectively. The total budget for recovery activities is 4000 units. For the objective functions of the lower level models, parameter q = 10. For the BPR function, parameters α = 0.5 and β = 2.
Three UEE types are considered: terrorist attack (u = 1), earthquake (u = 2), and freeze (u = 3). The influence parameters are set as ch1 = 3, ch2 = 5, and ch3 = 4 for terrorist attack, earthquake, and freeze, respectively. Table 5 describes the characteristics of each UEE type. Simple illustrations of each UEE type are provided in Figure 6. Detailed information about how the UEE scenarios of different UEE types affect the infrastructures of the PHCLTS are shown in Appendix A.
Following Chen and Miller-Hooks [4], the Monte Carlo simulation is employed to simulate UEE occurrence and obtain the resilience level of the PHCLTS under each UEE type:
(1)
For each UEE scenario, the Monte Carlo simulation is applied, a UEE subtype is randomly chosen and a set of numbers is randomly generated using the preset probability distribution function based on the particular UEE subtype following the setting information in Appendix A. This numerical set represents the decreasing capacity level of each affected arc. Thus, a scenario in which a UEE harms the test PHCLTS is created.
(2)
The bi-level programming model is solved with corresponding parameters to implement an efficient decision of recovery activities and obtain the resilience of the test PHCLTS under a particular UEE scenario, called point resilience [12].
(3)
The mean value of point resilience in all UEE scenarios is calculated, and it is the resilience level of the test PHCLTS under this UEE type.
For the PSO algorithm, the swarm size n = 40, the learning factors c1= 2 and c2 = 2, the maximum inertia weight wMax = 0.9, the minimum inertia weight wMin = 0.4, and the tolerance ε   = 0.01. Six UEE scenarios are randomly selected to test the convergence of the algorithm, including the SO and UE lower level models. These UEE scenarios converge on stable values when the number of iterations is approximately 20 (Figure 7). In consideration of the balance of calculation accuracy and calculation time to find an approximate global optimal solution in reasonable time, the maximum number of iterations, iterMax, is set to 50.
The capacities of arcs in the PHCLTS may significantly reduce after a UEE scenario, and the freight demand of hinterland points may become unsatisfiable if the government does not prepare an adequate budget to implement sufficient recovery activities in time. Under this circumstance, containers will be stranded in the port, and the unsatisfied containers will increase as time goes on, which will result in heavy economic and social losses. The total transportation time of containers post-disaster can be considered infinite, which is the denominator of Equation (1) for calculating the resilience level. Therefore, if the freight demand of hinterland points cannot be met post-disaster, the point resilience of the test PHCLTS is defined as 0. The government should increase the budget and implement additional recovery activities in this circumstance.
For each UEE type, 100 random scenarios are generated. Table 6 shows the resilience levels and average total transportation time. The UEE scenarios with point resilience equal to 0 are not counted in the average total transportation time. Suppose the occurrence probability of each UEE type is equal, the resilience levels of PHCLTS are 0.8723 and 0.8737 with the lower level SO and UE models, respectively. Under all three UEE types, the PHCLTS resilience levels with the lower level UE model are slightly higher, whereas the average total transportation time of the test PHCLTS with the lower level SO model is shorter. The main reason for the higher resilience levels with the lower level UE model is that the initial total transportation time with the lower level UE model ( T t o t a l s i n g l e m o d e U E = 50,813), which is the numerator of Equation (1) for calculating the resilience level, is greater than the initial total transportation time with the lower level SO model ( T t o t a l s i n g l e m o d e S O = 50,358). The resilience level of the test PHCLTS involving an earthquake UEE is lowest, which indicates that an earthquake is the most destructive to the test PHCLTS.

5.2. Comparison of PHCLTS with and without Recovery Activities

On the basis of the lower level SO model, the resilience of the test PHCLTS with and without recovery activities is compared. A specific terrorist attack scenario is regarded as an example in which arcs 6, 14, and 16 are affected. Assume that no recovery activity will be taken, the demands of warehouses 4 and 5 cannot be met. In accordance with the previous definition in Section 5.1, the point resilience of the test PHCLTS is 0 under this circumstance. As illustrated in Table 7, when the government implements recovery activities, the point resilience of the test PHCLTS steadily increases as the recovery budget increases.
The implication is that the implementation of recovery activities has a positive impact on the point resilience of the test PHCLTS. For all scenarios of each UEE type, the number of scenarios with point resilience increases by 0–0.2, 0.2–0.4, 0.4–0.6, 0.6–0.8, and 0.8–1 compared to the corresponding point resilience of the test PHCLTS without recovery activities, are shown in Figure 8. For all three UEE types, the point resilience of the test PHCLTS with recovery activities has been effectively improved. Taking the earthquake as an example (Figure 8b), after the implementation of recovery activities with a budget of 2000 units, the point resilience of 56 and 34 scenarios increases 0.4–0.6 and 0.6–0.8, respectively; After the implementation of recovery activities with a budget of 6000 units, the point resilience of 95 scenarios increases 0.8–1, and there is indication of a significant improvement on resilience of the test PHCLTS. Overall, appropriate recovery activities are crucial to the resilience of the PHCLTS.
The implementation of appropriate recovery activities can improve the resilience of PHCLTSs effectively, reducing the negative influence caused by UEEs. The resilience of PHCLTSs will be at a low level if insufficient recovery activities are implemented. Consequently, the total transportation time of containers will increase, and even the transportation demands at the hinterlands will not be met. Improvement of the resilience of PHCLTSs through appropriate recovery activities is crucial for the economy and society.

5.3. Analysis of the Influence of Maximum Budget for Recovery Activities

For each UEE type, the maximum budget for recovery activities is changed to analyze its influence on the resilience of the PHCLTS. Figure 9 illustrates the resilience levels of the test PHCLTS on the basis of the lower level SO model. The resilience levels under different UEE types increase with budget growth, and this condition is particularly evident in the UEE type of freeze (Figure 9c). Under all the three UEE types, the improved speed of resilience levels slows down with the increase in budget. This result is consistent with “resilience follows the law of diminishing returns,” which is mentioned by Nair et al. [12].
For each UEE type, the variances in point resilience of the test PHCLTS decrease with the increase in recovery budgets (Table 8). This case is particularly evident for the UEE type of freeze. That is, the randomness of the negative effect on the test PHCLTS is gradually reduced as the budget increases. The implementation of immediate recovery activities helps the test PHCLTS absorb disturbances efficiently. Different UEE types will affect the test PHCLTS in various ways. The government should prepare an appropriate budget in accordance with possible future UEE types and improve the resilience of the PHCLTS by implementing appropriate recovery activities.
In one case, suppose the possibility of a freeze is predicted to be particularly high, the government should prepare a budget of 4000 units if it wants to achieve a resilience level above 0.8 (Figure 9c). In another case, suppose the possibility of an earthquake is predicted to be particularly high, the government should raise the budget to 5000 units if it wants to reach a resilience level above 0.8 (Figure 9b). The proposed bi-level programming models can help the government select a reasonable budget level to obtain a satisfactory resilience level in the face of UEEs.

5.4. Comparison between the Lower Level SO and UE Models

The influence of PHCLTS user behaviors on the resilience of the PHCLTS is discussed. Two different UEE scenarios are taken as examples (Budget = 4000). UEE scenario 1 is another terrorist attack scenario in which arcs 1, 13, 14, and 15 are affected. UEE scenario 2 is an earthquake scenario in which arc 10 is at the epicenter, and arcs 6–9 and 11–14 are also affected. The decision made by the government about recovery levels of the damaged arc capacity by implementing recovery activities and the corresponding decisions of the truck carriers about freight volume arrangement in the recovered PHCLTS with different lower level SO and UE models are shown in Table 9.
The following findings are obtained:
(1)
PHCLTS user behaviors will influence the recovery decisions of the government, thereby affecting the point resilience of the test PHCLTS. In UEE scenario 1, the recovery decisions of the government are different in the lower level SO and UE models. In UEE scenario 2, the recovery decisions of the government are the same in the lower level SO and UE models, but the freight volume decisions made by the truck carriers are different, which lead to a different point resilience of the test PHCLTS.
(2)
The point resilience of the test PHCLTS with lower level SO and UE models has no absolute advantage. For each UEE type, most of the point resilience of the test PHCLTS with the lower level UE model is higher (Table 10).
(3)
For all UEE scenarios, the total transportation time with the lower level SO model is less than or equal to the corresponding value with the lower level UE model (Table 10). This result shows that the truck carriers that cooperate fully to reach an SO assignment perform better in transportation time consumption than those that make decisions independently to reach a UE assignment when facing UEEs.
In sum, the behaviors of PHCLTS users will affect the resilience of a PHCLTS. The government’s optimal decisions about the recovery activities may be different, while the behavior modes of truck carriers also differ. When studying the resilience problem of PHCLTSs, the behaviors of PHCLTS users should be considered to achieve the goal of improving the resilience of PHCLTSs.

5.5. Comparison with Single-Level Programming Models

Several factors are included in the bi-level programming models. By comparison, owing to the limitations of single-level programming models, studying all of them at the same time is difficult. Two different single-level programming models based on the SO model are established. In single-level programming model 1, the topology of PHCLTS infrastructures and the goal of the government are considered, and the goal of truck carriers is ignored (see the mathematical expression in Appendix B). In single-level programming model 2, only the transportation capacity of the PHCLTS is focused on to solve its maximum flow problem to reach the maximum value of total transportation capacity, and the demands of shippers in hinterlands and behaviors of truck carriers are ignored (see the mathematical expression of in Appendix C). These two single-level programming models are solved by the TOMLAB/CONOPT solver, firstly, by determining the government’s decision about recovery activities by solving the single-level programming models and then by determining the actual freight volume on each arc according to the SO principle.
Another specific earthquake scenario is regarded as an example in which arc 16 is at the epicenter, and arcs 6 and 10–15 are also affected. This scenario includes the same recovery budget of 4000 units, the decision made by the government about recovery levels of damaged arc capacity by implementing recovery activities, the corresponding decisions of the truck carriers about freight volume arrangement in the recovered PHCLTS, and the point resilience of the test PHCLTS obtained in different models as shown in Table 11.
The point resilience obtained using the bi-level programming model is better than those obtained using the two single-level programming models. When a single-level programming model is employed to solve this problem, the recovery activities taken by the government cannot correctly predict the behaviors of PHCLTS users and their influence on the PHCLTS resilience. This condition will lead to ineffectively improved PHCLTS resilience. Taking the result of single-level programming model 1 in Table 11 as an example, the damaged capacity of arc 16 is recovered by 72%. However, the actual freight volume on arc 16 is 0; i.e., to reach an SO assignment, the truck carriers do not arrange any container to transport through arc 16, and the recovery resources spent on arc 16 are meaningless.
For all UEE scenarios involving an earthquake, the point resilience of the test PHCLTS obtained using the bi-level programming model is the best (Figure 10). For the results of single-level programming model 2, the point resilience of the test PHCLTS under more than 80% of UEE scenarios is 0. The reason is that the container demands of shippers are ignored in single-level programming model 2. Consequently, the transportation process of containers cannot be completed in the recovered PHCLTS via this model. When studying the resilience problem of PHCLTSs, not only the goals and behaviors of truck carriers should be considered, but also the demands of shippers, who are another transportation PHCLTS user.
To sum up, single-level programming models cannot help the government implement the optimum recovery decisions to restore PHCLTS capacity. Bi-level programming models can support the government to formulate efficient decisions to improve the resilience of a PHCLTS in the face of UEEs.

5.6. Analysis of the Influence of Demand in Hinterland Points

From a long-term perspective, with economic development, daily transportation demands will increase. For the UEE type involving an earthquake, the resilience levels of the test PHCLTS under different demand levels with a recovery budget of 4000 units are shown in Figure 11.
With the continuous growth of the demand level, the decreasing rate of the resilience level of the test PHCLTS shows an accelerating trend. The increasing demand makes the freight flows in the PHCLTS considerably saturated, i.e., the PHCLTS redundancy, which refers to the remaining capacity that can be provided by the test PHCLTS aside from completing the transportation demand [14], decreases continuously. Once a UEE scenario occurs, less redundancy capacity can be employed to deal with it. Therefore, the government can adjust the recovery budget and the decisions of recovery activities in time by using the proposed bi-level programming model in accordance with the changes in demand level. When the daily transportation demand level has increased by 30%, the government should raise the budget to 6000 units to achieve a resilience level no less than the resilience level (0.77) under the original demand level (Figure 12).

6. Managerial Insights and Conclusions

This study focuses on improving the time resilience of PHCLTSs. Managerial insights are obtained from numerical studies as follows:
(1)
After the occurrence of UEEs, the implementation of appropriate recovery activities can effectively improve the resilience of PHCLTSs and reduce the negative influence caused by UEEs. When insufficient recovery activities are implemented, the resilience of PHCLTSs will be at a low level. Consequently, the total transportation time of containers will increase, and even the transportation demand in the hinterlands cannot be met, leading to negative effects on the economy and society. Appropriate recovery activities are crucial to reduce the negative effects after a UEE occurrence.
(2)
The randomness of the negative effect on the test PHCLTS is gradually reduced as the budget increases. The resilience of PHCLTSs follows the law of diminishing returns. The government should prepare an appropriate budget in accordance with possible future UEE types.
(3)
PHCLTS user behaviors will influence the recovery decisions of the government, thereby affecting the resilience of PHCLTSs. Although facing the same UEE situation, the behavior patterns of PHCLTS users are different, and the government’s decision on the recovery activities also differs. The behaviors of PHCLTS users should be considered when studying this issue to improve the resilience of PHCLTSs.
(4)
From a long-term perspective, with the continuous growth in demand level, the resilience level of a PHCLTS shows a decreasing trend because the redundancy in the PHCLTS decreases. The government should adjust the recovery budget and the decisions on recovery activities in accordance with the changes in demand level to achieve a satisfactory resilience level of the PHCLTS.
(5)
The government should correctly predict the behaviors of PHCLTS users and their influence on the resilience of PHCLTSs when making decisions about recovery activities after UEEs. Otherwise, the resilience of PHCLTSs will be ineffectively improved, and resources will not be well utilized. The proposed bi-level programming models in this study can support the government in selecting a reasonable budget level and implementing efficient recovery decisions to improve the resilience of PHCLTSs in the face of UEEs.
In this study, a medium- and long-term perspective, which focuses on the time resilience of logistics transportation systems under multiple transportation periods after UEEs, is adopted. A bi-level programming approach is applied to analyze the resilience problem of PHCLTSs. The topological structure of a PHCLTS infrastructure, the goals and behaviors of the PHCLTS planner and users, and their interactions are involved in the bi-level programming model. In consideration of the characteristics of the bi-level programming model, an algorithm combining the PSO algorithm and a traditional optimization algorithm realized using the TOMLAB/CONOPT solver is developed to solve the model.
The proposed bi-level programming models considering the hierarchical relationship and interaction among different participants can support the government in selecting a reasonable budget level and implementing efficient recovery decisions to improve the resilience of PHCLTSs in the face of UEEs. This study analyzes how to improve the resilience of PHCLTSs under UEE situations. It expands the prearranged planning about UEEs and provides a decision-making reference and scientific theoretical support for efficient response to UEEs.
Several implications are available for future research. First, only single mode PHCLTSs were studied, which can be extended to multimodal PHCLTSs. Second, the pre-disaster activities can be studied to obtain a PHCLTS with higher resilience. Third, algorithms for the bi-level programming model can be further developed to improve solving efficiency and validity.

Author Contributions

Conceptualization, S.G. and N.L.; methodology, S.G. and N.L.; software, S.G.; writing—original draft preparation, S.G.; writing—review and editing, N.L.; supervision, N.L.; funding acquisition, N.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 71471162), NSFC-Zhejiang Joint Fund for the Integration of Industrialization and Informatization (Grant No. U1509221), Zhejiang University Social Science Special Fund for Risk Identification and Security System of Port Maritime Logistics under the ‘One Belt One Road’ (Grant No. 102000-541903/111), and Major Project of National Social Science Foundation of China (Grant No. 20&ZD129).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data supporting the reported results in the present study will be available on request from the corresponding author or the first author.

Conflicts of Interest

No potential conflict of interest was reported by the authors.

Appendix A. Settings of Different UEE Types

In the numerical studies, three subtypes with the same occurrence probability are set for each UEE type. The information about how the UEE scenarios of different UEE types affect the infrastructures of PHCLTS are shown in Table A1 (Terrorist attack), Table A2 (Earthquake), and Table A3 (Freeze).
Table A1. Influence of terrorist attack on PHCLTS.
Table A1. Influence of terrorist attack on PHCLTS.
Detail a
Subtype 1 of terrorist attack3 of 16 arcs in the PHCLTS are randomly chosen, and their capacity drops to 0.
Subtype 2 of terrorist attack4 of 16 arcs in the PHCLTS are randomly chosen, and their capacity drops to 0.
Subtype 3 of terrorist attack5 of 16 arcs in the PHCLTS are randomly chosen, and their capacity drops to 0.
a Each arc has an equal probability of being selected.
Table A2. Influence of earthquake on PHCLTS.
Table A2. Influence of earthquake on PHCLTS.
Detail
Subtype 1 of earthquakeArc 9 is located at the epicenter, and the decrease percentage of its capacity is a uniform random number in the interval [0.9, 1]. Arcs 2, 5–8, 10–13 are affected by the earthquake, and the decrease percentages of their capacity are uniform random numbers in the interval [0.6, 0.8].
Subtype 2 of earthquakeArc 10 is located at the epicenter, and the decrease percentage of its capacity is a uniform random number in the interval [0.9, 1]. Arcs 6–8, 11–14 are affected by the earthquake, and the decrease percentages of their capacity are uniform random numbers in the interval [0.6, 0.8].
Subtype 3 of earthquakeArc 16 is located at the epicenter, and the decrease percentage of its capacity is a uniform random number in the interval [0.9, 1]. Arcs 6, 10–15 are affected by the earthquake, and the decrease percentages of their capacity are uniform random numbers in the interval [0.6, 0.8].
Table A3. Influence of freeze on PHCLTS.
Table A3. Influence of freeze on PHCLTS.
Detail
Subtype 1 of freeze1 warehouse in the PHCLTS is randomly chosen, and the decrease percentages of capacity of the arcs connected to this warehouse are uniform random numbers in the interval [0.7, 1].
Subtype 2 of freeze2 adjacent warehouses in the PHCLTS are randomly chosen, and the decrease percentages of capacity of the arcs connected to these warehouses are uniform random numbers in the interval [0.7, 1].
Subtype 3 of freeze 3 adjacent warehouses in the PHCLTS are randomly chosen, and the decrease percentages of capacity of the arcs connected to these warehouses are uniform random numbers in the interval [0.7, 1].

Appendix B. Single-Level Programming Model 1

In single-level programming model 1, the topology of the PHCLTS infrastructures and the goal of the government are considered. The model is described as follows:
Single-level programming model 1
min   T t o t a l s i n g l e m o d e ( ξ ) = a A h a ( ξ ) · t a ( ξ )
a b a · l a ( ξ ) · ( c a c a u ( ξ ) ) c a B
0 l a ( ξ ) 1 ,   if   a A u
l a ( ξ ) = 0 ,   if   a A n u
a A w h a ( ξ ) = D w
a ( : , k ) A h a ( ξ ) a ( k , : ) A h a ( ξ ) = 0
h a ( ξ ) c a ( ξ )
h a ( ξ ) 0
a A ,   k K
The objective function (A1) is the same as the objective function (3) in Section 3.3.1, which seeks to minimize total transportation time post-disaster to improve the resilience of the PHCLTS under each UEE scenario. The meaning of Constraints (A2)–(A8) is the same as the Constraints (4)–(6) and (10)–(13) in Section 3.3.1 and Section 3.3.2, respectively.

Appendix C. Single-Level Programming Model 2

Single-level programming model 2 only focuses on the transportation capacity of the PHCLTS. It is a maximum flow problem of the PHCLTS. The model is described as follows:
Single-level programming model 2
max a A h a ( ξ )
a b a · l a ( ξ ) · ( c a c a u ( ξ ) ) c a B
0 l a ( ξ ) 1 ,   if   a A u
l a ( ξ ) = 0 ,   if   a A n u
a ( : , k ) A h a ( ξ ) a ( k , : ) A h a ( ξ ) = 0
h a ( ξ ) c a ( ξ )
h a ( ξ ) 0
a A ,   k K
Objective function (A9) seeks to maximize the total transportation capacity of the PHCLTS. The meaning of Constraints (A10)–(A15) is the same as the Constraints (4)–(6) and (11)–(13) in Section 3.3.1 and Section 3.3.2, respectively.

References

  1. World Trade Organization. World Trade Report. 2019. Available online: https://www.wto.org/english/res_e/publications_e/wtr19_e.htm (accessed on 16 December 2020).
  2. Mansouri, M.; Sauser, B.; Boardman, J. Applications of systems thinking for resilience study in Maritime Transportation System of Systems. In Proceedings of the 3rd Annual IEEE International Systems Conference, Vancouver, BC, Canada, 23–26 March 2009; pp. 211–217. [Google Scholar]
  3. Chen, H.; Lam, J.S.L.; Liu, N. Strategic investment in enhancing port–hinterland container transportation network resilience: A network game theory approach. Transp. Res. Part B Methodol. 2018, 111, 83–112. [Google Scholar] [CrossRef]
  4. Chen, L.; Miller-Hooks, E. Resilience: An Indicator of Recovery Capability in Intermodal Freight Transport. Transp. Sci. 2012, 46, 109–123. [Google Scholar] [CrossRef]
  5. Ta, C.; Goodchild, A.V.; Pitera, K. Structuring a Definition of Resilience for the Freight Transportation System. Transp. Res. Rec. J. Transp. Res. Board 2009, 2097, 19–25. [Google Scholar] [CrossRef]
  6. Bard, J.F. Practical Bilevel Optimization: Algorithms and Applications; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1998. [Google Scholar]
  7. Talley, W.K.; Ng, M. Hinterland transport chains: A behavioral examination approach. Transp. Res. Part E Logist. Transp. Rev. 2018, 113, 94–98. [Google Scholar] [CrossRef]
  8. Bruneau, M.; Chang, S.E.; Eguchi, R.T.; Lee, G.C.; O’Rourke, T.D.; Reinhorn, A.M.; Shinozuka, M.; Tierney, K.; Wallace, W.A.; Von Winterfeldt, D. A Framework to Quantitatively Assess and Enhance the Seismic Resilience of Communities. Earthq. Spectra 2003, 19, 733–752. [Google Scholar] [CrossRef] [Green Version]
  9. Murray-Tuite, P.M. A Comparison of Transportation Network Resilience under Simulated System Optimum and User Equilibrium Conditions. In Proceedings of the Winter Simulation Conference, Monterey, CA, USA, 3–6 December 2006. [Google Scholar]
  10. Omer, M.; Mostashari, A.; Nilchiani, R.; Mansouri, M. A framework for assessing resiliency of maritime transportation systems. Marit. Policy Manage. 2012, 39, 1–19. [Google Scholar] [CrossRef]
  11. Faturechi, R.; Miller-Hooks, E. Travel time resilience of roadway networks under disaster. Transp. Res. Part B Methodol. 2014, 70, 47–64. [Google Scholar] [CrossRef]
  12. Nair, R.; Avetisyan, H.; Miller-Hooks, E. Resilience Framework for Ports and Other Intermodal Components. Transp. Res. Rec. J. Transp. Res. Board 2010, 2166, 54–65. [Google Scholar] [CrossRef]
  13. Miller-Hooks, E.; Zhang, X.; Faturechi, R. Measuring and maximizing resilience of freight transportation networks. Comput. Oper. Res. 2012, 39, 1633–1643. [Google Scholar] [CrossRef]
  14. Chen, H.; Cullinane, K.; Liu, N. Developing a model for measuring the resilience of a port-hinterland container transportation network. Transp. Res. Part E Logist. Transp. Rev. 2017, 97, 282–301. [Google Scholar] [CrossRef]
  15. Friesz, T.L.; Harker, P.T. Freight network equilibrium: A review of the state of the art. In Analytical Studies in Transport Economics; Cambridge University Press: Cambridge, UK, 1985; pp. 161–206. [Google Scholar]
  16. Jones, D.A.; Farkas, J.L.; Bernstein, O.; Davis, C.E.; Turk, A.; Turnquist, M.A.; Nozick, L.K.; Levine, B.; Rawls, C.G.; Ostrowski, S.D.; et al. import/export container flow modeling and disruption analysis. Res. Transp. Econ. 2011, 32, 3–14. [Google Scholar] [CrossRef]
  17. Meng, Q.; Wang, X.C. Intermodal hub-and-spoke network design: Incorporating multiple stakeholders and multi-type containers. Transp. Res. Part B Methodol. 2011, 45, 724–742. [Google Scholar] [CrossRef]
  18. Corman, F.; Viti, F.; Negenborn, R.R. Equilibrium models in multimodal container transport systems. Flex. Serv. Manuf. J. 2017, 29, 125–153. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, B.; Yao, T.; Friesz, T.L.; Hongcheng, L. Urban Freight Transportation Planning: A Dynamic Stackelberg Game-Theoretic Approach. 2012. Available online: https://arxiv.org/abs/1211.3950 (accessed on 28 September 2019).
  20. Jiang, J.; Zhang, D.; Meng, Q.; Liu, Y. Regional multimodal logistics network design considering demand uncertainty and CO2 emission reduction target: A system-optimization approach. J. Clean. Prod. 2020, 248, 119304. [Google Scholar] [CrossRef]
  21. Yu, S.; Jiang, Y. Network design and delivery scheme optimisation under integrated air-rail freight transportation. Int. J. Logist. Res. Appl. 2021, 1–17. [Google Scholar] [CrossRef]
  22. Moreno-Quintero, E.; Fowkes, T.; Watling, D. Modelling Planner-Carrier Interactions in Road Freight Transport: Optimization of Road Maintenance Costs Via Overloading Control. Transp. Res. Part E Logist. Transp. Rev. 2013, 50, 68–83. [Google Scholar] [CrossRef] [Green Version]
  23. Lee, H.; Song, Y.; Choo, S.; Chung, K.-Y.; Lee, K.-D. Bi-level optimization programming for the shipper-carrier network problem. Clust. Comput. 2014, 17, 805–816. [Google Scholar] [CrossRef]
  24. Qiu, R.; Xu, J.; Ke, R.; Zeng, Z.; Wang, Y. Carbon pricing initiatives-based bi-level pollution routing problem. Eur. J. Oper. Res. 2020, 286, 203–217. [Google Scholar] [CrossRef]
  25. Li, S.; Liang, Y.; Wang, Z.; Zhang, D. An Optimization Model of a Sustainable City Logistics Network Design Based on Goal Programming. Sustainability 2021, 13, 7418. [Google Scholar] [CrossRef]
  26. De Jong, G.; Gunn, H.; Walker, W. National and International Freight Transport Models: An Overview and Ideas for Future Development. Transp. Rev. 2004, 24, 103–124. [Google Scholar] [CrossRef] [Green Version]
  27. Sheffi, Y. Urban Transportation Networks: Equilibrium Analysis with Mathematical Programming Methods; Prentice-Hall, Inc.: Hoboken, NJ, USA, 1985. [Google Scholar]
  28. Grey, A. The generalised cost dilemma. Transportation 1978, 7, 261–280. [Google Scholar] [CrossRef]
  29. Jeroslow, R.G. The polynomial hierarchy and a simple model for competitive analysis. Math. Program. 1985, 32, 146–164. [Google Scholar] [CrossRef]
  30. Assadipour, G.; Ke, G.Y.; Verma, M. A toll-based bi-level programming approach to managing hazardous materials shipments over an intermodal transportation network. Transp. Res. Part D Transp. Environ. 2016, 47, 208–221. [Google Scholar] [CrossRef]
  31. Soares, I.; Alves, M.J.; Antunes, C.H. Designing time-of-use tariffs in electricity retail markets using a bi-level model–Estimating bounds when the lower level problem cannot be exactly solved. Omega 2020, 93, 102027. [Google Scholar] [CrossRef]
  32. Zhang, S.; Chen, N.; She, N.; Li, K. Location optimization of a competitive distribution center for urban cold chain logistics in terms of low-carbon emissions. Comput. Ind. Eng. 2021, 154, 107120. [Google Scholar] [CrossRef]
  33. Kennedy, J.; Eberhart, R. Particle Swarm Optimization. In Proceedings of the ICNN’95-International Conference on Neural Networks, Perth, WA, Australia, 27 November–1 December 1995; pp. 1942–1948. [Google Scholar]
  34. Shi, Y. A Modified Particle Swarm Optimizer. In Proceedings of the IEEE Icec Conference, Anchorage, AK, USA, 4–9 May 1998. [Google Scholar]
  35. Zhao, Z.; Gu, X.; Li, T. Particle Swarm Optimization for Bi-level Programming Problem. Syst. Eng. Theory Pract. 2007, 27, 92–98. (In Chinese) [Google Scholar]
  36. Holmström, K.; Göran, A.O.; Edvall, M.M. User’s Guide for TOMLAB/CONOPT. 2007. Available online: https://tomopt.com/tomlab/products/conopt/ (accessed on 1 March 2021).
  37. Boyd, S.; Vandenberghe, L. Convex Optimization; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  38. López, L.F.M.; Blas, N.G.; Albert, A.A. Multidimensional knapsack problem optimization using a binary particle swarm model with genetic operations. Soft Comput. 2018, 22, 2567–2582. [Google Scholar] [CrossRef]
  39. Alves, M.J.; Antunes, C.H.; Carrasqueira, P. A Hybrid Genetic Algorithm for the Interaction of Electricity Retailers with Demand Response. In Applications of Evolutionary Computation 2016; Squillero, G., Burelli, P., Eds.; Springer: Berlin/Heidelberg, Germany, 2016; pp. 459–474. [Google Scholar]
Figure 1. Infrastructure structure of a PHCLTS.
Figure 1. Infrastructure structure of a PHCLTS.
Sustainability 14 00180 g001
Figure 2. Decision framework of the bi-level programming model.
Figure 2. Decision framework of the bi-level programming model.
Sustainability 14 00180 g002
Figure 3. Solution procedure of the algorithm.
Figure 3. Solution procedure of the algorithm.
Sustainability 14 00180 g003
Figure 4. Illustration of particle encoding and iterative process. (a) A particle and the relationship with the solutions of upper and lower level models. (b) Iterative and repair process. (c) Evaluation process.
Figure 4. Illustration of particle encoding and iterative process. (a) A particle and the relationship with the solutions of upper and lower level models. (b) Iterative and repair process. (c) Evaluation process.
Sustainability 14 00180 g004
Figure 5. Test PHCLTS.
Figure 5. Test PHCLTS.
Sustainability 14 00180 g005
Figure 6. Illustrations of each UEE type. (a) Terrorist attack. (b) Earthquake (c) Freeze.
Figure 6. Illustrations of each UEE type. (a) Terrorist attack. (b) Earthquake (c) Freeze.
Sustainability 14 00180 g006aSustainability 14 00180 g006b
Figure 7. Convergence of the algorithm in different UEE scenarios.
Figure 7. Convergence of the algorithm in different UEE scenarios.
Sustainability 14 00180 g007
Figure 8. Improvement in point resilience of the test PHCLTS with recovery activities. (a) Terrorist attack. (b) Earthquake. (c) Freeze.
Figure 8. Improvement in point resilience of the test PHCLTS with recovery activities. (a) Terrorist attack. (b) Earthquake. (c) Freeze.
Sustainability 14 00180 g008aSustainability 14 00180 g008b
Figure 9. Resilience levels of PHCLTS under different recovery budgets. (a) Terrorist attack. (b) Earthquake (c) Freeze.
Figure 9. Resilience levels of PHCLTS under different recovery budgets. (a) Terrorist attack. (b) Earthquake (c) Freeze.
Sustainability 14 00180 g009
Figure 10. Comparison with single-level programming models (Earthquake).
Figure 10. Comparison with single-level programming models (Earthquake).
Sustainability 14 00180 g010
Figure 11. Resilience levels of PHCLTS under different demand levels.
Figure 11. Resilience levels of PHCLTS under different demand levels.
Sustainability 14 00180 g011
Figure 12. Resilience levels under different recovery budgets (30% increase in demand level).
Figure 12. Resilience levels under different recovery budgets (30% increase in demand level).
Sustainability 14 00180 g012
Table 1. Notation description for the proposed models.
Table 1. Notation description for the proposed models.
NotationsMeaning
Sets
W Set   of   all   warehouses   in   hinterland ,   w W
KSet of roadway intersections, i.e., set of nodes in a PHCLTS, except for the port and warehouses in hinterland, k K
A Set   of   all   arcs   in   a   PHCLTS ,   a A
AwSet of arcs connected to warehouse w, AwA
AuSet of arcs affected by UEEs, AuA
AnuSet of arcs unaffected by UEEs, AnuA
U Set   of   UEE   types ,   u U
Indices
spA port; only one port exists in a PHCLTS
wA warehouse in hinterland, which represents the demand area
kA roadway intersection
aAn arc, which represents one road of the PHCLTS. a(sp, k1) denotes an arc with origin sp and destination k1
ξ A UEE scenario
Deterministic parameters
T t o t a l s i n g l e m o d e Total transportation time of the PHCLTS in each transportation period pre-disaster
h a Freight volume on arc a pre-disaster
t a Transportation time on arc a pre-disaster
ba Cost of recovering 100% of the original capacity c a on arc a post-disaster. If the capacity of arc a recovered is x % · c a , the recovery cost on arc a will be x % · b a .
BMaximum budget for recovery activities
c a Original transportation capacity of arc a
t a 0 Free-flow travel time on arc a pre-disaster, which is the transportation time on arc a without road impedance
chuInfluence parameter of UEE type u on the free-flow travel time of arcs
G C t o t a l s i n g l e m o d e Total generalized transportation costs of truck carriers in each transportation period pre-disaster. G C t o t a l s i n g l e m o d e S O refers to the total generalized transportation costs pre-disaster following the SO principle; G C t o t a l s i n g l e m o d e U E denotes the total generalized transportation costs pre-disaster following the UE principle
raTransportation cost on arc a
qValue of time
α and βParameters of “Bureau of Public Roads function”, which capture the influence of freight flow on transportation time of arc a
DwDaily transportation demand for warehouse w
UEE related variables (under UEE scenario ξ)
T t o t a l s i n g l e m o d e ( ξ ) Total transportation time of the PHCLTS in each transportation period post-disaster and after recovery
t a ( ξ ) Transportation time on arc a post-disaster and after recovery
c a u ( ξ ) Capacity of arc a post-disaster and before recovery
c a ( ξ ) Capacity of arc a post-disaster and after recovery
t a r 0 ( ξ ) Free-flow travel time on arc a post-disaster and after recovery
G C t o t a l s i n g l e m o d e ( ξ ) Total generalized transportation costs of truck carriers in each transportation period post-disaster and after recovery. G C t o t a l s i n g l e m o d e S O ( ξ ) refers to the total generalized transportation costs post-disaster and after recovery following the SO principle; G C t o t a l s i n g l e m o d e U E ( ξ ) denotes the total generalized transportation costs post-disaster and after recovery following the UE principle
Decision variables
In the upper level model
l a ( ξ ) Recovery levels of damaged capacity that the government implements on arc a under UEE scenario ξ ( ( 0 l a ( ξ ) 1 : la = 1 means that the capacity of arc a is recovered to its original capacity, and la = 0 means that no recovery activity is implemented on arc a)
In the lower level model
h a ( ξ ) Freight volume on arc a post-disaster and after recovery under UEE scenario ξ
Table 2. Notation description for the PSO algorithm.
Table 2. Notation description for the PSO algorithm.
NotationsMeaning
nSwarm size, i.e., the number of particles
Pi Position   of   particle   i ,   i [ 1 ,   n ]
vi Velocity   of   particle   i ,   i [ 1 ,   n ]
pbestParticle’s best position
gbestWhole swarm’s best position
c1 and c2Learning factors
r1 and r2Uniform random numbers in the interval [0, 1]
wInertia weight
wMaxMaximum inertia weight
wMinMinimum inertia weight
iterNumber of iterations
iterMaxMaximum number of iterations
ε Tolerance
Table 3. Evaluation criteria for different particles.
Table 3. Evaluation criteria for different particles.
Value of Fitness Function
f i t n e s s ( P i , H i ) < f i t n e s s ( P j , H j ) f i t n e s s ( P i , H i ) f i t n e s s ( P j , H j )
(Pi, Hi)/(Pj, Hj)
whether or not to be
a feasible solution of
the bi-level programming model
Feasible/Feasible(Pi, Hi) is better(Pj, Hj) is better
Feasible/Infeasible(Pi, Hi) is better(Pi, Hi) is better
Infeasible/Feasible(Pj, Hj) is better(Pj, Hj) is better
Infeasible/Infeasible(Pi, Hi) is better(Pj, Hj) is better
Table 4. Arc parameters of the test PHCLTS.
Table 4. Arc parameters of the test PHCLTS.
Arc IDOriginDestination c a t a 0 raba
1PortWarehouse 11200111502400
2PortNode 116005603200
3Node 1Warehouse 16006801200
4Node 1Warehouse 210005702000
5Node 1Warehouse 39006601800
6PortNode 224005604800
7Node 2Warehouse 19007701800
8Node 2Warehouse 26006601200
9Node 2Warehouse 36005701200
10Node 2Warehouse 46006601200
11Node 2Warehouse 58007801600
12PortNode 318006403600
13Node 3Warehouse 310006702000
14Node 3Warehouse 410005802000
15Node 3Warehouse 58007901600
16PortWarehouse 51300131402600
Table 5. Characteristics of UEEs.
Table 5. Characteristics of UEEs.
UEE TypeDetail
Terrorist attackSome arcs in the PHCLTS are randomly chosen, and their capacity drops to 0.
EarthquakeAn area is randomly chosen, and the capacity of the arcs in this area decreases.
Freeze aSome nodes are randomly chosen, and the capacity of the arcs connected to these nodes decreases.
a is snow and ice weather wherein a mass of snow and ice covers the roads and negatively influences transportation.
Table 6. Resilience levels and average total transportation time of the test PHCLTS.
Table 6. Resilience levels and average total transportation time of the test PHCLTS.
Resilience Level (Average Total Transportation Time)
Terrorist AttackEarthquakeFreeze
Lower level model 1—SO0.9164 (54,213)0.7761 (65,514)0.9244 (55,162)
Lower level model 2—UE0.9174 (54,635)0.7787 (65,921)0.9249 (55,652)
Table 7. Comparison of PHCLTS with and without recovery activities.
Table 7. Comparison of PHCLTS with and without recovery activities.
BudgetRecovery Levels of Damaged Arc CapacityTOTAL
Transportation Time
Point
Resilience
Arc 6Arc 14Arc 16
With recovery activities 60001.00 0.60 053,7410.9370
40000.83 0058,0440.8676
20000.42 0068,5640.7345
Without recovery activity-----0
Table 8. Variances of the point resilience under different recovery budgets.
Table 8. Variances of the point resilience under different recovery budgets.
Recovery Budget
20003000400050006000
Terrorist attack0.10270.06650.02210.02060.0027
Earthquake0.03740.00750.00620.00380.0017
Freeze0.17220.09870.00910.00290.0012
Table 9. Comparison between the lower level SO and UE models.
Table 9. Comparison between the lower level SO and UE models.
Lower Level ModelRecovery Levels of
Damaged Arc Capacity a
(Decision of the Government)
Freight Volume on Each Arc
(Decisions of the Truck Carriers)
Total
Transportation Time
(Upper Level)
Point
Resilience
Total
Generalized Costs
(Lower Level)
UEE scenario 1SO modelArc 1Arc 2Arc 3Arc 4Arc 1Arc 2Arc 3Arc 452,5180.95891,058,137
1---5121001208328
Arc 5Arc 6Arc 7Arc 8Arc 5Arc 6Arc 7Arc 8
----4651516280172
Arc 9Arc 10Arc 11Arc 12Arc 9Arc 10Arc 11Arc 12
----335243486257
Arc 13Arc 14Arc 15Arc 16Arc 13Arc 14Arc 15Arc 16
00.80-02570714
UE modelArc 1Arc 2Arc 3Arc 4Arc 1Arc 2Arc 3Arc 453,1260.95651,051,428
0---01133404287
Arc 5Arc 6Arc 7Arc 8Arc 5Arc 6Arc 7Arc 8
----4421580595213
Arc 9Arc 10Arc 11Arc 12Arc 9Arc 10Arc 11Arc 12
----35876338858
Arc 13Arc 14Arc 15Arc 16Arc 13Arc 14Arc 15Arc 16
011-0425433429
UEE scenario 2SO modelArc 1Arc 2Arc 3Arc 4Arc 1Arc 2Arc 3Arc 456,6740.88861,098,601
----7801520220500
Arc 5Arc 6Arc 7Arc 8Arc 5Arc 6Arc 7Arc 8
-000800000
Arc 9Arc 10Arc 11Arc 12Arc 9Arc 10Arc 11Arc 12
0001000994
Arc 13Arc 14Arc 15Arc 16Arc 13Arc 14Arc 15Arc 16
01--0500494706
UE modelArc 1Arc 2Arc 3Arc 4Arc 1Arc 2Arc 3Arc 456,8810.89331,099,634
----7281572272500
Arc 5Arc 6Arc 7Arc 8Arc 5Arc 6Arc 7Arc 8
-000800000
Arc 9Arc 10Arc 11Arc 12Arc 9Arc 10Arc 11Arc 12
00010001046
Arc 13Arc 14Arc 15Arc 16Arc 13Arc 14Arc 15Arc 16
01--0500546654
a “-” means this arc is not affected by UEE scenario, and no recovery activity will be taken on it.
Table 10. Point resilience and total transportation time with lower level SO and UE models.
Table 10. Point resilience and total transportation time with lower level SO and UE models.
Higher Point Resilience Equal Point
Resilience
Lower Total Transportation TimeEqual Total
Transportation Time
Lower Level
SO Model
Lower Level
UE Model
Lower Level
SO Model
Lower Level
UE Model
Terrorist attack a22%76%096%02%
Earthquake11%89%0100%00
Freeze9%79%12% b100%00
a Two percent of terrorist attack scenarios with point resilience equal to 0 are not included. b The point resilience with lower level SO and UE models is equal to 1 under these 12% of freeze scenarios; i.e., the test PHCLTS is recovered to the original status.
Table 11. Comparison with single-level programming models.
Table 11. Comparison with single-level programming models.
Recovery Levels of
Damaged Arc Capacity a
(Decision of the Government)
Freight Volume on Each Arc
(Decisions of the Truck Carriers)
Total
Transportation Time
(Upper Level)
Point
Resilience
Bi-level
programming model
Arc 1Arc 2Arc 3Arc 4Arc 1Arc 2Arc 3Arc 467,1470.7500
----7801520220500
Arc 5Arc 6Arc 7Arc 8Arc 5Arc 6Arc 7Arc 8
-0.29--80050000
Arc 9Arc 10Arc 11Arc 12Arc 9Arc 10Arc 11Arc 12
-100050000
Arc 13Arc 14Arc 15Arc 16Arc 13Arc 14Arc 15Arc 16
00010001200
Single-level programming model 1 Arc 1Arc 2Arc 3Arc 4Arc 1Arc 2Arc 3Arc 487,6360.5746
----7801520220500
Arc 5Arc 6Arc 7Arc 8Arc 5Arc 6Arc 7Arc 8
-0.47--800130000
Arc 9Arc 10Arc 11Arc 12Arc 9Arc 10Arc 11Arc 12
-1110500800400
Arc 13Arc 14Arc 15Arc 16Arc 13Arc 14Arc 15Arc 16
0000.72004000
Single-level programming model 2 Arc 1Arc 2Arc 3Arc 4Arc 1Arc 2Arc 3Arc 4119,7000.4207
----7311346269500
Arc 5Arc 6Arc 7Arc 8Arc 5Arc 6Arc 7Arc 8
-0.64--57775000
Arc 9Arc 10Arc 11Arc 12Arc 9Arc 10Arc 11Arc 12
-000.24223230297379
Arc 13Arc 14Arc 15Arc 16Arc 13Arc 14Arc 15Arc 16
0000.630270109794
a “-” means this arc is not affected by UEE scenario, and no recovery activity will be taken on it.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gao, S.; Liu, N. Improving the Resilience of Port–Hinterland Container Logistics Transportation Systems: A Bi-Level Programming Approach. Sustainability 2022, 14, 180. https://doi.org/10.3390/su14010180

AMA Style

Gao S, Liu N. Improving the Resilience of Port–Hinterland Container Logistics Transportation Systems: A Bi-Level Programming Approach. Sustainability. 2022; 14(1):180. https://doi.org/10.3390/su14010180

Chicago/Turabian Style

Gao, Song, and Nan Liu. 2022. "Improving the Resilience of Port–Hinterland Container Logistics Transportation Systems: A Bi-Level Programming Approach" Sustainability 14, no. 1: 180. https://doi.org/10.3390/su14010180

APA Style

Gao, S., & Liu, N. (2022). Improving the Resilience of Port–Hinterland Container Logistics Transportation Systems: A Bi-Level Programming Approach. Sustainability, 14(1), 180. https://doi.org/10.3390/su14010180

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