Optimization Applications in Maritime Logistics and Operations
Optimization Applications in Maritime Logistics and Operations
By reading and using the thesis, the reader understands and agrees to the following terms:
1. The reader will abide by the rules and legal ordinances governing copyright regarding the
use of the thesis.
2. The reader will use the thesis for the purpose of research or private study only and not for
distribution or further reproduction or any other purpose.
3. The reader agrees to indemnify and hold the University harmless from and against any loss,
damage, cost, liability or expenses arising from copyright infringement or unauthorized
usage.
IMPORTANT
If you have reasons to believe that any materials in this thesis are deemed not suitable to be
distributed in this form, or a copyright owner having difficulty with the material being included in
our database, please contact [email protected] providing details. The Library will look into
your claim and consider taking remedial action upon receipt of the written requests.
Pao Yue-kong Library, The Hong Kong Polytechnic University, Hung Hom, Kowloon, Hong Kong
http://www.lib.polyu.edu.hk
OPTIMIZATION APPLICATIONS IN
MARITIME LOGISTICS AND OPERATIONS
WANG KAI
PhD
2019
The Hong Kong Polytechnic University
Wang Kai
July 2019
Certificate of Originality
I hereby declare that this thesis is my own work and that, to the best of my knowledge
and belief, it reproduces no material previously published or written, nor material that
has been accepted for the award of any other degree or diploma, except where due
acknowledgment has been made in the text.
_________________________ (Signed)
I would like to thank my supervisor Dr Shuaian Wang for his full support all the
time and thank my thesis board members, Professor Miao Song, Professor Xiangtong
Qi, Professor Yanzhi Li, and Dr Shuaian Wang for their time and valuable feedback.
I would also like to thank the financial support from the Hong Kong Ph.D. Fellowship
Scheme 2016/2017 and thank the research funding support from the Department of
Logistics and Maritime Studies, The Hong Kong Polytechnic University.
Figure 2.1: Four foldable containers and a standard container (Shintani et al.,
2010) .............................................................................................................. 4
Figure 2.2: A transpacific shipping service route operated by CMA CGM
(CMA CGM, 2017) ........................................................................................ 9
Figure 2.3: A sub-network for standard containers without considering the
short-term container leasing ......................................................................... 15
Figure 2.4: The sub-network defined for the preliminary NF model ......................... 18
Figure 2.5: A sub-network for foldable containers .................................................... 20
Figure 2.6: Connections between standard containers and foldable containers in
TNF model ................................................................................................... 20
Figure 2.7: An example of the pivot cycle W ............................................................ 24
Figure 2.8: Adding dummy nodes to the network ..................................................... 26
Figure 2.9: Three selected shipping services routes operated by CMA CGM
(CMA CGM, 2017) ...................................................................................... 32
Figure 2.10: Effect of the long-term leasing cost on the foldable container
usage ............................................................................................................ 38
Figure 2.11: Effect of the folding and unfolding cost on foldable container
usage ............................................................................................................ 39
Figure 2.12: Long-term leasing cost dependent sensitivity to the folding and
unfolding cost............................................................................................... 41
Figure 3.1: A sequence of container ships in a string ................................................ 50
Figure 3.2: The shipping routes involved in three case studies (CMA CGM,
2017) ............................................................................................................ 60
Figure 4.1: Itinerary of 7 Day Western Caribbean of Carnival (Carnival Cruise
Line, 2016) ................................................................................................... 77
Figure 4.2: Utility distribution for one day ................................................................ 78
Figure 4.3: Utility distributions and time spent at two ports ..................................... 79
Figure 4.4: Statistic on trips for the Carnival Vista (Cruise Ship Schedule,
2016) ............................................................................................................ 84
Figure 4.5: City ports in “14 Night Singapore to Fremantle Cruise” (Google
Map, 2018a) ................................................................................................. 92
Figure 4.6: City ports around the Caribbean Sea (Google Map, 2018b) ................... 95
Figure 4.7: Comparisons between the designed and actual itineraries ...................... 96
Figure 5.1: Concavity of the profit function ............................................................ 115
Figure 5.2: Locations of the port cities .................................................................... 129
Figure 5.3: The robustness test on 𝒑 ........................................................................ 131
This thesis includes four essays with a focus on maritime logistics and operations. It
studies four optimization application problems in the field, where the first two problems are
related to container ships that transport cargo and the other two problems are related to cruise
ships that ship cruise passengers. Different operations research methods are adapted to provide
effective solution approaches for optimization problems. Extensive experiments are conducted
to draw managerial insights on the problems.
1.1 BACKGROUND
In the container shipping industry, container ships are used to transport cargo among
seaports in the world, which contributes to the majority of shipping volume in the global trade.
The cost-effectiveness of sea transportation is the major competitiveness compared with other
transportation modes, such as road and air transportation. Thus, shipping companies endeavor
to optimize their operations such that they can reduce their operations costs. This is especially
important now since the industry is experiencing a transformation to an environmentally
friendly industry (Xia et al., 2019), which request more investments from the shipping
companies.
In the cruise shipping industry, cruise ships serve the purposes of providing pleasure
voyages. Once cruise passengers are on aboard, the voyages and the activities by the cruise
ships bring utility experience to the passengers, for which the cruise ships earn profits. The
transportation is not the major purpose of the cruise shipping since the cruise passengers
normally will return to the same port where they embark. The cruise industry is a very
promising industry and it keeps booming in recent years, especially for the Asia market (Wang
et al., 2016; Wang et al., 2017). However, the studies on this industry are rare in the literature,
but the operations of the industry are very urgent to be investigated.
The container shipping industry and the cruise shipping industry share many similarities
since they all belong to the maritime logistics and operations. They all follow some pre-
determined schedules to traverse a set of fixed ports of call, and they need to dwell at berths
of ports to conduct some on-shore activities. The major operational cost of the container and
cruise ships is the fuel consumption cost. However, they have several remarkable differences.
For instance, the container shipping aims to transport the cargo from its origin port to its
destination port, but the cruise shipping aims to provide voyage service for the cruise
passengers on aboard.
Chapter 1: Introduction 1
In this thesis, four operation optimization problems will be addressed and solved for the
container shipping and the cruise shipping, from which some similarities and differences
between the two industries can be sensed. The first two problems aim to optimize the
operations for container ships such that the operational costs can be reduced. The other two
problems aim to improve the operations for cruise ships such that the operating profits can be
increased. Different operations research methods are adopted to solve the problems, and
extensive numerical experiments are conducted to explore the hidden insights behind the
problems.
The remainder of this thesis is organized as follows. Chapter 2 and 3 are related to the
operations of container ships, and Chapter 4 and 5 are related to the operations of cruise ships.
In specific, Chapter 2 presents a study of ship type decision considering empty container
repositioning and foldable containers. Chapter 3 studies the optimal reefer slot conversion for
container freight transportation. Chapter 4 considers the problem of cruise itinerary schedule
design for a cruise ship. Chapter 5 addresses the cruise service planning considering berth
availability and decreasing marginal profit. Chapter 6 concludes the thesis.
Chapter 1: Introduction 2
Chapter 2: Container Ship Type Decision
This chapter addresses a ship fleet type decision problem with considering empty
container repositioning and foldable containers. This decision problem determines the capacity
of container ships deployed in a given shipping route. In reality, the ship fleet deployment
decision is usually made empirically according to laden container transportation and does not
consider the possible empty container repositioning. Meanwhile, a novel mode ‘foldable
container’ has shown its economic and logistical viability in recent years. This study hence
also considers the use of foldable containers and aims to find under what conditions, a shipping
liner needs to use the foldable containers in its liner shipping services. To solve the problem,
we formulate a network flow model with a revised network simplex algorithm, based on which
an exact solution approach is designed to determine the optimal ship type. A mixed-integer
programming model is also formulated for the problem. Numerical experiments based on real-
world voyages are conducted to find some managerial implications on the ship fleet
deployment and the foldable container usage.
2.1 INTRODUCTION
A shipping liner normally operates weekly-serviced ship routes with fixed schedules to
transport containers. Given a shipping route, a shipping liner deploys a fleet of container ships
for the operation over a planning horizon, e.g., six months. One of the critical decisions for the
shipping liner on the fleet is ship type decision, which determines the capacity of container
ships of the fleet deployed on the shipping route. Empirically, the shipping liner deploys a
suitable fleet type of ships on each route based on the laden container transportation over the
planning horizon, which guarantees that the deployed ships have the capacity to accommodate
all the laden containers in all the voyages. Under this circumstance, the shipping liner would
not deploy a ship fleet with a larger capacity as it increases the fixed operating cost for
maintaining the fleet. It is reasonable for the shipping liner to make such decision only
considering the laden container transportation. However, if we further consider the empty
container repositioning on the route, the ship type decision can be more complicated.
The empty container repositioning originates from the imbalance of container flow between
different regions in liner shipping routes. Take the trans-Pacific trade lane for example:
according to UNCTAD (2016), in 2015, the annual container flow from Asia to North America
(i.e., the eastbound) was around 15.8 million twenty-foot equivalent units (TEUs), and the
container flow in the opposite westbound direction was 7.4 million TEUs, which generated
Figure 2.1: Four foldable containers and a standard container (Shintani et al., 2010)
Storing empty containers in the depots and repositioning empty containers among the ports
inevitably incur storage cost and repositioning cost for the shipping liner, respectively (Lee
and Yu, 2012). To reduce costs, the usage of foldable containers is an effective method. The
idea of foldable containers is not so new, and several container companies have developed
foldable containers, such as Fallpac AB and Holland Container Innovation. Those foldable
containers have equivalent storage capacity and size as standard containers and a foldable
container only occupies one-quarter storage space of a standard container in folded status, as
shown in Figure 2.1 (See Appendix A for the specification comparison between a standard
container and a foldable container). After becoming empty, the foldable containers will be in
folded status for the storage in the depots or for the repositioning to other ports. As four
foldable empty containers in the folded status equal one standard empty container, it saves
75% storage space by using foldable containers, which leads to significant decrease in the
repositioning cost and the storage cost. However, using the foldable containers could incur
There have been numerous studies related to the ship fleet deployment and empty container
repositioning problems. For the ship fleet deployment, to the best of our knowledge, Perakis
and Jaramillo (1991) was the first study to address it, in which they built integer linear
programming models for the problem. Thereafter, there were generally two types of studies
on the problem. One type assumes that the container shipment demand is deterministic
(Gelareh and Meng, 2010; Brouer et al., 2013; Plum et al., 2013). For instance, Gelareh and
Meng (2010) developed a mixed integer nonlinear programming model for a short-term fleet
deployment problem, in which the optimal vessel speeds for different vessel types on different
routes are considered. The other type of studies relaxed the deterministic demand assumption
and treated the container shipment demand in a stochastic manner (Meng and Wang, 2012;
Ng, 2014; Ng, 2015). Meng and Wang (2012) addressed a practical ship fleet problem under
the background of week-dependent container shipment demand. Their study generated
practical container routes considering transit time constraints by using space-time network
approach. For more works on the ship fleet deployment, one can refer to Ng (2016), in which
it elaborated a class of fleet deployment models.
For the empty container repositioning or empty container allocation, Crainic et al. (1993)
introduced two dynamic formulations for empty container allocation in single and multi-
commodity cases, which provided a general framework to formulate this class of problems.
Cheung and Chen (1998) considered a dynamic empty container allocation problem, which
helped to determine the number of containers leased to fulfill the demands of customers over
time. The management of importing and exporting empty containers in a port was analyzed by
Li et al. (2004) based on the multi-stage inventory theory, and Markov decision processes were
proposed for the problem. Li et al. (2007) extended the previous study to a multi-port
application with a proposed heuristic algorithm for the problem. The empty container
repositioning problem for general shipping service routes was formulated by Song and Dong
Our problem focuses on ship type decision considering the empty container repositioning
and foldable containers for a given shipping service route. A fleet of container ships with
certain capacity is to be determined for weekly serving the ports along the shipping route. The
incoming container ship transports weekly laden containers originating from the ports to
destination ports. In each port, the empty container availability is critical for the shipping liner
to meet the laden container transportation. Generally, there are three empty container supplies
to fulfill laden container transportation. The economic supply is the arriving containers with
fully loaded goods from the incoming ship. Once these fully loaded containers arrive at the
ports and are delivered to consignees for unloading, the containers become empty and are
stored in depots for the laden consignment. The second supply is to reposition the empty
containers from surplus ports to deficit ports along the shipping route, which occupy the
capacity of the container ships among voyage legs. However, if the stored and repositioned
empty containers cannot meet the weekly laden container transportation, the shipping liner
must lease empty containers from container companies by short-term container leasing, which
is the third supply for the empty containers. The short-term container leasing is on an O-D pair
basis (See Section 2.3.2), which is different from the long-term container leasing that is on a
planning horizon basis (See Section 2.3.1).
Figure 2.2: A transpacific shipping service route operated by CMA CGM (CMA CGM,
2017)
In this study, the given shipping route has a fixed port rotation, shown by an illustrative
example in Figure 2.2. The itinerary of this route forms a loop: we can arbitrarily deem that
this itinerary starts at Shanghai and ends at Shanghai. Let 𝑝 ∈ 𝑃 represent the index of the
ports on a round trip for the route. Then, for this route, we can define Shanghai as Port 1,
Ningbo as Port 2, Pusan as Port 3, Los Angeles as Port 4 and Oakland as Port 5. Based on the
given route, the shipping liner has a set of O–D (Origin-Destination) port pairs 𝐷. The laden
2.3.3 Assumptions
Before addressing the model for our problem, we clarify some underlying assumptions:
(i) In the beginning, all the ships depart from the first port in the shipping route, and they
depart one by one with one-week interval (i.e., the first ship departs in the first week; the
second ship departs in the second week, and so on).
Assumption (i) indicates that when a shipping liner starts to operate a weekly shipping
service, it will mass a ship fleet at a homeport (i.e., the first port), and dispatch ships with a
one-week interval to form the weekly service pattern.
In essence, our problem has a latent network structure when the capacity of the ship fleet is
determined. Given the ship capacity, we can transfer the problem to a nonstandard minimum
cost flow problem by formulating a network flow (NF) model, rather than a mathematical
programming (MP) model. An MP model normally needs some commercial solvers to
optimize, such as CPLEX and Gurobi, which are not desirable for many shipping liner
companies. However, several network algorithms can easily solve an NF model to optimality,
for example, the network simplex algorithm. Therefore, in this section, we construct an NF
2.4.1 Notations
Indices and sets:
𝑝: Index of ports,
𝑃: Set of all the ports in the shipping service route,
(𝑜, 𝑑): Index of O–D port pairs, where 𝑜 ∈ 𝑃, 𝑑 ∈ 𝑃,
𝐷: Set of O–D port pairs and 𝐷 ≔ {(𝑜, 𝑑) ∈ 𝑃 × 𝑃|𝑜 ≠ 𝑑},
𝑣: Index of ship types,
𝑉: Set of ship types for the shipping service route,
𝑒: Index of round trips,
𝐸: Set of all the round trips in the planning horizon,
Input parameters:
𝑁: Number of ships deployed in the shipping service route,
𝑀: Number of foldable containers that can be folded into one standard container,
𝐶 𝑣 : Fixed operation cost in the planning horizon when the ships in Type 𝑣 are deployed in
the shipping service route, 𝑣 ∈ 𝑉,
𝐾 𝑣 : Container capacity of the ships in Type 𝑣 with the unit of TEUs, 𝑣 ∈ 𝑉,
𝑤𝑝 : Number of weeks that are needed for the devanning process in Port 𝑝, 𝑝 ∈ 𝑃,
𝑑𝑜𝑑,𝑒 : Number of laden container transportation consignments from Port 𝑜 to Port 𝑑 that
should be transported by the 𝑒 𝑡ℎ round trip, (𝑜, 𝑑) ∈ 𝐷, 𝑒 ∈ 𝐸,
𝑠𝑝𝑆 : Unit weekly storage cost of a standard container in Port 𝑝, 𝑝 ∈ 𝑃,
as the number of empty containers deficit or surplus in Port 𝑝 when visited by the 𝑒 𝑡ℎ round
trip, calculated by Eq. (2.1) and Eq. (2.2).
∀𝑝 ∈ 𝑃, 𝑒 ∈ 𝐸/{1}, (2.2)
If 𝑄𝑝,𝑒 > 0 (𝑄𝑝,𝑒 ≤ 0), it indicates that there are 𝑄𝑝,𝑒 (−𝑄𝑝,𝑒 ) number of empty containers
surplus (deficit) in Port 𝑝 when visited by the 𝑒 𝑡ℎ round trip. The long-term containers leasing
generates empty containers, and the empty container repositioning induces the empty container
flow between the deficit ports (i.e., the demand ports) and the surplus ports (i.e., the supply
ports). Henceforth, we can construct a flow network for the empty containers.
Figure 2.3: A sub-network for standard containers without considering the short-term
container leasing
The sub-network contains |𝑃| + 2 ∙ |𝐸| ∙ |𝑃| + 2 nodes. Among them, the |𝑃| nodes define
the flow conversion for the long-term container leasing in Port 𝑝 at the beginning of the
planning horizon, labelled as 𝐿𝑡𝑒𝑟𝑚 𝑆 _𝑝. The 2 ∙ |𝐸| ∙ |𝑃| nodes are categorized into |𝐸| ∙ |𝑃|
groups of two kinds of nodes. One kind of nodes in a group denotes the flow conservation in
the ship for the 𝑒 𝑡ℎ round trip after visiting Port 𝑝, labeled as 𝑆ℎ𝑖𝑝 𝑆 _𝑝_𝑒. The other kind of
nodes denotes the flow conservation in Port 𝑝 after visited by the 𝑒 𝑡ℎ round trip, labeled as
𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒. Two additional dummy nodes represent the source node (labelled as 𝑆𝑜𝑢𝑟𝑐𝑒) and
the sink node (labelled as 𝑆𝑖𝑛𝑘) respectively. To facilitate the understanding of the sub-
The arc from Node 𝑆ℎ𝑖𝑝 𝑆 _𝑝_𝑒 to 𝑆ℎ𝑖𝑝 𝑆 _(𝑝 + 1)_𝑒 represents the voyage leg from Port 𝑝 to
Port 𝑝 + 1 by the ship for the 𝑒 𝑡ℎ round trip. The amount of the flow on the arc denotes the
number of empty containers carried on the ship in the voyage leg. Here note that: (i) The
voyage leg from Port |𝑃| to Port 1 by the same ship for the 𝑒 𝑡ℎ round trip is the arc from
Node 𝑆ℎ𝑖𝑝_|𝑃|_𝑒 to 𝑆ℎ𝑖𝑝_1_(𝑒 + 𝑁) . (ii) The cost for the arc is zero, as we have
decomposed the repositioning cost. (iii) As the flows of laden containers are pre-determined,
the remaining capacity in each voyage leg (on each arc) is confirmed, denoted as 𝐾𝑝,𝑒 (the
remaining capacity in the voyage leg from Port 𝑝 to Port 𝑝 + 1 by the ship for the 𝑒 𝑡ℎ round
trip), which can be calculated as follows.
𝐾𝑝,𝑒 = 𝐾 − ∑𝑜≤𝑝,𝑑>𝑝,(𝑜,𝑑)∈𝐷 𝑑𝑜𝑑,𝑒 − ∑𝑑≥𝑝+1,𝑜>𝑑,(𝑜,𝑑)∈𝐷 𝑑𝑜𝑑,𝑒−𝑁 ∀𝑝 ∈ 𝑃, 𝑒 ∈ 𝐸, (2.3)
The arc from Node 𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 to 𝑆ℎ𝑖𝑝 𝑆 __𝑝_𝑒 (resp., Node 𝑆ℎ𝑖𝑝 𝑆 _𝑝_𝑒 to 𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 )
represents the empty container loading process from Port 𝑝 to the ship (resp., unloading
process from the ship to Port 𝑝) for the 𝑒 𝑡ℎ round trip. The amount of the flow on the arc
denotes the number of empty containers loaded to the ship (resp., unloaded to the port), and
the cost for the arc is 𝑎𝑝𝑆 (resp., 𝑏𝑝𝑆 ), i.e., the container loading cost at Port 𝑝 (resp., the
container unloading cost at Port 𝑝). The capacity on the arc is ∞.
The arc from Node 𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 to 𝑃𝑜𝑟𝑡 𝑆 _𝑝_(𝑒 + 1) represents the empty container
inventory (i.e., the empty containers left) in Port 𝑝 after visited by the 𝑒 𝑡ℎ round trip. The
amount of the flow on the arc denotes the number of empty containers left in the port after
visited by the 𝑒 𝑡ℎ round trip (i.e., the inventory level 𝛿𝑝,𝑒
𝑆
). The cost for the arc is 𝑠𝑝 . The
capacity on the arc is ∞.
With the two additional dummy nodes, i.e., Node 𝑆𝑜𝑢𝑟𝑐𝑒 and Node 𝑆𝑖𝑛𝑘 , we firstly
construct the dummy arcs from Node 𝑆𝑜𝑢𝑟𝑐𝑒 to all the nodes 𝐿𝑡𝑒𝑟𝑚 𝑆 _𝑝 with zero cost
nodes 𝐵𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 to Node 𝑆𝑖𝑛𝑘 such that 𝑄𝑝,𝑒 ≤ 0). The costs for all the arcs are zero. The
capacity on the arc from Node 𝑆𝑜𝑢𝑟𝑐𝑒 to 𝐵𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 is 𝑄𝑝,𝑒 (resp., from Node 𝐵𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒
to 𝑆𝑖𝑛𝑘 is −𝑄𝑝,𝑒 ). (iv) If 𝑄𝑝,𝑒 > 0 ( 𝑄𝑝,𝑒 ≤ 0), we construct the dummy arc from Node
𝐵𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 to 𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 (resp., from Node 𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 to 𝐵𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒) with arc capacity as
∞ and arc cost as zero. (v) For two nodes 𝐵𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒, there exist an origin node 𝐵𝑃𝑜𝑟𝑡 𝑆 _𝑜_𝑒
It is worthwhile to mention that in this subsection, we refer the short-term container leasing
to the Master Lease Agreement mentioned in Section 2.3.2, in which empty containers are
leased and the cost is charged on an O-D port pair basis. In real-world operations, the shipping
liner company and container leasing company could have a lease agreement that is on a round
trip basis. Under the agreement, the shipping liner company can lease a certain number of
empty containers at Port 𝑝, and this number of empty containers must be returned at Port 𝑝
after several weeks (or round trips). If this is the agreement applied between some shipping
liner companies and container leasing containers, we can make the following modifications to
the network. Firstly, we define 𝑔𝑝 as the number of weeks that the shipping liner company
needs to return leased empty containers at Port 𝑝, and define 𝑞𝑝 as the unit cost of leasing
empty containers for 𝑔𝑝 weeks at Port 𝑝. Then, for modifying the network, we connect the
that although some empty containers (say 𝜉 empty containers) are leased by the 𝑒 𝑡ℎ round trip
and are returned by the (𝑒 + 𝑞𝑝 )𝑡ℎ round trip at Port 𝑝, the empty container flow in the
network is in the opposite direction. This is due to that the empty container supply by the
(𝑒 + 𝑞𝑝 )𝑡ℎ round trip (resp., the 𝑒 𝑡ℎ round trip) at Port 𝑝 will decrease by (resp., increase by)
𝜉 empty containers for the returning process (resp., the leasing process).
Until now, we have constructed a whole sub-network for the problem when only considering
using standard containers, which leads to the preliminary NF model (𝑷𝑵𝑭 model). For 𝑷𝑵𝑭
model, we need to fulfill Node 𝑆𝑖𝑛𝑘 with total empty container demand ∑𝑝∈𝑃,𝑒∈𝐸(−𝑄𝑝,𝑒 )+ by
originating from Node 𝑆𝑜𝑢𝑟𝑐𝑒. The goal is to minimize the total cost through the sub-network,
which is a standard minimum cost flow problem.
2.4.3 The NF model for both standard containers and foldable containers
In this subsection, we incorporate the foldable containers into 𝑷𝑵𝑭 model for the NF model
that solves our problem given the ship capacity 𝐾 (𝑻𝑵𝑭 model). Firstly, we construct a sub-
network for foldable containers. Except from the short-term container leasing, the foldable
containers also have the long-term container leasing and the empty container repositioning.
Henceforth, the sub-network for foldable containers is similar to the sub-network shown for
the standard containers in Figure 2.3. Figure 2.5 shows an example for the sub-network for
foldable containers, in which the nodes 𝑆ℎ𝑖𝑝𝐹 _𝑝_𝑒, 𝑃𝑜𝑟𝑡𝐹 _𝑝_𝑒 and 𝐿𝑡𝑒𝑟𝑚𝐹 _𝑝 correspond to
the nodes 𝑆ℎ𝑖𝑝 𝑆 _𝑝_𝑒, 𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 and 𝐿𝑡𝑒𝑟𝑚 𝑆 _𝑝 in Figure 2.3, respectively. Note that the
costs on the arcs from Node 𝑃𝑜𝑟𝑡𝐹 _𝑝_𝑒 to 𝑆ℎ𝑖𝑝𝐹 __𝑝_𝑒 and from Node 𝑆ℎ𝑖𝑝𝐹 _𝑝_𝑒 to
𝑃𝑜𝑟𝑡𝐹 _𝑝_𝑒 are slightly different from that of standard containers, as there are additional
folding and unfolding costs (𝐴𝑜 and 𝐵𝑑 ) for foldable containers.
To embed the sub-network for the foldable containers into 𝑷𝑵𝑭 model defined in Section
2.4.2, we need to build some connections as shown in Figure 2.6: (i) we add a pairwise dummy
node for nodes 𝐵𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 and 𝑃𝑜𝑟𝑡𝐹 _𝑝_𝑒, denoted as 𝑇𝑃𝑜𝑟𝑡_𝑝_𝑒. Here, Node 𝑇𝑃𝑜𝑟𝑡_𝑝_𝑒
accumulates the empty container surplus or deficit for standard containers and foldable
containers; (ii) we disconnect all dummy arcs from Node 𝑆𝑜𝑢𝑟𝑐𝑒 to all other nodes or from
all other nodes to Node 𝑆𝑖𝑛𝑘; (iii) we construct the dummy arcs from Node 𝑆𝑜𝑢𝑟𝑐𝑒 to all the
nodes 𝑇𝑃𝑜𝑟𝑡_𝑝_𝑒 such that 𝑄𝑝,𝑒 > 0 (resp., from all the nodes 𝑇𝑃𝑜𝑟𝑡_𝑝_𝑒 to Node 𝑆𝑖𝑛𝑘 such
that 𝑄𝑝,𝑒 ≤ 0). The costs for all the arcs are zero. The capacity on the arc from Node 𝑆𝑜𝑢𝑟𝑐𝑒
to 𝑇𝑃𝑜𝑟𝑡_𝑝_𝑒 is 𝑄𝑝,𝑒 (resp., from Node 𝐵𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 to 𝑆𝑖𝑛𝑘 is −𝑄𝑝,𝑒 ); (iv) if 𝑄𝑝,𝑒 > 0, we
construct the dummy arcs from Node 𝑇𝑃𝑜𝑟𝑡_𝑝_𝑒 to 𝐵𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 and Node 𝑇𝑃𝑜𝑟𝑡_𝑝_𝑒 to
𝐵𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 to 𝑇𝑃𝑜𝑟𝑡_𝑝_𝑒 and Node 𝑃𝑜𝑟𝑡𝐹 _𝑝_𝑒 to 𝑇𝑃𝑜𝑟𝑡_𝑝_𝑒 to accumulate empty
container flow from the sub-network of standard containers and foldable containers,
respectively. Note that 𝑇𝑃𝑜𝑟𝑡_𝑝_𝑒 serve as the node to allocate the empty container demand
to standard containers by node 𝐵𝑃𝑜𝑟𝑡 𝑆 _𝑝_𝑒 and foldable containers by node 𝑃𝑜𝑟𝑡𝐹 _𝑝_𝑒.
Figure 2.6: Connections between standard containers and foldable containers in TNF model
As we have constructed the whole flow network for our problem, it is worthwhile to show
the total number of nodes involved in the network by referring Figure 2.3 to Figure 2.6. Figure
for the capacity restriction sharing. This characteristic leads to the trouble that we cannot
directly transfer our problem given the ship capacity to a standard minimum cost flow problem.
Alternatively, to avoid the trouble, we introduce two dynamic capacity restrictions for the two
𝑆 𝐹 𝑆 𝐹
parallel arcs, denoted as 𝜇𝑝,𝑒 and 𝜇𝑝,𝑒 , respectively. Based on 𝜇𝑝,𝑒 and 𝜇𝑝,𝑒 , we separate the
sharing capacity restriction 𝐾𝑝,𝑒 to two individual ones such that:
𝑆 𝑆
0 ≤ 𝑥𝑝,𝑒 ≤ 𝜇𝑝,𝑒 ∀𝑝 ∈ 𝑃, 𝑒 ∈ 𝐸, (2.5)
𝐹 𝐹
0 ≤ 𝑥𝑝,𝑒 ≤ 𝜇𝑝,𝑒 ∀𝑝 ∈ 𝑃, 𝑒 ∈ 𝐸. (2.6)
Meanwhile, the two individual capacity restrictions must hold that:
𝐹
𝜇𝑝,𝑒
𝑆
𝜇𝑝,𝑒 + = 𝐾𝑝,𝑒 ∀𝑝 ∈ 𝑃, 𝑒 ∈ 𝐸. (2.7)
𝑀
𝑆 𝐹
Note 𝜇𝑝,𝑒 and 𝜇𝑝,𝑒 keep dynamically changed in the algorithm developed in the next section,
𝑆 𝐹
so long as Eq. (2.7) holds. Here, we can initialize any values for 𝜇𝑝,𝑒 and 𝜇𝑝,𝑒 without
considering that the values will lead to no feasible solutions for the 𝑻𝑵𝑭 model. This is due
to that the short-term container leasing of standard containers can always treat the empty
The solution approach for our problem is an iterative procedure. In this section, we firstly
design a revised network simplex algorithm to solve the 𝑻𝑵𝑭 model given a ship capacity.
Then, for the solution approach, we initialize at 𝐾 = 𝐾 1 that is the capacity of Type 1 ship
(assuming to be the smallest ship type) for running the algorithm. Based on the information
from the previous running of the algorithm, we select another ship capacity 𝐾 𝑣 , given which
it again invokes the algorithm to derive the minimum flow cost for empty containers, denoted
as 𝜎 𝑣 . The summation of the fixed operation cost 𝐶 𝑣 and 𝜎 𝑣 is the total minimum cost when
using ships in Type 𝑣 for the fleet. By comparing those costs, we can obtain the optimal
solution for our problem as well as the total optimal cost. In this section, we also formulate a
mixed-integer linear programming model (MILP) for the problem, which can be solved
directly by CPLEX solver for verifying the optimality of the proposed solution approach.
nodes: it first assigns a node potential 𝜋(1) = 0 for the root node 1. Then by holding,
𝜋
𝑐𝑖𝑗 = 𝑐𝑖𝑗 − 𝜋(𝑖) + 𝜋(𝑗) = 0 ∀(𝑖, 𝑗) ∈ 𝑻, (2.8)
𝜋
it obtains all the nodes potentials 𝜋(𝑖) in the network, based on which 𝑐𝑖𝑗 of the non-tree arcs
𝜋 𝜋
is derived by 𝑐𝑖𝑗 = 𝑐𝑖𝑗 − 𝜋(𝑖) + 𝜋(𝑗), ∀(𝑖, 𝑗) ∈ 𝑳, (𝑖, 𝑗) ∈ 𝑼 . Note that 𝑐𝑖𝑗 has a similar
meaning with the reduced cost in the primal simplex algorithm for the linear programming
problem, which indicates the cost changed in the objective if we increase one more unit flow
on the arc (𝑖, 𝑗).
(iii) The network simplex algorithm generally maintains a feasible spanning tree structure
in each iteration and moves from one structure to another one until reaching the optimality. (a)
Given a (𝑻, 𝑳, 𝑼), the algorithm checks the optimality conditions. If all the conditions are
satisfied, the algorithm stops, otherwise selects the most violation arc (by finding the
𝜋 𝜋
maximum |𝑐𝑖𝑗 | ) from the arc set 𝒱 = {(𝑖, 𝑗) ∈ 𝑳 ∪ 𝑼: 𝑖𝑓 (𝑖, 𝑗) ∈ 𝑳, 𝑐𝑖𝑗 < 0; 𝑖𝑓 (𝑖, 𝑗) ∈
𝜋
𝑼, 𝑐𝑖𝑗 > 0}. (b) It adds the violation arc (called the entering arc) into the current spanning tree
structure, by which it obtains a negative cost cycle. That means increasing the flow on the
forward direction of the cycle will decrease the current objective of the total flow cost. (c) It
augments the maximum possible flow on the negative cost cycle until the flow of one arc in
the cycle reaches zero or its capacity restriction (called the leaving arc). (d) It replaces the
leaving arc with the entering arc for a new feasible spanning tree structure (𝑻′, 𝑳′, 𝑼′). The
algorithm repeats the above procedure (a)-(d) until finding the optimal spanning tree structure
(𝑻∗ , 𝑳∗ , 𝑼∗ ).
the arcs following the cycle direction belong to a forward arc set 𝑾, and the arcs following
the opposite cycle direction belong to a backward arc set 𝑾. Figure 2.7 shows an example for
the cycle, in which arc (3, 4) is the entering arc, arcs (5, 0), (0, 1), (2, 3) belong to the set 𝑾,
and arcs (5, 4), (2, 1) belong to the set 𝑾.
Supposing that 𝑥𝑖𝑗 is the existing flow on arc (𝑖, 𝑗), the maximum flow change on each arc
in the cycle is denoted by 𝛿𝑖𝑗 , calculated by:
𝒢 𝐹 ⊂ 𝒜. Supposing that the current allocated capacity for two parallel arcs are 𝜇𝑝,𝑒
𝑆 𝐹
and 𝜇𝑝,𝑒
satisfying Eq. (2.5)-(2.7), we will dynamically change the capacity allocation in the revised
pivot operation to reach the optimality and guarantee the feasibility. Here, we suppose that
𝑆 𝐹
𝑥𝑝,𝑒 and 𝑥𝑝,𝑒 are existing flow on the two parallel arcs.
In the revised pivot operation, we need to check each pair of parallel arcs on the cycle 𝑾.
Basically, we need to distinguish three situations for one of two parallel arcs, i.e., (i) the arc
belongs to the forward arc set 𝑾; (ii) the arc belongs to the backward arc set 𝑾; (iii) the arc
does not belong to the cycle 𝑾. Then, for an integration of two parallel arcs, there are nine
scenarios, for each of which we have different ways to change the capacity restrictions for
the parallel arcs.
𝑆 𝐹
Scenario 1: if 𝜖𝑝,𝑒 ∈ 𝑾 and 𝜖𝑝,𝑒 ∈ 𝑾, this is the scenario that needs the most attention, as
we need to care about their capacity restrictions simultaneously based on Eq. (2.8). To make
𝑆 𝐹 𝑆
Proposition 2.1: if 𝜖𝑝,𝑒 ∈ 𝑾 and 𝜖𝑝,𝑒 ∈ 𝑾, the arc 𝜖𝑝,𝑒 cannot be the leaving arc.
𝑆 𝐹 𝐹
Scenario 3: if 𝜖𝑝,𝑒 ∈ 𝑾 and 𝜖𝑝,𝑒 ∉ 𝑾, as the arc 𝜖𝑝,𝑒 is not in the cycle, its flow will not
𝐹 𝐹 𝑆 ∗ ∗
change during the revised pivot operation. Then, we adjust 𝜇𝑝,𝑒 = 𝑥𝑝,𝑒 and 𝜇𝑝,𝑒 = 𝐾𝑝,𝑒 −
𝐹 ∗ 𝐹
𝜇𝑝,𝑒 𝑥𝑝,𝑒 𝑆
= 𝐾𝑝,𝑒 − such that the capacity restriction on arc 𝜖𝑝,𝑒 is maximum.
𝑀 𝑀
𝑆 𝐹
Scenario 4: if 𝜖𝑝,𝑒 ∈ 𝑾 and 𝜖𝑝,𝑒 ∈ 𝑾, the operation is similar to Scenario 2, which is not
repeated here. We can also obtain a proposition that is similar to Proposition 2.1:
𝑆 𝐹 𝐹
Proposition 2.2: if 𝜖𝑝,𝑒 ∈ 𝑾 and 𝜖𝑝,𝑒 ∈ 𝑾, the arc 𝜖𝑝,𝑒 cannot be the leaving arc.
𝑆 𝐹
Scenario 5: if 𝜖𝑝,𝑒 ∈ 𝑾 and 𝜖𝑝,𝑒 ∈ 𝑾, we do not change the capacity restrictions.
𝑆 𝐹
Scenario 6: if 𝜖𝑝,𝑒 ∈ 𝑾 and 𝜖𝑝,𝑒 ∉ 𝑾, we do not change the capacity restrictions.
optimality. Note that the revised network simplex algorithm is applicable to any minimum cost
flow problems with sharing a capacity restriction among arcs.
Conditions. Starting from the ship capacity 𝐾 1 , if we increase the ship capacity to 𝐾 2 , the
fixed operation cost increases by 𝐶 2 − 𝐶 1 . However, the minimum flow cost may decrease as
we relax the ship capacity from 𝐾 1 to 𝐾 2 . Therefore, when increasing ship capacity, there is a
trade-off between the increasing of fixed operation cost and the decreasing of minimum flow
cost. As the increasing of the fixed operation cost is pre-determined (𝐶 2 − 𝐶 1 ), we need to
judge whether the decreasing of minimum flow cost can compensate for the increasing
operation cost 𝐶 2 − 𝐶 1 .
For notational convenience, we denote ∆1,2 = 𝜎 1 − 𝜎 2 as the decreasing of minimum flow
𝜋
cost by changing Type 1 ship to Type 2 ship. Here, we can use the reduced costs 𝑐𝑖𝑗 at ship
𝜋
capacity 𝐾 = 𝐾 1 to derive an upper bound for ∆1,2 . Recall that 𝑐𝑖𝑗 ≤ 0, ∀(𝑖, 𝑗) ∈ 𝑼 actually
means if the capacity of arc (𝑖, 𝑗) increases by one unit, the minimum flow cost has the
𝜋
potential to decrease by |𝑐𝑖𝑗 |. Considering a minimum cost flow problem is equivalent to
𝜋 𝜋
min{∑(𝑖,𝑗)∈𝑳 𝑐𝑖𝑗 𝑥𝑖𝑗 − ∑(𝑖,𝑗)∈𝑼 |𝑐𝑖𝑗 | 𝑥𝑖𝑗 } (Ahuja et al., 1993), ∆1,2 has the upper bound
𝜋
∑(𝑖,𝑗)∈𝑼|𝑐𝑖𝑗 | (𝐾 2 − 𝐾 1 ) if we do not have the parallel arcs and all arcs (𝑖, 𝑗) ∈ 𝑼 are restricted
by the ship capacity. However, under the situation that our problem has parallel arcs that
sharing capacity restrictions, the increased ship capacity 𝐾 2 − 𝐾 1 may be used by
repositioning standard containers or foldable containers. Thus, we propose Procedure 2.1 to
𝑈
derive the upper bound for ∆1,2 (denoted by ∆1,2 ), in which 𝑐𝜖𝜋𝑆 and 𝑐𝜖𝜋𝐹 denote the reduced
𝑝,𝑒 𝑝,𝑒
costs for the parallel arcs. Note that if we increase one unit ship capacity, the capacity on the
parallel arcs for foldable containers can increase by 𝑀 units. After obtaining the upper bound
𝑈 𝑈
∆1,2 , we can compare it with 𝐶 2 − 𝐶 1 . If ∆1,2 ≤ 𝐶 2 − 𝐶 1 , it is unnecessary to further invoke
the revised network simplex algorithm to derive 𝜎 2 for Type 2 ship, as the ∆1,2 is impossible
to compensate 𝐶 2 − 𝐶 1 . We will repeat the Procedure 1 by increasing ship capacity to 𝐾 𝑣 until
𝑈
we find a ship type 𝑣 such that ∆1,𝑣 ≥ 𝐶 𝑣 − 𝐶 1 , and then invoke the revised network simplex
To facilitate the understanding, we show an example to exclude ship types. Supposing that
we have invoked the algorithm to derive the reduced costs for the ship capacity at 𝐾 = 𝐾 1, we
𝑈 𝑈
derive ∆1,2 = 120 by Procedure 2.1. Then, if 𝐶 2 − 𝐶 1 = 200 ≥ ∆1,2 , it is unnecessary to
invoke the algorithm for the ship capacity at 𝐾 = 𝐾 2. In the next, we replace 𝐾 2 with 𝐾 3 in
𝑈
Procedure 2.1. If we have ∆1,3 = 300 and 𝐶 3 − 𝐶 1 = 250, we need to invoke the algorithm
Decision variables:
𝜀𝑣 : binary, set to one if the ships in Type 𝑣 are deployed, otherwise zero, 𝑣 ∈ 𝑉;
𝜁𝑝𝑆 : integer, number of long-term leasing standard containers in Port 𝑝, 𝑝 ∈ 𝑃;
Mathematical model:
[𝑴𝟐] min ∑𝑣∈𝑉 𝐶𝑣 𝜀𝑣 + ∑𝑝∈𝑃(𝐿𝑆𝑝 𝜁𝑝𝑆 + 𝐿𝐹𝑝 𝜁𝑝𝐹 ) + ∑𝑒∈𝐸 ∑(𝑜,𝑑)∈𝐷 𝑙𝑜𝑑 𝛾𝑜𝑑,𝑒 +
𝑆 𝑆 𝐹 𝐹
∑𝑒∈𝐸 ∑(𝑜,𝑑)∈𝐷(𝑟𝑜𝑑 𝛽𝑜𝑑,𝑒 + 𝑟𝑜𝑑 𝛽𝑜𝑑,𝑒 ) + ∑𝑒∈𝐸 ∑𝑝∈𝑃(𝑠𝑝𝑆 𝛿𝑝,𝑒
𝑆
+ 𝑠𝑝𝐹 𝛿𝑝,𝑒
𝐹
) (2.14)
In this section, based on three real-world shipping service routes, we conduct extensive
computational experiments to find insights on the ship type decision and the foldable container
usage by using our proposed solution approach. We run the experiments by a PC equipped
with 3.30GHz of Intel Core i5 CPU and 16GB of RAM. For all the test instances in the
experiments, the planning horizon is half a year, i.e., 26 weeks. As the shipping service is on
a weekly basis, the total round trips are 26 (i.e., |𝐸|=26).
Figure 2.9: Three selected shipping services routes operated by CMA CGM (CMA CGM,
2017)
According to Word Shipping Council (2013), we estimate that the imbalance ratio of laden
container flow for Asia-North America trade route is about 1.99 (i.e., Eastbound / Westbound),
By referring to Moon and Hong (2016), and Konings (2005), we set all the relevant costs
for foldable and standard empty containers, which are shown in Table 2.1. Here, we assume
that all the costs are the same for different ports, and four foldable containers can be folded as
one standard container. The short-term leasing cost is charged based on the travel time between
origin ports and destination ports. According to Moon and Hong (2016), the unit short-term
leasing cost is set as US$170/week, for example, if the travel time between an origin port and
a destination port is two weeks, the short-term leasing cost is US$340. For the three selected
Table 2.3: Comparing the proposed method with the MILP model by CPLEX solver
Table 2.4: Comparing the proposed model and the traditional model without empty
container allocation
Note: (i) 𝑍1(million US$) and 𝑍2(million US$) represent the objective values of the two models with
the unit of million US$. (ii) “Ship fleet” shows the capacity (in TEUs) of the selected ship type. (iii)
“Gap” shows the difference on the total costs by the two models, which is calculated by (𝑍2 − 𝑍1)/𝑍1.
The comparison between the proposed model and the traditional model is listed in Table
2.4, where “Ship fleet” shows the capacity of the selected ship fleet types by the two models,
and “Gap” shows the difference in the total costs by the two models. In the majority of the
instances in Table 2.4, the selected ship fleet types are different by the two models. More
importantly, the ship fleet capacity selected by the proposed model is no less than the ship
fleet capacity selected by the traditional model, which is consistent with our previous
foldable containers
Note: (i) 𝑍1(million US$) and 𝑍2(million US$) represent the objective values with the unit of million
US$. (ii) “Ship fleet” shows the capacity (in TEUs) of the selected ship type. (ii) “Container fleet”
shows the total number of empty containers used for the whole planning horizon, including the total
number of the empty containers owned initially and the long-term leasing empty containers. (iii) “Gap”
shows the difference on the total costs by the two models, which is calculated by (𝑍2 − 𝑍1)/𝑍1.
In Table 2.5, the container fleet (i.e., the total number of containers) using foldable
containers is no less than the container fleet that does not use foldable containers. This result
suggests that using foldable containers motivates the shipping liner to enlarge its container
fleet, which makes it more powerful to handle the laden container demands. Meanwhile, there
is nearly no advantage on cost reduction for the route LIBERTY2 by using foldable containers
as the total cost gaps are less than 0.10%, and using foldable containers has the biggest impact
on the route AANAANLCMA among the three routes. The results are in line with the results
of Table 2.4. However, the impacts of using foldable containers on the total cost for all the
three routes are not so significant as all the total cost gaps are less than 0.70%, which implies
that under the current cost settings, the shipping liner does not have a strong incentive on cost
reduction to use the foldable containers. This may be the reason why the foldable containers
are still not prevalent among the shipping services. To find strong incentives for using the
foldable containers for the shipping liner, we will analyze the effects of cost settings on the
foldable container usage in Section 2.6.5. Meanwhile, Table 2.5 shows the foldable container
usage could affect the ship type decision. In two instances of the route BOHAI (Instance
the usage of foldable containers, where 𝜁𝑝𝑆 and 𝜁𝑝𝐹 are the number of standard containers and
foldable containers used for the long-term container leasing. The bigger the ratio is the more
foldable containers are leased for the usage of the planning horizon.
Figure 2.10: Effect of the long-term leasing cost on the foldable container usage
Figure 2.10 illustrates the effect of the long-term leasing cost on the foldable container
usage, where the y-axis shows the ratio 𝜌, and x-axis indicates the long-term leasing cost of a
foldable container. In Figure 2.10, we keep the long-term leasing cost of a standard container
unchanged (i.e., US$480), and increase the long-term leasing cost of a foldable container from
US528 to US$1200. As can be seen, all three curves for the three shipping routes descend fast
when the long-term leasing cost increases, and a formal proof (See Appendix B) can verify the
non-increasing trend of using foldable containers. Under the current cost setting, i.e., the long-
Table 2.6: Total cost saving after using foldable containers under different long-term leasing
cost
long-term leasing
528 576 624 672 720 768 816
cost
Total cost saving 27.51% 22.26% 17.04% 12.49% 8.56% 5.44% 3.19%
long-term leasing
864 912 960 1008 1056 1104 1152
cost
Total cost saving 1.83% 0.86% 0.26% 0.11% 0.02% 0.00% 0.00%
Figure 2.11: Effect of the folding and unfolding cost on foldable container usage
percentage of foldable containers used under US$20 folding and unfolding cost (the baseline
setting), and 𝜌1 is the percentage under other costs, such as US$10 and US$30. In the
experiments, we set the long-term leasing cost to US$960, US$860, US$760, US$660 and
US$560, under each of which, we change the folding and unfolding cost to detect the
sensitivity.
Figure 2.12 shows the relationship between the sensitivity to the folding and unfolding cost
and the long-term leasing cost, where the y-axis shows the sensitivity ratio 𝜎, and x-axis
Figure 2.12: Long-term leasing cost dependent sensitivity to the folding and unfolding cost
An intuitive explanation for the phenomenon is that when the long-term leasing cost drops
to a certain level (say US$860), using foldable containers becomes nearly the same cost-
effective as using standard containers. At that level, the changes in the folding and unfolding
cost may have evident effectiveness towards the foldable container usage. However, if the
long-term leasing cost further reduces to a low level (say US$560), using foldable containers
is much more cost-effective than using standard containers (cf. Figure 2.10, approximately
97.8% foldable container usage). Henceforth, decreasing the folding and unfolding cost may
have a limited incentive to increase the foldable container usage. Based on the phenomenon,
we can obtain some managerial insights for shipping lines when popularizing the usage of
foldable containers. (i) If container leasing companies charge a high price for the long-term
leasing of foldable containers, it may not be economic to use foldable containers considering
2.6.7 Sensitivity analysis on the number of weeks for the devanning process
In Section 2.4.1, there is an input parameter 𝑤𝑝 showing the number of weeks allowed for
the devanning process. As discussed in Section 2.3.3, this parameter indicates the required
time for consignees to return empty containers. In essential, the parameter constructs a trade-
off between shipping liner companies and customers (or consignees). If 𝑤𝑝 becomes larger,
the customers will have more flexibility to deal with the cargo carried in containers, but this
would increase the opportunity cost for shipping liner companies, as they need empty
containers as soon as possible to fulfill next transportation consignments. Here, to investigate
the opportunity cost, we conduct a sensitivity analysis on the parameter 𝑤𝑝 , the results of
which are given in Table 2.7.
Table 2.7: Sensitivity analysis on the number of weeks allowed for returning empty
containers
No. of weeks 0 1 2 3 4 5 6 7
Total cost 70.93 75.53 85.48 96.13 107.35 118.80 129.98 140.81
Slope 4.60 9.95 10.65 11.22 11.45 11.18 10.83 NaN
Note: (i) “Total cost” with a unit of million US$ denotes the total cost on average of ten randomly
generated instances. (ii) “Slope” shows the increasing rate at a specific number of weeks, for
example, at “0” weeks, the rate is (75.53 − 70.93)/(1 − 0) = 4.60. (iii) “0” week for the devanning
process is unrealistic in the real-world operations. The “0” week setting only serves as a benchmark in
the sensitivity analysis.
In the table, when the number of weeks allowed for the devanning process increases, the
total cost grows significantly, which suggests the high opportunity cost for allowing more
devanning time. Meanwhile, the slope of total cost increases before reaching “4” weeks and
decreases after “4” weeks, which reveals that the total cost is increasing convex in the number
of weeks first and then becomes increasing concave. More importantly, given “0” weeks as
the benchmark, we can see that from “0” week to “1” week, the total cost increases by 4.60
million US$, which is far less than other increasing rates (e.g., 9.95 million US$ from “1”
week to “2” weeks). Here, we can derive a managerial insight for shipping liner companies:
Allowing one week for the devanning process is a better choice for shipping companies, which
2.7 CONCLUSION
This paper makes an explorative study on the ship type decision considering the empty
container repositioning and the foldable containers. Different from traditional research works
on the ship fleet deployment, our study incorporates both the laden container transportation
and the empty container repositioning into ship type decision in order to achieve the global
optimum for a shipping service route over a whole planning horizon. Meanwhile, as
researchers have shown the economic and logistical viability of foldable containers, the
problem also considers the use of foldable containers, which aims to find under what
conditions, the shipping liner needs to use the foldable containers in its liner shipping services.
In this study, we find that given the ship type with a certain capacity, the problem transfers
to a nonstandard minimum cost flow. Henceforth, we build a network flow model for the
problem by constructing a network. When considering standard containers and foldable
containers, trouble arises in the network construction that is some parallel arcs share the same
capacity restriction. To overcome this trouble, we design a revised network simplex algorithm
that changes the standard pivot operation. The algorithm is applicable to any minimum cost
flow problem with sharing capacity restrictions. Based on the algorithm, we develop a solution
approach by using reduced costs for excluding some ship type, which can find the optimal ship
type in the end. By using the solution approach, we conducted extensive numerical
experiments to find some managerial implications on the ship fleet deployment and the
foldable container usage.
Some useful managerial implications of this study are summarized from three perspectives.
(i) Ship type decision: when deciding the ship type deployed in a ship fleet, only involving
laden container transportation leads to sub-optimal solutions, as the possible empty container
repositioning affects the decision. After including foldable containers, a smaller ship type can
be expected to deploy, because foldable containers have the storage space advantage in the
empty container repositioning. (ii) Foldable container usage: under the current cost setting, it
is not cost-effective for shipping lines to use foldable containers, as the long-term leasing cost
is high. However, if the long-term leasing cost cuts down, using foldable containers is
encouraged, as foldable container usage is highly dependent on the long-term leasing cost. The
foldable container usage’s sensitivity to the folding and unfolding cost depends on the long-
2.8 REFERENCES
Ahuja, R. K., Magnanti, T. L., Orlin, J. B., 1993. Network Flows: Theory, Algorithms, and
Applications, New Jersey, USA.
APL, 2017. The perfect fit for all your different needs: your APL equipment guide,
https://www.apl.com/wps/wcm/connect/659ca9b6-7eb3-4c16-82de-
8de04abab7e9/Equipment+guide.pdf?MOD=AJPERES. Accessed on 10 August 2017.
Bell, M.G.H., Liu, X., Rioult, J., Angeloudis, P., 2013. A cost-based maritime container
assignment model. Transportation Research Part B 58, 58–70.
Brouer, B.D., Alvarez, J.F., Plum, C.E., Pisinger, D., Sigurd, M.M., 2013. A base integer
programming model and benchmark suite for liner-shipping network design.
Transportation Science 48(2), 281–312.
Cheung, R.K., Chen, C.Y., 1998. A two-stage stochastic network model and solution methods
for the dynamic empty container allocation problem. Transportation Science 32(2), 142–
162.
Container Auction, 2017. Shipping container lease agreements.
https://containerauction.com/read-news/shipping-container-lease-agreements. Accessed
on 29 June 2016.
COSCO, 2017. Aust port charges: import container detention,
http://www.fivestarshipping.com.au/australian-port-charges/container-detention/.
Accessed on 12 August 2017.
CMA CGM, 2017. Lines and Services, https://www.cma-cgm.com/products-services/line-
services. Accessed on 12 August 2017.
Crainic, T.G., Gendreau, M., Dejax, P., 1993. Dynamic and stochastic models for the
allocation of empty containers. Operations Research 41, 102–126.
Gelareh, S., Meng, Q., 2010. A novel modeling approach for the fleet deployment problem
within a short-term planning horizon. Transportation Research Part E 46(1), 76–89.
Hanh, L. D. 2003. The logistics of empty cargo containers in the Southern California region.
Final Report, Long Beach, CA.
Konings, R., 2005. Foldable containers to reduce the costs of empty transport? A cost-benefit
analysis from a chain and multi-actor perspective. Maritime Economics and Logistics
7(3), 223–249.
Lee, C.-Y., Yu, M.Z., 2012. Inbound container storage price competition between the
container terminal and a remote container yard. Flexible Services and Manufacturing
24(3), 320–348.
Li, J.A., Liu, K., Leung, S.C.H., Lai, K.K., 2004. Empty container management in a port with
long-run average criterion. Mathematical and Computer Modelling 40, 85–100.
Li, J.A., Leung, S.C.H., Wu, Y., Liu, K., 2007. Allocation of empty containers between multi-
ports. European Journal of Operational Research 182, 400–412.
Meng, Q., Wang, T.S., Wang, S., 2012. Short-term liner ship fleet planning with container
transshipment and uncertain container shipment demand. European Journal of
Operational Research 223, 96–105.
Meng, Q., Wang, S., 2011. Liner shipping service network design with empty container
repositioning. Transportation Research Part E 47(5), 695–708.
Meng, Q., Wang, S., 2012. Liner ship fleet deployment with week-dependent container
shipment demand. European Journal of Operational Research 222, 241–252.
Moon, I. K., Do, Ngoc A. D., Konings, R., 2013. Foldable and standard containers in empty
container repositioning. Transportation Research Part E 49, 107-124.
Moon, I. K., Hong, H., 2016. Repositioning of empty containers using both standard and
foldable containers. Maritime Economics & Logistics 18, 61-77.
Myung, Y. S., Moon, I. K., 2014. A network flow model for the optimal allocation of both
foldable and standard containers. Operations Research Letters 42(6-7), 484-488.
Ng, M.W., 2014. Distribution-free vessel deployment for liner shipping. European Journal of
Operational Research 238(3), 858–862.
This chapter addresses the reefer slot conversion problem for container freight transportation.
Given a fleet of container ships of varying capacity, a cost-efficient approach for improving
fleet utilization and reducing the number of delayed containers is to optimize the sequence of
container ships in a given string, a problem which belongs to the ship-deployment class. A
string sequence with ‘uniformly’ distributed ship capacity is more likely to accommodate a
random container shipment demand. The number of one’s total ship slots acts as a gauge of
the capacity of the container ships. Meanwhile, in reality, there are actually two types of ship
slots: dry slots and reefer slots. A dry slot only accommodates a dry container, while a reefer
slot can accommodate either a dry or a reefer container. The numbers of dry and reefer slots
for ships in a string are different. Therefore, in this study we propose a model that considers
both dry and reefer slots and use it to elucidate the optimal ship-deployment sequence. The
objective is to minimize the delay of dry and reefer containers when the demand is uncertain.
Furthermore, based on the optimal sequence deduced, the study also investigates the need to
convert some dry slots to reefer slots for the container ships.
3.1 INTRODUCTION
In a liner shipping network, the liner shipping company normally operates weekly-serviced
ship routes with fixed schedules. That is, for a given shipping route, there is a fleet of container
ships deployed such that the ports along the shipping route are visited on a weekly basis. In
order to provide weekly services for all the ports, the number of ships deployed in the fleet
should be the number of weeks needed to complete the shipping route (Bell et al., 2011; Meng
and Wang, 2012; Lin and Tsai, 2014; Wang, 2017). For example, if traversing a route takes
two weeks, the number of ships deployed should be two. When a container ship visits a port,
the containers that have arrived at the port during the past week would then be loaded onto the
ship to send them onwards to their destination ports. However, the ship’s capacity, as measured
by the number of ‘twenty-foot equivalent units’ (TEUs), is limited. Thus, some containers
cannot be loaded onto the ship, and will therefore be delayed for one week.
In a container ship, ‘‘slots’’ are used to accommodate containers. The number of slots in
a ship reflects the ship’s capacity. Generally, there are two types of slots for a container ship,
dry slots, and reefer slots. Reefer slots are equipped with an electrical outlet so that reefer
containers can be accommodated which have an integrated cooling unit. Reefer containers
cannot be placed in dry slots (due to lack of electricity supply) while dry containers can be
1 ARCHIMIDIS 7,943
2 CMA CGM NABUCCO 8,488
3 CSCL EAST CHINA SEA 10,036
4 CSCL SOUTH CHINA SEA 10,036
5 NAVARINO 8,533
6 XIN DA YANG ZHOU 8,533
In maritime studies, a great deal of research effort has been devoted to container ship fleet
deployment problems. This effort is mainly focused on determining the ship type for a given
shipping service route (Gelareh and Meng, 2010; Meng et al., 2012; Song and Dong, 2013;
Ng, 2014). In such studies, it is assumed that the container ships are categorized into different
types. Moreover, all the ships belonging to a particular category are homogenous, i.e. they
have the same capacity. In addition, the ship fleet deployed for a given route only consists of
one type of ship. However, in reality, the ships in a fleet may not all have the same capacity
(Lin and Liu, 2011; Du et al., 2017). For instance, consider the example listed in Table 3.1
(relating to the China–USA Yangtze Service route operated by France’s CMA CGM Group,
Figure 3.1 depicts the sequence of container ships corresponding to a given string of the
weekly-serviced Yangtze Service shipping route. The sequence shown is 1→4→6→5→2→3
(which is equivalent to 4→6→5→2→3→1, and the other sequences obtained by cyclically
permuting the original string, as the string of ships forms a loop). Without losing generality,
we assume that the sequence is written so that the first ship in the sequence is the ship with the
smallest capacity (e.g. ship 1 for the Yangtze Service route). Given such a string, we derive in
this section some equations that can be used to calculate the weekly profit of the shipping route.
Note that all the equations derived here are considered to be on a weekly basis. In the following
subsections, the variables and parameters used in the equations will be elaborated upon in
detail. We also outline some practical container delay and rejection rules to be applied when
the container ships visit the ports.
When a container ship visits a port, the containers accumulated in the past week need to be
loaded into the available slots in the container ship for shipment to their destination ports. If
the slots in the container ship are insufficient to stow all the accumulated containers, then some
containers m be delayed for one or more weeks. Such a delay incurs additional costs for the
shipping line. The main costs are: (i) the cost of storing the delayed container in the yard space
of the port, (ii) the cost to customer satisfaction (who, presumably, will not be happy about the
delay), and (iii) the cost incurred supplying electricity to a delayed container if it is of the
reefer variety.
In the following week, when the next container ship visits the port, the containers that
have been delayed will have priority when it comes to being loaded onto the ship. If some
delayed containers still cannot be transported due to the limited availability of slots in the ship,
then those delayed containers must be transported using slots from other shipping lines or by
transshipment (Hasheminia and Jiang, 2017; Jiang et al., 2017). Such containers subsequently
transported by other shipping lines are referred to as ‘‘rejected containers’’. Container
rejection is the last thing that the shipping liner wants as it leads to two consequences for the
shipping line: (i) it has to pay the freight rates for the shipment provided by the other shipping
line, and (ii) it loses the goodwill of the customers affected. Here, one-week delays reflect the
service level offered by the shipping liner with respect to guaranteeing the transit time for
transportation of the container. We can adjust it to two or more weeks delay depending on the
policy of the shipping liner. Meanwhile, container rejection does not mean the shipping line
will lose profit. The marginal profit associated with the container transportation will decrease,
but can still be positive after some additional incurred costs are deducted.
Based on the above discussion, the decision variables and parameters required are as
follows.
Decision variables:
Dvd1 Dvr1 : The number of dry (reefer) containers that were delayed from ship v – 1 in
Dvd Dvr : The number of dry (reefer) containers that are delayed from ship v in the
current week.
D
d
D
r
: The average number of dry (reefer) containers that are delayed.
Rvd Rvr : The number of dry (reefer) containers that are rejected from ship v in the
current week.
pnd pnr : The mass probability that n dry (reefer) containers need to be transported in the
current week.
qd qr : The expected number of dry (reefer) containers that need to be transported in the
current week.
d
r
: The realization of new dry (reefer) container demand in the current week.
Evdr : The capacity of reefer slots that are available to dry containers in the current week.
Nd Nr : The maximum number of dry (reefer) containers in the current week.
d
mass function of the random variable is assumed to be known based on historical data:
d d
Pr( n) pnd , n 0,1,..., N d . The expectation value (E) of is q d : E( ) n 0 npnd .
d Nd
r
The symbols , N , pn and q are correspondingly defined for reefer containers. Notice that
r r r
we consider the container demands on the ‘‘hit-haul’’ leg (i.e. the leg on which the highest
number of containers is carried) of the long-haul liner service routes rather than all the legs in
the shipping route. For instance, for Asia–Europe service routes, the leg after the last port in
Asia is normally the hit-haul leg, as the container demands from Asian ports to European ports
d r
are much higher than in the opposite direction. Thus, we only have and to denote the
demands on the hit-haul leg for dry and reefer containers, respectively.
Suppose that the V ships are in the sequence given by the string 1 2 v V . The dry
(reefer) container capacity of the ship v is Evd ( Evr ). In a particular week, when the ship v
on the values of Dv 1 and Dv 1 (which were determined one week ago) and the values of
d r d
and (which are just observed), the shipping line needs to determine the number of dry
r
d r
(reefer) containers to postpone. Denote this quantity by Dv ( Dv ) and the number they need
d r
to reject by Rv ( Rv ). We analyze the decisions as follows:
Dvr1 r Evr , and Dvd1 d Dvr1 r Evd Evr , then all of the containers will be
(iii) We assume that reefer containers have higher priority because they bring in more profit
than dry containers. If Dv 1 Ev , then not all reefer containers can be transported
r r r
immediately, and we allow some reefers to be postponed to the next week. Note that in
r
the following week, ship v 1 with reefer capacity of Ev 1 will arrive. If
Dvr1 r Evr Evr1 and if all of the Dvr1 r Evr reefers are postponed, then some of
them still cannot be transported in the following week. In this study, we assume that a
container (dry or reefer) can only be postponed by one week. Hence, if Dv 1 Ev
r r r
Dvr Dvr1 r Evr and Rvr 0 ; if Dvr1 r Evr and Dvr1 r Evr Evr1 , then Evr1
(iv) We now analyze the dry containers. We allow dry containers to occupy reefer slots that
are not used in the current week, but we do not allow dry containers to reserve reefer slots
in the following week. The available reefer slots in the current week are
Evdr : max 0, Evr ( Dvr1 r )
. If Dv 1 Ev Ev and Dv 1 Ev Ev Ev 1 ,
d d d dr d d d dr d
Rvd 0 ; if Dvd1 d Evd Evdr and Dvd1 d Evd Evdr Evd1 , then Evd1 dry
The shipping line aims to maximize their profit per week. To this end, they try to transport as
many containers as possible. However, due to the randomness of the demand, it may occur
that in some weeks the demand is very high and not all the containers can be transported. In
g d qd R
d
g q R .
r r r
(3.1)
The container shipment also incurs some costs. First, if a dry (reefer) container is rejected,
the cost associated with loss of goodwill will be denoted by c d ,1 ( c r ,1 ). Second, reefer
containers need electricity to maintain the desired temperature during transportation, and the
cost of the electricity/fuel when transporting a reefer container is c r ,2 . Third, if a dry (reefer)
container is postponed, it has to be stored in the container yard for one week, and so we
represent the corresponding storage cost of using the yard space by c d ,3 ( c r ,3 ). Fourth, the
additional electricity cost when storing a delayed reefer container is denoted by c r ,4 . Fifth, if
a dry (reefer) container is delayed, the cost of customer dissatisfaction is denoted by c d ,5 ( c r ,5 ).
d r
Denote the average number of delayed dry (reefer) containers per week by D ( D ).
Then, the expected cost per week is given by the expression
d r
c d ,1 R c r ,1 R c r ,2 q r R r
c d ,3
c d ,5 D c r ,3 c r ,4 c r ,5 D . (3.2)
d r
g d qd Rd g r qr Rr
c d ,1 R d c r ,1 R r c r ,2 q r R r c d ,3 c d ,5 D d c r ,3 c r ,4 c r ,5 D r
(3.3)
g q g c
d d r r ,2
q r
g d c d ,1 R d g r c r ,1 c r ,2 R r c d ,3 c d ,5 D d c r ,3 c r ,4 c r ,5 D r .
Note that the first term in the right-hand side of Eq. (3.3) is a constant and is independent of
the sequence of the ships in the string.
In this section, we use the profit equation derived above to optimize the sequence of ships in
the string in order to maximize the weekly profit. Note that sequence optimization is a tactical
level decision, which is unchanged over the weeks. There are two major reasons why the
sequence may need to be re-optimized:
(ii) The demand probability changes. In other words, the probability mass function changes
to a different one. For example, new competitors may enter the shipping market which
will affect the demand pattern for existing shipping liners.
Given the capacities of the ships in the fleet and demand probability mass function, we
design a simulation-based approach for the optimization process. Given V ships, any two of
which have different capacity, there will be S (V 1)! different possible ship sequences (as
the sequence forms a loop, it does not matter which ship is chosen to be the first one). Let the
sequence be denoted by s (1, 2, , S ) and the corresponding weekly profit be P ( s ) , as
d r d r
determined by Eq. (3.3). (With the obvious notation that D ( s ) , D ( s ) , R ( s ) , and R ( s ) are
d r d r
used in place of D , D , R , and R , respectively.) Mathematically, therefore, the weekly
profit of sequence s, P ( s ) , is given by:
P ( s ) g d q d g r c r ,2 q r
(3.4)
g d c d ,1 R d ( s ) g r c r ,1 c r ,2 R r ( s ) c d ,3 c d ,5 D d ( s ) c r ,3 c r ,4 c r ,5 D r ( s ) .
d r d r
We use Monte-Carlo simulations to calculate D ( s ) , D ( s ) , R ( s ) , and R ( s ) in
Eq. (3.4). The essence of the approach is as follows. (i) Define the number of weeks T we wish
to simulate (e.g. T 100, 000 weeks). (ii) Randomly generate the demand for each week using
the probability mass function. The demand for dry (reefer) containers in the week t are
denoted by td ( tr ). (iii) Simulate the decision-making process from weeks 1 to T. For each
week t, record the number of dry (reefer) containers postponed Dtd ( s) ( Dtr ( s) ) and rejected
d r d r
Rtd ( s) ( Rtr ( s ) ). (iv) The quantities D ( s ) , D ( s ) , R ( s ) , and R ( s ) can be subsequently
estimated using
1 T d 1 T 1 T 1 T
Dt ( s ), D ( s ) Dtr ( s ), R ( s ) Rtd ( s ), R ( s ) Rtr ( s ) . (3.5)
d r d r
D (s)
T t 1 T t 1 T t 1 T t 1
By substituting the results from Eq. (3.5) into Eq. (3.4), we can calculate P ( s ) , and
therefore choose the sequence with the maximum profit:
s * arg max P ( s ) . (3.6)
s 1,2 S
To find the optimal sequence, we apply the two-stage simulation method proposed by Nelson
et al. (2001). In the first stage, a given number of weeks of the process are simulated for each
Step 0: (i) Select a value for to represent a practically significant difference. For instance,
could be set at 1000 US$/week, meaning that we are indifferent to two solutions if
their expected difference in weekly profit is less than 1000 US$. This also implies that
if we find a sequence that is not the optimal one, but has a profit less than that of the
optimal one by at most , then we also consider it to be an optimal solution. (ii) Select
that the chance that the found solution is the optimal one is at least 90% (as mentioned
Step 1: Select a confidence level 1 0 for the first stage. For example, if 0 is 5%, it means
that the probability that the optimal sequence is not removed in the first stage is at
least 95%.
Step 2: Select a value for the grouping parameter U, e.g., 30. To appreciate the significance
of this parameter, consider the following discussion of the pertinent random variable.
The average profit of sequence s over U V weeks is given by PUV ( s ) . The central
parameter T1 that is related to the sample size of the first stage (e.g. T1 might be 50).
PUV ,1 ( s), PUV ,2 ( s), , PUV ,T1 ( s) . Now compute the sample mean
1 T1
PUV ( s ) (1) PUV ,l ( s ) ,
T1 l 1
(3.7)
where the superscript ‘‘(1)’’ means the first stage and sample variance
which is the 1 0
1 ( S 1)
quantile of the t distribution with T1 1 degrees of
freedom. Let
1
1 2
1 (3.10)
Ws , s ' : Var PUV ( s ) Var PUV ( s ') , s 1, 2, , S, s ' 1, 2, , S, s s'
T1 T1
Remove all of the s ' 1, 2, , S if there exists an s 1, 2, , S with s s ' such that
1 5% it means that, if the set I contains the optimal sequence, then the chance
that the optimal sequence will be identified in the second stage is at least 97.5%. For
h 2
T2 s max T1 , Var PUV ( s ) (3.12)
where h : h(1 1 , T1 , S ) is Rinott’s constant and rounds up to the next largest
integer.
Step 5: Simulate the process for each sequence s I for more (T2 s T1 )UV weeks. Then,
Our exploratory experiments show that the above algorithm is time-consuming and not
able to find the optimal sequence using a ‘‘reasonable’’ amount of computation time. (All
For instance, if PUV ,l ( s ) is 5,343,000 US$ and 1000 , the weekly profit becomes 5,343
where the units used are 1,000 US$. Such a step allows us to deal with the above problem and
accelerate the algorithm. Table 3.2 shows the computation times required for different values
of .
Table 3.2: Comparison of computation times for different values of the parameter 𝜷
As can be seen from Table 3.2, increasing the value of accelerates the algorithm.
However, we cannot increase as much as we like in order to simply accelerate the algorithm,
as there are repercussions. When 1000 , we can only claim that the selected sequence is
the optimal one to an accuracy of one unit, i.e. 1,000 US$. However, this will not necessarily
be the optimal sequence to the desired accuracy of 1 US$. As shown in the table, different
values lead to different sequences being selected. Therefore, we need to make a comprise
between speed and accuracy. In this study, we select 10 . This should be sufficiently close
In this section, we use the revised two-stage simulation algorithm to optimize the sequence of
a string of some particular cases of interest. First of all, according to EPRI (2010) and Ting
and Tzing (2004), we have to fine-tune the input parameters required. To this end, the cost
coefficients which we shall use were collated and are presented in Table 3.3. Here, we take
c d ,3 c r ,3 , the storage cost incurred to keep a dry (reefer) container in the yard for one week,
as an example of how we fine-tuned the costs. According to EPRI (2010), the storage space in
a yard is charged at a rate of US$0.21 per square foot per day. Thus, a twenty-foot standard
container corresponds to a weekly storage cost of US$30.
Table 3.3: The relevant input costs used in the case studies
We focus in this work on three shipping routes operated by CMA CGM: ‘‘Northwest
Express’’, ‘‘South Atlantic Express’’, and ‘‘Europe Pakistan India Consortium 1’’. In the
interest of brevity, we denote the three shipping routes as ‘‘1N’’, ‘‘2S’’, and ‘‘3E’’,
respectively. Figure 3.2 indicates the port rotations of these three shipping routes.
1N
3E
Figure 3.2: The shipping routes involved in three case studies (CMA CGM, 2017)
Some key information about the shipping routes and ships deployed are given in
Table 3.4 — further details can be found by referring to CMA CGM (2017). During string
optimization, the number of ships deployed is another key element that can affect the scale of
the problem. This is because the number of ships determines the number of possible sequences
to consider. Therefore, we chose three shipping routes that deploy different numbers of ships
(determined by the rotation time, considering the weekly service frequency). We assume that
the dry- and reefer-container demands follow a uniform distribution with their ranges as given
in Table 3.4. The last column in this table shows the variance of the total container demand
for both dry and reefer containers.
Table 3.4: Information on the shipping routes and ships deployed
Route Ship no. Ship name Capacity (TEU) Dry slots Reefer slots
COSCO
1 9,469 8,769 700
GUANGZHOU
2 COSCO INDONESIA 8,501 8,501 0
1N 3 COSCO JAPAN 8,501 7,801 700
4 COSCO PACFIC 10,020 9,220 800
5 COSCO PHLIPPINES 8,501 7,801 700
6 COSCO PUSAN 9,572 8,872 700
1 E.R. LONDON 6,008 5,208 800
2 MSC KATYAYNI 5,711 5,177 534
3 MSC KRYSTAL 5,782 5,222 560
2S 4 MSC MARGARITA 5,770 5,138 632
5 MSC ORIANE 5,782 5,082 700
NORTHERN
6 6,732 6,232 500
MAJESTIC
7 RIO BARROW 5,551 5,001 550
1 APL CHARLESTON 9,336 8,322 1,014
CMA CGM
2 9,130 8,530 600
AMAZON
CMA CGM
3 9,130 7,630 1,500
URUGUAY
3E 4 MSC ALBANY 8,762 7,762 1,000
5 MSC ALGHERO 8,827 8,827 0
6 MSC SILVANA 8,400 7,700 700
7 MSC TOMOKO 8,400 7,700 700
8 UASC AL KHOR 9,400 8,900 500
Based on the information presented on the shipping routes, we implemented the revised
two-stage simulation algorithm to find the optimal sequence to employ. The results are shown
in Table 3.6. The table shows that the string optimization process changes all of the original
sequences used in the shipping routes. The last column shows the difference or gap in profit
between the original and optimized sequences (which are all positive). As tactical-level
decisions, these increases in profit can be claimed to be significant considering that the string
optimization procedure is inexpensive, i.e. the approach is cost effective. The results prove the
effectiveness of using the proposed algorithm for string optimization. One remarkable feature
of the table is that the string optimization procedure produced a much larger profit increase for
The next natural question to ask is whether a shipping line should convert some of the dry slots
of a ship into reefer slots. Before addressing this question, we first consider the feasibility of
carrying out such conversions.
Reefer containers have an integral refrigeration unit with a water-cooling system to keep
the cargo cold. This refrigeration system needs an external power supply when the container
is stored in a ship. Therefore, a reefer slot has to be equipped with an electrical outlet which is
connected to the power system of the ship. Thus, an electrical outlet or plug needs to be
installed to convert a dry slot to a reefer slot. Technically, the conversion process is not much
trouble. Container ships usually have independent power sub-distribution panels that supply
power connections for refrigerated containers. When there is a need for more electrical outlets
to create a reefer slot, therefore, a ship can simply group several electrical outlets and supply
them with electricity using one power cable connected to a power sub-distribution panel (DNV
GL SE, 2015).
The slot conversion problem is addressed in the following way. Suppose that the cost of
converting one dry to one reefer slot is ĉ (US$/week). We need to determine how many dry
slots we should convert. Note that ĉ is an average cost per week. The cost of slot conversion
is mainly that associated with installing the electrical outlet or electrical plug which is used to
power the reefer container. According to EPRI (2010), the installation cost is approximately
1,250 US$ per electrical outlet. As we measure the profit on a weekly basis, we transfer the
installation cost to a depreciation cost of 48 US$ per week (i.e. cˆ 48 ). This assumes that
such a tactical level decision (i.e. decision to implement string optimization and slot
Input: V ships whose identities (IDs) are labeled 1, 2, ,V . The initial number of reefer slots
Step 0: Find the optimal sequence using the revised two-stage simulation algorithm. The
u th ship in the sequence ( u 1, 2,...,V ). Thus, we have found PUV ( ) , that is, the
optimal weekly profit derived using the revised two-stage simulation algorithm.
Step 1: If the reefer slots of all of the V ships have reached the limit R, then stop the process.
Step 3: If the reefer slots on the ship with ID (u ') have reached the limit R, then set
Step 4: Temporarily convert one dry slot on the ship with ID (u ') into a reefer slot. The
resulting new sequence is denoted by ' , which is the same as sequence except
that ship '(u ') has one more reefer slot (and, of course, one fewer dry slot) than the
(i) If PUV ( ') PUV ( ) cˆ , permanently convert the dry slot into a reefer one and
set ' (i.e. permanently convert a dry slot into a reefer slot on the ship
(u ') ). Go to Step 1.
Step 5: If u ' V , go to Step 4. Otherwise, this means that just converting dry slots without
changing the sequence is not economically viable. Hence, we need to check what
Step 6: If the reefer slots of all of the V ships have reached the limit R , stop.
Step 8: If the reefer slots of all of the V ships have reached the limit R , set v ' v ' 1 and
go to Step 10.
Step 9: Temporarily, convert one dry slot on the ship with ID v ' into a reefer slot. Use the
'' (i.e. permanently convert a dry slot into a reefer one on the ship whose
better limit to slot conversion can be derived by using information about the electrical loads
in the container ships and electricity usage of the reefer containers. For instance, EPRI (2010)
estimates that an electric reefer container needs 2.8875 kW per hour on average. The container
ship Hanjin Paris has a generator installed with a capacity of 7,600 kW and a load factor (%
of capacity) of 63%. This indicates that the ship has 2,812 kW of electrical power available
(Khersonsky et al., 2007). Some detailed information on electricity demand in container ships
can be found in the work of Zis et al. (2014). In the Hanjin Paris case, the container ship can
be equipped with an additional 973 reefer slots (at most) considering the generator’s
limitations (i.e. 2812/2.8875 slots can be supplied). Normally, the larger the capacity of the
ship, the greater the engine power available from the ship to support reefer slots (Zis et al.,
2013).
String optimization is the first step carried out in our research problem. We now want to use
the slot conversion algorithm (which embeds the revised two-stage optimization algorithm to
optimize strings) to conduct a further investigation of the three case studies given in Section
3.3.2 (i.e. using the shipping routes and information given in Tables 3.4 and 3.5).
Using the same input parameters, we ran the slot conversion algorithm to obtain the
number of reefer slots to convert. The results are shown in Table 3.7. As can be seen, all the
ships involved have some dry slots that are converted to reefer slots, apart from one (namely,
the CMA CGM URUGUAY in shipping route 3E which already has the maximum number of
reefer slots permitted in the conversion algorithm). This verifies that the benefits gained by
the greater flexibility of reefer slots outweigh the conversion costs incurred to change the
existing slot configurations of the ships.
The last column in Table 3.7 shows the ‘‘optimal increment’’ that the optimal change in
a number of reefer slots represents. This quantity corresponds to the ratio of the optimal
In the above cases, we compare weekly profits after slot conversion and after string
optimization. The profit improvement is not so clear cut, as we optimized the ship fleet using
string optimization in the first place. In this subsection, we conduct two further case studies
based on two shipping routes with very large shipping fleets and extremely high demand
variances. Our intention is to provide further motivation for integrating string optimization
with slot conversion.
The two routes selected are the French Asia Line 1 and Columbus JAX (both belonging
to CMA CGM) which involve 12 and 17 ships, respectively. CMA CGM (2017) gives detailed
information on the two shipping routes but some of the basic information of interest is shown
in Table 3.9. Compared with the previously considered routes (Table 3.4), these routes have
many more ships involved and their variance in demand is much larger.
We first calculated the weekly profits of the original sequences used in the shipping
routes (Table 3.10). Then, we used the slot conversion algorithm (with string optimization) to
optimize the sequence of ships and convert some dry slots to reefer slots. Note that the slot
conversion limit was increased to 2,500 in these two case studies as the deployed ships have
larger capacities. Table 3.10 shows the effectiveness of the proposed approach. Compared to
the original sequences, the weekly profits of the optimized fleets were improved by 9.06%
(French Asia Line 1) and 7.90% (Columbus JAX). These profit increases are very significant
considering the cost-efficiency of the methods used. Note also that the two routes have larger
dry-to-reefer conversion rates (8.63% and 7.12%) compared to those found in the three
previous cases. This result verifies the previous finding that slot conversion is more critical for
shipping routes with large fleets of ships and high demand variances.
Table 3.10: Results obtained using string optimization and slot conversion
3.5 CONCLUSION
This paper presents an improved algorithm to search for the optimal number of reefer slots to
have on a container ship. It is assumed that all the relevant parameters (e.g. freight rates,
storage costs, etc.) are already known. We first used a revised two-stage simulation approach
to optimize the sequence of ships deployed. Based on this, we then formulated a slot
conversion algorithm to determine the optimal slot configurations of the ships, which embeds
the two-stage simulation algorithm for string optimization.
In this study, we used several real shipping routes operated by CMA CGM to highlight
the effectiveness of our approach. Our results reveal that the algorithm is highly efficient and
can help shipping liners to significantly improve their profits. However, there are also several
issues that are worth studying further in future work:
(i) When converting the dry slots to reefer slots, draft and load capacities are not taken into
consideration. This would be of great use when making ship stowage plans.
3.6 REFERENCES
Arduino, G., Carrillo Murillo, D., & Parola, F., 2015. Refrigerated container versus bulk:
evidence from the banana cold chain. Maritime Policy & Management 42(3), 228-245.
Bell, M. G., Liu, X., Angeloudis, P., Fonzone, A., & Hosseinloo, S. H., 2011. A frequency-
based maritime container assignment model. Transportation Research Part B:
Methodological 45(8), 1152-1161.
Cheaitou, A., & Cariou, P., 2012. Liner shipping service optimisation with reefer containers
capacity: an application to northern Europe–South America trade. Maritime Policy &
Management 39(6), 589-602.
Chen, J., & Yahalom, S., 2013. Container slot co-allocation planning with joint fleet agreement
in a round voyage for liner shipping. The Journal of Navigation 66(4), 589-603.
Cool Logistics, 2014. The Hanjin crisis: what about the reefer cargoes?
http://coollogisticsresources.com/global/the-hanjin-crisis-what-about-the-reefer-cargoes.
Accessed 18 January 2017.
CMA CGM, 2016. Line details: Yangtze Service. https://www.cma-cgm.com/products-
services/line-services/flyer/YANGTZE. Accessed 30 June 2016.
CMA CGM, 2017. Products & Services: Search Line Services. http://www.cma-
cgm.com/products-services/line-services. Accessed 20 January 2017.
DNV GL SE, 2015. Rules for Classification and Construction: Chapter 3 Electrical
Installations. Hamburg, Germany. http://rules.dnvgl.com/docs/pdf/gl/maritimerules/gl_i-
1-3_e.pdf. Accessed 29 April 2017.
Du, Y., Meng, Q., & Wang, S., 2017. Mathematically calculating the transit time of cargo
through a liner shipping network with various trans-shipment policies. Maritime Policy
& Management 44(2), 248-270.
EPRI., 2010. Electric Refrigerated Container Racks: Technical Analysis.
http://www.epri.com/abstracts/Pages/ProductAbstract.aspx?ProductId=0000000000010
19926. Accessed 18 January 2017.
This chapter addresses the cruise itinerary schedule design problem for a cruise ship.
This problem determines the optimal sequence of a given set of ports of call (a port of call is
an intermediate stop in a cruise itinerary) and the arrival and departure times at each port of
call for maximizing the monetary value of the utility at ports of call minus the fuel cost. To
solve the problem, in view of the practical observations that most cruise itineraries do not have
many ports of call, we first enumerate all sequences of ports of call and then optimize the
arrival and departure times at each port of call by developing a dynamic programming
approach. To improve the computational efficiency, we propose effective bounds on the
monetary value of each sequence of ports of call, eliminating non-optimal sequences without
invoking the dynamic programming algorithm. Extensive computational experiments are
conducted and the results show that, first, using the bounds on the profit of each sequence of
ports of call considerably improves the computational efficiency; second, the total profit of the
cruise itinerary is sensitive to the fuel price and hence an accurate estimation of the fuel price
is highly desirable; third, the optimal sequence of ports of call is not necessarily the sequence
with the shortest voyage distance, especially when the ports do not have a naturally
geographical sequence.
4.1 INTRODUCTION
A cruise itinerary is a cruise route operated by a cruise company: A cruise ship picks up
cruise passengers at an embarkation port, calls at a set of ports of call for cruise passengers to
visit the port cities, and returns to a disembarkation port where cruise passengers get off the
cruise ship. The ship that is deployed on the itinerary, the embarkation port, the sequence of
ports of call, the disembarkation port, and the time schedule are all pre-determined in a cruise
itinerary. Cruise ships are different from other ships such as tankers, bulk carriers and
containerships in that transportation is not the purpose of cruise ships.
The cruising industry has maintained a steady increase in supply for the past 20 years. In
2014, the number of cruise passengers reached a total of 22.04 million and the global cruise
industry generated revenues of 37.1 billion U.S. dollars (Statista, 2015). Meantime, the world
cruise fleet had 296 ships (Cruise Industry News, 2015) with a total of 482,000 lower berths 1
(Statista, 2015). The Caribbean and the Mediterranean areas are the most important cruising
1It is often considered that one cabin has two beds (two lower berths) when calculating the capacity of cruise ships.
Any extra beds in a cabin are referred to as “upper berths”. The actual average number of beds per cabin in a
cruise ship is usually higher than two.
Figure 4.1: Itinerary of 7 Day Western Caribbean of Carnival (Carnival Cruise Line, 2016)
In the CISD problem, the following information is required as inputs: (i) The departure time
when a cruise ship leaves the home port and its arrival (return) time at the home port for
termination; (ii) The utility distribution at each port of call for the cruise passengers to
experience (for instance, the utility for the cruise passengers to spend time at a port city at 3:00
am is marginal); (iii) The relationship between bunker consumption and speed on each leg (a
leg is the voyage from one port to the next port). Then, based on these inputs, we make two
critical decisions: (i) The sequence of ports of call for the cruise ship to visit; and (ii) The
arrival and departure time at each port of call. The sequence displays the order list of ports, in
which ports of call must be visited by the cruise ship one by one. The arrival and departure
times confirm the staying time that the ship spends at each port of call and the voyage time on
each leg. The objective of the CISD problem is to maximize the total monetary value from the
utilities brought to cruise passengers at port cities minus the bunker cost of the cruise ship.
4.3.1 Departure time from the home port and the return time
We assume that the cruise ship departs from the home port and returns to the home port at
a pre-determined time. This assumption does not restrict the model but simply aims to simplify
the notation. Without loss of generality, we define that the cruise ship departs at Time 0 and
returns at Time 𝑇. Hence, there are a total of 𝑇 + 1 time points to complete the cruise, denoted
by set 𝕋 = {0, ⋯ , 𝑇}. Here, one time period could be set as one hour, as using one hour in the
schedule for cruise itineraries is precise enough (our model can also handle other time periods,
e.g., half an hour). Note that when we mention “at time 𝑡 ∈ 𝕋” we refer to the time at the end
of the 𝑡th time period (or equivalently, at the beginning of the (𝑡 + 1)th time period).
We make three comments on the utility shown in Figure 4.2. (i) The utility mentioned here
actually refers to the extra utility by spending time at port cities compared with spending time
at sea on cruise ships (as cruise passengers also have a lot of fun at sea). (ii) Different ports
have different utility distributions. For instance, a world-renowned city like Rome should have
high utilities; cities in which people tend to go to bed and wake up early, have different profiles
from those in which people stay late at nights in bars. (iii) It is possible that the ship stays at a
port when the utility is zero. For instance, in Figure 4.3, when two ports are very close, e.g,
two hours’ sailing, it is possible that the ship stays in Port 𝑗 overnight, when there is no utility,
and leaves the port at 8:00 am on the next day.
To facilitate the solution approach development, we define 𝑔𝑖 (𝑡) as the utility at Port 𝑖 in
Time period 𝑡 ∈ {1,2, ⋯ , 𝑇} . 𝑔𝑖 (𝑡) can be derived based on the daily distribution of the utility.
Considering the example in Figure 4.2, if the cruise ship leaves the home port at 4:00 pm
(which, by our definition, is Time 0), then 6:00-7:00 am of the next day corresponds to Time
period 15 and we thus have 𝑔𝑖 (15) = 2 as the utility for the hour from 6:00 am to 7:00 am is
two. Similarly, we have 𝑔𝑖 (15 + 24) = 2 and 𝑔𝑖 (15 + 48) = 2 in which 24 means one day
and 48 means two days.
To synthesize the total amount of utility that the cruise passengers will experience with the
bunker cost in the objective function, we denote 𝑝𝑢 as the monetary value for the cruise
The cruise ship can decelerate (or accelerate) if its sailing speed exceeds (or is lower than) the
optimal speed in order to save fuel.
Given 𝑑 and 𝜏, the average speed is 𝑑/𝜏. If 𝑑/𝜏 ≥ 𝑣 ∗ , the ship should sail at a constant
speed that is equal to 𝑑/𝜏 for saving fuel consumption. Otherwise, the ship should sail at its
optimal speed to the destination and then wait at the destination (Wang et al., 2013). Thus,
given 𝑑 and 𝜏, the minimum fuel consumption, denoted by function 𝐹(𝑑, 𝜏), is denoted as:
𝐹̃ (𝑑, 𝜏) , 𝑖𝑓 𝑑/𝜏 ≥ 𝑣 ∗
𝐹(𝑑, 𝜏) = { (4.4)
[𝑘 + 𝑘 ′ ∙ (𝑣 ∗ )𝑠 ] ∙ (𝑑⁄𝑣 ∗ ) , 𝑜𝑡ℎ𝑒𝑟𝑤𝑖𝑠𝑒
In the above model 𝐹0, Objective (4.5) maximizes the monetary value of the total utilities
minus the bunker costs. Constraints (4.6) and (4.7) guarantee that the cruise ship visits each
port of call exactly once. Constraints (4.8) and (4.9) ensure that each cruise starts at the home
port and goes back to the home port finally. Constraints (4.10) states that ports of call in the
same country must be visited without interruption for the sake of the visa restrictions.
Constraints (4.11) show the relationship between the departure time at one port and the arrival
time at the next visited port. Constraints (4.12) ensures that the cruise ship dwells in each port
for at least a certain period of time (i.e., 𝑚𝑖 ). Constraints (4.13) ensure that the real voyage
time between two ports should be more than the minimal voyage time, in which the cruise ship
sails at the maximum speed. Constraint (4.14) and Constraint (4.15) guarantee that the cruise
ship departs from the home port and returns to it at specified times. Constraints (4.16) indicate
that there is no partial connection between two ports. Constraints (4.17) and Constraints (4.18)
ensure that the cruise ship must arrive at one port and leave from it during its opening hours.
Figure 4.4: Statistic on trips for the Carnival Vista (Cruise Ship Schedule, 2016)
We now consider a special sequence 𝑠 that visits the ports of call according to their IDs, i.e.,
1, ⋯ 𝑗, 𝑗 + 1 ⋯ 𝑁, in which 1 and 𝑁 refer to the home ports and 𝑗 is the 𝑗 𝑡ℎ port visited on the
cruise itinerary. The cruise ship leaves Port 1 at Time 0 (i.e., 𝑏1 = 0), and returns to Port 𝑁 at
Time 𝑇 (i.e., 𝑎𝑁 = 𝑇). We need to determine the optimal arrival and departure times at each
Port 2, ⋯ 𝑗, 𝑗 + 1 ⋯ 𝑁 − 1 . The purpose of examining this special case is for notational
convenience; any sequence can be addressed using the same DP approach.
To apply DP, we firstly define 𝑈𝑗 (𝑡) as the maximum possible profit (to be determined)
from Time 𝑡 to 𝑇, if the cruise ship arrives at Port 𝑗 at Time 𝑡, and 𝑉𝑗 (𝑡) as the maximum
possible profit (to be determined) from Time 𝑡 to 𝑇, if the cruise ship leaves Port 𝑗 at Time 𝑡.
𝑝
Then, we define 𝜏𝑗𝑣 as the voyage time from the 𝑗 𝑡ℎ port to the (𝑗 + 1)𝑡ℎ port, and 𝜏𝑗 as the
𝑝
staying time at the 𝑗 𝑡ℎ port. The 𝜏𝑗𝑣 and 𝜏𝑗 are two decision variables in the DP algorithm.
Third, since the cruise ship has a maximum sailing speed 𝑣 𝑚𝑎𝑥 and the minimum time spent
at Port 𝑗 is 𝑚𝑗 , the earliest possible arrival time at Port 𝑗 can be computed by assuming the
cruise ship sails at its highest speed and spends the least time at previous ports of call:
𝑗−1 𝑑
∑𝑘=1 ⌈ 𝑘,𝑘+1
𝑚𝑎𝑥 ⌉ ,𝑗 = 2
𝑎𝑗𝑚𝑖𝑛 = { 𝑗−1 𝑑𝑣 𝑗−1
(4.20)
∑𝑘=1 ⌈ 𝑣𝑘,𝑘+1
𝑚𝑎𝑥 ⌉ + ∑𝑘=2 𝑚𝑘 , 𝑗 = {3,4, ⋯ 𝑁 − 1}
Define
𝑎𝑗𝑚𝑎𝑥 : = 𝑏𝑗𝑚𝑎𝑥 − 𝑚𝑗 , 𝑏𝑗𝑚𝑖𝑛 ≔ 𝑎𝑗𝑚𝑖𝑛 + 𝑚𝑗 ∀𝑗 ∈ 𝑃𝑐 (4.22)
range (i.e., ⌈𝑑𝑗,𝑗+1 /𝑣 𝑚𝑎𝑥 ⌉ ≤ 𝜏𝑗𝑣 ≤ 𝑇 − 𝑡). Here, we notice that the range may be infeasible
when time 𝑡 is large. For instance, when 𝑡 = 𝑇, the upper bound of the range is zero, which is
less than the lower bound of the range. To avoid these infeasible cases, we set 𝑉𝑗 (𝑡) to −∞
directly without calculation when the upper bound is less than the lower bound.
Eq. (4.25) is used to maximize the profit from Time 𝑡 to 𝑇 given that the cruise ship arrives
at Port 𝑗 at Time 𝑡. Similar to Eq. (4.24), the possible arrival time set at port 𝑗 is also to be
defined before calculation, which is the intersection of 𝒯𝑗 ′ and [𝑎𝑗𝑚𝑖𝑛 , 𝑎𝑗𝑚𝑎𝑥 ]. The possible
arrival time range (i.e., [𝑎𝑗𝑚𝑖𝑛 , 𝑎𝑗𝑚𝑎𝑥 ]) based on the speed of the ship and minimal staying hours
at the ports of call is also achieved by Eq. (4.20), (4.21) and (4.22). For all possible arrival
time spent at Port 𝑗 is 𝛿𝑗 . Since the minimum time spent at Port 𝑗 is 𝑚𝑗 , to have a feasible port
time solution, the time spent at Port 𝑗 must be between 𝑚𝑗 and 𝑇 − 𝑇̂ − ∑𝑖∈𝑃𝑐 \{𝑗} 𝑚𝑖 . The
model for obtaining 𝑈𝐵𝑝 (𝑇 − 𝑇̂) is:
𝑇−𝑇−∑ 𝑚 ̂
[𝑭𝟏] 𝑈𝐵𝑝 (𝑇 − 𝑇̂) = 𝑚𝑎𝑥 ∑𝑗∈𝑃𝑐 ∑𝛿 =𝑚 𝑖∈𝑃𝑐\{𝑗} 𝑖 𝐺𝑗 (𝛿𝑗 )𝑧𝑗𝛿𝑗 (4.31)
𝑗 𝑗
subject to:
𝑇−𝑇̂−∑𝑖∈𝑃𝑐\{𝑗} 𝑚𝑖
∑𝑗∈𝑃𝑐 ∑𝛿 𝛿𝑗 𝑧𝑗𝛿𝑗 = 𝑇 − 𝑇̂ (4.33)
𝑗 =𝑚𝑗
Based on above-mentioned definitions and formulations, the algorithm for improving the
enumeration can be designed (denoted as Algorithm 4.2) and is shown as:
Algorithm 4.2: Improving the enumeration:
Define the incumbent best profit: 𝑈𝐵𝑏𝑒𝑠𝑡 = −∞;
Calculate 𝑈𝐵𝑝 (𝑇 − 𝑇̂) for all possible values of 𝑇 − 𝑇̂;
Enumerate each sequence: sequence 𝑠 ∈ 𝑆;
For 𝑠 = 1 𝑡𝑜 𝑠 = (𝑁 − 2)!
If the sequence 𝑠 violates the visa restriction Then
Continue;
End If
If 𝑇 − ∑𝑁−1 𝑁−1
𝑖=1 𝑚𝑖 < ∑𝑗=1 ⌈𝑑𝑗(𝑠),(𝑗+1)(𝑠) /𝑣
𝑚𝑎𝑥
⌉ Then
Continue;
End If
Calculate the upper bound on the profit of the sequence 𝑠 by Eq. (4.35)
If 𝑈𝐵(𝑠) ≤ 𝑈𝐵𝑏𝑒𝑠𝑡 Then
Continue;
Else
Use DP (i.e., Algorithm 4.1) to find the maximum profit of the sequence 𝑠,
denoted by 𝑈𝐵𝑠𝐷𝑃 .
If 𝑈𝐵𝑏𝑒𝑠𝑡 < 𝑈𝐵𝑠𝐷𝑃 Then
Set 𝑈𝐵𝑏𝑒𝑠𝑡 = 𝑈𝐵𝑠𝐷𝑃 and record the sequence 𝑠;
End If
End If
End For
When the procedure for Algorithm 4.2 is finished, the optimal profit for the CISD problem is
achieved, which is 𝑈𝐵𝑏𝑒𝑠𝑡 , and the best sequence 𝑠 ∗ is recorded. The details of arrival times
and departure times can also be checked in the results of the 𝐷𝑃 for the sequence 𝑠 ∗ .
In order to validate the effectiveness of the proposed model 𝐹1 and the efficiency of the
developed solution method, we conduct numerical experiments by using a PC (Intel Core i5,
2.1G Hz; Memory, 4G). The solution method is implemented by Matlab R2013b. The integer
linear program 𝑀1 is solved by CPLEX12.1 with concert technology of C# (VS2008).
Table 4.2: Computational efficiency with and without using enumeration improving
Enumeration No enumeration
Instance Time ratio
improving method improving
# of 𝑇1
Cycle Instance
ports of 𝑍0 𝑇0 𝑍1 𝑇1
time id 𝑇0
call
7 3_7 2.115 4 2.115 5 1.25
3 8 3_8 3.196 5 3.196 8 1.60
9 3_9 3.690 8 3.690 20 2.50
9 5_9 2.440 10 2.440 14 1.40
5 10 5_10 3.374 13 3.374 25 1.92
11 5_11 4.096 15 4.096 43 2.87
11 7_11 2.646 22 2.646 62 2.82
7 12 7_12 3.573 53 3.573 198 3.74
13 7_13 4.404 85 4.404 552 6.49
13 9_13 3.803 236 3.803 2035 8.62
9 14 9_14 4.697 374 4.697 3961 10.59
15 9_15 5.413 642 5.413 9753 15.19
Average: 4.92
Note: (i) “ 𝑍0” and “ 𝑍1 ” columns list the optimal profits obtained by two methods with the unit of one
million US dollars. (ii) “ 𝑇0 ” and “ 𝑇1 ” columns show the CPU time (seconds) to solve the problem.
As can be seen in Table 4.2, both two methods obtain optimal results. This is verified by
the same optimal profits shown in the column “ 𝑍0 ” and the column “ 𝑍1 ”. However, the
method without using the enumeration improving is extremely time-consuming. The average
ratio of CPU time between this enumeration method and the method with the enumeration
improving is on average 4.92, which implies the proposed method saves approximately 80%
of the computational time on average. More importantly, the time ratio between two methods
4.5.4 Sensitivity analysis for a real case with a natural geographical sequence
As we have already tested the efficiency of the proposed method, we will use this method
to conduct sensitivity analysis in terms of three inputs, which are fuel price for the ship,
minimal staying hours of ports of call and opening hours of ports of call. Here, we take a
popular cruise line from the Royal Caribbean, “14 Night Singapore to Fremantle Cruise”, as
an example, and suppose that we are helping the Carnival Cruise Line to design a similar one-
way cruise. The nine ports involved in this cruise are shown in Figure 4.5, and it can easily be
seen that this itinerary has a naturally geographical sequence, under which we may design a
good visiting sequence by direct observation.
Figure 4.5: City ports in “14 Night Singapore to Fremantle Cruise” (Google Map, 2018a)
For this cruise, assume that we deploy the Carnival Vista, which is the biggest cruise ship
of the Carnival Cruise Line. It has 133,500 deadweight tons. Same as the “14 Night Singapore
to Fremantle Cruise” operated by the Royal Caribbean International, in this cruise, the cruise
ship departs from the port of Singapore at 5:00 pm on Day 1 and arrives at the port of
Fermantle, Australia at 7:00 am on Day 15. Singapore and Fermantle are the home ports. The
Figure 4.6: City ports around the Caribbean Sea (Google Map, 2018b)
As the bunker cost is one of the major concerns for cruise companies, we would like to
further conduct analysis on the fuel price to see how important an optimal port sequence is
needed when the fuel price fluctuates. In order to conduct such analysis, we select 16 real
cruise services operated by Royal Caribbean International in the Caribbean Sea area (Royal
Caribbean International, 2016). Before conducting the analysis, we first redesign the itinerary
schedules for the 16 real cruise service. The comparisons in terms of port sequence and port
staying hours between the 16 designed itineraries and the 16 actual itineraries are listed in
Figure 4.7. Here, the y-axis shows the total port staying hours for each itinerary, and the x-
axis indicates the index of each cruise service. Symbol “Y” (Symbol “N”) above pairs of bars
indicates that the designed itinerary and the actual itinerary have the same port sequence
(different port sequences) for the corresponding cruise service.
For further analysis of the fuel price, based on the above-mentioned 16 cruise services, the
optimal sequence under the baseline setting of the fuel price is obtained by using the proposed
method at first. Then, the profits under other fuel price settings of this optimal sequence (i.e.,
Baseline itinerary) are calculated, which are defined as Baseline itinerary profit. The optimal
sequences and the optimal profits under other fuel price settings are also re-optimized by using
the proposed method, which is compared with the Baseline itinerary and Baseline itinerary
profit. The results of the comparisons are reported in Table 4.7.
From Table 4.7, we can observe that when the fuel price fluctuates, the optimal sequences
for the itineraries are highly likely to change. Especially for large degrees of fluctuation, such
as the 0.4-differentiation and the 1.6-differentiation on the fuel price, the majority of optimal
sequences in the 16 cruise services are different from the optimal sequences obtained under
the baseline setting of the fuel price. However, when the fuel price fluctuates, if the optimal
sequence obtained in the baseline setting is fixed for the itinerary, there is significant profit
loss based on the average deviation in the table. In extreme cases, for instance, with 0.4-
differentiation on the fuel price, the average profit loss is 4.52%. Therefore, it is critical for
the cruise companies to re-optimize the visiting sequence when the fuel price fluctuates
significantly in order to achieve a higher profit.
Number of instances
having the same Average
Differentiation on the Average Baseline Average
optimal port optimal
fuel price itinerary profit deviation
sequence as the profit
baseline setting
0.4 3 5.022 4.805 4.52%
0.6 6 4.894 4.772 2.56%
0.8 8 4.779 4.721 1.23%
Baseline setting 16 4.663 4.663 0.00%
1.2 10 4.538 4.490 1.07%
1.4 6 4.424 4.317 2.48%
1.6 5 4.322 4.159 3.92%
Note: (i) “Average optimal profit” column list the average profit among the 16 cruise services under
each fuel price setting. (ii) “Average optimality gap” column shows the gap between the average
optimal profit and the average Baseline itinerary profit.
Based on the 16 selected cruise services from Royal Caribbean International, we also
conduct some sensitivity analysis on the utility distributions of the ports of call. For each port
of call in one instance, its utility distribution among hours in one day may have different values
of mean and standard deviation (SD). The SD indicates the degree of variation of the utility
among 24 hours. Here, we change the SD of the utility distribution for each port of call by
The
Minimum fuel cost The shortest path
Instance group proposed
heuristic heuristic
method
# of
Average Average Average Average Average
ports Cycle time
profit profit deviation profit deviation
of call
9 2.352 2.225 5.40% 2.325 1.15%
5 10 3.249 3.044 6.31% 3.184 2.00%
11 3.894 3.551 8.81% 3.764 3.34%
11 2.501 2.352 5.96% 2.457 1.76%
7 12 3.297 3.084 6.46% 3.194 3.12%
13 4.235 3.814 9.94% 4.074 3.80%
13 3.712 3.478 6.30% 3.604 2.91%
9 14 4.494 4.121 8.30% 4.284 4.67%
15 5.108 4.578 10.38% 4.818 5.68%
Average: 7.54% 3.16%
Note: (i) In each instance group, there are 10 instances, which are randomly generated based on major
ports or islands in the Caribbean Sea area. (ii) The average profit of 10 instances in one group is
calculated for the proposed method and the two heuristics. (iii) “Average deviation” shows the
deviations between the average profit by the proposed method and the average profits by the heuristics.
Attributes Levels
4.7 CONCLUSION
This paper addresses the cruise itinerary schedule design problem that determines the
optimal sequence of a given set of ports of call and the arrival and departure times at each port
to maximize the monetary value of the utility minus the fuel cost. In view of the practical
observations that most cruise itineraries do not have many ports of call, we first enumerate all
sequences of ports of call and then optimize the arrival and departure times at each port of call
by developing a dynamic programming approach. To improve the computational efficiency,
we propose effective bounds on the profit of each sequence of ports of call, eliminating non-
optimal sequences without invoking the dynamic programming algorithm. The computational
experiments show that, first, the proposed bounds on the profit of each sequence of ports of
call can considerably improve the computational efficiency; second, the total profit of the
cruise itinerary is sensitive to the fuel price and hence, it is acceptable to use the shortest
voyage distance method to design the schedule when the ports of call have a naturally
geographical distance; in contrast, determining the sequence of ports of call solely by
minimizing the overall voyage distance frequently leads to a significant reduction in the total
profit when the ports do not have a naturally geographical sequence.
Given that cruise itineraries have fixed sequences of ports of call and fixed schedules,
optimization-based itineraries planning tools should be able to increase the profit or save the
cost for cruise shipping companies and improve the service quality for cruise passengers.
Compared with other areas of transportation such as a truck, rail, and air (we note that
4.8 REFERENCES
Applegate, D.L., Bixby, R.E., Chvatal, V., Cook, W.J., 2011. The Traveling Salesman
Problem: A Computational Study. Princeton University Press, New Jersey.
Cormen, T.H., Leiserson, C.E., Rivest, R.L., Stein, C., 2009. Introduction to Algorithms, Third
Edition. The MIT Press, London, England.
Cordeau, J. F., Laporte, G., Mercier, A., 2001. A unified tabu search heuristic for vehicle
routing problems with time windows, Journal of the Operational Research Society 52(8),
928–936.
Carnival Cruise Line, 2016. 7 day Western Caribbean from Miami, FL.
http://www.carnival.com/itinerary/7-day-western-caribbean-cruise/miami/glory/7-
days/wek. Accessed 4 April 2016.
Ding, M., Park, Y.H., Bradlow, E.T., 2009. Barter markets for conjoint analysis, Management
Science 55(6), 1003–1017.
Du, Y., Chen, Q., Quan, X., Long, L., Fung, R.Y.K., 2011. Berth allocation considering fuel
consumption and vessel emissions, Transportation Research Part E 47, 1021–1037.
Du, Y., Chen, Q., Lam, J. S. L., Xu, Y., Cao, J. X., 2015. Modeling the impacts of tides and
the virtual arrival policy in berth allocation, Transportation Science 49(4), 939-956.
Feillet, D., Dejax, P., Gendreau, M., 2005. Traveling salesman problems with profits,
Transportation Science 39(2), 188–205.
Fransoo, J.C., Lee, C.Y., 2013. The critical role of ocean container transport in global supply
chain performance, Production and Operations Management 22(2), 253–268.
Gilbride, T.J., Lenk, P.J., Brazell, J.D., 2008. Market share constraints and the loss function
in choice-based conjoint analysis, Marketing Science 27(6), 995–1011.
Green, P.E., Krieger, A.M., 1996. Individualized hybrid models for conjoint analysis,
Management Science 42(6), 850–867.
Hongmin, L., Woonghee, T.H., 2011. Pricing multiple products with the multinomial logit and
nested logit models: concavity and implications, Manufacturing & Service Operations
Management 13(4), 549–563.
Kimes, S.E., 1989. Yield management: a tool for capacity-considered service firms, Journal
of Operations Management 8, 348–363.
Koppelman, F.S., Goldren, G.M., Parker, R.A., 2008. Schedule delay impacts on air-travel
itinerary demand, Transportation Research Part B: Methodological 42, 263–273.
Kuhfeld, W., 2010. Marketing research methods in SAS. Cary: SAS Institute Inc..
Li, B., 2011. The multinomial logit model revisited: A semi-parametric approach in discrete
choice analysis, Transportation Research Part B: Methodological 45(3), 461–473.
Maddah, B., Moussawi-Haidar, L., El-Taha, M., Rida, H., 2010. Dynamic cruise ship revenue
management, European Journal of Operational Research 207(1), 445–455.
Meng, Q., Wang, S., Andersson, H., Thun, K., 2014. Containership routing and scheduling in
liner shipping: overview and future research directions, Transportation Science 48 (2),
265–280.
Ng, M.W., 2014. Distribution-free vessel deployment for liner shipping, European Journal of
Operational Research 238(3), 858–862.
Ng, M.W., 2015. Container vessel fleet deployment for liner shipping with stochastic
dependencies in shipping demand, Transportation Research Part B 74, 79–87.
Paat, R., Huseyin, T., 2012. Robust assortment optimization in revenue management under the
multinomial logit choice model, Operation Research 60(4), 865–882.
Pang, K.W., Liu, J., 2014. An integrated model for ship routing with transshipment and berth
allocation, IIE Transactions 46(12), 1357–1370.
Qi, X.T., Song, D.P., 2012. Minimizing fuel emissions by optimizing vessel schedules in liner
shipping with uncertain port times, Transportation Research Part E 48(4), 863–880.
Rodrigue, J., Notteboom, T., 2013. The geography of cruises: Itineraries, not destinations,
Applied Geography 38, 31–42.
Royal Caribbean International, 2016. Cruise itineraries sailing to the Caribbean area.
http://www.royalcaribbean.com/cruises?destinationRegionCode_CARIB=true.
Accessed 22 Aug 2016.
Song, D.P. and Dong, J.X. 2011. Effectiveness of an empty container repositioning policy with
flexible destination ports, Transport Policy 18 (1), 92–101.
Song, D.P. and Dong, J.X. 2012. Cargo routing and empty container repositioning in multiple
shipping service routes, Transportation Research Part B 46(10), 1556–1575.
Song, D.P. and Dong, J.X. 2013. Long-haul liner service route design with ship deployment
and empty container repositioning, Transportation Research Part B 55, 188–211.
Soriani, S., Bertazzon, S., Cesare, F.D., Rech, G., 2009. Cruising in the Mediterranean:
structural aspects and evolutionary trends, Maritime Policy & Management 36(3), 235–
251.
Statista, 2015. Cruise industry—Statista Dossier: Facts and statistics on the cruise industry.
http://www.statista.com/study/11547/cruise-line-industry-statista-dossier/. Accessed 28
Aug 2015.
Talluri, K., van Ryzin, J., 2004. The Theory and Practice of Revenue Management, Kluwer
Academic, Boston.
Toubia, O., Hauser, J., Garcia, R., 2007. Probabilistic polyhedral methods for adaptive choice-
based conjoint analysis: theory and application, Marketing Science 26(5), 596–610.
Toth, P., Vigo, D., 2001. The Vehicle Routing Problem. Society for Industrial and Applied
Mathematics, Philadelphia.
Vermeulen, B., Goos, P., & Vandebroek, M., 2008. Models and optimal designs for conjoint
choice experiments including a no-choice option. International Journal of Research in
Marketing 25(2), 94–103.
Veronneau, S., Roy, J., 2009. Global service supply chains: An empirical study of current
practices and challenges of a cruise line corporation, Tourism Management 30, 128–139.
Veronneau, S., Roy, J., Beaulieu, M., 2015. Cruise ship suppliers: A field study of the supplier
relationship characteristics in a service supply chain, Tourism Management Perspectives
16, 76–84.
Wang, K., Wang, S., Zhen, L., Qu, X., 2016. Cruise shipping review: Operations planning and
research planning. Maritime Business Review 1(2), 133-148.
Wang, S., Meng, Q., 2012b. Liner ship route schedule design with sea contingency time and
port time uncertainty, Transportation Research Part B 46(5), 615–633.
Wang. S., Meng. Q., Liu. Z., 2013. Bunker consumption optimization methods in shipping: A
critical review and extensions, Transportation Research Part E: Logistics and
Transportation Review 53, 49–62.
Wang. S., Wang. H., Meng. Q., 2015. Itinerary provision and pricing in container liner
shipping revenue management, Transportation Research Part E: Logistics and
Transportation Review 77, 135–146.
Xie, H., Kerstetter, D.L., Mattila, A.S., 2012. The attributes of a cruise ship that influence the
decision making of cruisers and potential cruisers, International Journal of Hospitality
Management 31, 152–159.
This chapter addresses a decision problem on planning cruise services for a cruise ship
so as to maximize the total profit during a planning horizon. The service is a sequence of ports
(harbor cities) the cruise ship visits. In this decision problem, the constraint about the
availability of berths at each port is taken into account. In reality, if a cruise service is executed
by the ship repeatedly for several times, the profit earned by the cruise service in each time
decreases gradually. This effect of decreasing marginal profit is also considered in this study.
We propose a nonlinear integer programming model to cater to the concavity of the function
for the profit of operating a cruise service repeatedly. To solve the nonlinear model, two
linearization methods are developed, one of which takes advantage of the concavity for a
tailored linearization. Some properties of the problem are also investigated and proved by
using dynamic programming (DP) and two commonly used heuristics. In particular, we prove
that if there is only one candidate cruise service, a greedy algorithm can derive the optimal
solution. Numerical experiments are conducted to validate the effectiveness of the proposed
models and the efficiency of the proposed linearization methods. In case some parameters
needed by the model are estimated inexactly, the proposed decision model demonstrates its
robustness and can still obtain a near-optimal plan, which is verified by experiments based on
extensive real cases.
5.1 INTRODUCTION
In the cruising industry, service planning is often independent among different cruise ships,
and planning problems of different cruise ships can be solved individually. According to
Rodrigue and Notteboom (2013), the cruise ship deployment focuses on a specific cruise ship
rather than a fleet of cruise ships. Cruise ships are often unique, even if the cruise ships have
the same capacity (in passengers), different cruise ships have significantly different onboard
activities, which are part of cruising experience. The same cruise route traversed by two cruise
ships constructs two different cruise services, as the onboard activities are different. Our
research is applicable to this situation. If two ships are very similar and serve the same region
and visit common ports of call, then the competition between the services provided by the
ships have to be accounted for.
This study assumes that when a cruise ship is repositioned to a new region, a home port and
a set of candidate cruise services are chosen in advance and addresses the Cruise Service
Planning (CSP) problem. This problem aims to determine how to plan cruise services for the
cruise ship to operate in the region for a period of time. In other words, over the period of time,
The cruise shipping related studies could belong to the area of tourism research as cruise
ships provide cruise passengers with tourism service. Meanwhile, it can be also sorted into the
area of maritime research as the cruise services are akin to container liner services. However,
the past research on cruise shipping is limited, the reason of which may include: (i) the
worldwide cruise ship tourism just accounts for about 2% of the world tourism market revenue,
thus the tourism-related researchers have not paid much attention to the cruise shipping related
studies (Gui and Russo, 2011); (ii) the maritime logistic related researchers mainly focus on
freight transportation (e.g., Bell et al., 2011; Meng and Wang, 2012; Song et al., 2015).
The majority of past research on the cruise shipping analyzed the cruising industry as a
tourism service supply chain. Gui and Russo (2011) constructed an analytic framework
connecting the global structure of cruise value chains to the regional land-based cruise
services. The demand side and the supply side of the cruise shipping at the worldwide level
were analyzed by Soriani et al. (2009). They also investigated the main characteristics of the
cruising in the Mediterranean region and examined the main cruising ports in the region. A
field study of a large Florida-based global cruise company’s practices in re-supplying ships
globally was conducted by Veronneau and Roy (2009), which makes them amongst the first
to take a comprehensive study of a specific service supply chain. Veronneau et al. (2015)
investigated the relationships between a major cruise line corporation and its suppliers by a
field study. Rodrigue and Notteboom (2013) focused on capacity deployment and itineraries
in two important markets: the Caribbean and Mediterranean. They found that these two market
areas interact with each other due to seasonal variations in demand.
Although those researchers devoted significant efforts into cruising shipping research, their
works are mainly descriptive and belong to empirical studies. The majority of existing related
works did not provide the cruise industry with quantitative analysis on cruise shipping, which
could be critical for some detailed problems. Maddah et al. (2010) conducted a quantitative
study on cruise shipping, in which a discrete-time dynamic capacity control model was built
to improve the profit of cruise ships. Their model could give cruise ship managers some
suggestions about which requests from customers should be accepted based on the remaining
cabin and lifeboat capacities and the type of requests. The research is adaptable for cruise
companies to improve profit, and it focuses on an operational level problem. As the cruising
industry has developed dramatically, more research efforts should be made on the problems in
strategic level or tactical level.
In this section, we provide a brief description of the CSP problem for a cruise ship
considering the decreasing marginal profit and the berth availability and formulate it as integer
programming models.
When planning services for the cruise ship in the planning horizon, we assume that all the
services must be finished before 𝑇. To ensure this assumption, we initially set input data 𝛿𝑟𝑡
in the following way: for each candidate service 𝑟 with rotation time 𝑠𝑟 , 𝛿𝑟𝑡 ≔ 0 if 𝑡 is bigger
than 𝑇 − 𝑠𝑟 + 1 . This setting can guarantee that all the services are finished within the
planning horizon.
For the interest of simplicity, firstly, this study assumes that the marginal profit of operating
Service 𝑟 is constant when the service is repeated several times. Under this assumption, an
integer programming model is formulated, denoted as model 𝐹𝑀1. Then, in order to be close
to the real situation, we change the profit setting to that the marginal profit of operating a
service decreases when the service is repeated several times. Based on this, a nonlinear integer
programming model is formulated, denoted as model 𝐹𝑀2.
In the model with decreasing marginal profit, we assume that when a service is repeated for
several times, the marginal profit of the service will decrease gradually. This consideration
makes the problem close to the reality as few potential passengers would choose a cruise
service when the cruise service has been repeated many times. To show the decreasing pattern
of the marginal profit in the model, we define a strictly increasing concave function 𝐺𝑟 (𝑥𝑟 ) as
Here, it is worthwhile to mention that there is a better way to derive the upper bound of the
number of repeats for Service 𝑟. We propose a greedy algorithm to find the maximum number
of repeats for Service 𝑟 in the planning horizon (note that the greedy algorithm is applied to
derive optimal 𝑥𝑟 in the numerical experiments in Section 5.5), which is a tighter upper bound
compared with Eq. (5.9). The details of the greedy algorithm and why the greedy algorithm
can find the optimal number of repeats for Service 𝑟 will be given in the proof of Proposition
5 in Section 5.4.3.
The first linear model:
The first linear model is defined as 𝐹𝑀2′ , which does not take advantage of the problem
structure for the linearization. In this linear model, 𝑔𝑟𝑥 ≔ 𝐺𝑟 (𝑥) − 𝐺𝑟 (𝑥 − 1) is used to define
the marginal profit of the 𝑥 𝑡ℎ repeat of Service 𝑟. Meanwhile, the binary decision variable
(i.e., 𝑧𝑟𝑡 ,) in 𝐹𝑀2 is changed to another binary decision variable denoted by 𝑧𝑟𝑥𝑡 , ∀𝑟 ∈ ℝ, 𝑥 ∈
{1, 2 ⋯ 𝑥𝑟 }, 𝑡 ∈ 𝕋. Here, 𝑧𝑟𝑥𝑡 equals one if and only if Service 𝑟 is operated for the 𝑥 𝑡ℎ time
starting from 0:00 am (Day 𝑡 ). Based on the above redefinition of some variables, the
formulation of 𝐹𝑀2′ is:
In this model, Objective (5.16) maximizes the sum of the total profits for operating services
by repeating different times. Constraints (5.17) are used to calculate the total profit of each
service in a linear manner. Here, notice that 𝐺𝑟 (0) ≔ 0, ∀𝑟 ∈ ℝ.
Both linear models provide the linearization for the nonlinear model (i.e., 𝐹𝑀2), which will
be further compared in Section 5.5 for the computation efficiency to solve the problem.
In this section, for the complexity analysis, two dynamic programming (DP) based pseudo-
polynomial algorithms are developed for the model 𝐹𝑀1 and the model 𝐹𝑀2 respectively.
Some properties on commonly used heuristic methods for the CSP problem are also
investigated and proved.
∀𝑡 ∈ 𝕋. (5.20)
To calculate 𝑈(𝑡) at each stage, we firstly enumerate all the services and select feasible
services by considering the availability of berths at the ports visited (𝛿𝑟𝑡 ). Then, we start the
feasible services one by one to derive the profits for the feasible services, and also calculate
into the planning horizon sequentially until it is impossible to place more “detailed services”.
Based on the above two commonly used heuristics, we can design two myopic heuristic
rules to solve the general case of the model 𝐹𝑀2 when considering the berth availability, in
which the decision is made on a daily basis from Day 1 to Day 𝑇.
Myopic Rule_1: For a specified Day 𝑡, we determine all the cruise services that can be
operated with considering the berth availability (i.e., ∀𝑟, 𝛿𝑟𝑡 = 1). Based on those cruise
services, we select the optimal cruise service 𝑟 ∗ with the maximal daily operating profit (i.e.,
𝑔𝑟∗ /𝑠𝑟∗ ). Then, the time is updated to Day 𝑡 + 𝑠𝑟∗ for the next selection. If no cruise service
can be operated, the time is updated to Day 𝑡 + 1 for the next selection.
Myopic Rule_2: For a specified Day 𝑡, we determine all the cruise services that can be
operated with considering the berth availability (i.e., ∀𝑟, 𝛿𝑟𝑡 = 1) at first. Among those cruise
services, we select the optimal cruise service 𝑟 ∗ with the maximal operating profit (i.e., 𝑔𝑟∗ )
to operate. Then, the time is updated to Day 𝑡 + 𝑠𝑟∗ for the next selection. If no cruise service
can be operated, the time is updated to Day 𝑡 + 1 for the next selection. Notice that if two
cruise services have the same daily operating profit or the same operating profit in the rules,
the priority will be given to the cruise service with shorter rotation time as it would occupy
few days in the planning horizon.
Here, Myopic Rule_1 (Myopic Rule_2) is designed by the unit-profit-first heuristic (the
operating-profit-first heuristic) as it selects the optimal cruise service 𝑟 ∗ with the maximal
For the CSP problem considering berth availability, the two commonly used heuristic rules
do not give the priority on the berth availability, which leads to tremendous profit loss. Thus,
in cruise shipping, the operation managers should keep well informed about the berth
availability from cruise terminals. Based on the information, the managers should make
schedules on the overall picture for a whole period rather than from day to day. In reality, it is
impossible to guarantee the sufficient berth availability in all cruise terminals, but we do
encourage that the cruise lines own some cruise terminals such that they have more flexibility
on berths to operate their cruise services.
Proposition 5.5: When there is only one candidate cruise service 𝑟 ′ , the solution obtained
by Myopic Rule_1 or Myopic Rule_2 is the optimal solution of the model 𝐹𝑀2.
Proof: When there is only one candidate cruise service 𝑟 ′ , the rotation time is constant, then
we can transfer the two myopic heuristic rules to a greedy algorithm for the problem, in which
the decision is made on a daily basis from Day 1 to Day 𝑇: for a specified Day 𝑡, if the cruise
service 𝑟 ′ can be operated with considering the berth availability, this cruise service is settled
for the operation. Then, the time is updated to Day 𝑡 + 𝑠𝑟′ for the next decision. Otherwise,
the time is updated to Day 𝑡 + 1 for the next decision.
As there is only and one candidate cruise service (the cruise service 𝑟 ′ ), the objective of the
model 𝐹𝑀2 aims to maximize the total profit earned by the one cruise service (maximize 𝑍 =
𝐺𝑟′ (𝑥𝑟′ )). Based on the concavity of the function 𝐺𝑟′ (𝑥𝑟′ ), as it is shown in Figure 5.1, the
In this section, in order to validate the effectiveness of the proposed models and efficiency
of solving the models, we conduct extensive numerical experiments by using a PC (Intel Core
i5, 2.3G Hz; Memory, 8G). The integer programs 𝐹𝑀2′ and 𝐹𝑀2′′ are solved by CPLEX12.5
with concert technology of C# (2012).
In Section 5.3.3, we have proposed two linear models 𝐹𝑀2′ and 𝐹𝑀2′′ for the nonlinear
model 𝐹𝑀2. Here, we test which linear model has a higher efficiency to derive solutions for
the problem. As we have verified that the model 𝐹𝑀2′ without Constraints (5.14) can be
solved faster, the comparison is conducted between this case of the model 𝐹𝑀2′ and the model
𝐹𝑀2′′ . Table 5.3 illustrates the comparison between the model 𝐹𝑀2′ and 𝐹𝑀2′′ . Both linear
models are valid for the linearization of the nonlinear model as the optimal solutions are
obtained in all instance groups. Whereas, the model 𝐹𝑀2′′ can be solved much faster than the
model 𝐹𝑀2′ by CPLEX. The ratio of CPU times between two models is 0.23 on average,
which shows the advantage of using the concavity of 𝐺𝑟 (𝑥𝑟 ) for the linearization. In a
technical perspective of CPLEX, the model 𝐹𝑀2′ spends too much CPU time on pre-solving
the problem, and the nodes explored in CPLEX for two models are more or less the same,
shown by “B&B nodes”.
Without the
Instance With the constraints Comparison
constraints
# of candidate Instance
𝑍𝑤 (US$) 𝑇𝑤 (s) 𝑍𝑜 (US$) 𝑇𝑜 (s) 𝑇𝑜 / 𝑇𝑤
service ID
2_20_1 1.156 35 1.156 16 0.46
2_20_2 1.093 47 1.093 18 0.38
20 2_20_3 1.136 29 1.136 15 0.52
2_20_4 1.084 33 1.084 22 0.67
2_20_5 1.217 24 1.217 11 0.46
2_40_1 1.312 218 1.312 29 0.13
2_40_2 1.311 56 1.311 20 0.36
40 2_40_3 1.313 90 1.313 38 0.42
2_40_4 1.349 47 1.349 21 0.45
2_40_5 1.304 124 1.304 43 0.35
2_80_1 1.392 236 1.392 39 0.17
2_80_2 1.378 504 1.378 57 0.11
80 2_80_3 1.446 378 1.446 38 0.10
2_80_4 1.405 870 1.405 67 0.08
2_80_5 1.407 573 1.407 88 0.15
Average: 0.32
Note: (i) “# of candidate service” column denotes the total number of candidate services. (ii) “ 𝑍𝑤 ”
and “ 𝑍𝑜 ” columns list the optimal profits under two cases with the unit of ten million US dollars. (iii)
“ 𝑇𝑤 ” and “ 𝑇0 ” columns show the CPU time (seconds) to solve the problem.
Table 5.3: Comparison between the two linear models
Compariso
Instance The first model The second model LP-relaxation
n
# of B&B B&B
Instanc 𝑍𝑓 𝑇𝑓 𝑇𝑠
candidat node 𝑍𝑠 ($) node 𝑇𝑠 /𝑇𝑓 𝑍𝑙 Gap
e id ($) (s) (s)
e service s s
3_20_1 1.182 19 1 1.182 4 1 0.21 1.185 0.23%
3_20_2 1.154 20 1 1.154 5 1 0.25 1.158 0.37%
20 3_20_3 1.127 13 1 1.127 5 1 0.38 1.131 0.36%
3_20_4 1.147 12 1 1.147 2 1 0.17 1.150 0.26%
3_20_5 1.167 18 1 1.167 5 1 0.28 1.173 0.55%
3_40_1 1.321 28 162 1.321 6 141 0.21 1.323 0.19%
3_40_2 1.308 21 41 1.308 7 60 0.33 1.312 0.31%
40 3_40_3 1.294 33 453 1.294 10 407 0.30 1.299 0.37%
3_40_4 1.313 28 79 1.313 7 83 0.25 1.317 0.32%
3_40_5 1.301 54 179 1.301 8 303 0.15 1.306 0.39%
3_80_1 1.398 134 593 1.398 17 537 0.13 1.400 0.13%
3_80_2 1.387 64 83 1.387 15 100 0.23 1.389 0.13%
80 3_80_3 1.412 85 154 1.412 14 294 0.16 1.414 0.17%
3_80_4 1.405 101 317 1.405 20 357 0.20 1.407 0.12%
3_80_5 1.394 88 177 1.394 16 326 0.18 1.396 0.16%
Average: 0.23 1.291 0.27%
Table 5.4: Comparison between the model 𝑭𝑴𝟐′′ and two myopic rules
For decision makers of a cruise ship, a challenge of implementing our model is to estimate
the marginal profit of operating a cruise service accurately. Usually, as the money spent by
35.0
30.0
25.0
20.0
15.0
10.0
5.0
0.0
C1 C2 C3 C4 C5 C6 C7 C8 C9 C10 C11 C12 C13
Figure 5.5: Upper bound 𝒙𝒓 obtained by the greedy algorithm and Eq. (5.9)
Figure 5.5 reports 𝑥𝑟 , upper bound of the number of repeats for Service 𝑟 in one planning
horizon. Eq. (5.9) exhibits a way to approximate 𝑥𝑟 . In Proposition 5, we proved that the
greedy algorithm can obtain the optimal 𝑥𝑟 . Here, we report 𝑥𝑟 obtained by the greedy
algorithm under the five berth availability scenarios and by Eq. (5.9) are given in Figure 5.5.
Note that the 𝑥𝑟 obtained by Eq. (5.9) is constant under different berth availability (given by
𝑇
Bar Eq. (5.9)), as it is derived by ⌊𝑠 ⌋, and the 𝑥𝑟 obtained by the greedy algorithm is different
𝑟
under different berth availability scenarios (given by five bars from BA1 to BA5). Therefore,
each candidate cruise service (indexed by C1 to C13) in Figure 5.5 contains six bars. As can
be seen, the upper bound 𝑥𝑟 obtained by Eq. (5.9) is much worse than that of the greedy
algorithm, especially when the berth availability is low (BA1). For example, for Cruise service
13, the 𝑥𝑟 obtained by Eq. (5.9) is more than ten times as large as that of the greedy algorithm
under the berth availability scenario BA1. Thus, to implement our proposed model, the greedy
algorithm should be applied to approximate 𝑥𝑟 .
Figure 5.6 illustrates the average number of repeats of each cruise service under each berth
availability scenario. In general, cruise services 4, 5 and 6 outperform other cruise services
with a higher average number of repeats. This is due to the fact that those cruise services have
comparatively higher marginal profits and shorter rotation times (cf. Table 5.5). Cruise service
11 is the least selected cruise service to be operated, especially when the berth availability is
low. The cruise service with the longest rotation time (cruise service 1) also performs badly in
BA1, but the performance improves when the berth availability increases and around 1.5
Figure 5.6: Number of repeats of cruise services under different berth availability scenarios
Table 5.6 shows the effects of different berth availability scenarios on major outputs of the
model. When the berth availability increases, operation days and total profit rise
simultaneously. However, the increase of operation days or total profit does not keep the pace
with the increase of the berth availability. For instance, from BA3 to BA5, the berth
availability grows by 20%, but the total profit increases by 9.8%. Average profit per day shares
the same trend as the total profit because the number of total days is constant. By comparison,
the average profit per operation day keeps nearly unchanged when the berth availability
fluctuates.
Table 5.6: Outputs of the model under different berth availability scenarios
Average Average
Total Operation Deviation Total Deviation
ID profit per profit per
days days 1 profit 2
day operation day
BA1 180 127.8 −18.1% 8,568,362 −19.1% 47,602 67,043
BA2 180 146.4 −6.2% 9,839,851 −7.1% 54,666 67,186
BA3 180 156.0 0.0% 10,588,128 0.0% 58,823 67,901
BA4 180 168.2 7.8% 11,363,633 7.3% 63,131 67,577
BA5 180 172.8 10.8% 11,622,052 9.8% 64,567 67,266
Note: (i) “Total days” shows the length of one planning horizon. (ii) “Operation days” indicates the
number of days that the cruise ship operates cruise services, i.e., the cruise ship is traveling. (iii)
“Deviation 1” lists the deviation of the operation days between the corresponding berth availability
and BA3. (iv) “Deviation 2” lists the deviation of the total profit between the corresponding berth
availability and BA3.
5.7 REFERENCES
Arrival Time, 2016. Arrival Time—Wusong Kou Cruise Terminal.
http://www.wskcruise.com/article/48.html. Accessed 22 February 2016.
Arthur, H., Ronald, T.W., 2000. Stochastic prediction in multinomial logit models.
Management Science 46(8), 1137–1144.
Bell, M.G.H., Liu, X., Rioult, J., Angeloudis, P., 2013. A cost-based maritime container
assignment model. Transportation Research Part B 58, 58–70.
Bell, M.G.H., Liu, X., Angeloudis, P., Fonzone, A., Hosseinloo, S.H., 2011. A frequency-
based maritime container assignment model. Transportation Research Part B 45 (8),
1152–1161.
Bretthauer, K., Shetty, B., 2002. The non-linear knapsack problem - algorithms and
applications. European Journal of Operational Research 138, 459–472.
Cruise Industry, 2015. Cruise industry—Statista Dossier: Facts and statistics on the cruise
industry. http://www.statista.com/study/11547/cruise-line-industry-statista-dossier/.
Accessed 26 December 2015.
Cruise route: Quantum of the Seas, 2016. Royal Caribbean itinerary 2016 cruises.
http://www.wskcruise.com/tours /. Accessed 22 February 2016.
Esteve-Perez, J., Garcia-Sanchez, A., 2014. Cruise market: Stakeholders and the role of ports
and tourist hinterland. Maritime Economics and Logistics 17(3), 1–18.
Fransoo, J.C., Lee, C.Y., 2013. The critical role of ocean container transport in global supply
chain performance. Production and Operations Management 22(2), 253–268.
Giallombardo, G., Moccia L., Salani, M., Vacca, I., 2010. Modeling and solving the tactical
berth allocation problem. Transportation Research Part B 44(2), 232–245.
Hongmin, L., Woonghee, T.H., 2011. Pricing multiple products with the multinomial logit and
nested logit models: concavity and implications. Manufacturing & Service Operations
Management 13(4), 549 –563.
Iris, C., Pacino, D., Ropke, S., and Larsen, A., 2015. Integrated berth allocation and quay crane
assignment problem: Set partitioning models and computational results. Transportation
Research Part E 81, 75–97.
Kameshwaran, S., Narahari, Y., 2009. Nonconvex Piecewise Linear Knapsack Problems.
European Journal of Operational Research 192, 56–68.
Li, B., 2011. The multinomial logit model revisited: A semi-parametric approach in discrete
choice analysis. Transportation Research Part B 45(3), 461 –473.
Maddah, B., Moussawi-Haidar, L., El-Taha, M., Rida, H., 2010. Dynamic cruise ship revenue
management. European Journal of Operational Research 207(1), 445–455.
Meisel, F., Bierwirth. C., 2009. Heuristics for the integration of crane productivity in the berth
allocation problem. Transportation Research Part E 45(1), 196–209.
Meng, Q., Wang, S., 2012. Liner ship fleet deployment with week-dependent container
shipment demand. European Journal of Operational Research 222(2), 241–252.
Meng, Q., Wang, S., Andersson, H., Thun, K., 2014. Containership routing and scheduling in
liner shipping: overview and future research directions. Transportation Science 48 (2),
265–280.
Meng, Q., Wang, T., Wang, S., 2012. Short-term liner ship fleet planning with container
transshipment and uncertain container shipment demand. European Journal of
Operational Research 223(1), 96–105.
Ng, M.W., 2014. Distribution-free vessel deployment for liner shipping. European Journal of
Operational Research 238(3), 858-862.
Ng, M.W., 2015. Container vessel fleet deployment for liner shipping with stochastic
dependencies in shipping demand. Transportation Research Part B 74, 79-87.
Pang, K.W., Liu, J., 2014. An integrated model for ship routing with transshipment and berth
allocation. IIE Transactions 46(12), 1357–1370.
Paat, R., Huseyin, T., 2012. Robust assortment optimization in revenue management under the
multinomial logit choice model. Operation Research 60(4), 865–882.
Poirriez, V., Yanev, N., Andonov, R., 2009. A hybrid algorithm for the unbounded knapsack
problem. Discrete Optimization 6, 110–124.
Rodrigue, J., Notteboom, T., 2013. The geography of cruises: Itineraries, not destinations.
Applied Geography 38, 31–42.
Soriani, S., Bertazzon, S., Cesare, F.D., Rech, G., 2009. Cruising in the Mediterranean:
structural aspects and evolutionary trends. Maritime Policy & Management 36(3), 235–
251.
Song, D.P., Dong, J., 2012. Cargo routing and empty container repositioning in multiple
shipping service routes. Transportation Research Part B 46 (10), 1556–1575.
Song, D.P., Dong, J.-X., 2013. Long-haul liner service route design with ship deployment and
empty container repositioning. Transportation Research Part B 55, 188–211.
Song, D. P., Li, D., Drake, P., 2015. Multi-objective optimization for planning liner shipping
service with uncertain port times. Transportation Research Part E 84, 1–22.
Statista, 2015. Passenger capacity of the global cruise industry from 2007 to 2017.
http://www.statista.com/statistics/301582/passenger-capacity-of-the-global-cruise-
industry/. Accessed 28 Aug 2015.
Tang, L., Sun, D., Liu, J., 2015. Integrated storage space allocation and ship scheduling
problem in bulk cargo terminals. IIE Transactions 48(5), 428-439.
Vacca, I., Salani, M., Bierlaire, M., 2012. An Exact Algorithm for the Integrated Planning of
Berth Allocation and Quay Crane Assignment. Transportation Science 47(2), 148–161.
Veronneau, S., Roy, J., 2009. Global service supply chains: An empirical study of current
practices and challenges of a cruise line corporation. Tourism Management 30, 128–139.
Veronneau, S., Roy, J., Beaulieu, M., 2015. Cruise ship suppliers: A field study of the supplier
relationship characteristics in a service supply chain. Tourism Management Perspectives
16, 76–84.
Wang, Y., Meng, Q., Du, Y., 2015. Liner container seasonal shipping revenue management.
Transportation Research Part B 82, 141–161.
Zhen, L., Chew, E.P., Lee, L.H., 2011. An integrated model for berth template and yard
template planning in transshipment hubs. Transportation Science 45(4), 483–504.
Zhen, L. 2015. Tactical berth allocation under uncertainty. European Journal of Operational
Research 247, 928–944.
This thesis investigates four problems in maritime logistics and operations, where the first
two problems are related to container ships that transport cargos and other two problems are
related to cruise ships that transport passengers. The first problem focuses on the ship type
decision considering empty container repositioning and foldable containers. In this study,
given the ship type with a certain capacity, the problem transfers to a nonstandard minimum
cost flow problem. Then, a network flow model for the problem is formulated. When
considering both standard containers and foldable containers, trouble arises in the network
construction that is some parallel arcs share the same capacity restriction. To overcome this
trouble, a revised network simplex algorithm that changes the standard pivot operation is
designed. Based on the algorithm, a solution approach is proposed to solve the optimization
problem. Some useful managerial implications of this study are obtained after conducting real-
case experiments, which mainly includes three perspectives, i.e., ship type decision, foldable
container usage, and container devanning time.
The second problem addresses the optimal reefer slot conversion for container freight
transportation. This study presents an algorithm to search for the optimal number of reefer
slots to have on a container ship. It is assumed that all the relevant parameters (e.g. freight
rates, storage costs, etc.) are already known. To optimize the sequence of ships deployed, a
revised two-stage simulation approach is proposed. Based on this, a slot conversion algorithm
that determines the optimal slot configurations of the ships is formulated, which embeds the
two-stage simulation algorithm for string optimization. In this study, to highlight the
effectiveness of our approach, several real shipping routes operated by CMA CGM are
investigated. The overall results reveal that the algorithm is highly efficient and can help
shipping liners to significantly improve their profits.
The third problem concerns the cruise itinerary schedule design. It aims to determine the
optimal sequence of a given set of ports of call and the arrival and departure times at each port
to maximize the monetary value of the utility minus the fuel cost. To solve the problem, it first
enumerates all sequences of ports of call and then optimizes the arrival and departure times at
each port of call by developing a dynamic programming approach. To improve the
computational efficiency, effective bounds on the profit of each sequence of ports of call are
proposed. The computational experiments show that, first, the proposed bounds on the profit
of each sequence of ports of call can considerably improve the computational efficiency.
Second, the total profit of the cruise itinerary is sensitive to the fuel price and hence, it is
acceptable to use the shortest voyage distance method to design the schedule when the ports
Table A.1 shows the major specifications of standard containers and foldable containers,
which are almost the same. The specifications of standard containers are collected from APL
(2017) and the specifications of foldable containers are collected from Holland Container
Innovations (2017).
Informal, we can prove the non-increasing trend by using a simplified mathematical model.
Assuming 𝑋 represents the vector for the number of standard containers in ports and 𝑌
𝑌𝑒 𝑇
represents the vector for the number of foldable containers in ports. Then, 𝜌 = 𝑋𝑒 𝑇 +𝑌𝑒 𝑇 shows
the percentage of foldable container usage, where 𝑒 = {1, … ,1} and 𝑒 𝑇 is the transposition of
𝑒. Given the defined vector variables 𝑋 and 𝑌, we can use the following simplified standard
model to represent the formulation of our problem.
𝑀𝑖𝑛 𝐶1 𝑋 + 𝐶2 𝑌
𝑠. 𝑡. 𝐴1 𝑋 + 𝐴2 𝑌 = 𝐵
𝑋, 𝑌 ≥ 0
where all coefficient matrixes or vectors (𝐶1 , 𝐶2 , 𝐴1 , 𝐴2 , 𝐵) are positive. In the next, we can
derive:
𝑋 = 𝐴1−1 𝐵 − 𝐴1−1 𝐴2 𝑌
Appendices 141
𝑀𝑖𝑛 𝐶1 𝐴1−1 𝐵 − 𝐶1 𝐴1−1 𝐴2 𝑌 + 𝐶2 𝑌
Based on which, if the cost coefficient 𝐶2 for foldable containers increase, 𝑌 will decrease. As
1 𝑋𝑒 𝑇 +𝑌𝑒 𝑇 𝑋𝑒 𝑇
a result, 𝑋 will increase. As we have 𝜌 = 𝑌𝑒 𝑇
= 𝑌𝑒 𝑇 + 1, the increasing of cost coefficient
1
𝐶2 will lead to the increasing of 𝜌, that is the decreasing of 𝜌. Such a proof verifies the result
Appendices 142
References
APL, 2017. The perfect fit for all your different needs: your APL equipment guide,
https://www.apl.com/wps/wcm/connect/659ca9b6-7eb3-4c16-82de-
8de04abab7e9/Equipment+guide.pdf?MOD=AJPERES. Accessed on 10 August 2017.
Wang, K., Wang, S., Zhen, L., & Qu, X., 2016. Cruise shipping review: operations planning
and research opportunities. Maritime Business Review 1(2), 133-148.
Wang, K., Wang, S., Zhen, L., & Qu, X., 2017. Cruise service planning considering berth
availability and decreasing marginal profit. Transportation Research Part B:
Methodological 95, 1-18.
Wang, K., Zhen, L., Wang, S., & Laporte, G., 2018. Column generation for the integrated
berth allocation, quay crane assignment, and yard assignment problem. Transportation
Science 52(4), 812-834.
Xia, J., Wang, K., & Wang, S., 2019. Drone scheduling to monitor vessels in emission control
areas. Transportation Research Part B: Methodological 119, 174-196.
Zhen, L., Wang, S., & Wang, K., 2016. Terminal allocation problem in a transshipment hub
considering bunker consumption. Naval Research Logistics 63(7), 529-548.
References 143