Next Article in Journal
A Versatile Aerial Manipulator Design and Realization of UAV Take-Off from a Rocking Unstable Surface
Previous Article in Journal
Recent Development and Future Perspectives of Quality Control and Assurance for the Deep Mixing Method
Previous Article in Special Issue
Container Volume Prediction Using Time-Series Decomposition with a Long Short-Term Memory Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Resource Capacity Requirement for Multi-Terminal Cooperation in Container Ports

1
Department of Industrial Systems and Management, National University of Singapore, 1 Engineering Drive 2, Singapore 117576, Singapore
2
Lee Kong Chian School of Business, Singapore Management University, 50 Stamford Rd., Singapore 178899, Singapore
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(19), 9156; https://doi.org/10.3390/app11199156
Submission received: 4 August 2021 / Revised: 25 September 2021 / Accepted: 28 September 2021 / Published: 1 October 2021

Abstract

:

Featured Application

This study examines the capacity requirement planning for resources utilized in multiterminal operations over time to understand the effect of levels of multiterminal cooperation on resource capacity requirement.

Abstract

Capacity sharing among neighboring terminals offer a means to meet increasing or unexpected demand for cargo-handling without additional capital investment. This study proposes a model for capacity requirement planning of major resources, such as quay cranes (QCs), storage space, and gate, in multiterminal port operations where demand is time dependent. A resource profile simulation is run to generate random events across the terminals and estimate the capacity requirement in the form of workload distributions on port resources over time-shifts. The effects on workload requirement, arising from multiterminal cooperation, are subsequently evaluated in consideration of different container flows among terminals. Experimental results suggest that higher transferring rate between terminals will reduce the QC intensity and storage space requirements but increase gate congestion. Variabilities in the QC intensity and storage space requirements also increase due to shorter stays and more movements in container inventory at the yard. The interaction effect between transferring and trans-shipment rates further shows that the average resource requirements for a terminal can be greatly reduced when the interterminal transferring of containers contributes positively to a more even workload redistribution across terminals. The most significant improvements occur when trans-shipment rate is 85% and transferring rate is 75% for QC intensity; trans-shipment rate is 90% and transferring rate is 60% for storage capacity; and trans-shipment rate is 80% and transferring rate is 75% for gate congestion.

1. Introduction

Container ports provide infrastructure and facilities that are essential for shipping companies to provide a seamless transportation service for cargo from diverse sources to destinations and improve the efficiency of a supply chain as a whole. Proper planning in port capacity, involving estimation of service capability of a system, was shown to alleviate instances of costly underutilization and congestions at any single time point. Taking it a step further, this study proposes a capacity requirement planning model that assesses the capability of the system to meet demands even over time-shifts where there is a sudden surge in demand.
A container port typically consists of several terminals. Each of them can be independently run by a different operator, or the same operator may manage more than one terminal in the port. For instance, PSA (https://www.singaporepsa.com) runs all the container terminals in the Port of Singapore while the terminals of the Port of Hong Kong are operated by several competitive firms [1]. Although the operators are running terminals as a single port entity, the terminals compete with one another to attract shipping liners for their own benefit. Since these terminals are geographically located nearby, they could potentially cooperate to reap greater scale of economics and enhance their market merit. Several container terminals at Port of Busan have the experience to share their wharfs and yards to improve the vessel service though different operators are running their own terminals [2], and the upcoming Tuas Maritime Hub of Singapore operates as a single terminal consisting of four fingers (subterminals) and cooperates with existing Pasir Panjang terminal [3].
One compelling reason for interterminal cooperation via resource sharing is that a standalone terminal may not be able to meet the increasing demand due to insufficient capacity or a mismatch of resource requirement and resource availability. The capacity requirement on the resources could be reduced by sharing these resources among neighboring terminals within the container port such that the capacity shortage in one terminal can be fulfilled with the unutilized capacity in another terminal. Interterminal cooperation also brings about other benefits such as higher bargaining power as a single port entity against other regional competing ports via the provision of large capacity in container handling service. The need for intelligent capacity requirement planning for multiterminal operations is probably even more critical for trans-shipment hub ports, as the ability to provide prompt and adequate handling services for large shipping lines is strongly associated with revenue.
Nonetheless, resource sharing practices among terminals present operational challenges for handling strategies in berth scheduling, storage space allocation, interterminal process planning, etc. Particularly, when a container terminal shares its resources with others, the container flows among the terminals rely on berth schedules that indicate expected arrival times of vessels [4,5]. These interterminal flows involve a subset of trans-shipment containers discharged from a vessel at a terminal but required to be loaded onto another vessel berthing at another terminal. Several studies were conducted in this area. Hendriks et al. [4] studied a multiterminal container port and proposed an interterminal container movement minimization problem of spreading a set of cyclically calling vessels over different terminals and building a berth schedule for the vessels to balance the quay crane workload over the terminals and over time-shifts. Taking the yard as a single resource representing a terminal, Lee et al. [5] investigated the allocation problem of terminals and yards instead of berths to manage interterminal container flows in a multiterminal trans-shipment hub. The survey study conducted by [6] saw that most of the studies on interterminal transportation (interterminal container transportation, which is covering any type of land and sea transportation moving containers and cargo between organizationally separated areas (e.g., container terminals, empty container depots) within a port) discovered their problems at Port of Rotterdam, where many terminals could cooperate within a port. While this survey also suggested research themes for the future, such as transport scheduling and vehicle routing to coordinate container flows between terminals, harnessing information technologies to facilitate data communication, collaborative resource sharing, green (sustainable and economical) logistics, and deployment models for real-world systems, container transportation is still noted as a key operation for a single container terminal. Mishra et al. [7] developed a semiopen queueing network model for interterminal transportation, whereby the service times at the terminal handling stations depend on the number of containers being loaded/unloaded. The proposed model described the uncertainties in the arrival times of containers, the handling times of containers at the different terminals, the travel times of vehicles between terminals, and the interaction between containers and vehicles. Gharehgozli et al. [8] compared various transportation modes for interterminal transportation and applied a game theoretical approach to determine a coalition minimizing the infrastructure cost for collaborative operations among multiple trucking operators. Heilig et al. [9] studied the multiobjective problem consisting of minimizing variable and fixed costs of vehicles as well as greenhouse gas (GHG) emissions for interterminal truck routing. Despite the reduction of empty trips of the truck being found to reduce emissions, the authors highlighted the cost implications arising from the conflict between the number of trucks and empty trips with the reduction in the number of trucks leading to an increase of empty trips, and vice versa. More recently, Jin and Kim [2] provided a time-space graphical model for collaborative operations among multiple trucking operations. Hu et al. [10] formulated the vehicle routing and container transport problem in an interterminal transportation network (including the hinterland railway terminal) and optimized the container delivery based on fixed fleet size. Hu et al. identified a tradeoff between the degree of interterminal transportation and the delay of rail transportation. Schepler et al. [11] formulated a coordination optimization model providing the berthing positions and time windows to serve vessels, the rail tracks, and time windows to serve trains, as well as the time windows for truck appointment under the interterminal transportation environment. The objective of the study is to minimize weighted turnaround time of vessels or trains.
From the existing literature, it is not difficult to see that much attention was paid to the operations planning and scheduling of container transportation among terminals in a port. However, while many scholars have supported the idea of interterminal cooperation, there is no study that examines the capacity requirement measured in terms of workload on the resources over time-shifts. An evaluation of workload information on the resources (i.e., the availability and requirement) serves three purposes. Firstly, the provision of workload information on critical resources is necessary to confirm the feasibility of the operations schedule design. In particular, knowledge on the changes in workload requirement over time-shifts boosts the effectiveness of the operations scheduling in achieving a desirable throughput volume at the terminal and port levels. Secondly, workload information supports long term planning of capacity investment in a terminal while taking care of the sporadic demands that may or may not be met with the available capacity of its counterparts. And most importantly, when dealing with multiterminal cooperation, the provision of workload information on the shared resources becomes particularly imperative to allow a fair allocation of workload and utilization of resources.
This paper investigates the efficiency gains from multiuser terminal (studies on port cooperation typically focus on the cooperation and competition among nearby ports or terminals belonging to different ports, and hence, the cooperation is conditionally recommended or discouraged [12,13,14,15,16,17]) co-operations where a single port entity serves calling vessels from different shipping lines and inbound, outbound and trans-shipment containers arrive at terminals at different rates. While some of these terminals may be operating under a single operator, it is interesting to note that the terminals basically function independently of one another with little coordination among themselves. Port of Singapore at Pasir Panjang is one example. Such governance arrangement is made on premise that interterminal competition with a port promotes port productivity as a whole. However, as competitions among peripheral ports intensify with the advancements in logistical system, the conception of having terminals to work together towards a common good of the port is being put forward. The paper proposes a design process for a system-wide capacity requirement planning when terminals within a port cooperate with one another. The logistics process of a container port can be seen as a typical batch process in production planning as containers arrived vessel-by-vessel. According to Tenhiälä [18], among the classical capacity methods (that includes rough-cut capacity planning, finite loading with capacity leveling and finite loading with optimization), capacity requirement planning is the most suitable for batch processes. By providing a detailed technique for ascertaining the feasibility of material plans, capacity requirement planning enables the computation of the utilization of individual resources through resource requirement distributions generated by the container arrival plans and help to generate more realizable operations schedules [19,20].
The analysis begins by first generating a set of resource profile that represents the specific timing of the projected cumulative workloads at individual resources. The resource profiles form a basis for capacity requirement planning, of which a resource profile simulation is developed to estimate the changes in the workload of resources with container arrivals and handling across the neighboring terminals over time-shifts. This study investigates the capacity requirement for the three representative resources, namely, QCs, storage yards, and gates. Results on key performance indicators such as QC intensity, the TEU (twenty equivalent unit) workload at the yard, and the queue length at the gate are recorded for varying inbound, outbound and trans-shipment container arrival rates, and compared over a range of transferring container volumes across terminals. This study enriches the applications of capacity requirement planning beyond the classical applications of production systems. In its consideration of interterminal cooperation and logistics, the study integrates probability models and discrete-event simulation to generate resource profiles and estimate the capacity requirement resulting from resources sharing.
The rest of this study is organized as follows: Section 2 designs the capacity requirement planning for multiterminal operations via simulating resource profiles; Section 3 conducts the simulation experiment; Section 4 draws the management implications from the simulation results; and Section 5 concludes this study.

2. Resource Profiles for Representing Workloads

According to Jacobs et al. [19] and Fogarty et al. [20], capacity requirement planning in container ports can make use of the planned container arrivals from various modes (i.e., vessels, rail, and trucks) of container transportation to calculate the workload for specific resources. To estimate the handling capability of the resources over time-shifts, capacity planners require information on the capacity requirement for different resource types to meet the workload. In particular, capacity planning should take into account the interterminal container flows and deployment of shared resources arising from the cooperation efforts among neighboring terminals.
When an operational task (for examples, unloading or loading container onto a vessel, retrieving a container from yard, or putting a container into storage etc.) is distributed to resources, the corresponding workload is formed at each resource. The workload distribution representing the capacity requirement on the resources at different time points are recorded in individual resource profiles that can be used to manage the operational tasks over time-shifts. Section 2.1 and Section 2.2 introduce resource profiles and cumulative workload distributions for a resource, Section 2.3 provides the resource profile simulation to examine capacity requirement on resources when events leading to cargo handling are triggered over time-shift in multiterminal operations. This study uses the terminology TEUs instead of containers to indicate workload for container arrivals and/or handling.

2.1. Resource Profiles

Three major resources in a multiterminal container system are the wharf, gate, and yard. Quay cranes (QCs) on the wharf are interfacing resources between vessels and a terminal; the gates are interfacing resources between the terminal and consignees, and the yards are a representative resource for container storage supporting the operations process across the unloading, loading, receiving, and delivering. A resource profile that considers the varying of the capacity requirements for each type of resource over time-shifts is drawn out. This time-phased resource profile provides a means of estimating the workload on a resource throughout the time-shifts, and thus, assess the capacity adequacy of the system in managing the sporadic demands.
Figure 1 represents the resource profiles of a yard and QCs graphically for two terminals when a vessel arrives at a terminal. The relevant notations are described as follows:
t o p e n
= Time for starting the window for outbound or trans-shipment TEUs to arrive at a terminal. It is typically set to two weeks before the expected vessel arrival.
t c l o s e
= Time for ending the window for outbound or trans-shipment TEUs to arrive at the terminal. It is typically set to one day before vessel arrival.
T a r r i v a l
= (Random variable) Time that the vessel arrives at the terminal. The ETA of a vessel is nominated when constructing a berth schedule. For purposes of simplicity, this study excludes the berth scheduling, and the vessel arrival is set to follows a probability distribution.
t u n l o a d
= Time for starting the unloading operations for the vessel. It is typically set to 2 h after vessel arrival.
T s w i t c h
= (Random variable) Time for switching operation modes. The duration between time when an operation mode is switched and the time when the unloading begins, i.e., T s w i t c h t u n l o a d , depends on the QC intensity requested by the vessel and the amount of TEUs for the unloading from the vessel.
T l o a d
= (Random variable) Time for TEUs to be loaded onto a vessel, transferred to another terminal (if any) or delivered to mainland. Typically, time window T l o a d T s w i t c h = 0 .
T e n d
= (Random variable) Time for completing the loading operation for a vessel. T e n d T l o a d also depends on the QC intensity requested by a vessel and the amount of TEUs for the loading operation of the vessel.
t d e p a r t
= Time for vessel departure. It is typically set to 2 h after the loading completion.
T t r a n s f e r
= (Random variable) Time for TEUs to be loaded onto trucks for transfer to another terminal for trans-shipment on another vessel. Trans-shipment TEUs typically leaves the terminal in 2–3 days after unloading.
T f i n i s h
= (Random variable) Time that TEUs left the terminal. T f i n i s h T t r a n s f e r is dependent on the travel time of a truck between the two terminals and the amount of TEUs scheduled to be transferred.
C a
= (Random variable) the amount of TEUs arrived in the terminal before the arrival of the mainline vessel. The actual arrivals of TEU quantity are uncertain but variations are usually minimal.
C b
= (Random variable) the amount of inbound and trans-shipment TEUs unloaded from the mainline vessel. There is uncertainty in the amount of TEUs for unloading with a small variance.
C c
= (Random variable) the amount of outbound and trans-shipment TEUs loaded to the mainline vessel to be transferred to another terminal or delivered to mainland. It is also uncertain with a small variance.
C d
= (Random variable) the amount of TEUs received from transferring operations. It depends on the berth schedule and the TEU information of the corresponding vessel. There is also a small level of uncertainty.
C a f
= (Random variable) C a for a feeder vessel.
C b f = (Random variable) the amount of inbound TEUs unloaded from the feeder vessel.
C c f
= (Random variable) the amount of outbound TEUs loaded to the feeder vessel.
Figure 1 depicts the container storage and retrieval activities associated with C a , C b , C c , and C d . The resource profiles for C a , C b , C c , and C d can be represented by cumulative workload distributions. The QC intensity refers to the number of QCs allocated to a vessel and is conditional upon the vessel category and the TEU quantity. The QC intensity affects the unloading and loading rates, and hence, can be used to represent the workload of QCs at the resource profile. Since the yard is a storage space for containers, the number of containers for the accommodation can be straightforwardly regarded as the workload. The truck queue length is a proxy for gate workload.
When a mainline vessel arrives, the outbound containers from mainland and the transferred outbound containers from other terminals were stacked in the yard in advance (TEU workload in C a ). The TEUs from the unloading operations of the vessel temporarily increases the TEU quantity in the yard (TEU workload in C b ) as much as the inbound and trans-shipment TEUs. The stacked outbound and transferred containers are loaded onto the vessel and the corresponding TEU quantity (TEU workload in C c ) decreases during the loading operation. A part of unloaded containers can be transferred to another terminal (TEU workload in C d ) for trans-shipment according to the predetermined berth schedules of the terminals.
When applying the resource profiles for feeder vessels, as depicted in Figure 2, the resource requirement is relatively simple compared to that of mainline vessels. Since the feeder vessels transport containers plying between trans-shipment and neighboring terminals, there are only unloading and loading operations at a terminal without the need of transferring to another. When a given amount of TEUs are unloaded from a feeder vessel, they are eventually loaded onto a mainline vessel resulting in the increase of TEU quantity in the yard (TEU workload in C b f ). Meanwhile, when a feeder vessel is scheduled to load TEUs from the terminal, a small proportion of total TEUs arrived earlier at the yard (TEU workload in C a f ) from a mainline vessel, and a certain amount of TEUs are loaded onto the feeder vessel (TEU workload in C c f ).
Comparing the unloading and loading processes, the speed of unloading speed is faster at the beginning of operations and subsequently slows down in preparation for the switch to loading operations. The unloading operations, especially for containers on the deck, can usually be done efficiently with a well-designed vehicle dispatching strategy. On the other hand, the loading process is relatively slower because the process needs to consider the loading plan (i.e., terminal-level stowage plan) and ensure the ship’s stability and safety. The loading process are also dependent upon the supporting operations of vehicles and yard cranes that need to be carefully sequenced and dispatched. Meanwhile, the hatch cover handling operations, involving the uncovering and covering of the hatch, increases the durations of both the unloading and loading process as the quay cranes switch the spreader to wired hooks.

2.2. Workload Distributions

The random variables C a , C b , C c , C d , C a f , C b f ,   and   C c f change the workload of the resources over time-shifts. The set of randomness of the random variables, described in the above subsection, also generate probability distributions for requesting workload for storage and retrieval operations as conceptually depicted in the Figure 3.
The cumulative conditional workload distribution for TEU arrivals before the mainline vessel arrival, associated with the random variable C a , is
W a r r i v a l t o p e n t t c l o s e |   T a r r i v a l = t a r r i v a l = i = t o p e n t E C a P C a | T i ,
where T a r r i v a l is a random variable with mild variation and P C a | T i is the cumulative distribution of outbound TEU arrivals from t o p e n . The cumulative workload distribution is conditioned by T a r r i v a l as t o p e n and t c l o s e are determined by the realization of T a r r i v a l . The cumulative conditional workload distribution of the unloading operation is
W u n l o a d T s w i t c h T a r r i v a l + t u n l o a d t   |   T a r r i v a l = t a r r i v a l = i = T a r r i v a l + t u n l o a d t E C b P C b | T i = m i n t i H Q U , E C b ,
where H Q U is the random variable for the QC gross productivity of a vessel during the unloading operation measured by the amount of TEUs handled per time unit. P C b | T i is replaced with the time-independent unloading rate. The gross productivity is dependent on the QC intensity, which in turn depends on the vessel size and the TEU volume. Hence, H Q U = h Q U I Q , where h Q U is the gross productivity distribution of a QC for a vessel for the unloading operation, which can be realized by the number of moves per time unit. For S being a random variable with a small deviation for the vessel category, I Q = F C b + C c , S is the QC intensity distribution for the vessel operation that varies with the amount of TEUs and the vessel size as summarized in Appendix A. Accordingly, the cumulative conditional workload distribution for the loading operation is given as,
W l o a d T e n d T l o a d t   |   T l o a d = t l o a d = i = T l o a d t E C c P C c | T i = m i n t i H Q L , E C c ,
where H Q L = h Q L I Q . h Q L is the gross productivity distribution for the loading operations. Note that the workload for the loading operation take a positive effect on reducing the storage space requirement but it requires handling effort to be performed. As for the transferring operations from a terminal to another, the cumulative conditional workload distribution is
W t r a n s f e r T f i n i s h T t r a n s f e r t   |   T t r a n s f e r = t t r a n s f e r = i = T t r a n s f e r t E C d P C d | T i ,
where P C d | T i represents the transferring productivity of trucks that transport TEUs out to other terminals including the traffic congestion and the gate process. T t r a n s f e r should correspond to the t o p e n for another vessel incoming to another terminal (Terminal 2 in Figure 1).
The same goes for the workload distributions C a f , C b f and C c f of the feeder vessels in Figure 2. The cumulative conditional workload distribution for the unloading operation of a feeder vessel is
W u n l o a d f T e n d T a r r i v a l + t u n l o a d t   |   T a r r i v a l = t a r r i v a l = i = T a r r i v a l + t u n l o a d t E C b f P C b f | T i = m i n t i H Q U ,   E C b f ,
where H Q U = h Q U I Q consisting of the gross productivity distribution of a QC for a feeder vessel during the unloading operation h Q U and the QC intensity distribution dependent on the amount of TEUs and the feeder vessel size I Q = F C b f + C c f , S . For a feeder vessel requiring only the loading operation, the cumulative conditional workload distribution for TEU arrivals for a feeder vessel is
W a r r i v a l f t o p e n t t c l o s e |   T a r r i v a l = t a r r i v a l = i = t o p e n t E C a f P C a f | T i
The cumulative conditional workload distribution for the loading operation of a feeder vessel is
W l o a d f T e n d T a r r i v a l + t l o a d t   |   T a r r i v a l = t a r r i v a l = i = T a r r i v a l + t l o a d t E C c f P C c f | T i = m i n t i H Q L , E C c f ,
where H Q L = h Q L I Q consisting of the gross productivity distribution of a QC for the loading operation h Q L .

2.3. Resource Profile Simulation

A resource profile simulation is developed to generate events imposing workload requirement for resources, such as vessel arrivals, truck transportation, loading and unloading operations, etc. presented in Figure 4, and simulate the workload cascading over the resources across terminals. The events are associated with the random variables provided in the previous section and each event triggers another. The resource profiles simulation estimates the workload distributions for resources to estimate the capacity requirement over time-shifts across the terminals.
  • E a i , x -Vessel arrival. This event is triggered when a vessel with type i arrives at the port and randomly assigned to any of one of the terminals x M with an equal probability, where M is the set of cooperative terminals at the port. The vessel type is characterized by the attributes such as the category, the TEU quantity, the required QC intensity. When the event is executed, the attributes of the vessel are randomly assigned by referring to Appendix A. The event triggers an immediate execution of E u , E t , and E s and schedule a new event E a i , x with an exponential interarrival rate λ . This event is associated with the workloads C a and C a f .
  • E t -TEU arrival. This event is triggered when a vessel arrives at the port. The arrivals of TEUs was recorded earlier than the time of executing E a by tracing back in time. The workload C a is distributed only in the time ranged from t o p e n to t c l o s e . Both t o p e n and t c l o s e are dependent on the generated T a r r i v e .
  • E u -Unloading. This event is triggered when the first vessel in the queue arrives at a randomly selected terminal. The event execution time is delayed until t u n l o a d as both the vessel and the QCs operations must be ready for the unloading operations (which generally takes 2 h) to be performed. The event allocates QCs to the vessel as requested (generated) in the event of E a . It is set that the required number of QCs is always available at a terminal. The time window for the unloading operation T s w i t c h t u n l o a d is conditionally determined by the generated H Q U . When the event is triggered, the amount of TEUs in the yard increases at the rate of H Q U . The workloads C b and C b f are associated with the event. The simulator sets unloaded TEUs to leave the terminal within 5 days following the unloading. The unloading event E u triggers the execution of the loading event E l through determining T s w i t c h .
  • E l -Loading. The event is triggered immediately upon the completion of the unloading operation on the vessel. T s w i t c h is used to execute the loading event E l but the event execution time is delayed until T l o a d . Similar to the unloading operation, the loading time window T e n d T l o a d is conditionally determined by the generated H Q L   while the TEUs at the yard also decreases accordingly. The workloads C c and C c f are adjusted with the event. The event triggers the execution of the vessel departure event E d though determining T e n d .
  • E d -Vessel departure. The execution time of this event is delayed until t d e p a r t , which is determined by the generated T e n d as the both the vessel and the QCs operations must be ready to be released. When executing this event, the vessel returns the assigned number of QCs, releases the quay space of the terminal, and triggers an immediate execution of E a .
  • E s -TEU Transfer. This event is triggered when the event of vessel arrival E a i , y where y M \ x is executed. The event generates C d over the time window T f i n i s h T t r a n s f e r for interterminal trans-shipment at each cooperating terminal, with the time window for trans-shipment TEUs being 2–3 days after unloading at a terminal. Hence, the TEU transferring activities are recorded earlier than the time of executing E a i , y . T t r a n s f e r is generated by the execution of E a i , y and T f i n i s h is dependent on the transferring productivity of trucks.
  • E t -Truck transport. This event is triggered when E a i , y determines T t r a n s f e r . An amount of TEUs at a terminal, C d , will need to be transported to terminal y from x by trucks. A generated single truck typically transports 1.5 TEUs at a time. When a truck transports transferring TEUs from a terminal to another, the process takes a number of time-units including the congestion on traveling paths. The congestion has the effect of reducing the speed of truck flows. Underwood formula is used to adjust the truck speed for the traffic congestion on a traveling path, taking into account parameters such as length of the travel path, free flow speed, and a predetermined amount of maximum number of trucks. The effective traveling speed u n in a traveling path between two terminals with n trucks is estimated by
    u n = u f e x p k n k 0 ,
    where   u f is the free speed and   k 0 is the traffic density at maximum flow. For L t representing the length of the traveling path, k n = n 1 / L t   is the traffic density in the traveling path with n trucks. The event triggers an immediate execution of E g with the delay of traveling and congestion of a truck.
  • E g -Gate check-in. The execution time of this event is delayed until the arrival of a truck, as the traveling time and delays caused by congestion should be included. Since the gate is a resource with limited capacity, the queue length will be a meaningful measure for gate workload. The gate is modeled as a single server machine and its service times are adjusted for multiple pass lanes. Whenever a truck arrived at a gate (i.e., the execution of E g ) , it joins the gate queue and waits for the check-in service. The truck queue length is affected by the workload distribution of Equation (4) and counted as the gate congestion measure in the simulation. The gate check-in event triggers an immediate event E c .
  • E c -Arrival of Transferred TEU. This event is triggered when a truck passes the gate, but the execution time is deferred until the truck is processed at the gate. It increases the TEU workload in the yard accordingly at most C d . This event generates T f i n i s h , and calls the event E s to transport TEUs continuously until cumulative probability for TEU workload, P C d | T i , becomes 1. Transferring time window T f i n i s h T t r a n s f e r essentially includes the two delay elements, namely, congestion and queueing that increase uncertainty.

3. Simulation Experiment

A set of simulation experiment is conducted to understand the capacity requirements for different volumes of trans-shipment ratios of containers relative to inbound and outbound containers at a terminal, as well as different transferring container volumes across terminals. The base-case setting is as follows: four equal-size of terminals located at a trans-shipment hub port are sharing the resources. Each terminal has 1000 m quay length and initial inventory at the yard is 20,000 TEUs. The distance between adjacent terminals is set to 500 m, and the terminals are sequentially laid in a line. The truck free speed in between terminal is set to 14 m/s. The traffic density at maximum flow ( k 0 ) in the Underwood formula is calculated by setting the effective speed under jammed traffic ( u N m a x ) to 0.01 m/s, where the maximum number of trucks ( N m a x ) for transferring activities between terminals is limited to 25. The gate service time is preset to follow U 5 ,   10 min. This service time measures the duration-of-stay of a truck in the boundary of gate system and includes the time spent on checking the driver identification, truck number, container number, and measuring the weight, temperature (for reefer), and size (for out-of-gauge) of the container and deciding a storage slot for the container. Additionally, when the container door is not facing outward, a reach stacker needs to turn the direction of the container.
To avoid the situation whereby insufficient resource capacity leads to a distortion in the actual capacity requirement reflected in the resource profile, the study assumed that the requested QC intensity of every vessel can be met. The QC intensity is determined by the vessel category (size) and the TEU quantity on the vessel described in Appendix A. The vessel and TEU categories in the experiment is drawn upon their respective pools based on a uniform distribution. The gross productivity of a QC is assumed to be normally distributed, N 28 ,   4 2 , hourly. t o p e n and t c l o s e are set to 2 weeks and 1 day, respectively, and the vessel arrival T a r r i v a l is set to follows exponential interarrival time with a mean of 2.05 hourly. t u n l o a d is set to 2 h after the vessel arrival. It is further assumed that the incoming vessels arrive at any terminal with equal probabilities. The probability that an arriving vessel is a mainline vessel is 90%, while the probability of receiving a feeder vessel is 10%. When the latter occurs, half of these vessels are for unloading and the other half for loading only (note: this contrasts with the mainline vessel which will first perform the unloading operations and then the loading operations). The terminals have equal probabilities to be selected when transferring the trans-shipment TEUs to them.
The simulation scenarios are represented by varying the rates of transshipment and outbound and inbound TEUs, as well as the rates for which the trans-shipment TEUs would remain in its arrival terminal or be transferred to other terminals. The simulation runs 10 replications for each setting for 200 days after the warming-up period. Since the warming-up period depends on experiment settings, the simulation results are sampled out after eliminating the initialization bias by shown in the output plots. The simulation model is implemented using C# programming language as part of the Visual Studio Community 2017 and run on a computer with Intel processor i5, 1.8 GHz with 16 GB RAM.

3.1. Simulation Results with Low Transferring

The proportion of TEUs for inbound and trans-shipment-inbound are set to 20% and 80%, respectively, and the same proportion applies to the TEUs for outbound and transshipment-outbound. The proportions of TEUs remaining in a terminal and transferring to the other terminals are set to 40% and 60%, respectively. All the other terminals have equal probabilities to take the transferred TEUs and the transferring is only limited to trans-shipment TEUs.
Based on the collected samples, X ¯ R Charts for QC intensity, the TEU workload at the yard, and the queue length at the gate can be drawn with the estimated control limits consisting of upper control limits (UCLs), central lines (CLs), and lower control limits (LCLs) for each set of output [21]. Specifically, the simulation plots the QC intensities currently serving all vessels at the terminal whenever a vessel leaves a terminal ( E d ). Figure 5 shows the simulation results of QC intensities over time-shifts. The mean requirement of QC intensity for vessel service is 8.12 and the upper limit is estimated as 9.78. It means that a terminal needs to prepare 9.78 QCs to meet the handling capacity requested by the incoming vessels, even though on average only 8.12 QCs are needed to meet the capacity request.
The storage space requirement in TEUs is plotted when any of these three events takes place: (i) a vessel arrives at the terminal ( E a ); (ii) the interterminal transferring is completed ( E c ), and (iii) a vessel departs ( E d ). Figure 6 shows that the mean storage capacity required to support the requested vessel service is estimated to be 117,109.64 TEUs, with 125,006.81 and 109,212.50 TEUs as the upper and lower control limits, respectively. The mean difference between maximum and minimum storage requirements is estimated as 13,686.60 TEUs.
Gate availability is critical for the smooth execution of interterminal transferring activities. A long truck queue is usually an indication of gate congestion. Figure 7 plots truck queue length at every instance when a truck reaches the gate ( E g ). The estimated queue length is 1.56 trucks, and the upper limit shows 2.32 trucks in the queue. Thus, it can be inferred that the gate needs to possess adequate capacity to clear an average queue length of 1.56 trucks to provide the handling service requested by vessels during the transferring activities. Note that the truck queue is estimated only for transferring containers. The results are summarized in Table 1.

3.2. Simulation Results with High Transferring

A parallel experiment is conducted to examine the effect of a high transferring TEUs rate. The proportions of TEUs remaining in a terminal and transferring to the other terminals are set to 25% and 75%, respectively. The proportion of TEUs for inbound, trans-shipment-inbound, outbound, and transshipment-outbound remains to be same as Section 3.1. The simulation results are drawn in Figure 8, Figure 9 and Figure 10 using X ¯ R Charts and summarized in Table 2.
The simulation results of QC intensities for all vessel arrivals over the run-time suggest that a terminal is expected to prepare 9.71 QCs if it wishes to meet the handling capacity requested by the incoming vessels entirely, but on average 7.92 QCs are sufficient to meet the capacity request (Figure 8). The QC intensity is generally lower than what was found previously in Section 3.1. It indicates the workload in the quayside operation will be reduced by 2.5% when the transferring rate increases from 60% to 75%.
Figure 9 shows that the mean storage capacity requirement to support the requested vessel service is estimated to be 107,120.78 TEUs, with 127,060.92 and 87,180.64 TEUs as UCL and LCL, respectively. Compared to the results in Section 3.1, the increase of transferring rate has helped reduce the storage capacity requirement by 8.5% on average. However, the gap between the two limits has also widened by 152.5% at the same time. This widened gap indicates increased variation of storage capacity requirement with the mean difference between maximum and minimum storage requirements estimated at 34,558.30 TEUs. Meanwhile, the estimated truck queue length is 1.69 trucks with 2.52 in UCL for transferring TEUs (Figure 10). The gate processing requirement is increased by 8.3% in terms of the mean queue length as the transferring TEUs increase.

4. Managerial Insights

The simulation results in Section 3 gives the capacity requirements for a range of setting of transferring TEUs rates across terminals when each terminal handles a combination of inbound, outbound, and trans-shipment TEUs. This section further discusses the managerial insights over various combinations of TEU workload requirement for multiterminal cooperation.

4.1. Analysis on Workload Significance

In the previous section, the simulation results are obtained under the experimental setting of 80% trans-shipment TEUs and the remaining 20% equally split into inward-bound and outward-bound TEUs. This forms the base model for comparison.
Further experimental analysis allows the summation of equal inbound and outbound TEU rates to be determined by the TEU rate for transshipment, i.e., the larger the percentage of trans-shipment TEUs, the lower the aggregate inbound and outbound TEUs. Similarly, the rate for TEU remaining in a terminal is also dependent upon the transferring rate to other terminals since both rates must sum to 1. As before, only trans-shipment TEUs can be transferred, and the other terminals have equal probabilities to take up the transferred TEUs. The purpose of this extended analysis is to understand the effects of the relative volumes of trans-shipment containers against those of the inbound and outbound containers and the volumes of transferring containers across terminals on the capacity requirement.
The four levels of the trans-shipment rates (i.e., Upper-90%, Higher-85%, Norm-80%, and Lower-75%) and the two levels of the interterminal container transferring (i.e., Higher-75% and Lower-60%) are simulated. Each combination has a sample of size 200. It means that 200 time-points of vessel departures are randomly selected from the simulation results and the relevant dataset (for the selected 200 time-points) for the QC intensity, the storage space requirement, and the gate congestion is taken. The summarized descriptive statistics are presented in Appendix B.
Figure 11a compares the mean QC intensities across the trans-shipment rates and the transferring rates. From the figure, the QC intensity is increased by 4.1% and 12% when transferring rates are at a high and low of 75% and 60%, respectively, when the trans-shipment rate is increased up to 80% from 75%. The increase in QC intensity is logical, as the shortened duration-of-stay of containers at yard increases higher handling request at the quayside. The smaller increment observed under the higher transferring rate is due to the higher support of distributing workload by means of transferring. Further increase to the trans-shipment TEUs rates (for example, 85%) leads to 5.6% and 4.1% reduction of QC intensity under higher and lower transferring rates, respectively. This reduction can be attributed to the possibility of transferring workload to other terminals, which in turn, translates into better chances for a more evenly distributing QC handing workload across terminals. Subsequent increase in the trans-shipment rate to 90% produces an increase of QC intensity by nearly net of 0.9 and 0.3%, respectively. At this level of trans-shipment rate, the workload distribution effect of transferring work is dominated by the even higher handling request for QCs.
Another important element concerns the uncertainties in resource requirements. As shown in Figure 11b, despite having a high level of the trans-shipment rate (i.e., 90%) that provides opportunities to reduce QC intensity through TEUs redistribution across terminals, higher uncertainty in QC intensity requirement is reported when there are active transferring activities among the terminals (i.e., transferring rate at 75%). On the other hand, the lower transferring rate (i.e., 60%) contributes to the monotonic decreases of QC intensity as the trans-shipment rate increases.
When the trans-shipment rates are reduced, the corresponding rates for outbound and inbound TUEs increase. This leads to longer expected duration-of-stays for TEUs at the yard and increasing storage space requirement (Figure 12a). When the transferring TEUs rate to other terminals is lowered from 75% to 60%, the terminal will require a 1.5% higher storage space capacity on average as more space is needed to accommodate containers at the yard. However, as the trans-shipment rate rises up to 90%, the storage space capacity associated with the lower transferring rate (60%) drops by 14.2%. While the shortened expected duration-of-stays associated with the high rate of trans-shipment TEUs (i.e., at 90%) leads to the reduced storage space capacity requirement, the higher rate of transferring TEUs increases the requirement for storage space capacity to put up incoming containers (thereby producing an offsetting effect for storage space capacity savings). The offsetting effect revisits when the trans-shipment rate is lowered to 75% leading to the storage space requirement increases to nearly 2.6%.
Figure 12b shows that the requirement uncertainty on storage space is somewhat unaffected when the trans-shipment rate decreases from 90% to 80%. The requirement uncertainty on storage space saw a sharp drop of 68.1% with further reductions in the trans-shipment rate from 80% to 75% under the higher transferring rate of 75%. As the trans-shipment rate falls, the expected duration-of-stays for TEUs at the yard lengthens and reduces the requirement uncertainty though the higher transferring rate (i.e., 75%) works in the opposite direction. Further reduction of trans-shipment rate (i.e., 75%) magnifies the uncertainty reduction of storage space requirement with longer duration-of-stays for TEUs. The requirement uncertainty on storage space is much variable for the trans-shipment rates when the transferring rate is lower (i.e., 60%). Overall, there is a tendency towards the declination of the requirement uncertainty when the trans-shipment rate decreases for the same reason as that of the higher transferring rate.
The truck queue length is generally shorter when the rate of trans-shipment TEUs and interterminal rates of transferring TEUs are low (Figure 13a). On the contrary, under the higher transferring rate, increasing the rate of trans-shipment TEUs results in increased traffic congestion on the path transporting containers between terminals. This leads to a reduction in arrival rate of trucks which, in turn, results in an average 2.4% reduction in queue length at the gate. When the trans-shipment rate is decreased to 75%, the traffic congestion is less effective in alleviating gate congestion. Instead, a reduced amount of trans-shipment TEUs transferred is shown to be more effective in reducing the truck arrival rate with the truck queue at the gate being shortened by 4.1%. The plot of gate congestion range in Figure 13b behaves very similarly to that of queue length in Figure 13a since the uncertainty on truck queue is attributed to the traffic congestion for trucks transporting containers between terminals, as well as the arrival rate of trucks at the gate. Variations of truck queues at the gates decrease with a shorter average queue length and increase with a longer queue of trucks.

4.2. Analysis on Interaction Effects

A two-factor ANOVA with replications is applied to examine the effect of trans-shipment and transferring rates on capacity requirement of QC intensity (Table 3), storage space requirement (Table 4), and gate congestion (Table 5). The p-values reject the null hypothesis in favor of the alternative hypothesis and conclude that the difference in trans-shipment rates has an effect on the QC intensity, storage space requirement, and gate congestion at the 95% significant level. The same applies for the transferring rates and the interaction effect of trans-shipment and transferring rates.
The ANOVA contrast determines whether there is a significant difference in levels of the QC intensity, the storage space requirement, and the gate congestion when the factor levels (i.e., transferring and trans-shipment rates) vary. The number of contrast coefficients is set to 2, and Bonferroni correction is applied when testing the significance between levels of a factor.
Table 6 suggests statistical evidence towards the association between trans-shipment rates and storage space requirement, as well as trans-shipment rates and gate congestion. On the other hand, there is no significant difference in QC intensity in the case for trans-shipment rates of 90% and 85%. As explained earlier in Figure 11, the higher levels of trans-shipment rate would trigger the higher volume of transferring containers, and hence the QC requirement is shared among the collaborating terminals. The rest of cases demonstrate that changes in trans-shipment rates have a significant effect on QC intensity, storage space requirement, and gate congestion.
The interaction terms consider if the effect of one predictor variable, in this case the transferring or trans-shipment rates, on the response variable, i.e., QC intensity, is different at different values of the other predictor variable. Unlike the results in Table 6, Table 7 shows a significant difference QC intensity between a 90% trans-shipment rate and an 85% trans-shipment rate, whereas an insignificance difference is observed between 90% and 75% of trans-shipment rates under a 75% transferring rate. When transferring rates among the terminals are high, the surge in trans-shipment rates smooths out the QC intensity across the terminals and produces a similar effect to the reduction in the trans-shipment rates that attenuates the QC intensity for trans-shipment TEUs. This finding is underpinned by nonexistence of insignificant interactions in the lower levels of transferring rate (for example, at 60%), whereas high transferring rate contributes the nonsignificant relationship between the two extreme trans-shipment rates.
Meanwhile, the differences in degree of gate congestions between the 80% and 75% trans-shipment rates with a 75% transferring rate, and between 90% and 85% trans-shipment rates with a 60% transferring rate are found to be insignificant. The estimates in Cohen’s d, which indicate the standardized mean differences between the two groups to be about 9% and 15%, provides further support to these observations (Table 8). As effect size is known to be independent of the sample size, the nonsignificances in the two cases concluded through p-value in Table 7 can rightfully be attributed to the close similarity between the two samples in terms of a comparison of percentiles [22]. This similarity can be explained by the fact that truck arrivals are controlled by the traffic conditions using Equation (8) for the high transferring rates.

5. Conclusions

When a container terminal collaborates with its neighboring terminals within a port or across ports in the form of resource sharing, the operations in the collaborating terminals inevitably become more complicated. To realize the gains from such collaborations, an operations management system that can make use of the available resources efficiently and effectively needs to be in place. This research studies three of the most important port resources, namely, QCs, yard, and gate, to manage the capacity requirements over terminals. A resource profile simulation, which provides a platform for simulating random components of resource profiles and estimating the workload on the resources over time-shifts, is developed. The estimated workloads decide capacity requirement on the resources as represented by the QC intensity, the storage space requirement, and gate congestion.
The experiment results provide X ¯ R charts for QC intensity, storage space requirement, and gate congestion. The significance tests are examined for the rates of trans-shipment containers compared to inbound and outbound containers, and the rates of transferring containers to other terminals. For different settings of trans-shipment and transferring rates, there was a statistically significant effect of the trans-shipment and transferring rates on the capacity requirement representing workload of the QCs, yard, and gate. The two-factor ANOVA results supported the statistical significance of the trans-shipment and transferring rates, as well as their interaction terms on capacity requirements. In particular, high trans-shipment rates are found to contribute to the reduced workloads for QCs, yard, and gate in terms of capacity requirement when the transferring rates take positive effect on distributing workloads among the terminals for the trans-shipment rates.
The realism of the discussed multiterminal capacity requirement planning could be increased with a study of fidelity that estimates the capacity requirements at the equipment-level (e.g., yard cranes, vehicles, etc.) and operational-level (e.g., single vs. dual cycling) while considering the financial commitment (e.g., operational and capital costs for resources).

Author Contributions

Conceptualization, B.K.L.; methodology, B.K.L. and J.M.W.L.; software, B.K.L.; validation, J.M.W.L.; formal analysis, B.K.L. and J.M.W.L.; investigation, B.K.L.; resources, J.M.W.L.; data curation, B.K.L.; writing—original draft preparation, B.K.L.; writing—review and editing, J.M.W.L.; visualization, B.K.L.; funding acquisition, J.M.W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Maritime and Port Authority of Singapore (MPA) Research Fellowship Grant.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Estimating QC Intensity

Port operators have their own operational practices to allocate the number of QCs to a vessel according to the amount of TEUs and the vessel category (size). Note that feeder vessels belong to the vessel category 1. Table Table A1, Table A2, Table A3 show the parameters of the triangular distributions generating the vessels and TEU quantity, respectively. The QC intensity is ranged from 1.0 to 6.0 based on the vessel and TEU categories that a generated vessel and its TEU quantity belong to, respectively. The collected data were from various port operators and consulted by an industry collaborator.
Table A1. Parameters applied to triangular distribution generating vessels.
Table A1. Parameters applied to triangular distribution generating vessels.
Vessel CategoryLength Overall (LOA) (m)
Triangular Distribution Parameters
Lower LimitModeUpper Limit
14590140
2140170190
3190230270
4270320400
Table A2. Parameters applied to triangular distribution generating TEU quantities on a vessel.
Table A2. Parameters applied to triangular distribution generating TEU quantities on a vessel.
TEU CategoryTEU Quantity
Triangular Distribution Parameters
Lower LimitModeUpper Limit
1150350500
25008501100
3110014501900
41900290015,000
Table A3. Parameters applied to triangular distribution generating QC intensity based on the generated vessels and TEU quantities.
Table A3. Parameters applied to triangular distribution generating QC intensity based on the generated vessels and TEU quantities.
Vessel CategoryTEU CategoryQC Intensity
Triangular Distribution Parameters
Lower LimitModeUpper Limit
111.01.062.0
121.01.583.0
131.42.142.9
142.72.702.7
211.01.262.9
221.01.863.1
231.22.313.7
241.62.563.6
311.01.442.7
321.02.003.7
331.32.624.7
341.53.344.9
411.01.453.3
421.02.084.3
431.42.884.5
442.13.986.7

Appendix B. Descriptive Statistics for Two-factor ANOVA

This appendix displays the sample statistics applied to two-factor ANOVA. The 200 time-points of vessel departures are randomly selected, and the corresponding dataset of QC intensity, storage space requirement, and gate congestion is collected from simulation results per setting. Note that a setting is a combination of the trans-shipment TEU rate and the transferring rate. Table A4, Table A5 and Table A6 shows the descriptive statistics of the collected samples.
Table A4. Sample statistics for QC intensity for collected 200 simulation results.
Table A4. Sample statistics for QC intensity for collected 200 simulation results.
TransshipmentMean Variance
TransferringHigher (75%) Lower (60%) Overall Higher (75%) Lower (60%) Overall
Upper (90%)7.437.877.650.060.150.15
Higher (85%)7.267.967.610.160.140.27
Norm (80%)8.448.098.270.220.240.26
Lower (75%)7.447.607.520.100.320.22
Overall7.647.887.760.350.250.31
Table A5. Sample statistics for storage space requirement for collected 200 simulation results.
Table A5. Sample statistics for storage space requirement for collected 200 simulation results.
TransshipmentMean Variance
TransferringHigher (75%) Lower (60%) Overall Higher (75%) Lower (60%) Overall
Upper (90%)72,733.7962,976.3167,855.051,047,840.23456,104.7024,611,871.06
Higher (85%)81,190.6491,660.9486,425.791,439,681.152,397,562.5229,389,265.57
Norm (80%)106,006.40116,520.41111,263.401,606,384.97805,290.1328,908,207.94
Lower (75%)135,075.54132,135.77133,605.652,270,371.301,097,727.513,845,800.66
Overall98,751.59100,823.3699,787.47591,666,169.87687,846,471.42640,429,945.84
Table A6. Sample statistics for gate congestion for collected 200 simulation results.
Table A6. Sample statistics for gate congestion for collected 200 simulation results.
TransshipmentMean Variance
TransferringHigher (75%) Lower (60%) Overall Higher (75%) Lower (60%) Overall
Upper (90%)1.731.551.640.070.050.07
Higher (85%)1.611.581.590.050.040.05
Norm (80%)1.211.421.310.020.050.04
Lower (75%)1.231.661.440.030.040.08
Overall1.441.551.500.090.050.08

References

  1. Ng, A.K.Y.; Liu, J.J. Port-Focal Logistics and Global Supply Chains; Palgrave Macmillan: New York, NY, USA, 2014; pp. 77–113. [Google Scholar]
  2. Jin, X.; Kim, K.H. Collaborative Inter-Terminal Transportation of Containers. Ind. Eng. Manag. Syst. 2018, 17, 407–416. [Google Scholar] [CrossRef]
  3. Maritime and Port Authority of Singapore. Sustaining a Future-Ready Maritime Singapore. Maritime and Port Authority of Singapore Annual Report 2018. Available online: https://www.mpa.gov.sg/web/portal/home/about-mpa/annual-report (accessed on 3 August 2021).
  4. Hendriks, M.P.M.; Armbruster, D.; Laumanns, M.; Lefeber, E.; Udding, J.T. Strategic allocation of cyclically calling vessels for multi-terminal container operators. Flex. Serv. Manuf. J. 2011, 24, 248–273. [Google Scholar] [CrossRef] [Green Version]
  5. Lee, D.-H.; Jin, J.G.; Chen, J.H. Terminal and yard allocation problem for a container transshipment hub with multiple terminals. Transp. Res. Part E Logist. Transp. Rev. 2012, 48, 516–528. [Google Scholar] [CrossRef]
  6. Heilig, L.; Voß, S. Inter-terminal transportation: An annotated bibliography and research agenda. Flex. Serv. Manuf. J. 2016, 29, 35–63. [Google Scholar] [CrossRef]
  7. Mishra, N.; Roy, D.; van Ommeren, J.-K. A Stochastic Model for Interterminal Container Transportation. Transp. Sci. 2017, 51, 67–87. [Google Scholar] [CrossRef]
  8. Gharehgozli, A.H.; de Koster, R.; Jansen, R. Collaborative solutions for inter terminal transport. Int. J. Prod. Res. 2016, 55, 6527–6546. [Google Scholar] [CrossRef]
  9. Heilig, L.; Lalla-Ruiz, E.; Voß, S. Multi-objective inter-terminal truck routing. Transp. Res. Part E Logist. Transp. Rev. 2017, 106, 178–202. [Google Scholar] [CrossRef]
  10. Hu, Q.; Corman, F.; Wiegmans, B.; Lodewijks, G. A tabu search algorithm to solve the integrated planning of container on an inter-terminal network connected with a hinterland rail network. Transp. Res. Part C Emerg. Technol. 2018, 91, 15–36. [Google Scholar] [CrossRef]
  11. Schepler, X.; Balev, S.; Michel, S.; Sanlaville, É. Global planning in a multi-terminal and multi-modal maritime container port. Transp. Res. Part E Logist. Transp. Rev. 2017, 100, 38–62. [Google Scholar] [CrossRef] [Green Version]
  12. Song, D.-W. Regional container port competition and co-operation: The case of Hong Kong and South China. J. Transp. Geogr. 2002, 10, 99–110. [Google Scholar] [CrossRef]
  13. Song, D.-W. Port co-opetition in concept and practice. Marit. Policy Manag. 2003, 30, 29–44. [Google Scholar] [CrossRef]
  14. Wang, K.; Ng, A.K.Y.; Lam, J.S.L.; Fu, X. Cooperation or competition? Factors and conditions affecting regional port governance in South China. Marit. Econ. Logist. 2012, 14, 386–408. [Google Scholar] [CrossRef]
  15. Kaselimi, E.N.; Notteboom, T.E.; De Borger, B. A game theoretical approach to competition between multi-user terminals: The impact of dedicated terminals. Marit. Policy Manag. 2011, 38, 395–414. [Google Scholar] [CrossRef]
  16. Yu, M.; Lee, C.-Y.; Wang, J.J. The regional port competition with different terminal competition intensity. Flex. Serv. Manuf. J. 2016, 29, 659–688. [Google Scholar] [CrossRef]
  17. Wang, K.; Zhang, A. Climate change, natural disasters and adaptation investments: Inter- and intra-port competition and cooperation. Transp. Res. Part B Methodol. 2018, 117, 158–189. [Google Scholar] [CrossRef]
  18. Tenhiälä, A. Contingency theory of capacity planning: The link between process types and planning methods. J. Oper. Manag. 2010, 29, 65–77. [Google Scholar] [CrossRef]
  19. Jacobs, F.R.; Berry, W.L.; Whybark, D.C.; Vollmann, T.E. Manufacturing Planning and Control for Supply Chain Management; McGraw-Hill: New York, NY, USA, 2011; pp. 215–316. [Google Scholar]
  20. Fogarty, D.W.; Blackstone, J.H.; Hoffmann, T.R. Production and Inventory Management, 2nd ed.; South-Western College Pub: Cincinnati, OH, USA, 1991. [Google Scholar]
  21. Montgomery, D.C. Statistical Quality Control, 7th ed.; Wiley: Hoboken, NJ, USA, 2012. [Google Scholar]
  22. Ellis, P.D. The Essential Guide to Effect Sizes: Statistical Power, Meta-Analysis, and the Interpretation of Research Results; Cambridge University Press: Cambridge, UK, 2010; pp. 3–44. [Google Scholar]
Figure 1. Simplified illustration of resource profiles for a mainline vessel.
Figure 1. Simplified illustration of resource profiles for a mainline vessel.
Applsci 11 09156 g001
Figure 2. Simplified illustration of resource profiles for unloading (left) and loading (right) operations for a feeder vessel.
Figure 2. Simplified illustration of resource profiles for unloading (left) and loading (right) operations for a feeder vessel.
Applsci 11 09156 g002
Figure 3. Illustrated TEU workload distributions for storage space requirement.
Figure 3. Illustrated TEU workload distributions for storage space requirement.
Applsci 11 09156 g003
Figure 4. Event entity diagram for resource profile simulation.
Figure 4. Event entity diagram for resource profile simulation.
Applsci 11 09156 g004
Figure 5. X ¯ R charts over time-shifts (x-axis) for QC intensity (y-axis) when trans-shipment and nontrans-shipment (inbound and outbound) TEUs are 80% and 20%, respectively, and 60% of arrived TEUs are transferred to other terminals.
Figure 5. X ¯ R charts over time-shifts (x-axis) for QC intensity (y-axis) when trans-shipment and nontrans-shipment (inbound and outbound) TEUs are 80% and 20%, respectively, and 60% of arrived TEUs are transferred to other terminals.
Applsci 11 09156 g005
Figure 6. X ¯ R charts over time-shifts (x-axis) for storage space requirement (y-axis) when trans-shipment and nontrans-shipment (inbound and outbound) TEUs are 80% and 20%, respectively, and 60% of arrived TEUs are transferred to other terminals.
Figure 6. X ¯ R charts over time-shifts (x-axis) for storage space requirement (y-axis) when trans-shipment and nontrans-shipment (inbound and outbound) TEUs are 80% and 20%, respectively, and 60% of arrived TEUs are transferred to other terminals.
Applsci 11 09156 g006
Figure 7. X ¯ R charts over time-shifts (x-axis) for gate congestion (truck queue) (y-axis) when trans-shipment and nontrans-shipment (inbound and outbound) TEUs are 80% and 20%, respectively, and 60% of arrived TEUs are transferred to other terminals.
Figure 7. X ¯ R charts over time-shifts (x-axis) for gate congestion (truck queue) (y-axis) when trans-shipment and nontrans-shipment (inbound and outbound) TEUs are 80% and 20%, respectively, and 60% of arrived TEUs are transferred to other terminals.
Applsci 11 09156 g007
Figure 8. X ¯ R charts over time-shifts (x-axis) for QC intensity (y-axis) when trans-shipment and nontrans-shipment (inbound and outbound) TEUs are 80% and 20%, respectively, and 75% of arrived TEUs are transferred to other terminals.
Figure 8. X ¯ R charts over time-shifts (x-axis) for QC intensity (y-axis) when trans-shipment and nontrans-shipment (inbound and outbound) TEUs are 80% and 20%, respectively, and 75% of arrived TEUs are transferred to other terminals.
Applsci 11 09156 g008
Figure 9. X ¯ R charts over time-shifts (x-axis) for storage space requirement (y-axis) when trans-shipment and nontrans-shipment (inbound and outbound) TEUs are 80% and 20%, respectively, and 75% of arrived TEUs are transferred to other terminals.
Figure 9. X ¯ R charts over time-shifts (x-axis) for storage space requirement (y-axis) when trans-shipment and nontrans-shipment (inbound and outbound) TEUs are 80% and 20%, respectively, and 75% of arrived TEUs are transferred to other terminals.
Applsci 11 09156 g009
Figure 10. X ¯ R charts over time-shifts (x-axis) for gate congestion (truck queue) (y-axis) when trans-shipment and nontrans-shipment (inbound and outbound) TEUs are 80% and 20%, respectively, and 75% of arrived TEUs are transferred to other terminals.
Figure 10. X ¯ R charts over time-shifts (x-axis) for gate congestion (truck queue) (y-axis) when trans-shipment and nontrans-shipment (inbound and outbound) TEUs are 80% and 20%, respectively, and 75% of arrived TEUs are transferred to other terminals.
Applsci 11 09156 g010
Figure 11. Changes of QC intensity requirement for trans-shipment and transferring levels.
Figure 11. Changes of QC intensity requirement for trans-shipment and transferring levels.
Applsci 11 09156 g011
Figure 12. Changes of storage space requirement for trans-shipment and transferring levels.
Figure 12. Changes of storage space requirement for trans-shipment and transferring levels.
Applsci 11 09156 g012
Figure 13. Changes of gate congestion requirement for trans-shipment and transferring levels.
Figure 13. Changes of gate congestion requirement for trans-shipment and transferring levels.
Applsci 11 09156 g013
Table 1. Control limits of X ¯ R charts for QC intensity, storage space requirement, and gate congestion corresponding to Figure 5, Figure 6 and Figure 7.
Table 1. Control limits of X ¯ R charts for QC intensity, storage space requirement, and gate congestion corresponding to Figure 5, Figure 6 and Figure 7.
X ¯ Charts R Charts
UCLCLLCLUCLCLLCL
QC intensity9.788.126.476.062.870
Storage space requirement125,006.81117,109.64109,212.5028,933.4813,686.600
Gate congestion2.321.560.792.801.320
Table 2. Control limits of X ¯ R charts for QC intensity, storage space requirement, and gate congestion corresponding to Figure 8, Figure 9 and Figure 10.
Table 2. Control limits of X ¯ R charts for QC intensity, storage space requirement, and gate congestion corresponding to Figure 8, Figure 9 and Figure 10.
X ¯ Charts R Charts
UCLCLLCLUCLCLLCL
QC intensity9.717.926.135.653.100
Storage space requirement127,060.92107,120.7887,180.6473,056.2534,558.300
Gate congestion2.521.690.863.031.430
Table 3. Two-factor ANOVA results for QC intensity in 5% significance level.
Table 3. Two-factor ANOVA results for QC intensity in 5% significance level.
Sum of Square (SS)Degree of Freedom (df)Mean Sum of Square (MS)F-Ratiop-ValueSignificance (Sig.)
Trans-shipment140.07346.69266.663.22 × 10−140yes
Transferring23.03123.03131.532.60 × 10−29yes
Interaction61.14320.38116.393.82 × 10−68yes
Error278.751,5920.18
Total502.991,5990.31
Table 4. Two-factor ANOVA results for storage space requirement in 5% significance level.
Table 4. Two-factor ANOVA results for storage space requirement in 5% significance level.
SSdfMSF-Ratiop-ValueSig.
Trans-shipment989,432,180,454.493329,810,726,818.16237,253.370yes
Transferring1,716,883,011.1511,716,883,011.151,235.069.21 × 10−201yes
Interaction30,685,348,400.53310,228,449,466.847,357.960yes
Error2,213,071,537.211,5921,390,120.31
Total1,024,047,483,403.381,599640,429,945.84
Table 5. Two-factor ANOVA results for gate congestion in 5% significance level.
Table 5. Two-factor ANOVA results for gate congestion in 5% significance level.
SSdfMSF-Ratiop-ValueSig.
Trans-shipment26.053.008.68195.784.58 × 10−108yes
Transferring4.621.004.62104.249.50 × 10−24yes
Interaction21.743.007.25163.442.35 × 10−92yes
Error70.601,592.000.04
Total123.011,599.000.08
Table 6. Significance test among trans-shipment levels for QC intensity, storage space requirement, and gate congestion using contrasts at 5% significance level.
Table 6. Significance test among trans-shipment levels for QC intensity, storage space requirement, and gate congestion using contrasts at 5% significance level.
Higher (85%)Norm (80%)Lower (75%)
QC intensityUpper (90%)noyesyes
Higher (85%)-yesyes
Norm (80%)--yes
Storage space requirementUpper (90%)yesyesyes
Higher (85%)-yesyes
Norm (80%)--yes
Gate congestionUpper (90%)yesyesyes
Higher (85%)-yesyes
Norm (80%)--yes
Table 7. Significance test between individual trans-shipment levels under transferring levels for QC intensity, storage space requirement, and gate congestion using contrasts in 5% significance level.
Table 7. Significance test between individual trans-shipment levels under transferring levels for QC intensity, storage space requirement, and gate congestion using contrasts in 5% significance level.
Higher (75%): TransferringLower (60%): Transferring
Higher (85%)Norm (80%)Lower (75%)Higher (85%)Norm (80%)Lower (75%)
QC intensityUpper (90%)yesyesnoyesyesyes
Higher (85%)-yesyes-yesyes
Norm (80%)--yes--yes
Storage space requirementUpper (90%)yesyesyesyesyesyes
Higher (85%) yesyes-yesyes
Norm (80%)--yes--yes
Gate congestionUpper (90%)yesyesyesnoyesyes
Higher (85%)-yesyes-yesyes
Norm (80%)--no--yes
Table 8. Cohen’s d and Pearson’s r between individual trans-shipment levels under transferring levels for QC intensity, storage space requirement, and gate congestion using contrasts in 5% significance level.
Table 8. Cohen’s d and Pearson’s r between individual trans-shipment levels under transferring levels for QC intensity, storage space requirement, and gate congestion using contrasts in 5% significance level.
Higher (75%): TransferringLower (60%): Transferring
Higher (85%)Norm (80%)Lower (75%)Higher (85%)Norm (80%)Lower (75%)
QC intensityUpper (90%)d = 040
r = 0.10
d = 2.43
r = 0.52
d = 0.03
r = 0.01
d = 0.22
r = 0.05
d = 0.52
r = 0.13
d = 0.65
r = 0.16
Higher (85%)-d = 2.83
r = 0.58
d = 0.42
r = 0.11
-d = 0.31
r = 0.08
d = 0.87
r = 0.21
Norm (80%)--d = 2.40
r = 0.52
--d = 1.18
r = 0.28
Storage space requirementUpper (90%)d = 7.17
r = 0.87
d = 28.22
r = 0.99
d = 52.88
r = 1.00
d = 24.33
r = 0.99
d = 45.41
r = 1.00
d = 58.66
r = 1.00
Higher (85%) d = 21.05
r = 0.98
d = 45.70
r = 1.00
-d = 21.08
r = 0.98
d = 34.33
r = 0.99
Norm (80%)--d = 24.66
r = 0.99
--d = 13.24
r = 0.96
Gate congestionUpper (90%)d = 0.56
r = 0.14
d = 2.46
r = 0.52
d = 2.36
r = 0.51
d = 0.15
r = 0.04
d = 0.61
r = 0.15
d = 0.54
r = 0.13
Higher (85%)-d = 1.90
r = 0.43
d = 1.81
r = 0.41
-d = 0.79
r = 0.19
d = 0.39
r = 0.10
Norm (80%)--d = 0.09
r = 0.02
--d = 1.15
r = 0.28
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, B.K.; Low, J.M.W. Resource Capacity Requirement for Multi-Terminal Cooperation in Container Ports. Appl. Sci. 2021, 11, 9156. https://doi.org/10.3390/app11199156

AMA Style

Lee BK, Low JMW. Resource Capacity Requirement for Multi-Terminal Cooperation in Container Ports. Applied Sciences. 2021; 11(19):9156. https://doi.org/10.3390/app11199156

Chicago/Turabian Style

Lee, Byung Kwon, and Joyce M. W. Low. 2021. "Resource Capacity Requirement for Multi-Terminal Cooperation in Container Ports" Applied Sciences 11, no. 19: 9156. https://doi.org/10.3390/app11199156

APA Style

Lee, B. K., & Low, J. M. W. (2021). Resource Capacity Requirement for Multi-Terminal Cooperation in Container Ports. Applied Sciences, 11(19), 9156. https://doi.org/10.3390/app11199156

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

Article Metrics

Back to TopTop