1. Introduction
Liner shipping handles over 80% of the transportation of finished products in world trade [
1]. Liner vessels operate container transportation services along fixed routes, calling at designated ports in a predetermined sequence, adhering to a published schedule, and charging relatively fixed freight rates. Typically, liner companies announce their schedules for various regions 3–6 months in advance [
2]. Shippers, such as global manufacturers, then plan their production and transportation schedules accordingly. However, recent years have seen a significant escalation in the frequency and severity of port congestion. For instance, the ports of Los Angeles and Long Beach in the United States have been grappling with prolonged and severe congestion since October 2020, with 109 vessels backlogged as of 9 January 2022 [
3]. This discrepancy between planned and actual vessel arrivals and departures disrupts the delivery timelines of maritime shipments and undermines the reliability of vessel schedules. Consequently, unreliable schedules compel shippers to maintain substantial safety stocks, hindering the implementation of efficient just-in-time production plans and posing challenges to the security and stability of global supply chains [
4].
Liner schedule punctuality reflects a liner company’s capability in route planning, port cooperation, and risk management, serving as a crucial factor for shippers when selecting maritime service providers. Notteboom [
5] investigated the primary reasons for low liner punctuality, finding that over 93.8% of schedule delays were related to port operations, with congestion causing 65.5% of the unexpected waiting time before cargo handling and port infrastructure capacity limitations accounting for 20.6%. Scholars [
6,
7] have extensively studied port congestion using methods like stochastic programming, robust optimization, and variational inequality model. However, since late 2023, events such as the significant container backlog at ports in Durban and Cape Town, drought-induced low water levels in the Panama Canal, and geopolitical conflicts in the Red Sea have led to frequent disruptions, maintaining the severity of global port congestion and undermining container liner schedule punctuality. The global mainline punctuality index released by the Shanghai Shipping Exchange in February 2024 stood at a mere 32.48%, plummeting by 16.16% compared to February 2023. Addressing the real challenges of port congestion while considering the uncertainties of vessel waiting and handling times in ports (where waiting time refers to the time ships spend at anchor after arrival, and handling time refers to the time spent handling containers by terminal cranes), as well as designing reliable liner schedules, presents a significant opportunity for liner companies to attract customers and seize market share.
On the one hand, port congestion primarily affects the waiting time for vessels before entering berths, known as the in-port waiting time. However, the existing literature often treats in-port waiting time as a uniform random variable [
8,
9,
10,
11]. In reality, vessel berthing patterns and transit times vary significantly across different ports. Queueing theory, recognized as a crucial analytical tool for congestion studies, has found applications in optimizing berth allocation and describing port performance, proving to be a highly effective method for quantifying port congestion. It allows for the estimation of critical parameters such as average vessel waiting time, queue length, and average number of vessels in port. Therefore, leveraging queueing theory models can effectively address the uncertainty of vessel in-port waiting time.
On the other hand, shipping companies and ports have also made collaborative efforts to address port congestion [
12,
13,
14]. Before opening a new route, shipping companies engage in information sharing with port operators, including vessel details (especially vessel length and draft) and the range of container handling volumes. Port operators, considering port throughput capacity and expected handling requirements, share available container handling time windows and average handling efficiency data with liner companies and offer multiple options for container handling efficiency at a higher-than-normal level, along with additional charges for this “VIP” service. Then, shipping companies predict in-port waiting times due to port congestion, determine vessel arrival and departure times at each port, and select suitable container handling efficiency based on a comprehensive evaluation of their schedules, available port operation time windows, and container handling plans, thus determining the vessel’s in-port operation time [
15]. This mechanism for selecting container handling efficiency is a win-win strategy, enhancing the reliability and stability of shipping company schedules, while also optimizing port revenue and reducing congestion issues without compromising port operator interests. However, this will also make the vessel handling time in port uncertain, thus increasing the difficulty of the design of the shipping schedule.
Therefore, this study addresses the liner schedule design in the context of port congestion. It utilizes queueing theory to predict and quantify vessel waiting times in port, integrating vessel schedules, cargo handling volumes, and available port operation time windows; proposing a container handling efficiency selection mechanism for inbound vessels; and determining vessel handling times in port. By jointly considering the uncertainty of vessel waiting and handling times, a mixed-integer nonlinear programming model is developed, aiming to minimize the total cost of shipping services. The model undergoes linearization and is solved using CPLEX to devise a reliable shipping schedule. Finally, an experiment analysis is conducted on a real-world Asia-to-Mediterranean shipping route. Our study makes the following contributions:
(1) This paper addresses the issue of liner schedule design, considering the uncertainties of vessel waiting and handling times caused by port congestion. It determines the vessel deployment quantities on the route, vessel speeds on segments, and arrival and departure times at each call port. This comprehensive framework for handling port congestion can be extended to other types of emergencies, such as disruption.
(2) In maritime logistics, the existing literature often treats vessel time in port as a uniform random variable, overlooking significant differences in vessel berthing and container handling across different ports. This paper treats vessel operations upon arrival as a queueing system and employs queueing theory models to predict and quantify vessel waiting times in port. In fact, this model can also be applied to terminal operations and truck-to-port transfers.
(3) We propose a port handling efficiency selection mechanism. Numerical simulations confirm that this information-sharing mechanism not only enhances the flexibility and reliability of liner shipping services but also increases profitability for both liner companies and ports. This mechanism offers a new approach to collaboration and information sharing in maritime logistics.
(4) A simulation analysis on a real Asia-to-Mediterranean liner shipping route demonstrates that extreme weather events and geopolitical conflicts can cause severe congestion at certain ports. Liner companies must adjust vessel schedule promptly in response. Additionally, such events can impact the marine fuel market, necessitating strategies such as increasing vessel operations and reducing vessel speeds under high fuel prices.
The remainder of this paper is organized as follows.
Section 2 presents related studies on liner schedule design, application of queueing theory in maritime operations, and port handling efficiency selection mechanism.
Section 3 develops a mixed-integer nonlinear programming model.
Section 4 linearizes the model and utilizes commercial solver CPLEX to solve.
Section 5 describes computational experiments, illustrates the result analysis, and derives managerial insights.
Section 6 provides conclusions.
3. Model Formulation
3.1. Problem Description
The current global port congestion significantly impacts liner schedule punctuality, necessitating consideration of port congestion uncertainty in schedule design. Vessels adhere to predetermined port call schedules, sequentially visiting port during each voyage, and they need to arrive at the port within designated time windows for operations. However, in the context of port congestion, vessel time in port is affected by waiting for berths and container handling, thus introducing significant uncertainty. To address the uncertainty of waiting for berths, we utilize queueing theory models to predict and quantify vessel berth times. By considering factors such as the daily vessel arrivals number, ; available berths number, ; and average number of vessels served per berth, , we calculate the queue service intensity, , at port . Subsequently, we determine the actual wait times, , for vessels at port by comparing the relationship of queue times, ; arrival times, ; and available service windows.
To address the uncertainty of handling times, we propose a container handling efficiency selection mechanism. This allows vessels to choose one of several container handling efficiency options, , at port based on their specific needs. We introduce a decision variable, , to represent vessel choices, determining container handling efficiency, , at port . Subsequently, considering container quantity, , and unit port service costs, , we calculate the total container handling fees that vessels need to pay at each port and determine vessel berth times, , accordingly.
After addressing vessel berth and handling uncertainties, vessels decide their sailing speeds, , based on segment distances, ; departure times, , from previous ports; and vessel speed ranges, , calculating sailing times, , for each segment. Fuel costs, determined by unit fuel costs, ; fuel coefficients, ; unit fuel consumption, ; and segment distances, , is incurred during voyages, in addition to and container inventory costs, , based on container quantities, , and sailing times, , per segment. If vessels fail to complete operations within specified time windows, we impose a penalty cost based on delay duration, , and unit penalty cost, , aiming to enhance liner service quality.
Last, shipping companies need to deploy vessels on the route, each incurring fixed operational costs, , per voyage to meet weekly service requirements. The number of vessels available for deployment has an upper bound, . In summary, this study establishes an optimization model minimizing total costs by considering vessel operational, fuel consumption, container handling, inventory, and delay penalty costs.
3.2. Assumptions
(1) Vessel berth time is primarily determined by waiting and container handling, without considering time spent on other activities [
26].
(2) The sequence and ports of call along the route are known, with liner companies offering weekly service frequencies [
43].
(3) The study focuses solely on optimizing schedules for a single route, excluding transshipment issues [
25].
(4) Similar vessel types are deployed on the route, sharing identical technical characteristics, and only the fuel consumption of main engines is considered [
25].
(5) The import process of vessel arrivals at ports follow a Poisson distribution [
44,
45], while berth occupancy times adhere to an exponential distribution [
46].
(6) It is assumed that vessels adhere to a first-come, first-served principle at ports. If vessels arrive before the available berth window, they must wait at anchor. Given predetermined port calls at a strategic level, liner companies exhibit patience in queueing without departing or diverting mid-route [
47].
3.3. Symbol Specification
Before formulating the model for this problem, we list the notation as follows.
Indices and sets:
: set of liner line call port or line segment; that is, the ith call port corresponds to the ith section, .
: set of container handling efficiency in port , .
Parameters:
: number of ports of call.
: fixed operating costs .
: unit inventory cost of containerized goods .
: unit penalty cost of vessel delay in port .
: unit fuel cost .
: vessel unit service cost when selecting container handling efficiency, , at port .
: container handling efficiency, , at port .
: minimum vessel speed .
: maximum vessel speed .
: distance of segment .
: start time of port ’s available time window .
: end time of port ’s available time window .
: correlation coefficients of fuel consumption function.
: number of containers transported by the vessel in segment .
: number of containers handled at port .
: maximum number of vessels that can be deployed on the route .
: average number of vessels arriving at port per day .
: average number of vessels serviced per day by a berth in port .
: number of berths available in port .
Intermediate variables:
: Unit fuel consumption of the vessel in segment .
: Arrival time of the vessel at port .
: Container handling time of the vessel at port .
: Departure time of the vessel from port .
: Sailing time of the vessel in segment .
: Delay time of the vessel at port .
: Queueing time of the vessel at port .
: Waiting time of the vessel at port .
: Queue service intensity of port .
Variables:
: speed of vessel on segment .
: 0–1 variable, 1 if the vessel selects the hth container handling efficiency, , in port ; otherwise, 0.
: Number of vessels deployed on the route .
Regarding the abbreviations used in this paper, TEU stands for Twenty-Foot Equivalent Unit, USD represents United States Dollar, and n mile denotes Nautical Mile. We use these abbreviations consistently throughout the text.
3.4. Mixed-Integer Nonlinear Programming Model
Based on the above description, the first liner schedule design model (M1) is established with the goal of minimizing the total costs of liner services as follows:
Equation (1) represents the objective function of the model, aiming to minimize the total costs of liner services. The first term pertains to vessel operating costs, which are contingent upon the number of vessels deployed along the route. The second term accounts for fuel costs, determined by vessel speed and distance traveled. The third term reflects port handling costs, influenced by container handling volumes and efficiency. The fourth term represents container inventory costs, primarily linked to the duration of containers on board. Lastly, the fifth term denotes penalties for vessel delays, contingent upon the duration of container delays.
3.4.1. Waiting Time Prediction Model Based on Queueing Theory
Initially, we quantify the queueing time for vessels waiting for berthing in port due to port congestion. When vessels randomly dock at any public berth at a port, their arrivals follow the
multi-server queueing system. As this study involves multiple time points related to port availability windows,
Figure 1 elucidates the interrelationships among these time points.
In
Figure 1, the vessel departs from port
at time
, after traveling for
units of time in segment
(from port
to port
) and arriving at the port
at time
. The start and end times of the available time window at port
are
and
, respectively. As the queueing time for the vessel at port
is uncertain and subject to variation within the green matrix, three scenarios can occur:
(1) The vessel completes queueing before the start of the available time window, , at port . In this case, the waiting time for berthing is .
(2) The vessel completes queueing precisely between the start time, , and the end time, , of the available time window at port (including these two time points). Here, the waiting time for berthing is .
(3) Vessels complete queueing after the end of the available time window at port . Despite being unable to avail port services at this point, there is still an opportunity cost for waiting, and the waiting time for berthing is . Certainly, this scenario corresponds to ; that is, the number of vessels exceeds the port’s service capacity for the day, theoretically resulting in prolonged queueing. For instance, since October 2023, Durban Port on the South African East Coast has experienced extreme weather conditions and equipment malfunctions by the port operator, Transnet. This has led to over 100,000 containers being stranded, with more than 100 container vessels becoming stuck, causing significant delays in liner schedules. Experts estimated that the backlog would not be cleared until February 2024. Opting to skip the port would be the optimal choice for liner companies during such circumstances. However, as these events are extreme emergencies beyond schedule design considerations, this study does not account for them. Hence, we assume that vessels always arrive before the end of the available time window at the port.
The waiting time prediction model based on queueing theory is as follows:
Formulas (2)–(7) quantify the vessel’s berth waiting time at port using queueing theory, with a progressive relationship between them. Equation (2) computes the queue service intensity, , at port . Furthermore, when the service intensity is , Equation (3) calculates the probability, , of no vessels arriving at port . Building upon Equations (3) and (4), we assess the probability, , of the vessel needing to queue at port when vessels arrive and . Equation (5) determines the queue length, , of the vessel waiting for service at port . Equation (6) computes the waiting time, , of the vessel at port . Additionally, Equation (7) calculates the waiting time of the vessel at port . Given the closed-loop nature of liner shipping, this waiting time is contingent upon the departure time from the last port of the previous loop and the sailing time to the first port of the current loop.
3.4.2. Liner Schedule Design Model
For liner enterprises, schedule design involves determining vessel arrival and departure times of vessels at various ports along established routes or shipping networks. Before finalizing these times, liner companies need to gather information about available working hours, container handling capacity, vessel turnaround times at ports, distances and weather conditions along route segments, cargo volume between ports, and constraints related to vessel attributes and quantity on the route. Based on this information, parameters are set according to principles of minimizing costs or optimizing service. With these parameters, round-trip voyage times can be calculated. This refers to the time that vessels spend traveling from the originating port, visiting all ports in sequence, and returning to the originating port, encompassing both sailing and port waiting and operation times. The calculation formulas and constraints of the model are provided below.
Formula (8) indicates that the vessel must select one container handling efficiency scheme at each port. Formula (9) calculates the sailing time for each route segment. Formula (10) computes the fuel consumption per unit distance based on the vessel sailing speed for each segment. Formula (11) ensures that the deployed vessel quantity maintains a weekly port service frequency. The left side represents the total time for a vessel to complete one voyage, including sailing, handling, and waiting times, while the right side is the product of the total hours in a week () and the vessel quantity, , needed to be deployed. Formula (12) calculates the departure time of the vessel from port based on arrival, waiting, and handling times. Formulas (13) and (14) determine the arrival times at the first port and subsequent ports. Formula (15) calculates the delay time for the vessel at each port. Formula (16) computes the handling time for the vessel at each port. Formula (17) restricts the deployed vessel quantity on the route in order to not exceed the maximum available vessel quantity. Formula (18) sets the constraint for vessel sailing speed.
4. Model Solving
The model’s constraints (9) and (10) involve reciprocal and power functions of vessel sailing speed, , making it a mixed-integer nonlinear programming problem. To facilitate solving, the model needs linearization.
For Formula (9), we replace vessel speed, , with its reciprocal, , to linearize the constraint.
Regarding the model, the liner speed is a continuous variable within a certain range, meaning the unit fuel consumption calculated using Formula (10) can be any rational number within that range. Various methods have been proposed in the liner transportation literature to address this nonlinearity.
(1) Enumerative method: assuming constant vessel speeds for route segments, minimizing the total service cost to solve vessel schedule.
(2) Discretization method: discretizing vessel speed reciprocals, estimating fuel consumption for each discretized reciprocal speed, and simplifying vessel schedule into a mixed-integer linear problem.
(3) Dynamic programming method: simplifying vessel schedule into a shortest-path problem in a spatiotemporal network, with time as the horizontal axis (usually in days) and ports as the vertical axis.
(4) Customized method: substituting nonlinear fuel consumption functions with approximate functions (e.g., sets of tangent or secant lines) to simplify vessel schedule into a mixed-integer linear problem.
(5) Second-order cone programming method: transforming the original mixed-integer nonlinear schedule model into a mixed-integer second-order cone programming model [
48].
Among these methods, dynamic programming, discretization, and customized methods usually more effectively approximate nonlinear fuel consumption functions [
49]. This paper focuses on liner schedule design, a static problem determined before the ship’s voyage. This static nature distinguishes it from dynamic scheduling problems that dynamic programming typically addresses. Therefore, dynamic programming is not suitable for our context. On the other hand, discretization can set appropriate speed precision for liner ships based on the shipping company’s expectations and is frequently used in schedule design and ship scheduling. Thus, we employed the discretization method.
According to this method, vessel speed reciprocals are discretized into a finite set of values,
. Let
be the reciprocal vessel speed at the discretization point,
. Then,
denotes the unit fuel consumption when using the reciprocal vessel speed value at the discretization point,
. The degree of discretization,
, increases the precision of approximating the fuel consumption function but also increases the number of variables in the model, potentially leading to longer solution times. Further discussion on this is provided in the
Section 5.2. We introduce a new parameter,
, which equals 1 if the vessel’s fuel consumption value for segment
is estimated using the discretization point,
, and it is 0 otherwise. The impact of different discretization precisions (i.e., speed selection ranges) on the unit fuel consumption is illustrated in
Figure 2.
Figure 2 shows the results for four different speed discretization scenarios. In each scenario, we selected one point with the same x-coordinate, where the reciprocal of the vessel speed,
, is 0.054750. The y-coordinates (unit fuel consumption
) vary due to the different ranges of speed options available. In the first scenario (k = 5), the speed selection range includes only five options, resulting in an approximated unit fuel consumption coefficient,
, of 0.173774. In contrast, in the fourth scenario (k = 20), with twenty speed options, the approximated unit fuel consumption coefficient,
, is 0.172868, yielding a precision difference of 0.524%. Although this difference in precision is not significantly large, it can substantially affect the final fuel costs. For instance, consider the example from
Section 5.1, where the total voyage distance is 20,948 n miles, the vessel speed is assumed to be 20 knots, and the fuel price is 300 USD/ton. The calculated fuel cost difference between the first and fourth scenarios is USD 125,688. Given that the precision of the discretization method significantly impacts the quality and computation time of the final results, we provide a more detailed analysis in
Section 5.4.
At this point, the original nonlinear model (M1) is converted to a mixed-integer linear programming model (M2):
Constraints:
Formulas (2), (3), (6)–(8), and (11)–(17).
Objective function (19) aims to minimize the total operational, fuel, port handling, container inventory, and vessel delay penalty costs. Formula (20) ensures selecting only one discretization point for estimating fuel consumption on each route segment. Formula (21) dictates the selection of discretization points for calculating vessel speed reciprocals on each route segment. Formula (22) computes vessel fuel consumption using the chosen discretization points on each route segment. Formula (23) calculates the sailing time for vessels on each route segment. Formula (24) imposes constraints on the vessel speed reciprocal values.
At this stage, M2 becomes a mixed-integer linear programming model, solvable using commercial solvers like CPLEX.
6. Conclusions
The uncertainty in vessel waiting and port handling times affects the stability of liner schedules. Current research typically treats vessel port times as random variables, failing to accurately predict port congestion based on observable factors, like vessel arrivals and port operations capacity. This study addresses the issue of liner schedule design under port congestion. We employ queueing theory models to describe vessel waiting time uncertainty and propose a container handling efficiency selection mechanism for arriving vessels to determine their port handling time. By jointly considering these two uncertainties, a robust liner schedule design model is established and solved using the CPLEX.
Numerical simulations on an Asia-to-Mediterranean liner route reveal that extreme weather events or geopolitical conflicts may cause severe port congestion, affecting vessel punctuality and requiring timely adjustments to vessel schedules. Additionally, such events impact the international maritime fuel market, prompting liner companies to consider strategies like increasing vessel operations and reducing vessel speed under high fuel prices. The container handling efficiency selection mechanism allows liner companies to flexibly design schedules, while balancing economic costs and service reliability.
Future studies can be conducted in the following areas. (1) Enhanced data accuracy: This study references Dulebenets’ study [
50] for simplifying vessel arrival data, employing queueing theory to predict port congestion and designing a robust liner schedule. However, real-world vessel arrivals are subject to uncertainties caused by events like COVID-19 or the Red Sea crisis, which can disrupt normal patterns. To further enhance the robustness of liner schedules, future work could involve using extensive port historical data to train queueing models for each port and utilizing automatic identification systems to track real-time vessel locations near ports, thus enhancing data accuracy. (2) Cooperative strategies among heterogeneous fleets: By coordinating schedules and operations, fleets can better manage uncertainties and optimize resource utilization, leading to a more resilient and efficient schedule. (3) Collaborative agreements between adjacent terminal operators: Such agreements can facilitate workload sharing and improve overall port efficiency, reducing congestion and delays. (4) Incorporation of carbon emission costs: Incorporating these costs into the scheduling model can help liner companies balance economic efficiency with environmental sustainability, promoting greener shipping practices.