Next Article in Journal
Third-Order Differential Subordinations Using Fractional Integral of Gaussian Hypergeometric Function
Next Article in Special Issue
Condition-Based Maintenance Optimization Method Using Performance Margin
Previous Article in Journal
Some Properties for Subordinations of Analytic Functions
Previous Article in Special Issue
A Game—Theoretic Model for a Stochastic Linear Quadratic Tracking Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Uncertain Programming Model for the Cross-Border Multimodal Container Transport System Based on Inland Ports

School of Traffic and Transportation, Beijing Jiaotong University, Beijing 100044, China
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(2), 132; https://doi.org/10.3390/axioms12020132
Submission received: 23 December 2022 / Revised: 24 January 2023 / Accepted: 26 January 2023 / Published: 28 January 2023
(This article belongs to the Special Issue Advances in Uncertain Optimization and Applications)

Abstract

:
The importance of inland ports in promoting current cross-border trade is increasingly recognized. In this work, we aim to design the entire network for the cross-border multimodal container transport system based on inland ports. Unlike previous studies, we consider strong uncertainty in cross-border transportation demand to be caused by a variety of realistic factors such as the global economic situation, trade policies among countries, and global epidemics, etc. To handle the demand uncertainty, we develop an uncertain programming model for the considered cross-border multimodal container transportation network design problem to minimize the expectation of the total costs, including carbon emissions, by imposing two types of chance constraints for capacity limitations. Under mild assumptions, we further convert the proposed uncertain model into its equivalent deterministic one, which can be solved by off-the-shelf solvers such as CPLEX, Gurobi, and Lingo. Finally, we illustrate the applicability of the proposed model by taking the Huaihai Economic Zone-Europe multimodal container transport system as a real-world case study. The computational results provide valuable suggestions and policy guidance regarding four issues: the inland port locations, the transportation route choices, the strategies for reducing the total cost, and the schemes for improving network performance against uncertain demand.

1. Introduction

Inland ports are the key nodes for inland cities to connect to foreign countries. In Europe, the inland ports are typically recognized as an extension of seaport functions in inland areas [1], whereas in China, the inland ports are more frequently recognized as independent logistics hubs, relying on multiple transportation modes with the China-Europe Railway Express (CRE) as the core to achieve direct channel connections with foreign hubs [2]. In contrast to conventional logistics parks, the inland ports in China are marked by larger handling capacities, more transportation modes, the ability to open CRE, customs functions in inland cities, and clear economies of scale [2]. According to the China Inland Port Development Report [3], the number of inland ports in China has risen to 203 by the end of 2021. From this, the function of inland ports in promoting international trade flow is gaining prominence.
With the implementation of the Belt and Road Initiative, cross-border trade in Chinese cities has seen rapid growth [4]. From 2013 to 2021, China’s cumulative import and export amount of products was CNY 262.3 trillion, with an average annual growth rate of 5.4% [4]. Likewise, the progress of the CRE is also outstanding [5]. The CRE has already operated 49,000 trains, transported 4,432,000 TEUs of goods, and reached 180 cities in 23 European countries by the end of 2021 [5]. With the logistics service network covering the entire territory of Asia and Europe, CRE is becoming an international public logistics product widely recognized by countries along the route. Therefore, the inland ports and CRE trains are both important components of China’s cross-border trade logistics network, especially under the Belt and Road Initiative [3,5,6].
Even though the number of Chinese inland ports and the volume of cross-border trade have grown rapidly in recent years, it should be noted that there are still issues with the network layout of the China-Europe multimodal container transport system. On the one hand, some city nodes are over-subsidizing the operation of CRE and adopting low-priced strategies to compete with maritime transport for cargo sources, which distorts the market order to a certain extent [6]. On the other hand, some regions in the network of China have demonstrated the problem of over-intensive inland ports, which has led to increased competition among inland ports and the waste of transport and train resources [7]. With this concern in mind, we concentrate on the network layout for the cross-border multimodal container transport system based on inland ports between China and Europe.
Generally speaking, the transportation demand between inland cities and foreign hubs is influenced by multiple factors, including the global economic situation, trade policies between countries, exchange rate fluctuations, global epidemics, sudden natural disasters, etc., all of which will bring uncertainty. Based on the aforementioned background, this study addresses the following research questions:
  • How to design the entire network for the cross-border multimodal container transport system based on inland ports under uncertain demand conditions?
  • How to identify the impact of different factors on the optimal network structure?
  • What strategies can we propose to improve network performance against uncertain demand?
In order to answer the above research questions, this paper introduces the cross-border multimodal container transportation network design problem. Given the importance of carbon neutral policies, this paper especially considers the cost of carbon emissions, which can be expressed as the product of the volume of carbon emissions and the amount of carbon tax, and analyses the situations when the carbon tax changes. It is worth noting that the assumption of deterministic demand often becomes restrictive in the real-world cross-border multimodal container transport system. By contrast, this paper employs the uncertain variables to describe the uncertainty associated with demand via the experiences and judgments of field experts. Furthermore, this paper proposes an uncertain programming model for the cross-border multimodal container transportation network design problem to minimize the expected total costs by imposing two types of chance constraints for capacity limitations. To solve the proposed model efficiently, this paper also adopts uncertainty theory to convert the uncertain model into an equivalent crisp one, which can be directly handled by off-the-shelf solvers such as CPLEX, Gurobi, and Lingo. Finally, this paper takes the Huaihai Economic Zone-Europe multimodal container transport system as a case study and presents several valuable suggestions and policy guidance with strong theoretical and practical implications.
This paper is structured as follows: Section 2 provides an overview of the current literature on inland ports, multimodal container transport systems, and uncertain theory. Section 3 develops the deterministic and uncertain programming models for the cross-border multimodal container transportation network design problem. Section 4 investigates a case study of the Huaihai Economic Zone-Europe multimodal container system. Section 5 analyzes the computational results and proposes policy guidance. Section 6 gives the main conclusions and makes recommendations for further research.

2. Literature Review

2.1. Inland Port and Its Role in Cross-Border Trade

The definition of the inland port is widely accepted as an inland intermodal terminal directly connected to seaport(s) with high-capacity transport mean(s), where customers can leave/pick up their standardized units as if directly to a seaport based on Leveque and Roso [8]. There have been a number of studies focusing on seaports. For example, Bernacki and Lis [9] investigated the future evolution of port systems, considering the case of Poland and the Rhine-Scheldt Delta region. Szaruga et al. [10] presented the synchronization of economic cycles of GDP and oil product cargo volumes in major Polish seaports. Kotowska et al. [11] noticed that inland shipping to serve the hinterland is a big challenge for seaport authorities. Focusing on inland ports, Witte et al. [12] presented a systematic and integrated review of inland port studies, covering 80 international peer-reviewed academic journal papers on inland port development between 1992 and 2017, suggesting that the inland port has become a widely accepted concept. Wiegmans et al. [1] focused attention on the characteristics of inland waterway ports in a European context, proving that the definition of inland ports has been enriched beyond just the extension of seaports. Wang et al. [13] investigated the impact of inland port development in China on the promotion of bilateral trade flows between China and South Korea, suggesting that inland port development will increase cross-border trade volume. Monios and Wang [14] studied the spatial and institutional characteristics of inland port development in China, providing inland port development suggestions in a new geographical context. Xie et al. [15] suggested that China’s foreign trade needs new momentum and that it is important to optimize the logistic network under the Belt and Road Initiative.
Despite years of research on inland ports, there is still a lack of solutions for the optimal layout of inland ports in the cross-border transport system. Additionally, there hasn’t been much research done on the over-concentration of inland ports. To fill the research gap, this study develops a mathematical model to design the entire network by locating the inland ports.

2.2. Multimodal Container Transport System

There have been a number of models focusing on multimodal container transport systems that consider different indicators [16,17,18,19,20,21,22]. For instance, Jiang et al. [16] established a theoretical basis for the related exploration of multimodal container transportation. Corman et al. [17] presented a study on the influence and sensitivity of different model parameters, in order to analyze the implications on strategic decisions, fostering a target modal share for freight transportation. In order to describe the optimal organization problem, an optimization model based on dynamic programming was presented and was satisfied with reality constraints [18]. Hryhorak et al. [19] analyzed the dynamics and structure of freight transport in Ukraine. Containerized commodity transportation schemes are highly efficient for the majority of transcontinental and long-distance deliveries optimizing costs, and time quality of transport operations based on exact forecasting of container turnover. Zhang et al. [20] introduced a modeling approach for the optimization of terminal networks, taking into account the costs of CO2 emissions and economies of terminal scale. Zehendner et al. [21] proposed a mixed integer linear programming model, based on a network flow representation of the terminal, to determine the number of appointments to accept per time slot and an allocation of internal resources minimizing service times of trains and barges simultaneously. Wei and Dong [22] proposed a cross-border logistics network based on inland ports and established a two-objective model in order to optimize freight cost and transport time at the same time.
The modeling of multimodal container transport systems has been the subject of a sizable body of research [16,17,18,19,20,21], but the cross-border multimodal container transport systems based on inland ports have received less attention. By combining the amount of carbon emissions with the carbon tax, this study also proposes a more precise method for estimating the cost of carbon emissions.

2.3. Application of Uncertainty Theory

Liu [23] creatively provided a self-contained, comprehensive, and up-to-date presentation of uncertainty theory. Liu [23,24,25] showed a more comprehensive understanding of uncertain theory, including uncertain programming, uncertain risk analysis, uncertain reliability analysis, uncertain process, uncertain calculus, uncertain differential equation, uncertain logic, uncertain entailment, and uncertain inference. With the development of uncertainty theory, there is a growing literature on the applications of uncertainty theory. For example, Gu and Zhu [26] investigated a new type of optimal control problem governed by a parabolic uncertain partial differential equation in which the expected value criterion was adopted as the objective function. Shen and Zhu [27] employed uncertain programming to deal with the job shop scheduling problem with uncertain processing time and cost. Zhang and Peng [28] proposed the concepts of the expected shortest route, α-shortest route, and distribution shortest route in the Chinese postman problem based on uncertainty theory. Ke et al. [29] proposed uncertain random multilevel programming for modeling decentralized decision-making problems with uncertain random parameters. Yang et al. [30] built a multi-period uncertain workforce planning model with job satisfaction level, with uncertainty in labor demands and operation costs.
However, since uncertainty theory is a relatively new methodology, there are not many applications for it in cross-border logistics. This work represents the first attempt to integrate uncertainty theory into the cross-border multimodal container transport system based on inland ports. By employing an uncertain programming approach to tackle the considered problem, the findings are more applicable to the practical cross-border logistics process.

3. Problem Formulations

In this study, we introduce the cross-border multimodal container transportation network design problem, which is to locate the inland ports and choose the transportation modes between nodes with minimum total costs. As displayed in Figure 1, we have classified the total costs into five types in this paper. The network optimization model for the case of deterministic demand is analyzed first. On this premise, the investigation into the uncertain demand case is conducted, and the problem is converted into a crisp model, leading to solutions for the inland port selection, transport route selection, and minimum cost for each case.
Based on the background described in Section 1, the network for the cross-border multimodal container transport system based on inland ports is illustrated in Figure 2. More specifically, the considered network contains five types of nodes: origin city, domestic logistics park, domestic inland port, domestic seaport, and foreign hub, and the inland port is formed by the expansion of the logistics park. The transportation mode from origin cities to logistics parks/inland ports, or seaports is road. The transportation modes between logistics parks/inland ports and seaports include road, railway, and inland waterways. The transportation mode between inland ports and foreign hubs is CRE. The transportation mode between seaports and foreign hubs is shipping. In this paper, we only consider the one-way transportation process from domestic cities to foreign hubs.

3.1. Definition of Symbols, Parameters, and Decision Variables

3.1.1. Symbols and Parameters

To formally characterize the problem, we first list all symbols and parameters used in the problem formulation, as shown in Table 1.

3.1.2. Decision Variables

In this problem, the decision variables concern the selection of inland ports and the transportation route from each origin city to each foreign hub. Considering different kinds of routes from origin city to foreign hub, we get six types of decision variables in total, as described below:
Y i k l j 1 = { 1 ,   if   route   i road k road l shipping j   is   selected 0 ,   otherwise Y i k l j 2 = { 1 ,   if   route   i road k railway l shipping j   is   selected 0 ,   otherwise Y i k l j 3 = { 1 ,   if   route   i road k inland   water l shipping j   is   selected 0 ,   otherwise Y i l j = { 1 ,   if   route   i road l shipping j   is   selected 0 ,   otherwise Y i k j = { 1 ,   if   route   i road k CRE j   is   selected 0 ,   otherwise X k = { 1 ,   if   node   k   is   converted   to   an   inland   port 0 ,   otherwise

3.2. Deterministic Model

3.2.1. Calculation of Carbon Emissions

There are a number of approaches that can be used to estimate carbon emissions. Essentially, carbon emissions from the transportation process are affected by two primary elements: the kind of vehicle used (which is also regarded as the transportation method in this study) and the kind of energy used to power it. In this research, based on the formula for unit transport CO2 emissions from the IPCC Guidelines for National Greenhouse Gas Inventories (2006) issued by the United Nations [31], we get the following formula to compute transportation carbon emissions in this study.
e m = E m H P λ
in which, em represents the carbon emissions per unit distance per TEU transported by transportation mode m (tonne). Em denotes the combined energy consumption of standard coal used per unit of transport workload for transport mode m (kg). H represents the calorific value of standard coal (GJ/kg). P represents the emission factor of CO2 (kg/GJ). λ denotes the carbon oxidation factor. According to the IPCC database, P is 25.8 if the fuel used is coking coal. In China, the calorific value of standard coal is 7000 kcal/kg, which is 29.3067 GJ/t in international units. The parameter λ is usually taken as 1.

3.2.2. Assumptions

Assumption 1. The transportation demand from the domestic cities to the foreign hubs is determined.
Assumption 2. For the identified domestic cities and foreign hubs, only one transportation route can be used to meet transportation demand, and for a predetermined route, only one transportation mode can be taken between two adjacent points, which is consistent with the actual transportation process.
Assumption 3. Only inland ports have the necessary conditions for the opening of CRE trains, i.e., they have customs functions, a large handling capacity, and multimodal transportation.
Assumption 4. For transfer costs, only the costs associated with changing the mode of transportation are considered, and the transfer cost differences caused by infrastructures at different nodes are not taken into account.

3.2.3. Mathematical Formulation with Deterministic Demand

The objective of the optimization model is to minimize the total cost. Based on the notations and assumptions made before, we build the objective functions as in Equation (2). The objective function minimizes the total cost, which is composed of five terms. C1 represents the transport cost considering different routes and different transportation modes. C2 represents the carbon emission cost, which corresponds to the volume of goods, the route and transportation mode option, and the carbon tax. C3 represents the transfer cost, including the transfer handling cost and container occupancy cost, between different transportation modes. C4 represents the cost of customs clearance and is related to whether the clearance point is an inland port or a seaport. C5 represents the average annual cost of inland port conversion.
min C t o t a l = C 1 ( q i j ) + C 2 ( q i j ) + C 3 ( q i j ) + C 4 ( q i j ) + C 5 C 1 ( q i j ) = i I j J k K l L m = 1 3 q i j [ Y i k l j m ( C i k 1 + C k l s + C l j 4 ) + Y i l j ( C i l 1 + C l j 4 ) + Y i k j ( C i k 1 + C k j 5 ) ] = i I j J k K l L m = 1 3 q i j [ Y i k l j m ( D i k 1 c 1 + D k l m c m + D l j 4 c 4 ) + Y i l j ( D i l 1 c 1 + D l j 4 c 5 ) + Y i k j ( D i k 1 c 1 + D k j 5 c 5 ) ] C 2 ( q i j ) = i I j J k K l L m = 1 3 q i j F [ Y i k l j m ( D i k 1 e 1 + D k l m e m + D l j 4 e 4 ) + Y i l j ( D i l 1 e 1 + D l j 4 e 5 ) + Y i k j ( D i k 1 e 1 + D k j 5 e 5 ) ] C 3 ( q i j ) = i I j J k K l L m = 1 3 q i j [ Y i k l j m ( W 1 m + T 1 m H + W m 4 + T m 4 H ) + Y i l j ( W 14 + T 14 H ) + Y i k j ( W 15 + T 15 H ) ] C 4 ( q i j ) = i I j J k K l L m = 1 3 q i j [ ( Y i k l j m + Y i l j ) G l + Y i k j G k ] C 5 = k K Z k X k
We next present some basic constraints that the network design must satisfy, as described below. Constraint (3) ensures that there is only one route from the origin city i to the foreign hub j. Constraint (4) ensures that only when the logistics park is upgraded to an inland port can a CRE train be opened at this node. Constraints (5) and (6) are capacity constraints for logistics parks, inland ports, and seaports. Constraint (7) ensures that the average annual total investment in inland ports is less than the investment limit. Constraints (8) and (9) are logical constraints on the decision variables.
k K l L Y i k l j 1 + Y i k l j 2 + Y i k l j 3 + Y i l j + Y i k j = 1 , i I , j J
sign ( i I j J Y i k j ) = X k , k K
i I j J l L q i j ( Y i k l j 1 + Y i k l j 2 + Y i k l j 3 + Y i k j ) Q k 1 + X k Q k 2 ,   k K
i I j J k K q i j ( Y i k l j 1 + Y i k l j 2 + Y i k l j 3 + Y i l j ) Q l , l L
k K Z k X k B
Y i k l j 1 { 0 , 1 } , Y i k l j 2 { 0 , 1 } , Y i k l j 3 { 0 , 1 } , Y i l j { 0 , 1 } , Y i k j { 0 , 1 }   i I , j J , k K , l L
X k { 0 , 1 } ,   k K

3.3. Uncertain Programming Model

3.3.1. Preliminaries of Uncertainty Theory

Uncertainty theory can offer an axiomatic system to deal with uncertain information when historical data is limited. It is different from probability theory; the experiences and judgments of field experts determine the possibility of an event happening. Moreover, compared with fuzzy set theory, uncertainty theory is more suitable for handling the transportation demand in this study because the possibility measure has no self-duality property [32]. A number of studies demonstrated that removing the self-duality property from mathematics could produce unexpected results when modeling uncertain transportation demand [23]. We now introduce several fundamental concepts and characteristics of uncertainty theory and uncertain programming in this section, which will be used throughout the rest of the work [33].
Let Γ be a nonempty set and be a σ-algebra over Γ. Each element Λ is called an event. A set function from to [0,1] is called an uncertain measure if it satisfies the following four axioms:
Axiom 1
 [23]. (normality axiom) { Γ } = 1 for the universal set Γ.
Axiom 2
[23]. (duality axiom) { Λ } + { Λ c } = 1 for any event Λ .
Axiom 3
[23]. (subadditivity axiom) For every countable sequence of events { Λ i } , we have
{ i = 1 Λ i } i = 1 { Λ i }
Axiom 4
[23]. (product axiom) Let ( Γ i , i , i ) be uncertainty spaces for i = 1,2,…. Then the product uncertain measure is an uncertain measure satisfying
{ i = 1 Λ i } = i = 1 i { Λ i }
where Λ i are arbitrarily chosen events from i for i = 1,2,…, respectively.
Definition 1
[23]. An uncertain variable ξ is a measurable function from an uncertainty space (Γ, , ) to the set of real numbers; that is, for any Borel set B of real numbers, the set
{ ξ B } = { γ Γ , ξ ( γ ) B }
is an event.
For a sequence of uncertain variables ξ 1 , ξ 2 , , ξ n and a measurable function f, Liu proved that
ξ = f ( ξ 1 , ξ 2 , , ξ n )
defined as ξ ( γ ) = f ( ξ 1 ( γ ) , ξ 2 ( γ ) , , ξ n ( γ ) ) , γ Γ , is also an uncertain variable.
Definition 2
[23]. The uncertainty distribution Ф of an uncertain variable ξ is defined by
Φ ( x ) = M { ξ x }
for any real number x.
Definition 3
[24]. An uncertain variable ξ is called zigzag if it has a zigzag uncertainty distribution
Φ ( x ) = { 0 ,                   if   x a x a 2 ( b a ) ,   if   a x b x + c 2 b 2 ( c b ) ,   if   b x c 1 ,                       if   x c
denoted by 𝒵 (a,b,c) where a,b,c are real numbers with a<b<c. Figure 3 shows the Zigzag uncertainty distribution.
The inverse uncertainty distribution of zigzag uncertain variable 𝒵 (a,b,c) is
Φ 1 ( α ) = { ( 1 2 α ) α + 2 a b ,   if   α < 0.5 ( 2 2 α ) b + ( 2 α 1 ) c ,   if   α 0.5
The zigzag uncertain variable ξ = 𝒵 (a,b,c) has an expected value
E [ ξ ] = a + 2 b + c 4
Theorem 1
[25]. Let ξ be an uncertain variable with regular uncertainty Ф; if the expected value exists, then
E [ ξ ] = 0 1 Φ 1 ( α ) d α
For ξ = f ( ξ 1 , ξ 2 , , ξ n ) , if ξ is strictly increasing functions with respect to ξ 1 , ξ 2 , , ξ n , we have
E [ ξ ] = 0 1 Φ 1 ( α ) d α = 0 1 f ( Φ 1 1 ( α ) , Φ 2 1 ( α ) , , Φ n 1 ( α ) ) d α
Theorem 2
[25]. Let ξ 1 , ξ 2 , , ξ n be independent uncertain variables with uncertainty distributions Ф12,…,Фn, respectively. If f ( x , ξ 1 , ξ 2 , , ξ n ) and g i ( x , ξ 1 , ξ 2 , , ξ n ) are strictly increasing functions with respect to ξ 1 , ξ 2 , , ξ n , for i = 1, 2,…, l, then the uncertain programming
{ min x E [ f ( x , ξ ) ] subject   to { g i ( x , ξ ) 0 } α i , i = 1 , 2 , , l .
is equivalent to the crisp mathematical programming
{ min x 0 1 f ( x , Φ 1 1 ( α ) , Φ 2 1 ( α ) , Φ n 1 ( α ) ) d α subject   to g i ( x , Φ 1 1 ( α i ) , Φ 2 1 ( α i ) , Φ n 1 ( α i ) ) 0 , i = 1 , 2 , , l .

3.3.2. Mathematical Formulation with Uncertain Demand

The deterministic model in the previous section neglects the uncertain factors in the real-world logistics process of the cross-border multimodal container transport system. However, under the influence of the evolving global economic situation and epidemics such as COVID-19, uncertainty has become increasingly apparent in recent years. Hence, we relax Assumption 1 by considering the uncertain transportation demand from domestic cities to foreign hubs. In the following discussions, the demands are described as uncertain variables with a zigzag distribution function. That is, the container transportation demand between city i and foreign hub j, which is mentioned as qij in Section 3.1, can be denoted by uncertain variable ξ i j in this section.
In addition, the total cost has been proven in Equation (2) to be related to container transportation demand between city i and foreign hub j. According to Definition 2, we use η to denote the total cost, and η = f ( ξ 11 , ξ 12 , ξ 1 j , ξ 21 , ξ 22 , ξ i j ) is also an uncertain variable. Moreover, let the uncertainty distribution of ξ i j be Фij, and the uncertainty distribution of η be Ψ . Given credibility confidence levels β1 and β2, the uncertain programming model for the problem considered in this study can be described as follows:
{ min x = E   [ f ( x , ξ ) ] subject   to { i I j J l L ξ i j ( Y i k l j 1 + Y i k l j 2 + Y i k l j 3 + Y i k j ) Q k 1 + X k Q k 2 } β 1 ,   k K { i I j J k K ξ i j ( Y i k l j 1 + Y i k l j 2 + Y i k l j 3 + Y i l j ) Q l } β 2 ,   l L   constrains   ( 3 )   to   ( 4 ) ,   constrains   ( 7 )   to   ( 9 )
It is clear that f ( x , ξ 11 , ξ 12 , ξ 1 j , ξ 21 , ξ 22 , ξ i j ) and g k ( x , ξ 11 , ξ 12 , ξ 1 j , ξ 21 , ξ 22 , ξ i j ) are strictly increasing functions with respect to ξ 11 , ξ 12 , ξ 1 j , ξ 21 , ξ 22 , ξ i j , for i = 1,2,…,r, j = 1,2,…,s. Therefore, according to Theorem 1 and Theorem 2, the uncertain program is equivalent to the crisp mathematical model:
{ min = C 1 ( E ( q i j ) ) + C 2 ( E ( q i j ) ) + C 3 ( E ( q i j ) ) + C 4 ( E ( q i j ) ) + C 5 where   C 1 ( ) C 4 ( )   and   C 5   are   explained   in   equation   ( 2 )   subject   to i I j J l L Φ i j 1 ( β 1 ) ( Y i k l j 1 + Y i k l j 2 + Y i k l j 3 + Y i k j ) Q k 1 + X k Q k 2 , k K i I j J k K Φ i j 1 ( β 2 ) ( Y i k l j 1 + Y i k l j 2 + Y i k l j 3 + Y i l j ) Q l , l L   constrains   ( 3 )   to   ( 4 ) ,   constrains   ( 7 )   to   ( 9 )
Apparently, the program (20) is a binary integer programming model. The independent variable in program (20) is Xk and Y i k l j 1 , Y i k l j 2 , Y i k l j 3 , Y i l j , Y i k j , and the objective is to find a network solution including nodes and routes to minimize the total cost under different scenarios. Program (20) can be directly solved by general optimizers that are suitable for binary integer programming models.

4. Case Study

To show the applicability of the proposed models, we conducted a case study concerning the Huaihai Economic Zone-Europe multimodal container transport system. All the experiments were implemented on a Lenovo Thinkpad laptop with an Intel(R) Core (TM) [email protected] GHz CPU and 16.00 GB RAM using Microsoft Windows 10 (64-bit). All optimization models were solved directly by applying Lingo 18.0 with default settings.

4.1. Description of the Huaihai Economic Zone—Europe Multimodal Container Transport System

4.1.1. Introduction of Huaihai Economic Zone

The Huaihai Economic Zone is situated in the eastern bridgehead area of the Asia-Europe Continental Bridge, bordered by the Zhongyuan Economic Zone to the west, the Yangtze River Delta Economic Zone to the south, and the Bohai Economic Zone to the north [34,35]. As illustrated in Figure 4, due to its prominent location at the eastern bridgehead of the Asia-Europe Continental Bridge, it has assumed the essential role of a relay station that serves the east-west convergence of China’s economy.
The concept of the Huaihai Economic Zone was first proposed in 1986, and the broad Huaihai Economic Zone includes 20 cities in four provinces: Jiangsu, Shandong, Henan, and Anhui. In 2018, the National Development and Reform Commission of China [34] clarified the exact scope of the Huaihai Economic Zone, which includes a total of ten cities, as listed in Figure 4.

4.1.2. Huaihai Economic Zone—Europe Multimodal Container Transport System

The Huaihai Economic Zone-Europe multimodal container transport system in this study contains four types, 21 nodes in total. For origin cities, we considered all cities in the Huaihai Economic Zone, as listed in Figure 4. For logistics parks, we selected six nodes in the Huaihai Economic Zone, Xuzhou logistics park (XZLP): Suzhou logistics park (SZLP), Yudong logistics park (YDLP), Linyi logistics park (LYLP), Yanzhou logistics park (YZLP), and Zaozhuang logistics park (ZZLP). For seaports, we selected Lianyungang port and Rizhao port. It is important to note that the Rizhao port, which is located in the broad Huaihai Economic Zone, has been included in the network in order to meet the diversity of node selection. For foreign hubs, we selected three hubs that have both seaport and inland port functions: Rotterdam, Hamburg, and Duisburg.

4.1.3. Data Collection

In this section, we specify the relevant constant data for the case study. According to [36], combined with Equation (1), Table 2 reports the transportation costs and carbon emissions under different transportation modes. Table 3 shows the average annual investment for inland port conversion and the annual freight handling capability before/after conversion. The data source for this table is based on development plans for each inland port [35]. It should be noted that the useful life of inland ports is calculated on a 50-year basis, and cross-border services account for approximately 20% of the overall investment in the conversion of inland ports.
Additionally, Table 3 also shows the capacity of domestic seaports and foreign hubs. Table 4 uses historical data provided by [38] as the basis for statistics and analysis of the relevant levels of demand to provide a rough measure of cross-border transportation demand between each city in the Huaihai Economic Zone and each foreign hub. For clarity of presentation, the transport distance between different nodes by different transport modes is relegated to Appendix A in this paper.

4.2. Computational Results for Deterministic Model

We first conduct a case study to validate the efficiency of the deterministic model. Figure 5 shows the network structure when the transportation cost of CRE is 2.0 or 2.8. It appears that the location of inland ports and the transportation route will change significantly when the transportation cost of CRE differs. More specifically, when the transportation cost of CRE is 2.0, four logistics parks (XZLP, SZLP, LYLP, and ZZLP) are converted into THE inland ports, and eight CRE routes are selected. However, when the transportation cost of CRE is 2.8, only XZLP and ZZLP are expanded into inland ports, and only four CRE routes are selected.
Note that the capacity of inland ports is also one of the most important factors influencing the layout of the network. The estimated capacity of the inland port in this paper refers to the capacity in the planning scenario. In practice, due to the phasing of inland port construction and subsequent expansion, the actual increased capacity of the logistics park after expansion into an inland port may also change. As can be seen in Figure 6, the proportion of CRE grows with inland port capacity, and the proportion of CRE decreases significantly when the actual added capacity at inland ports is less than the estimated added capacity.
Figure 7 suggests that the optimal total cost of the system increases as the CRE cost rises, while the total cost under no investment restrictions is lower when the CRE cost is less than 2.4. The total cost with no investment restrictions would also be lower when the CRE cost is less than 2.4. Moreover, as the transportation cost of CRE decreases, the proportion of CRE carried tends to increase progressively. The increase in the proportion of CRE with no investment restrictions is particularly noticeable when the CRE cost is less than 2.4. Figure 8 reveals that the optimal total cost of the system is also affected by the carbon emissions of CRE and the value of the carbon tax. There is evidence to indicate that the rise in carbon emissions of CRE and carbon tax gives rise to the total cost.

4.3. Computational Results for Uncertain Model

We then implement the uncertain programming model by adopting the uncertainty demand in Appendix B. It can be observed from Figure 9 that the optimal total cost rises when β1 and β2 rise. That is to say, in order to satisfy greater uncertainty, the total cost needs to be sacrificed. This result confirmed for us that the optimization process under uncertainty needed to be reconsidered, and we therefore carried out a further detailed analysis of the situation relevant when the parameters are varied under uncertainty.
It can be seen from Figure 10 that the total carbon emissions are influenced synergistically by the parameters β1 and β2. Note that the two extreme values occur when β1 = 0.95, β2 = 0.7, and β1 = 0.95, β2 = 0.80. It is obvious that the carbon emission from shipping and inland water transportation is much lower than other transportation modes. Generally speaking, the proportion of shipping rises when β1 increases and β2 decreases, and when β1 = 0.95, in order to meet the capacity constraints for inland ports, more goods are transported by shipping instead of CRE, which can reduce the total carbon emission. However, since the objective of the model is to minimize the expected value of the total cost of the system, sometimes the route with a lower cost and higher carbon emission will also be chosen, which explains why in the scenario β1 = 0.95, β2 = 0.75, the carbon emission is higher than in the scenario β1 = 0.95, β2 = 0.8.
Figure 11 depicts the changes in the volume of goods transported by CRE when β1, β2, and the capacity of inland ports change. Firstly, it can be found that when the capacity of the inland port is deterministic, the larger β1 is, the less volume is generally transported by CRE. In other words, when the constraint on the logistics park and inland port capacity is tightened, shipping will be more likely to be chosen as the transport mode linking the Huaihai Economic Zone to Europe. However, it can be noted that there is a set of anomalous data in the results: the proportion of CRE is higher when β1 = 0.95 and proportion = 75% than when β1 = 0.90 and ratio = 75%. This is because the inland port option when β1 = 0.90 and ratio = 75% cannot meet the capacity constraint when β1 = 0.95 and ratio = 75%. By replacing Yanzhou inland port with Linyi inland port, the overall capacity of the inland port increased, and the optimal routes under this scenario, by increasing the relative total cost, gave rise to the volume of goods transported by CRE. Moreover, the evidence in Figure 11 also suggests that when other factors are determined, the value of β2 does not have a major impact on the choice of inland ports and the proportion of CRE, and only changes to several individual routes will be made.
Apparently, according to Figure 12, the proportion of goods transported by CRE under uncertain demand shows a significant decrease compared with the situation under deterministic demand. In addition, when β1 is determined, the increase in inland port capacity will increase the proportion of goods transported by CRE.
Table 5 shows the inland port locations, CRE routes, and shipping routes of different plans. The optimal network plan for the deterministic demand case, called Plan A, is not feasible under an uncertain demand situation. Under plan A, the estimated volume of goods transported through the Xuzhou inland port is 728,500 TEU, which is beyond the capacity of the inland port. To better describe the actual problem, four optimal network plans were made for the solution under uncertainty, taking into account the different values of the two parameters, as shown in Figure 13.
As illustrated in Table 5 and Figure 13, the optimal network plan changes a lot in different scenarios. However, the inland port locations do not change a lot, except for Plan C. This may imply that when other factors are determined, the inland port locations will not change significantly due to the uncertainty of demand. Moreover, compared with Plan A, Plans B to E all add a shipping route from Rizhao Port to Rotterdam. It can be concluded that shipping is more stable than CRE under uncertain demand.

5. Discussions and Policy Guidance

Our research provides the following discussions and policy guidance according to the results of the numerical experiments:

5.1. Inland Port Selections

  • The selection plan for inland ports is robust. Our results demonstrated that the selected plan of inland ports is robust against uncertain demand. Table 5 and Figure 13 show that four of the five plans choose XZIP, SZIP, YZIP, and ZZIP. Additionally, investment in inland ports is a long-term process with multiple investment risks, and the findings of this paper provide some support for investment decisions in inland ports. That is, when the investment and construction of an inland port is a choice made after scientific analysis, the inland port will have a certain degree of robustness and will be able to cope with the situation under changing demand.
  • There should not be too many inland ports in a certain region. According to the results, none of the plans converts all six logistics parks into inland ports. It is reasonable to assume that the number of inland ports in a certain region should be limited. Otherwise, although some inland ports are built, they are not selected in the optimal network plan, which means the capital resources for construction will be tied up while the total cost of the whole system will increase.
  • The construction of inland ports has priority. As indicated in Figure 11, goods transported through XZIP and ZZIP by CRE account for 75% of the total goods transported by CRE, which suggests that the construction of these two inland ports has a high impact on the optimal solution for the whole network and should be first considered to be built when the investment amount is not sufficient to build all four inland ports. The conclusion is consistent with the real situation that XZIP and ZZIP are built now, especially for XZIP, which has run over 1000 CRE trains in 2021 [3,5].

5.2. Transportation Routes Choices

  • CRE offers advantages in reducing total cost. As illustrated in Figure 5 and Figure 6, the rise in the proportion of goods transported by CRE decreases the total cost. We can assume that, although the transport cost per unit transported by shipping is lower than by CRE, the distance by shipping is longer, and the transfer time is increasing, it is likely that in most of the cross-border transport processes, CRE is more cost-effective than shipping. In the actual transport process, it will generally take 30–48 days for goods to be transported from China to Europe by shipping, but only 20–25 days by CRE [5]. The result is also consistent with the fact that CRE will save nearly 8–20% of the total cost compared with shipping [5].
  • Shipping is robust against uncertain demand. According to Table 5 and Figure 13, it is likely to be assumed that shipping can satisfy a higher credibility confidence level against uncertain demand. A plausible explanation for this phenomenon can be attributed to the fact that seaports generally have a greater handling capacity than inland ports.
  • Road is the main transportation mode for short distances within the country. Combining Figure 5 and Figure 13 and Table 5, there is evidence to prove that road is still the main mode of transportation within the Huaihai Economic Zone. According to Table A1 and Table A2, the distances between nodes within the Huaihai Economic Zone are less than 500 km. Under all scenarios, the road is generally the primary mode of transport in the case we study.

5.3. Strategies for Reducing the Total Cost

  • Reduce transport costs for CRE. Results in Figure 6 and Figure 7 support the opinion that the total cost increases as the transport cost of CRE increases, and the lower transportation cost of CRE gives rise to the proportion of goods transported by CRE. In the actual operation process, we can reduce the transport cost of CRE by increasing the full load rate, innovating the organization of CRE trains, and making technical innovations to the CRE carriers. Moreover, the marginal effect is larger when the cost of the CRE is 2.1 as well as 2.6. Thus, when the transport cost of CRE is a little over this value, we can also consider using subsidies to achieve marginal benefits.
  • Reduce carbon emissions from CRE.Figure 8 proves the view that the lower carbon emissions of CRE and lower carbon tax help to reduce the total costs. The carbon tax in China is CNY 54.22 on average in 2021 [39], and when the carbon tax is CNY 600 (which is similar to the average carbon tax in Europe), the slope is large, which means that the carbon emissions of CRE influence the total cost a lot in this scenario. The evidence suggests reducing the carbon emissions of CRE by, for example, improving energy conversion rates to adapt to an increasing carbon tax in the future.
  • Expand investment limits under suitable conditions.Figure 7 shows that when the transportation cost of CRE is less than 2.4, it is advisable to expand investment limits. The primary way to raise the investment limit is to optimize the structure of the investors. The current main investors in the inland port are mainly the local government, but in the future, the role of logistics real estate developers in inland port investment can be fully exploited.

5.4. Schemes for Improving Network Performance against Uncertain Demand

  • Increase the capacity of inland ports.Figure 11 and Figure 12 provide evidence that, under an uncertain demand situation, the change in β1 has a greater impact on the network optimal plan compared with the change in β2. And Figure 11 also demonstrates that when inland port capacity is increased to 125%, the network plan is more stable under different changes of β1, the inland port location plan is not changing, and the route selection does not change a lot. In order to get a higher credibility confidence level, the inland port capacity should be increased as much as possible.
  • Add suitable shipping routes to the network.Figure 13 and Table 5 support the idea that by adding a new shipping route to the network, the network’s robustness under uncertain demand is improved. It can be assumed that the shipping mode can help the whole system to withstand the risk that the logistics network is not able to meet transport demand due to uncertainty.

6. Conclusions

In this paper, we investigated a cross-border multimodal container transport system based on inland ports. We formulated both deterministic and uncertain models for the cross-border multimodal container transportation network design problem. Based on uncertainty theory, we converted the uncertain model into a crisp model. We provided a case study of the Huaihai Economic Zone-Europe multimodal container transport system. The results concluded the suggestions and policy guidance over four issues: the inland port locations, the transportation route choices, the strategies to reduce the total cost, and the schemes to improve network performance against uncertain demand.
The novelty of this work can be summarized as follows:
  • It is distinct from the aforementioned work in that we focus on the important role of inland ports in the cross-border container transportation system, particularly considering CRE as a significant transportation mode.
  • It is the first work in which the uncertain transportation demand, carbon emission, and customs clearance cost are jointly considered to determine the inland port location as well as the cross-border transportation routes and modes. Based on this, the obtained results can provide strong practical guidance.
  • This work presents an integrated uncertain programming model by combining the expected value model and a chance-constrained program to formulate the cross-border multimodal container transportation network design problem under uncertain transportation demand.
It should be noted that the carbon emissions of different modes of transportation are changing greatly due to the application of new transportation carriers such as new energy vehicles and LNG ships. However, the above-mentioned situation is not considered in this paper. Hence, in future work, we will address the cost of carbon emissions in a more precise way. Furthermore, in this paper, we have found that the objectives of minimizing total costs and minimizing carbon emissions are not always consistent. Therefore, we will obtain the optimal plan under multiple objectives in the future.

Author Contributions

All authors contributed to the study conception and design. Methodology and software, J.M. and K.Y.; formal analysis and investigation, J.M., K.Y. and L.J.; data curation and writing—original draft preparation, J.M.; writing—review and editing, X.W. and K.Y; supervision, X.W.; funding acquisition, X.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Major Scientific and Technological Innovation Project of Shandong Province (No. 2019JZZY020715).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Transport Distance between Different Nodes by Different Modes

Table A1. Transport distance from origin city to domestic logistics parks/domestic seaports by road.
Table A1. Transport distance from origin city to domestic logistics parks/domestic seaports by road.
Unit: kmXZLPSZLPYDLPLYLPYZLPZZLPLYG PORTRZ PORT
Xuzhou2010016020018075230308
Lianyungang22530040015026022520130
Suqian125200270160300190170250
Suzhou8010190300290180320400
Huaibei6050160260260150287370
Shangqiu17012015370240250394470
Zaozhuang9017025014412515240265
Jining16025523020540152340301
Linyi21027536020190150125160
Heze280220130310140270500405
Table A2. Transport distance from domestic logistics park to domestic seaports by different modes.
Table A2. Transport distance from domestic logistics park to domestic seaports by different modes.
Unit: kmRoadRailwayInland Waterway
LYG PortRZ PortLYG PortRZ PortLYG PortRZ Port
XZLP250320185300500700
SZLP320400300350400800
YDLP400480350420-1-
LYLP200170180150--
YZLP310270400295450400
ZZLP250270350300500400
1 The two nodes cannot be linked by inland waterway
Table A3. Transport distance from seaports to foreign hubs by shipping.
Table A3. Transport distance from seaports to foreign hubs by shipping.
Unit: kmRotterdamHamburgDuisburg
LYG Port20,82220,22421,851
RZ Port21,52021,39620,056
Table A4. Transport distance from inland ports to foreign hubs by CRE.
Table A4. Transport distance from inland ports to foreign hubs by CRE.
Unit: kmRotterdamHamburgDuisburg
XZLP10,40011,68311,000
SZLP12,01012,20010,855
YDLP12,90012,35012,860
LYLP12,30011,80013,400
YZLP12,10012,00011,050
ZZLP12,78010,43010,890

Appendix B. Cross-Border Transportation Uncertain Demand between Different Origin Cities and Foreign Hubs

Table A5. Cross-border transportation uncertain demand between origin cities and Rotterdam.
Table A5. Cross-border transportation uncertain demand between origin cities and Rotterdam.
Unit: 10,000TEURotterdam (a) 1Rotterdam (b) 1Rotterdam (c) 1
Xuzhou52040
Lianyungang2530
Suqian2510
Suzhou0.5520
Huaibei0.5520
Shangqiu135
Zaozhuang2520
Jining3520
Linyi51020
Heze2515
1 (a) refers to the minimum demand between nodes; (b) refers to the most likely demand between nodes; (c) refers to the maximum demand between nodes
Table A6. Cross-border transportation uncertain demand between origin cities and Hamburg.
Table A6. Cross-border transportation uncertain demand between origin cities and Hamburg.
Unit: 10,000TEUHamburg (a)Hamburg (b)Hamburg (c)
Xuzhou3510
Lianyungang358
Suqian135
Suzhou125
Huaibei124
Shangqiu245
Zaozhuang358
Jining51015
Linyi81020
Heze257
Table A7. Cross-border transportation uncertain demand between origin cities and Duisburg.
Table A7. Cross-border transportation uncertain demand between origin cities and Duisburg.
Unit: 10,000TEUDuisburg (a)Duisburg (b)Duisburg (c)
Xuzhou3510
Lianyungang125
Suqian0.523
Suzhou134
Huaibei157
Shangqiu123
Zaozhuang2510
Jining2510
Linyi3510
Heze81015

References

  1. Wiegmans, B.; Witte, P.; Spit, T. Characteristics of European Inland Ports: A Statistical Analysis of Inland Waterway Port Development in Dutch Municipalities. Transp. Res. Part A Policy Pract. 2015, 78, 566–577. [Google Scholar] [CrossRef]
  2. Zeng, Q.; Maloni, M.J.; Paul, J.A.; Yang, Z. Dry Port Development in China. Transp. J. 2013, 52, 234. [Google Scholar] [CrossRef]
  3. Association of Development District in China (Inland Port Branch). China Inland Port Development Report; China Fortune Press: Beijing, China, 2021; ISBN 9787504775092. [Google Scholar]
  4. Tao, K.; Chao, Y. Rethinking Port Role as Transport Corridor under Symbiosis Theory-Case Study of China-Europe Trade Transportation. Int. J. Sustain. Dev. World Policy 2019, 8, 51–61. [Google Scholar] [CrossRef] [Green Version]
  5. China National Railway Group Company. China-Europe Railway Express Development Report, 2021; China National Railway Group Company: Beijing, China, 2021. Available online: https://www.ndrc.gov.cn/fzggw/jgsj/kfs/sjdt/202208/P020220818311703111697.pdf (accessed on 1 November 2021).
  6. Gong, X.; Li, Z.-C. Determination of Subsidy and Emission Control Coverage under Competition and Cooperation of China-Europe Railway Express and Liner Shipping. Transp. Policy 2022, 125, 323–335. [Google Scholar] [CrossRef]
  7. Zheng, S.; Zhang, Q.; van Blokland, W.B.; Negenborn, R.R. The Development Modes of Inland Ports: Theoretical Models and the Chinese Cases. Marit. Policy Manag. 2021, 48, 583–605. [Google Scholar] [CrossRef]
  8. Roso, V.; Woxenius, J.; Lumsden, K. The Dry Port Concept: Connecting Container Seaports with the Hinterland. J. Transp. Geogr. 2009, 17, 338–345. [Google Scholar] [CrossRef]
  9. Bernacki, D.; Lis, C. Investigating the Future Dynamics of Multi-Port Systems: The Case of Poland and the Rhine–Scheldt Delta Region. Energies 2022, 15, 6614. [Google Scholar] [CrossRef]
  10. Szaruga, E.; Kłos-Adamkiewicz, Z.; Gozdek, A.; Załoga, E. Linkages between Energy Delivery and Economic Growth from the Point of View of Sustainable Development and Seaports. Energies 2021, 14, 4255. [Google Scholar] [CrossRef]
  11. Kotowska, I.; Mańkowska, M.; Pluciński, M. Inland Shipping to Serve the Hinterland: The Challenge for Seaport Authorities. Sustainability 2018, 10, 3468. [Google Scholar] [CrossRef] [Green Version]
  12. Witte, P.; Wiegmans, B.; Ng, A.K.Y. A Critical Review on the Evolution and Development of Inland Port Research. J. Transp. Geogr. 2019, 74, 53–61. [Google Scholar] [CrossRef]
  13. Wang, C.; Chu, W.; Kim, C.Y. The Impact of Logistics Infrastructure Development in China on the Promotion of Sino-Korea Trade: The Case of Inland Port under the Belt and Road Initiative. J. Korea Trade 2020, 24, 68–82. [Google Scholar] [CrossRef]
  14. Monios, J.; Wang, Y. Spatial and Institutional Characteristics of Inland Port Development in China. GeoJournal 2013, 78, 897–913. [Google Scholar] [CrossRef] [Green Version]
  15. Xie, J.; Sun, Y.; Huo, X. Dry Port-Seaport Logistics Network Construction under the Belt and Road Initiative: A Case of Shandong Province in China. J. Syst. Sci. Syst. Eng. 2021, 30, 178–197. [Google Scholar] [CrossRef]
  16. Jiang, J.; Lu, J. Research on Optimum Combination of Transportation Modes in the Container Multimodal Transportation System. In Logistics: The Emerging Frontiers of Transportation and Development in China, Proceedings of the Eighth International Conference of Chinese Logistics and Transportation Professionals (ICCLTP), Chengdu, China, 8–10 October 2009; American Society of Civil Engineering: Reston, VA, USA, 2009. [Google Scholar] [CrossRef]
  17. Corman, F.; Viti, F.; Negenborn, R.R. Equilibrium Models in Multimodal Container Transport Systems. Flex. Serv. Manuf. J. 2015, 29, 125–153. [Google Scholar] [CrossRef] [Green Version]
  18. Hao, C.; Yue, Y. Optimization on Combination of Transport Routes and Modes on Dynamic Programming for a Container Multimodal Transport System. Procedia Eng. 2016, 137, 382–390. [Google Scholar] [CrossRef] [Green Version]
  19. Hryhorak, M.; Lyakh, O.; Sokolova, O.; Chornogor, N.; Mykhailichenko, I. Multimodal Freight Transportation as a Direction of Ensuring Sustainable Development of the Transport System of Ukraine. IOP Conf. Ser. Earth Environ. Sci. 2021, 915, 012024. [Google Scholar] [CrossRef]
  20. Zhang, M.; Wiegmans, B.; Tavasszy, L. Optimization of Multimodal Networks Including Environmental Costs: A Model and Findings for Transport Policy. Comput. Ind. 2013, 64, 136–145. [Google Scholar] [CrossRef]
  21. Zehendner, E.; Feillet, D. Benefits of a Truck Appointment System on the Service Quality of Inland Transport Modes at a Multimodal Container Terminal. Eur. J. Oper. Res. 2014, 235, 461–469. [Google Scholar] [CrossRef]
  22. Wei, H.; Dong, M. Import-Export Freight Organization and Optimization in the Dry-Port-Based Cross-Border Logistics Network under the Belt and Road Initiative. Comput. Ind. Eng. 2019, 130, 472–484. [Google Scholar] [CrossRef]
  23. Liu, B. Uncertainty Theory; Springer: New York, NY, USA, 2016; ISBN 9783662499887. [Google Scholar]
  24. Liu, B. Uncertainty Theory: A Branch of Mathematics for Modelling Human Uncertainty; Springer Science & Business Media: Cham, Switzerland, 2011; ISBN 9783642139581. [Google Scholar]
  25. Liu, B.; Kacprzyk, J. Theory and Practice of Uncertain Programming; Springer: Berlin/Heidelberg, Germany, 2009; ISBN 9783540894841. [Google Scholar]
  26. Gu, Y.; Zhu, Y. Optimal Control for Parabolic Uncertain System Based on Wavelet Transformation. Axioms 2022, 11, 453. [Google Scholar] [CrossRef]
  27. Shen, J.; Zhu, Y. Chance-Constrained Model for Uncertain Job Shop Scheduling Problem. Soft Comput. 2015, 20, 2383–2391. [Google Scholar] [CrossRef]
  28. Zhang, B.; Peng, J. Uncertain Programming Model for Chinese Postman Problem with Uncertain Weights. Ind. Eng. Manag. Syst. 2012, 11, 18–25. [Google Scholar] [CrossRef] [Green Version]
  29. Ke, H.; Su, T.; Ni, Y. Uncertain Random Multilevel Programming with Application to Production Control Problem. Soft Comput. 2014, 19, 1739–1746. [Google Scholar] [CrossRef]
  30. Yang, G.; Tang, W.; Zhao, R. An Uncertain Workforce Planning Problem with Job Satisfaction. Int. J. Mach. Learn. Cybern. 2016, 8, 1681–1693. [Google Scholar] [CrossRef]
  31. IPCC Publications—IPCC-TFI. Available online: https://www.ipcc-nggip.iges.or.jp/public/2006gl/index.html (accessed on 1 April 2007).
  32. Zhang, W.; Wang, X.; Yang, K. Incentive Contract Design for the Water-Rail-Road Intermodal Transportation with Travel Time Uncertainty: A Stackelberg Game Approach. Entropy 2019, 21, 161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Gao, Y.; Yang, L.; Li, S. Uncertain Models on Railway Transportation Planning Problem. Appl. Math. Model. 2016, 40, 4921–4934. [Google Scholar] [CrossRef]
  34. Tao, J.; Lu, Y.; Ge, D.; Dong, P.; Gong, X.; Ma, X. The Spatial Pattern of Agricultural Ecosystem Services from the Production-Living-Ecology Perspective: A Case Study of the Huaihai Economic Zone, China. Land Use Policy 2022, 122, 106355. [Google Scholar] [CrossRef]
  35. China National Development and Reform Commission. Huaihe Ecological and Economic Belt Development Plan. Available online: https://www.ndrc.gov.cn/xxgk/zcfb/ghwb/201811/t20181107_962252.html?code=&state=123 (accessed on 15 November 2018).
  36. Zhang, H. Research on Optimization of China-EU Container Multimodal Transportation Route Based on Robust Optimization; Dalian Maritime University: Dalian, China, 2020. [Google Scholar]
  37. Zhang, X.; Zhang, H.; Yuan, X.; Hao, Y. Optimization of low-carbon multimodal transport routes under double uncertainty. J. Beijing Jiaotong Univ. 2022, 21, 113–120. [Google Scholar]
  38. Jiangsu Development and Reform Commission. Huaihai International Inland Port Develop Plan (2021–2025). Available online: http://doc.jiangsu.gov.cn/zcq/newsinfo.html?id=17455 (accessed on 19 November 2021).
  39. Tan, Y. An Overview of China’s Carbon Market in Its First Year of Operation. Available online: https://chinadialogue.net/zh/3/75074/ (accessed on 1 February 2022).
Figure 1. A schematic diagram of the stated work.
Figure 1. A schematic diagram of the stated work.
Axioms 12 00132 g001
Figure 2. The network studied in this paper.
Figure 2. The network studied in this paper.
Axioms 12 00132 g002
Figure 3. Zigzag uncertainty distribution.
Figure 3. Zigzag uncertainty distribution.
Axioms 12 00132 g003
Figure 4. Location of Huaihai Economic Zone.
Figure 4. Location of Huaihai Economic Zone.
Axioms 12 00132 g004
Figure 5. (a) Network layout when transportation cost of CRE is 2.0; (b) network layout when transportation cost of CRE is 2.8.
Figure 5. (a) Network layout when transportation cost of CRE is 2.0; (b) network layout when transportation cost of CRE is 2.8.
Axioms 12 00132 g005
Figure 6. Changes in the proportion of CRE and total optimal costs when actual added capacity are changed.
Figure 6. Changes in the proportion of CRE and total optimal costs when actual added capacity are changed.
Axioms 12 00132 g006
Figure 7. Changes in the proportion of CRE and total optimization costs when CRE transport costs and investment limits are changed.
Figure 7. Changes in the proportion of CRE and total optimization costs when CRE transport costs and investment limits are changed.
Axioms 12 00132 g007
Figure 8. Changes in optimal costs when the carbon emissions of CRE and carbon tax change.
Figure 8. Changes in optimal costs when the carbon emissions of CRE and carbon tax change.
Axioms 12 00132 g008
Figure 9. Optimal total cost versus values of β1 and β2.
Figure 9. Optimal total cost versus values of β1 and β2.
Axioms 12 00132 g009
Figure 10. Carbon emissions versus values of β1 and β2.
Figure 10. Carbon emissions versus values of β1 and β2.
Axioms 12 00132 g010
Figure 11. Changes in volume of goods transported by CRE when β1, β2, and capacity of inland ports changes.
Figure 11. Changes in volume of goods transported by CRE when β1, β2, and capacity of inland ports changes.
Axioms 12 00132 g011
Figure 12. Changes in proportion of goods transported by CRE when β1, β2, and capacity of inland ports changes.
Figure 12. Changes in proportion of goods transported by CRE when β1, β2, and capacity of inland ports changes.
Axioms 12 00132 g012
Figure 13. Optimal network plans on different scenarios.
Figure 13. Optimal network plans on different scenarios.
Axioms 12 00132 g013
Table 1. Symbols and parameters used in this problem.
Table 1. Symbols and parameters used in this problem.
NotationsDetailed Definition
IThe set of all origin cities.
JThe set of all foreign hubs.
KThe set of all domestic logistics parks.
LThe set of all domestic seaports.
C a b m The transport costs per TEU cargo from node a to node b by transport mode m (yuan/TEU).
D a b m The transport distance from node a to node b by transport mode m (km).
mThe transport mode: m = 1 refers to road, m = 2 refers to railway, m = 3 refers to inland waterway, m = 4 refers to shipping, m = 5 refers to CRE.
cmThe transport costs per kilometer per TEU cargo by transport mode m (yuan/km·TEU).
cmThe carbon emissions per kilometer per TEU cargo by transport mode m (tonne/km·TEU).
Wm1m2The handling costs per TEU cargo converted from transport mode m1 to m2 (yuan/TEU).
Tm1m2The waiting time of cargo converted from transport mode m1 to m2 (day).
HThe container occupancy cost per day per TEU cargo (yuan/TEU·day).
GlThe customs clearance costs per TEU cargo from seaport l (yuan/TEU).
GkThe customs clearance costs per TEU cargo from inland port k (yuan/TEU).
FThe carbon tax per tonne of carbon emissions (yuan/tonne).
ZkThe average annual construction cost of converting node k to an inland port (10,000 yuan).
BThe average annual investment limit in the conversion of inland ports (10,000 yuan).
QLThe maximum annual handling capacity of the domestic seaport l (10,000 TEU).
Q k 1 The maximum annual handling capacity when node k is not expanded into inland port (10,000 TEU).
Q k 2 The maximum annual handling capacity when node k is expanded into inland port (10,000 TEU).
qijThe annual transport demand from origin city i to foreign hub j (10,000 TEU).
Table 2. Transportation costs and carbon emissions under different transportation modes.
Table 2. Transportation costs and carbon emissions under different transportation modes.
Transportation ModeRoadRailwayInland Water TransportationShippingCRE
Transportation cost 1 (yuan/TEU·km)102.71.01.52.5
Carbon emission 1
(tonne/TEU·km)
1.77 × 10−39.0 × 10−43.0 × 10−43.0 × 10−48.0 × 10−4
1 Data based on [31,36,37].
Table 3. Average annual investment for conversion and capacity before/after conversion.
Table 3. Average annual investment for conversion and capacity before/after conversion.
NodeXZLPSZLPYDLPLYLPYZLPZZLPLYG
Port
RZ
Port
Rotterdam
Port
Hamburg
Port
Duisburg
Port
Origin
Capacity
(10,000TEU)
105510510100100300300300
Capacity
Added by
Conversion
(10,000TEU)
50510201530-----
Average annual investment for conversion
(10,000 yuan)
20005001000100010001500-----
Table 4. Cross-border transportation demand for different cities in Huaihai Economic Zone.
Table 4. Cross-border transportation demand for different cities in Huaihai Economic Zone.
Unit: 10,000TEURotterdamHamburgDuisburg
Xuzhou2055
Lianyungang552
Suqian532
Suzhou523
Huaibei525
Shangqiu342
Zaozhuang555
Jining5105
Linyi10105
Heze5510
Table 5. Computational results in Plan A to Plan E.
Table 5. Computational results in Plan A to Plan E.
Inland Port
Locations
CRE RoutesShipping Routes
Plan A:
Deterministic demand situation
XZ, SZ,
YZ, ZZ
XZ inland port–RotterdamRZ port–Duisburg
XZ inland port–DuisburgLYG port–Rotterdam
SZ inland port–Duisburg
YZ inland port–DuisburgLYG port–Hamburg
ZZ inland port–Hamburg
Plan B:
Uncertain demand situation, β1 = 0.8, β2 = 0.8
XZ, SZ,
YZ, ZZ
XZ inland port–RotterdamRZ port–Duisburg
SZ inland port–DuisburgRZ port–Rotterdam
YZ inland port–Rotterdam
YZ inland port–DuisburgLYG port–Hamburg
ZZ inland port–DuisburgLYG port–Rotterdam
ZZ inland port–Hamburg
Plan C:
Uncertain demand situation, β1 = 0.9, β2 = 0.8
XZ, YZ, ZZXZ inland port–RotterdamRZ port–Duisburg
XZ inland port–DuisburgRZ port–Rotterdam
XZ inland port–Hamburg
YZ inland port–DuisburgLYG port–Hamburg
ZZ inland port–DuisburgLYG port–Rotterdam
ZZ inland port–Hamburg
Plan D:
Uncertain demand situation, β1 = 0.8, β2 = 0.9
XZ, SZ,
YZ, ZZ
XZ inland port–RotterdamRZ port–Duisburg
SZ inland port–DuisburgRZ port–Rotterdam
YZ inland port–Hamburg
YZ inland port–DuisburgLYG port–Hamburg
ZZ inland port–HamburgLYG port–Rotterdam
Plan E:
Uncertain demand situation, β1 = 0.9, β2 = 0.9
XZ, SZ,
YZ, ZZ
XZ inland port–RotterdamRZ port–Duisburg
XZ inland port–DuisburgRZ port–Rotterdam
SZ inland port–Duisburg
YZ inland port–DuisburgLYG port–Hamburg
ZZ inland port–HamburgLYG port–Rotterdam
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

Ma, J.; Wang, X.; Yang, K.; Jiang, L. Uncertain Programming Model for the Cross-Border Multimodal Container Transport System Based on Inland Ports. Axioms 2023, 12, 132. https://doi.org/10.3390/axioms12020132

AMA Style

Ma J, Wang X, Yang K, Jiang L. Uncertain Programming Model for the Cross-Border Multimodal Container Transport System Based on Inland Ports. Axioms. 2023; 12(2):132. https://doi.org/10.3390/axioms12020132

Chicago/Turabian Style

Ma, Junchi, Xifu Wang, Kai Yang, and Lijun Jiang. 2023. "Uncertain Programming Model for the Cross-Border Multimodal Container Transport System Based on Inland Ports" Axioms 12, no. 2: 132. https://doi.org/10.3390/axioms12020132

APA Style

Ma, J., Wang, X., Yang, K., & Jiang, L. (2023). Uncertain Programming Model for the Cross-Border Multimodal Container Transport System Based on Inland Ports. Axioms, 12(2), 132. https://doi.org/10.3390/axioms12020132

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