Next Article in Journal
Paleogeomorphology Restoration of Post-Rift Basin: Volcanic Activity and Differential Subsidence Influence in Xihu Sag, East China Sea
Next Article in Special Issue
Redefinition for Fundamental Characteristics of Waterway Traffic Stream Considering Vessel Displacement
Previous Article in Journal
Underwater Acoustic Signal Detection against the Background of Non-Stationary Sea Noise
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization of Integrated Tugboat–Berth–Quay Crane Scheduling in Container Ports Considering Uncertainty in Vessel Arrival Times and Berthing Preferences

1
School of Navigation, Jimei University, Xiamen 361021, China
2
Fujian Shipping Research Institute, Xiamen 361021, China
*
Authors to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2024, 12(9), 1541; https://doi.org/10.3390/jmse12091541
Submission received: 6 August 2024 / Revised: 23 August 2024 / Accepted: 25 August 2024 / Published: 4 September 2024
(This article belongs to the Special Issue Resilience and Capacity of Waterway Transportation)

Abstract

:
Influenced by the dynamics of supply and demand, the demand for maritime transport has been increasing annually, putting significant pressure on container ports. To alleviate this pressure, a new mixed-integer programming model for the integrated scheduling of tugboats, berths, and quay cranes has been established. This model considers the uncertainties in vessel arrival times, vessel berthing preferences, time-varying quay crane availability, and the constraint that quay cranes cannot cross each other. The objective is to minimize the total costs including fuel consumption during port stays, delays and waiting times for berthing and departure, berthing deviation costs, tugboat assistance costs, and quay crane handling costs. To obtain high-quality solutions, an adaptive large neighborhood search (ALNS) algorithm was employed to solve the model. The algorithm incorporated five destruction operators and five repair operators that were specifically designed to enhance the solution accuracy and efficiency for the integrated scheduling problem. Several case studies of varying scales, based on a port in China, were used to validate the effectiveness of the proposed model and algorithm. The experimental results demonstrate the model’s validity and show that the ALNS algorithm designed for the integrated scheduling problem outperformed CPLEX and other algorithms in terms of the accuracy and efficiency. Finally, a sensitivity analysis of the key parameters provides recommendations for the integrated scheduling of tugboats, berths, and quay cranes, offering valuable insights for port operations.

1. Introduction

With the rapid development of the global economy and the continuous increase in international trade volume, ports, as crucial nodes in the logistics and transportation system, play an essential role in the efficiency and service quality. In port operations, the scheduling of berths, quay cranes, and tugboats is a critical factor affecting port operational efficiency. Therefore, optimizing the scheduling of these resources to improve port throughput capacity and service levels has become a focal point of interest for both academics and practitioners.
Berth allocation mainly involves the mooring arrangement of vessels at the port. An optimized berth allocation can effectively reduce the vessel waiting times and improve the overall utilization rate of the port. Quay crane scheduling relates to the loading and unloading efficiency of cargo between the vessel and the dock, which is essential to ensure the quick turnaround of port operations. Tugboat scheduling affects not only the speed of vessel berthing and departure, but also the safety of port operations. In-depth research on the scheduling of port resources has been conducted by some scholars, but most of these studies are limited to the optimization of a single resource, lacking systematic research on the integrated scheduling of berths, quay cranes, and tugboats. Single-resource optimization cannot meet the needs of modern port efficient operations; integrated scheduling is a key approach to enhancing the overall port efficiency. Therefore, research on the integrated scheduling of berths, quay cranes, and tugboats has significant theoretical and practical implications in addressing bottlenecks in actual port operations. This work seeks to fill this research gap by constructing a comprehensive integrated scheduling model and using optimization algorithms to improve the collaborative utilization efficiency of port resources and provide theoretical support and practical guidance for port operations, facilitating the development toward intelligent and refined port management.
To improve the efficiency of port quayside operations, a mixed-integer programming model for the integrated scheduling of tugboats, berths, and quay cranes is proposed in this paper. The objective was to minimize the total costs including fuel consumption costs during the vessel’s stay in port, penalties for waiting and delayed departure, berth deviation penalties, tugboat fuel costs, and quay crane loading and unloading costs. Subsequently, the integrated scheduling model was linearized, converting it into a linear mixed-integer programming model. The CPLEX solver and the adaptive large neighborhood search (ALNS) algorithm were then employed to solve it. The effectiveness of the model and method was verified in the experimental results section. The contributions of this paper are as follows:
(1)
A linear integer programming model for the integrated scheduling of tugboats, berths, and quay cranes was constructed, which can yield the optimal berth allocation scheme using CPLEX for small-scale instances.
(2)
An adaptive large neighborhood search algorithm tailored to the characteristics of the integrated scheduling problem of tugboats, berths, and quay cranes was developed. For large-scale instances, the adaptive large neighborhood search algorithm can obtain near-optimal solutions in a relatively short time.
(3)
Practical constraints such as the uncertainty of vessel arrival times, vessel berthing preferences, the non-overlapping of quay cranes, and partially time-varying constraints of quay cranes were considered, thereby enhancing the applicability of the integrated scheduling scheme.
(4)
A continuous berth allocation layout was implemented, improving the utilization of the port quay space.
(5)
The model’s validity and the algorithm’s efficiency were verified through data experiments.
The rest of this paper is organized as follows. Section 2 reviews the related literature on port berth, quay crane, and tugboat scheduling, summarizing the deficiencies of the existing research. Section 3 constructs the linear integer programming model for the integrated scheduling of tugboats, berths, and quay cranes. Section 4 designs the solution method for the model. Section 5 presents a case analysis to verify the effectiveness of the model and algorithm. Section 6 summarizes the paper.

2. Literature Review

The literature related to the scheduling of tugboats, berths, and quay cranes is summarized in this section. It provides a summary of the relevant studies in four aspects: berth allocation, integrated scheduling of berths and quay cranes, tugboat scheduling, and integrated scheduling of berths, tugboats, and quay cranes.

2.1. Research on Berth Allocation

Berth allocation can be categorized into three types: discrete, hybrid, and continuous. In the research on discrete berth allocation problems, Imai et al. [1] aimed to minimize the waiting and handling times for all vessels by constructing a nonlinear integer programming model and designing a Lagrangian relaxation-based heuristic algorithm to solve it. Prencipe et al. [2] focused on minimizing the in-port time of vessels by developing a linear mixed-integer programming model and designing a swarm optimization algorithm to solve it. Iris et al. [3] developed a set-packing model for the strategic berth allocation problem at large container ports and demonstrated the model’s excellence through the analysis of numerous examples of varying scales.
In hybrid berth allocation, a large vessel can occupy multiple berths simultaneously if there is sufficient space, and several small vessels can also dock at one berth concurrently. Umang et al. [4] constructed a mixed-integer linear programming model with the objective of minimizing the in-port time of all vessels and designed a squeaky wheel optimization algorithm for its solution. Kordić et al. [5] aimed to minimize the total in-port costs of all vessels by developing a linear mixed-integer programming model and designing an exact algorithm based on the precipitation method to solve it. Venturini et al. [6] were the first to propose the concept of the multi-port berth allocation problem, considering both vessel speed optimization and emissions. They developed a MILP model for this problem and demonstrated the model’s effectiveness and applicability through extensive data experiments.
For continuous berth allocation problems, Ernst et al. [7] studied berth allocation in export ports and established a two-stage linear mixed-integer programming model, introducing two efficient inequalities to enhance the solvability of the model. Cheimanoff et al. [8] also researched berth allocation in export ports, creating a mixed-integer linear programming model, designing a variable neighborhood search algorithm to solve it, and employing machine learning methods to adjust the algorithm parameters. Martin-Iradi et al. [9] investigated continuous berth allocation problems in multi-port scenarios, developing a mixed-integer linear programming model for this problem.

2.2. Integrated Scheduling of Berths and Quay Cranes

The integrated scheduling of berths and quay cranes also requires the consideration of the division of the quay line, with related research primarily focusing on discrete and continuous types. Xiang et al. [10] studied the discrete berth and quay crane integrated scheduling problem, addressing uncertainties through the reservation of time slots. Wang et al. [11] researched the continuous berth and quay crane integrated scheduling problem under various uncertainties. In the continuous division, the quay line is divided into multiple smaller berths, and vessels need to occupy several berths when docking, thereby improving the utilization of the quay line.
Research on the integrated scheduling of berths and quay cranes involves numerous variables and constraints, making it challenging to solve all variables simultaneously. Therefore, related research often adopts a phased approach. Ji et al. [12] solved the berth and quay crane integrated scheduling problem in two stages: the first stage determines the berthing position, berthing time, and the number of quay cranes assigned to the vessel, while the second stage allocates specific quay cranes to the corresponding vessels. Liu et al. [13] used a rolling horizon strategy to solve the integrated scheduling problem of sea-rail intermodal ports in two stages: the first stage determines the berth occupied by the vessel and the number of quay cranes assigned to each vessel, and the second stage uses a heuristic algorithm to solve the time variables. Chargui et al. [14] employed a decomposition algorithm to divide the problem into a master problem and multiple scenario sub-problems: the first stage solves the deterministic master problem, and the second stage decides whether to solve the sub-problems based on the results of the master problem.
The research on the integrated scheduling of berths and quay cranes must consider time-invariant or time-variant constraints on the number of quay cranes. The time-invariant constraint means that the number of quay cranes assigned to the vessel remains unchanged during its berthed period; the time-variant constraint means that the number of quay cranes assigned to the vessel can change during its berthed period. Yu et al. [15] studied the integrated scheduling of berths and quay cranes under time-invariant constraints, aiming to balance the economic and environmental benefits. Li et al. [16] studied the integrated berth and quay crane scheduling problem by considering quay crane maintenance issues and employed a time-varying strategy for the number of quay cranes. Iris et al. [17] compared the integrated scheduling of berths and quay cranes under time-variant and time-invariant strategies and found that the scheduling efficiency of time-variant strategies was higher than that of time-invariant strategies due to the greater flexibility and efficiency of fully time-variant strategies, which could improve the utilization of quay cranes. When considering time-invariant constraints, the scheduling efficiency is generally lower than when considering time-variant constraints. Time-variant strategies can be fully time-variant or partially time-variant. The modeling and solving process for fully time-variant strategies is more complex, while the complexity of partially time-variant strategies is not significantly higher than that of time-invariant strategies.
When vessels enter and leave port channels, it is necessary to consider whether the channel and berth depths meet the vessel’s draft requirements. Due to tidal factors, the depths of the port berths and channels change periodically, significantly affecting vessel entry, departure, and berthing. Malekahmadi et al. [18] studied the impact of tidal factors on the berth depth, emphasizing that if the berth depth does not meet the vessel’s draft requirement, the vessel cannot dock. Krimi et al. [19] studied the integrated scheduling of berths and quay cranes in export ports, where vessels are unaffected by tides when entering the port, but must take advantage of the tides when departing.
Uncertainty factors are inevitable in port operations, and research on the integrated scheduling of berths and quay cranes should consider these uncertainties. Iris and Lam et al. [20] addressed planned uncertainties by setting buffer times and using different buffer times to handle different levels of uncertainty. Liu et al. [13] incorporated uncertainty sets into the model to address uncertainties related to vessels and quay cranes and introduced control parameters to manage the degree of uncertainty in the model.
The integrated scheduling problem of berths and quay cranes is an NP-hard problem. The type of model determines the solution method. For a mixed-integer linear programming (MILP) model, heuristic algorithms and CPLEX can be used; for a nonlinear mixed-integer programming model, only heuristic algorithms can solve the model. He et al. [21] considered the fully time-variant strategy for the number of quay cranes, establishing a mixed-integer nonlinear programming model and solving it with a memetic algorithm. Wang et al. [22] established a MILP model considering both single-rate and piece-rate port carbon tax policies, employing equivalent and relaxed constraints to simplify the model, and solving it with CPLEX. Cheimanoff et al. [23] developed a MILP model for the integrated scheduling of berths and quay cranes, employing CPLEX and a variable neighborhood search algorithm to solve the model. Iris et al. [24] developed a MILP model for the integrated berth and quay crane scheduling problem. They improved the model’s solving speed by introducing valid inequalities and designed an adaptive large neighborhood search (ALNS) algorithm based on the problem characteristics to solve the model. Finally, data experiments demonstrated the efficiency of the model improvements and the effectiveness of the algorithm designed in their study.

2.3. Tugboat Scheduling

The tugboat scheduling problem can be divided into static and dynamic tugboat scheduling problems. The static tugboat scheduling problem assumes that all vessels have already arrived at the port before the scheduling begins. Xu et al. [25] studied the static tugboat scheduling problem by considering the vessel relocation, establishing a three-stage model including vessel berthing, shifting, and departure, and designed a hybrid simulated annealing-ant colony optimization algorithm to solve the model.
The dynamic tugboat scheduling problem allows vessels to arrive at the port at any time during the planning period, with all vessel arrival times known. Wei et al. [26] pointed out that static scheduling does not conform to practical scenarios and proposed a mixed-integer programming model where vessels can arrive at different times during the planning period, designing a customized heuristic algorithm to solve the model. Sun et al. [27] studied the dynamic tugboat scheduling problem in large-scale ports by considering the tugboat operation uncertainties and tugboat transit area restrictions. A mixed-integer programming model was developed, and a reverse operation-based genetic algorithm was designed to solve the model. The algorithm’s effectiveness was validated with data from Zhoushan Port. Li et al. [28] considered the uncertainty of vessel arrival times, established a stability model with the objective of minimizing the sum of vessel waiting costs and tugboat fuel consumption, and then solved the model using a reinforcement learning algorithm. Zhong et al. [29] considered the impact of tidal factors on tugboat assistance for vessel entry, created a multi-objective mixed-integer linear programming model aiming to minimize tugboat handling time and fuel consumption, and designed an NSGA-II algorithm to solve the model. The algorithm and model’s effectiveness were validated through experiments using data from Guangzhou Port. Kang et al. [30] studied the tugboat scheduling problem considering uncertainties in vessel arrival times and tugboat towing processes, established a mixed-integer linear programming model for tugboat scheduling, and designed a self-organizing map (SOM) algorithm to solve the model.
Wei et al. [31] noted the path point selection patterns in port scheduling practices and developed a mixed-integer linear programming model incorporating these patterns, designing a branch-and-bound algorithm with six different effective inequalities to solve the model. The algorithm and model’s effectiveness were validated through experiments with data from Singapore Port. Jia et al. [32] considered the vessel berthing plans and different vessel entry paths with the objectives of minimizing vessel delays, tugboat operation costs, and unserved vessels. A mixed-integer linear programming model was established, and a heuristic method based on Lagrangian relaxation and Benders decomposition was designed to solve the model. The heuristic’s effectiveness was validated with data from Shanghai. Wang et al. [33] aimed to improve tugboat utilization by developing a mixed-integer programming model for tugboat scheduling by considering multiple vessel entry path modes. An adaptive large neighborhood search algorithm incorporating a feasible path detection procedure was designed to solve the model, and the model’s effectiveness was validated through case studies.

2.4. Integrated Scheduling of Tugboats, Berths, and Quay Cranes

There is limited research on the integrated scheduling of tugboats, berths, and quay cranes. To the best of our knowledge, Yang et al. [34] are the only researchers who have studied this problem. The integrated scheduling of tugboats, berths, and quay cranes, considering the economic benefits, was investigated. A nonlinear mixed-integer programming model for the integrated scheduling problem was developed, and an adaptive satin bower bird optimizer algorithm integrating chaotic mapping and quantum computing strategies was designed to solve the model. In the tugboat scheduling part of their study, each berthing and unberthing operation required a certain number of tugboats with specific horsepower to be allocated to the vessel. The time required for berthing and unberthing varied depending on the number of tugboats assigned. For berth allocation, the study focused on discrete berth allocation layouts. In quay crane scheduling, a time-invariant strategy for the number of quay cranes was employed, ensuring that the quay cranes did not cross each other. Our study also investigated the integrated scheduling of tugboats, berths, and quay cranes. Table 1 compares our study with the study by Yang et al. The key differences include accounting for the uncertainty of vessel arrival times, utilizing a continuous berth layout, implementing a time-variant strategy for quay crane allocation, developing a linear mixed-integer programming model, and employing an adaptive large neighborhood search algorithm to solve the model.

3. Model

3.1. Problem Description

As shown in Figure 1, the integrated scheduling problem of tugboats, berths, and quay cranes studied in this paper can be described as follows. After arriving at the port, vessels enter the anchorage area. If the number of available tugboats meets the required number for assisting a ship in entering the port, the corresponding tugboats will be allocated to facilitate its entry. If the available tugboats cannot meet the number requirements, the vessel waits at the anchorage until sufficient tugboats are available. Similarly, if there are not enough continuous berths available on the quay line, the vessel also needs to wait at the anchorage until sufficient berths are free to enter the port. Once berthed, quay cranes are immediately allocated to the vessel for loading and unloading operations. If the vessel’s berthing position deviates from its preferred berth, additional container truck transportation costs are incurred, referred to as berth deviation penalties in this paper. Since all quay cranes move on the same track, they cannot cross each other. If the number of available tugboats meets the required number for assisting a ship in leaving the port, the corresponding tugboats will be allocated to facilitate its departure. If the available tugboats do not meet the number requirements, the vessel waits at the berth until enough tugboats are available for departure. Time-variant constraints on the number of quay cranes were considered in this work. During the vessel’s waiting period after the completion of loading and unloading operations, the quay cranes allocated to this vessel can move to other berths to handle other newly berthed vessels. This study primarily addresses the uncertainty in vessel arrival times by implementing a buffer period. After a vessel completes its loading and unloading operations, the berth it occupied must remain vacant for a designated buffer time before accommodating a new incoming vessel. As illustrated in Figure 1, after entering the port, Ship 3 can only dock in the green-marked mooring area near the shoreline and is prohibited from docking in the yellow-marked non-mooring area because Ship 5 has just vacated that zone, which must remain unavailable for other vessels until the buffer time has expired. This buffer period allows for adjustments in the loading and unloading schedules when there are fluctuations in vessel arrival times. The larger the buffer interval, the higher the stability of the model.

3.2. Model Assumptions

(1) The port water depth meets the draft requirements of all vessels; (2) vessels will not be relocated, and each vessel can berth only once; (3) the berthing of vessels must meet physical conditions; (4) each vessel has a preferred berth; (5) the safe distance between vessels is included in the vessel length; (6) the number of quay cranes allocated to each vessel has both upper and lower limits; (7) vessel deviation from the preferred berth will incur additional container truck scheduling costs; (8) vessels start loading and unloading services immediately after berthing; (9) the number of tugboats allocated to a ship should meet the required number necessary for the ship’s entry or departure from the port; and (10) each tugboat provides the same amount of horsepower.

3.3. Model Parameters

The parameters and variables of the integrated scheduling model for tugboats, berths, and quay cranes constructed in this paper are shown in Table 2.

3.4. Scheduling Model

3.4.1. Objective Function

F = F 1 + F 2 + F 3 + F 4 + F 5
F 1 = C 1 v V D T v A T v
F 2 = C 2 v V R v + S 1 v A T v + S 2 v E H T v
F 3 = C 3 v V D S v
F 4 = C 4 v V s S s E T v s S A v s + s L T v s S E v s
F 5 = C 5 v V q Q q A Q v q H T v
In the objective function, F1 represents the fuel consumption costs of all vessels during their stay in the port, F2 represents the total penalties for the delayed departure and waiting time of vessels, F3 represents the total penalties for the berth deviation of vessels, F4 represents the costs incurred by tugboats assisting vessels in entering and leaving the port, and F5 represents the costs incurred by quay cranes loading and unloading vessels

3.4.2. Tugboat-Related Constraints

S 1 v A T v , v V
B T v = S 1 v + s S S A v s E T v s , v V
s S S A v s = 1 , v V
t R T v t + M 1 R T v t S 1 v , v V , t T
t T R T v t = B T v S 1 v , v V
t + 1 R T v t B T v , v V , t T
s S s S A v s S N v , v V
S 2 v E H T v , v V
D T v = S 2 v + s S S E v s L T v s , v V
s S S E v s = 1 , v V
t C T v t + M 1 C T v t S 2 v , v V , t T
t T C T v t = D T v S 2 v , v V
t + 1 C T v t D T v , v V , t T
s S s S E v s S N v , v V
v V s S s S A v s R T v t + v V s S s S E v s C T v t S , t T
Constraint (7) ensures that vessels can only enter the port and berth after arriving. Constraint (8) ensures that the berthing time with tugboat assistance is related to the number of tugboats assigned to the vessel. Constraint (9) ensures that the number of tugboats assigned to a vessel remains constant during the tugboat-assisted berthing process. Constraints (10) to (12) ensure the continuity of the vessel’s entry and berthing times. Constraint (13) ensures that the number of tugboats allocated to incoming ships satisfies the tugboat requirements for these ships. Constraint (14) ensures that vessels can only depart after completing loading and unloading operations. Constraint (15) ensures that the unberthing time with tugboat assistance is related to the number of tugboats assigned to the vessel. Constraint (16) ensures that the number of tugboats assigned to a vessel remains constant during the tugboat-assisted unberthing process. Constraints (17) to (19) ensure the continuity of the vessel’s unberthing and departure times. Constraint (19) ensures that the number of tugboats allocated to departing ships satisfies the tugboat requirements for these ships. Constraint (20) ensures that the number of tugboats assisting vessels in entering and leaving the port at any given time does not exceed the maximum number of tugboats available at the port.

3.4.3. Vessel-Related Constraints

S 2 v E H T v , v V
E H T v > = B T v + H T v , v V
D S v P v P B v , v V
D S v P B v P v , v V
S 2 v + G A P B T v + 1 X v v , v , v V , v v
P v + L v P v + 1 Y v v , v , v V , v v
X v v + X v v + Y v v + Y v v 1 , v , v V , v v
P v + L v B , v V
R v 0 , v V
R v D T v L D T v , v V
D T v H 1
Constraint (22) ensures that vessels can only depart after completing loading and unloading operations. Constraint (23) ensures that vessels receive loading and unloading services immediately after berthing and defines the relationship between the vessel’s berthing time, loading and unloading time, and the completion time of these operations. Constraints (24) to (25) ensure that the berth deviation distance for vessels is positive. Constraints (27) to (28) ensure that the quay space and time occupied by any two berthed vessels do not overlap. Constraint (29) ensures that vessels can only berth within the quay area. Constraints (30) to (31) ensure that the delayed departure time for vessels that did not experience a delayed departure is zero, and that it is positive for vessels that did experience a delayed departure. Constraint (32) ensures that the vessel’s departure time is within the planning period.

3.4.4. Quay Crane-Related Constraints

q Q q A Q v q H T v F T v , v V
q Q A Q v q = 1 , v V
L Q v R Q v , v V
q Q q A Q v q = R Q v L Q v + 1 , v V
q Q q A Q v q Q U v , v V
q Q q A Q v q Q L v , v V
R Q v Q 1 , v V
E H T v B T v + M 1 Z v v , v , v V , v v
R Q v L Q v + M 1 Y v v + M Z v v + M Z v v 1 , v , v V , v v
Constraint (33) ensures that vessels can complete loading and unloading tasks within the given time and conditions. Constraint (34) ensures that the number of quay cranes assigned to a vessel is unique. Constraints (35) to (36) ensure that the quay crane numbers assigned to the same vessel are consecutive. Constraints (37) to (38) ensure that the number of quay cranes assigned to a vessel is within the available range of quay cranes for that vessel. Constraint (39) ensures that the total number of quay cranes assigned to all vessels does not exceed the total number of quay cranes. Constraints (40) to (41) ensure that the quay cranes assigned to a vessel can transfer along the track to handle newly berthed vessels after completing the loading and unloading operations for the current vessel, even if the vessel needs to wait at the berth for sufficient tugboats to become available for departure. This ensures that the number of quay cranes assigned to a vessel meets the partial time-variant constraints.
X v v , Y v v , Z v v = 0 , 1 , v , v V , v v
A Q v q = 0 , 1 , v V , q Q
A T 1 v s , A T 2 v s = 0 , 1 , v V , s S
A T 1 v t , A T 2 v t = 0 , 1 , v V , t T
S 1 v , S 2 v , B T v , D T v , H T v , L Q v , R Q v , E H T v , P v , S T I v , E T I v , D S v , R v = 0 , 1 , 2 , , v V
Constraints (42) to (46) ensure the type and range of values for the decision variables.

3.5. Model Linearization

In Equations (6) and (42), AQvqHTv is nonlinear. To linearize it, we need to introduce an auxiliary variable AU1vq and use AU1vq to replace AQvqHTv. To ensure equivalence, additional constraints (47) to (51) need to be introduced.
A U 1 v q H T v M 1 A Q v q , v V , q Q
A U 1 v q H T v , v V , q Q
A U 1 v q M A Q v q , v V , q Q
A U 1 v q M A Q v q , v V , q Q
A U 1 v q = 0 , 1 , 2 , v V , q Q
In Equation (21), SAvsRTvt and SEvsCTvt need to be linearized by introducing auxiliary variables AU2vst and AU3vst to replace SAvsRTvt and SEvsCTvt, respectively. To ensure equivalence, additional constraints (52) to (58) need to be introduced.
A U 2 v s t S A v s , v V , s S , t T
A U 2 v s t R T v t , v V , s S , t T
A U 2 v s t S A v s + R T v t 1 , v V , s S , t T
A U 3 v s t S E v s , v V , s S , t T
A U 3 v s t C T v t , v V , s S , t T
A U 3 v s t S E v s + C T v t 1 , v V , s S , t T
A U 2 v s t , A U 3 v s t = 0 , 1 , 2 , , v V , s S , t T

4. Methodology

A mixed-integer linear programming (MILP) model for the integrated scheduling of tugboats, berths, and quay cranes was constructed in this paper. Therefore, the model can be solved using the CPLEX solver. However, due to the high complexity of the problem and the large number of variables, the solution time of CPLEX is significantly long, and the solution time increases exponentially with the size of the problem instance. This makes it challenging to apply to practical port scheduling. Hence, a heuristic algorithm is required to obtain a near-optimal integrated scheduling solution in a shorter time. To achieve high-quality integrated scheduling solutions within a shorter period, the adaptive large neighborhood search (ALNS) algorithm was employed to solve the model in this paper.

4.1. Solution Encoding

The solution encoding uses a multi-layer real number encoding format. As shown in Figure 2, the solution encoding for a set of 10 arriving vessels is illustrated.

4.2. Initial Solution Generation

Based on information such as the vessels, quay cranes, and tugboats, an initial feasible solution was generated. The process for generating the initial solution is illustrated in Figure 3.

4.3. Destruction Operators

4.3.1. Random Destruction Operator

Randomly select σ vessels and remove their related information from the solution. The IDs of the removed vessels are recorded in IDdel, and the solution after removing σ vessels is denoted as Xre.

4.3.2. Temporal Destruction Operator

Calculate the temporal penalty F2 for all vessels, and sequentially remove the vessels with the highest F2 values until σ vessels have been removed. The IDs of the removed vessels are recorded in IDdel, and the solution after removing σ vessels is denoted as Xre.

4.3.3. Berth Deviation Destruction Operator

Calculate F3 for all vessels, and sequentially remove the vessels with the highest F3 values until σ vessels have been removed. The IDs of the removed vessels are recorded in IDdel, and the solution after removing σ vessels is denoted as Xre.

4.3.4. Quay Crane Margin Destruction Operator

Define the quay crane margin for vessel v as ABQv = QUv q Q q A Q v q H T v . Calculate the quay crane margin for all vessels, and sequentially remove pairs of vessels with the highest and lowest ABQv values until σ vessels have been removed. The IDs of the removed vessels are recorded in IDdel, and the solution after removing σ vessels is denoted as Xre.

4.3.5. Tugboat Margin Destruction Operator

Define the tugboat margin for vessel v as ABTv = 2|S| − s S s E T v s + s L T v s . Calculate the tugboat margin for all vessels, and sequentially remove pairs of vessels with the highest and lowest ABTv values until σ vessels have been removed. The IDs of the removed vessels are recorded in IDdel, and the solution after removing σ vessels is denoted as Xre.

4.4. Repair Operators

4.4.1. Random Repair Operator

Randomly select a vessel from IDdel one by one and insert its schedule into Xre without conflicting with the schedules of the remaining vessels. The insertion process for a vessel is as follows: randomly allocate a certain number of tugboats to assist its entry into the port, randomly allocate a certain number of quay cranes to it, calculate the vessel’s loading and unloading time in the port, then randomly select a berth for it, and randomly allocate a certain number of tugboats to assist its departure from the port. This process continues until all vessels in IDdel are inserted, completing the loading and unloading operations.

4.4.2. Berth Priority Repair Operator

Sort the vessels in IDdel by the degree of deviation from their preferred berths in descending order, and insert them into the solution in that order. The insertion process is as follows: randomly allocate a certain number of tugboats to assist the vessel’s entry into the port, randomly allocate a certain number of quay cranes to the vessel, calculate the vessel’s loading and unloading time in the port, and then assign the vessel to the berth closest to its preferred berth that can be processed the soonest. Ensure that there are no conflicts with the schedules of the remaining vessels during insertion, and randomly allocate a certain number of tugboats to assist with the vessel’s departure from the port.

4.4.3. Time Priority Repair Operator

Sort the vessels in IDdel by their temporal penalty F2 in descending order, and insert them into the solution in that order. The insertion process is as follows: determine the earliest possible entry time for the vessel, allocate the maximum available number of tugboats at that time to assist with the vessel’s entry into the port, allocate the maximum available number of quay cranes at that time to the vessel, calculate the vessel’s loading and unloading time in the port, and then assign the vessel to the berth that can be processed the soonest, ensuring no conflicts with other vessels. Allocate the maximum available number of tugboats at that time to assist with the vessel’s departure from the port.

4.4.4. Quay Crane Margin Priority Repair Operator

Sort the vessels in IDdel by their corresponding ABQ values in descending order, and insert them into the solution in that order. The insertion process is as follows: first, randomly allocate a certain number of tugboats to assist with the vessel’s entry into the port, randomly assign the vessel to one of the available berths at that time, then allocate the maximum number of quay cranes available at that berth at that time (not exceeding the range of quay cranes available to that vessel), calculate the vessel’s port processing time and departure time, and randomly allocate a certain number of tugboats to assist with the vessel’s departure from the port.

4.4.5. Tugboat Margin Priority Repair Operator

Sort the vessels in IDdel by their corresponding ABT values in descending order and insert them into the solution in that order. The insertion process is as follows: allocate the maximum available number of tugboats at that time to assist the vessel’s entry into the port, randomly assign the vessel to one of the available berths at that time, then randomly allocate the quay cranes available at that berth at that time (within the vessel’s allowable range), calculate the vessel’s port processing time and departure time, and allocate the maximum available number of tugboats at that time to assist with the vessel’s departure from the port.

4.5. Operator Manager

The operator manager in the ALNS algorithm is one of its core components, responsible for selecting and managing the destruction and repair operators. The design of the operator manager allows the ALNS algorithm to dynamically adjust the frequency and priority of operator usage based on past performance, thereby improving the efficiency and quality of the solution. The following details the key functions and operational mechanisms of the operator manager:
(1) Roulette wheel selection method: The roulette wheel selection method is used by the operator manager to randomly select destruction and repair operators based on their weights. Each operator’s weight reflects its historical performance success, with higher-weighted operators having a greater probability of being selected in the future.
(2) Adaptive update: After each iteration, the performance of the operators is evaluated, and their weights are adjusted accordingly based on the improvement or deterioration of the solution. For example, if an operator generates a higher quality solution, its weight is increased; if the result does not improve or worsens, its weight is decreased.
(3) Performance evaluation: Performance evaluation is based on the extent of improvement in the solution generated by the operator. This typically involves comparing the quality of the solution before and after the operator’s application such as cost reduction, service level improvement, etc.
(4) Weight adjustment: An exponential smoothing update rule is used for weight adjustment to ensure that the operator’s weight reflects its recent performance. This helps the algorithm avoid the pitfall of continually selecting operators that are no longer effective.

4.6. Adaptive Large Neighborhood Search Algorithm Process

The process of the adaptive large neighborhood search (ALNS) algorithm is shown in Algorithm 1.
Algorithm 1. Algorithm execution process.
ALNS
Input Parameters:
T: Initial temperature
T0: Termination temperature
qt: Cooling rate
I: Maximum number of iterations
Methods: Set of destruction and repair methods
Weights: Initial weights corresponding to the methods
Start
1. solution ← GenerateInitialSolution ()
2. Set temperature T_initial to T
3. Set cooling rate to r
4. Initialize methods and weights
5. Repeat until termination condition is met:
  6. improv ← False
  7. For i: = 1 to I Do
    8. destroy_method ← SelectMethod (Methods, Weights)
    9. partial_solution ← Destroy (solution, destroy_method)
    10. repair_method ← SelectMethod (Methods, Weights)
    11. new_solution ← Repair (partial_solution, repair_method)
    12. If Evaluate(new_solution) better than Evaluate(solution) or AcceptByProbability():
      13. solution ← new_solution
      14. UpdateWeights (Methods, Weights, Positive)
      15. improv ← True
    16. Else:
      17. UpdateWeights (Methods, Weights, Negative)
    18. If not improv:
      19. TT * qt
    20. If T < T0:
      21. Break
End
Output Solution

4.7. Specific Quay Crane Scheduling Solution

After determining a feasible joint allocation plan, specific quay cranes are assigned to each vessel based on this plan, and the specific quay crane scheduling scheme is determined. Since the feasible berth and quay crane allocation plan has already been obtained, only the variables LQv and RQv remain to be solved. An integer programming method can be used to solve this, convert the solved variables into constants, and remove the constraints related to these already solved variables. The transformed model consists of Equations (1) and (35) to (41). Using CPLEX, the detailed quay crane scheduling scheme can be obtained within 0.1 s.

4.8. Rescheduling of Uncertain Arriving Vessels

When it is determined that the arrival time of certain vessels is delayed, it is necessary to reschedule the corresponding uncertain arriving vessels based on the original scheduling plan. In the rescheduling plan, the schedule of the rescheduled vessels must not conflict with the schedules of other vessels. This insertion process is the same as the solution repair process in the ALNS. The IDs of the vessels that need to be rescheduled are stored in IDdel. Therefore, during rescheduling, it is only necessary to call the repair operator in the ALNS algorithm, based on the modified arrival times of the uncertain arriving vessels, to obtain the corresponding rescheduling plan.

5. Case Study Analysis

A port in Xiamen was studied in this paper, where the container handling quay is 1200 m long, equipped with twelve quay cranes and five tugboats. There are three types of incoming vessel: small, medium, and large. The relevant parameters for these three types of vessels are shown in Table 3. The minimum number of tugboats required for different types of ships to enter or leave the port varies. Additionally, the more tugboats allocated to a ship for its entry or departure, the shorter the time needed. Table 4 illustrates the relationship between the time required for various types of ships to enter or leave the port and the number of tugboats assigned to them. For instance, if five tugboats are allocated to a large ship for its entry, the time required for this operation is 1 h. In this paper, time and space were discretized, with the minimum spatial unit being 50 m and the minimum time unit being 1 h. Therefore, the quay was divided into 24 small berths. Additionally, four different scales of examples with 5, 10, 20, and 40 incoming vessels were generated. The planning periods for these four different scales were 60, 120, 180, and 300 h, respectively. The experiments were conducted on a computer equipped with an Intel(R) Core (TM) i5-10400 CPU @ 2.90 GHz processor and 32 GB of RAM, utilizing example programs within the Python framework PyCharm version 2023 and CPLEX solver version 12.6.

5.1. Model Validation

To validate the effectiveness of the integrated scheduling model for tugboats, berths, and quay cranes constructed in this paper, a small-scale example with five vessels was randomly generated, and the GAP value was set to 1. The specific parameters of the example are shown in Table 5. The example data were imported into CPLEX for solving, and the integrated scheduling solution obtained by CPLEX is shown in Figure 4. This solution satisfies the constraints of tugboats, berths, and quay cranes, and the berthing time intervals between vessels with similar arrival times and the same berth meet the GAP reservation length. Therefore, the MILP model constructed in this paper is effective.
To validate the effectiveness of the partial time-variant constraint on the number of quay cranes used in this paper, a time-invariant MILP model was also constructed. In the time-invariant model, the quay cranes allocated to a vessel can only provide loading and unloading services for other vessels after the vessel departs. Five different small-scale examples with five incoming vessels were generated, with the GAP value set to 1. The example data were imported into the model for solving, and the results are shown in Table 6. In this table, OBJ represents the objective function value and CT represents the solving time. It can be observed that the objective function values of the integrated scheduling solutions obtained by considering the partial time-variant quay crane strategy were all less than or equal to those obtained by considering the time-invariant strategy, and the solving times were similar. This proves that the partial time-variant strategy for quay cranes considered in this paper is feasible and effective.

5.2. Algorithm Validation

5.2.1. Algorithm Parameter Tuning

The outer loop iteration counts I in the adaptive large neighborhood search (ALNS) algorithm was set to 10, based on related research. The inner loop iteration count is determined by T, T0, and qt. Let T0 and qt be 1 and 0.99, respectively. By using different scale examples, we can find the initial temperature T that results in a low objective function value and short runtime. Then, randomly generate four different scale examples with five, ten, twenty, and forty vessels, set different T values for solving, and calculate the objective function value OBJ and computation time CT, which are the averages of five runs under the corresponding conditions. Perform a quadratic fit of the results, as shown in Figure 5. The upper and lower limits of the curve on the vertical axis represent the maximum and minimum values obtained by the algorithm. Analysis indicates that the optimal T values for the four different scales of examples with 5, 10, 20, and 40 vessels should be 2000, 1000, 1000, and 500, respectively.

5.2.2. Algorithm Comparison Analysis

To verify the effectiveness of the ALNS algorithm designed in this paper, four different scales of examples were randomly generated based on the port’s vessel arrival patterns, with five examples generated for each scale. These examples of small-scale were solved using both ALNS and CPLEX, and the results are shown in Table 7. The analysis indicates that for the small-scale instances, CPLEX can obtain the optimal tugboat–berth–quay crane integrated scheduling solution within 5 h, while ALNS can achieve a feasible solution with an average deviation of 0.72% from the optimal solution within 2 min, proving the effectiveness of the proposed algorithm. Moreover, as the size of the test cases gradually increased, the CPLEX solving time escalated significantly. When the example scale increased to more than nine vessels, CPLEX failed to find the integrated scheduling solution within 5 h, while ALNS was less affected by the increase in scale and could obtain a feasible integrated scheduling solution within 3 min. In summary, the ALNS algorithm designed in this paper is effective, has high solving efficiency, and is more suitable for solving large-scale tugboat–berth–quay crane integrated scheduling problems for larger numbers of incoming vessels.
To validate the superiority of the proposed algorithm, four different scales of examples were randomly generated, with five examples generated for each scale. These examples were solved using ALNS, large neighborhood search (LNS), variable neighborhood search (VNS), and genetic algorithm (GA), respectively. The explanations regarding VNS, LNS, and GA are presented in Appendix A. The specific GA operation process can be found in Appendix B. The results are shown in Table 8. Analysis indicates that the solution quality of ALNS was higher than that of other neighborhood search algorithms (LNS and VNS) and the swarm intelligence algorithm (GA). Additionally, neighborhood search heuristics are less affected by the increase in example scale, whereas swarm intelligence heuristics experience a significant increase in solving time due to the complexity of the problem and the limitation on the number of solutions. Therefore, compared to other algorithms, the ALNS has more advantages in solving the integrated scheduling of tugboats, berths, and quay cranes.

5.3. Sensitivity Analysis

To validate the impact of uncertainty factors and berth preference on integrated scheduling, different buffer time intervals (GAP) and berth deviation control coefficients (C3) were set. A case with 20 arriving vessels was used as the research object, and the objective function value (OBJ) as well as F1, F2, F3, F4, and F5 under different GAP and C3 conditions were calculated. The results are shown in Figure 6. Analysis showed that as C3 increased, OBJ also increased. The increase in OBJ was mainly related to the increase in F1, F2, and F3. As C3 increased, the rate of increase in F3 gradually decreased, indicating that the berth deviation of all vessels became smaller as C3 increased. Meanwhile, the combined increase in F1 and F2 showed an increasing trend, indicating that an increase in C3 led to longer vessel waiting times in the port, delayed departure times, and increased time in the port. It is necessary to balance the berth deviation costs and time-related costs based on the actual operating conditions of the port and select an appropriate C3 value.
An increase in the vessel buffer time GAP led to an increase in F1, F2 and OBJ, while the values of F3, F4, and F5 remained relatively unchanged, resulting in an increased OBJ value. Reserving a certain buffer time interval causes vessels with similar arrival times and berths to wait for a period before entering the port for loading and unloading services, thereby delaying the berthing and departure times and increasing the waiting and delayed departure costs. When determining the situation of uncertain arriving vessels, it is necessary to reschedule the uncertain vessels based on the original scheduling plan. After rescheduling, the scheduling cost will inevitably increase. The higher the stability of the original schedule, the smaller the increase in cost compared to the original plan after rescheduling. The increase in rescheduling cost is denoted as Rc, where Rc = OBJ*OBJ, with OBJ* being the total cost after rescheduling and OBJ being the total cost of the original schedule. To verify the stability of the model under different GAP values, this study introduced various scenarios of uncertain vessel arrivals. The key metrics for measuring these uncertain scenarios were the number of delayed arriving vessels, denoted as θ, and the vessel arrival delay time, denoted as β. Figure 7 shows the Rc values obtained by the ALNS algorithm under different θ, β, and GAP conditions for a scenario with a fleet of 20 vessels. The analysis of Figure 7 revealed that the larger the GAP value, the smaller the RC value of the obtained plan, indicating higher model stability. Therefore, an appropriate GAP value should be selected based on the uncertainty fluctuation of vessel arrival times at the port, rather than arbitrarily choosing a large buffer time interval to cope with the uncertainty of vessel arrival times.

6. Conclusions

To address the increasing challenges of maritime demand at container ports and considering the limited scheduling resources, this study investigated the integrated scheduling of port resources such as tugboats, berths, and quay cranes to provide efficient scheduling solutions for ports. The research focused on integrated scheduling under conditions of uncertain vessel arrival times, vessel berthing preferences, and the constraint that quay cranes cannot cross each other. A mixed-integer linear programming (MILP) model was developed with the objective of minimizing the total costs including vessel fuel consumption, vessel waiting and departure delay costs, berthing deviation costs, quay crane handling costs, and tugboat fuel costs. An adaptive large neighborhood search (ALNS) algorithm was designed with destruction and repair operators specific to the characteristics of the integrated scheduling problem, enhancing the search efficiency and solution quality of the algorithm. The conclusions drawn from the data experiments are as follows:
(1)
The constructed MILP model is effective, providing optimal integrated scheduling solutions for tugboats, berths, and quay cranes in a short time for small-scale instances. Considering certain time-varying strategies can yield better scheduling solutions.
(2)
Compared to CPLEX, the ALNS algorithm can find optimal or near-optimal solutions in a shorter time for small-scale instances. For larger-scale instances where CPLEX could not find the optimal solution, ALNS still provided high-quality scheduling solutions efficiently. Compared to other algorithms such as LNS, VNS, and GA, ALNS delivers higher-quality solutions.
(3)
The scheduling cost increases as the buffer time interval (GAP) between vessels increases. The GAP value should be set based on the uncertainty level of vessel arrival times.
(4)
A larger C3 value reduces the degree of berthing deviation, ensuring that all vessels are serviced near their preferred berths, but it also increases the scheduling cost.
Future research on port landside scheduling can extend from the integration of tugboats, berths, and quay cranes to include other components such as yard cranes, achieving integrated port scheduling decisions. In terms of algorithm research, more efficient algorithms or machine learning methods can be used for scheduling studies.

Author Contributions

Conceptualization, L.C., J.Z. and Q.Y.; Methodology, L.C. and J.Z.; Software, J.Z. and Q.Y.; Validation J.Z. and X.C.; Formal analysis J.Z., Q.Y. and X.C.; Investigation, L.C. and J.Z.; Resources, L.C. and J.Z.; Data curation, J.Z., Q.Y. and X.C.; Writing—original draft preparation, L.C. and J.Z.; Writing—review and editing, L.C., J.Z. and Q.Y.; Visualization, J.Z., Q.Y. and X.C.; Supervision, L.C. and Q.Y.; Project administration, L.C.; Funding acquisition, L.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Social Science Fund of China [grant number 23&ZD138].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data, models, or code generated or used during the study are available from the corresponding author by request.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

MethodExplanation
VNSThe variable neighborhood search (VNS) is a metaheuristic algorithm based on local search, designed to solve combinatorial optimization problems. It systematically changes the neighborhood structure of solutions to escape local optima, thereby exploring a broader solution space. The core idea of VNS is to dynamically adjust and switch between different neighborhood structures during the search process, gradually finding better so-lutions.
LNSThe large neighborhood search (LNS) is a metaheuristic algorithm used to solve complex optimization problems. It explores a broader solution space by perturbing a larger portion of the solution in each iteration to find better solutions. The core idea of LNS is to selectively destroy and reconstruct the current solution to escape local optima and search for the global optimum.
GAThe genetic algorithm (GA) is an optimization algorithm based on the principles of natural selection and ge-netics. It simulates the process of biological evolution, iteratively optimizing candidate solutions through op-erations such as selection, crossover, and mutation. Genetic algorithms are particularly well-suited for prob-lems with large search spaces, complex objective functions, or those that are difficult to differentiate. By con-tinuously iterating, they can find approximate optimal solutions to the problem.

Appendix B

(1) Solution encoding: As shown in Figure A1, the solution encoding adopts a four-layer encoding structure.
Figure A1. Solution Encoding Method.
Figure A1. Solution Encoding Method.
Jmse 12 01541 g0a1
(2) Crossover operation: Figure A2 illustrates the method of solution crossover.
Figure A2. Solution crossover.
Figure A2. Solution crossover.
Jmse 12 01541 g0a2
(3) Mutation operation: Figure A3 illustrates the mutation operation of the solution.
Figure A3. Solution mutation.
Figure A3. Solution mutation.
Jmse 12 01541 g0a3
(4) Solution repair operation: Figure A4 illustrates the solution repair operation. When calculating the objective function value, the solution must first undergo repair. The repair process involves using a heuristic designed based on the FCFS (first-come, first-served) principle to repair the solution. Once the repair is completed, the objective function value can be calculated.
Figure A4. Solution repair.
Figure A4. Solution repair.
Jmse 12 01541 g0a4

References

  1. Imai, A.; Nishimura, E.; Papadimitriou, S. The Dynamic Berth Allocation Problem for a Container Port. Transp. Res. Part. B Methodol. 2001, 35, 401–417. [Google Scholar] [CrossRef]
  2. Prencipe, L.P.; Marinelli, M. A Novel Mathematical Formulation for Solving the Dynamic and Discrete Berth Allocation Problem by Using the Bee Colony Optimisation Algorithm. Appl. Intell. 2021, 51, 4127–4142. [Google Scholar] [CrossRef]
  3. Iris, Ç.; Lalla-Ruiz, E.; Lam, J.S.L.; Voß, S. Mathematical Programming Formulations for the Strategic Berth Template Problem. Comput. Ind. Eng. 2018, 124, 167–179. [Google Scholar] [CrossRef]
  4. Umang, N.; Bierlaire, M.; Vacca, I. Exact and Heuristic Methods to Solve the Berth Allocation Problem in Bulk Ports. Transp. Res. Part E Logist. Transp. Rev. 2013, 54, 14–31. [Google Scholar] [CrossRef]
  5. Kordić, S.; Davidović, T.; Kovač, N.; Dragović, B. Combinatorial Approach to Exactly Solving Discrete and Hybrid Berth Allocation Problem. Appl. Math. Model. 2016, 40, 8952–8973. [Google Scholar] [CrossRef]
  6. Venturini, G.; Iris, Ç.; Kontovas, C.A.; Larsen, A. The Multi-Port Berth Allocation Problem with Speed Optimization and Emission Considerations. Transp. Res. Part D Transp. Environ. 2017, 54, 142–159. [Google Scholar] [CrossRef]
  7. Ernst, A.T.; Oğuz, C.; Singh, G.; Taherkhani, G. Mathematical Models for the Berth Allocation Problem in Dry Bulk Terminals. J. Sched. 2017, 20, 459–473. [Google Scholar] [CrossRef]
  8. Cheimanoff, N.; Fontane, F.; Kitri, M.N.; Tchernev, N. A Reduced VNS Based Approach for the Dynamic Continuous Berth Allocation Problem in Bulk Terminals with Tidal Constraints. Expert Syst. Appl. 2021, 168, 114215. [Google Scholar] [CrossRef]
  9. Martin-Iradi, B.; Pacino, D.; Ropke, S. An Adaptive Large Neighborhood Search Heuristic for the Multi-Port Continuous Berth Allocation Problem. Eur. J. Oper. Res. 2024, 316, 152–167. [Google Scholar] [CrossRef]
  10. Xiang, X.; Liu, C. An Almost Robust Optimization Model for Integrated Berth Allocation and Quay Crane Assignment Problem. Omega 2021, 104, 102455. [Google Scholar] [CrossRef]
  11. Wang, Z.; Cheng, J.; Hu, H. A Proactive-Reactive-Based Approach for Continuous Berth Allocation and Quay Crane Assignment Problems with Hybrid Uncertainty. J. Mar. Sci. Eng. 2024, 12, 182. [Google Scholar] [CrossRef]
  12. Ji, B.; Huang, H.; Yu, S.S. An Enhanced NSGA-II for Solving Berth Allocation and Quay Crane Assignment Problem With Stochastic Arrival Times. IEEE Trans. Intell. Transp. Syst. 2023, 24, 459–473. [Google Scholar] [CrossRef]
  13. Liu, W.; Zhu, X.; Wang, L.; Li, S. Rolling Horizon Based Robust Optimization Method of Quayside Operations in Maritime Container Ports. Ocean Eng. 2022, 256, 111505. [Google Scholar] [CrossRef]
  14. Chargui, K.; Zouadi, T.; Sreedharan, V.R.; El Fallahi, A.; Reghioui, M. A Novel Robust Exact Decomposition Algorithm for Berth and Quay Crane Allocation and Scheduling Problem Considering Uncertainty and Energy Efficiency. Omega 2023, 118, 102868. [Google Scholar] [CrossRef]
  15. Yu, J.; Tang, G.; Voß, S.; Song, X. Berth Allocation and Quay Crane Assignment Considering the Adoption of Different Green Technologies. Transp. Res. Part E Logist. Transp. Rev. 2023, 176, 103185. [Google Scholar] [CrossRef]
  16. Li, Y.; Chu, F.; Zheng, F.; Liu, M. A Bi-Objective Optimization for Integrated Berth Allocation and Quay Crane Assignment With Preventive Maintenance Activities. IEEE Trans. Intell. Transport. Syst. 2022, 23, 2938–2955. [Google Scholar] [CrossRef]
  17. Iris, Ç.; Pacino, D.; Ropke, S.; Larsen, A. Integrated Berth Allocation and Quay Crane Assignment Problem: Set Partitioning Models and Computational Results. Transp. Res. Part E Logist. Transp. Rev. 2015, 81, 75–97. [Google Scholar] [CrossRef]
  18. Malekahmadi, A.; Alinaghian, M.; Hejazi, S.R.; Assl Saidipour, M.A. Integrated Continuous Berth Allocation and Quay Crane Assignment and Scheduling Problem with Time-Dependent Physical Constraints in Container Terminals. Comput. Ind. Eng. 2020, 147, 106672. [Google Scholar] [CrossRef]
  19. Krimi, I.; Benmansour, R.; El Cadi, A.A.; Deshayes, L.; Duvivier, D.; Elhachemi, N. A Rolling Horizon Approach for the Integrated Multi-Quays Berth Allocation and Crane Assignment Problem for Bulk Ports. Int. J. Ind. Eng. Comput. 2019, 10, 577–591. [Google Scholar] [CrossRef]
  20. Iris, Ç.; Lam, J.S.L. Recoverable Robustness in Weekly Berth and Quay Crane Planning. Transp. Res. Part B Methodol. 2019, 122, 365–389. [Google Scholar] [CrossRef]
  21. He, J.; Huang, Y.; Yan, W.; Wang, S. Integrated Internal Truck, Yard Crane and Quay Crane Scheduling in a Container Terminal Considering Energy Consumption. Expert Syst. Appl. 2015, 42, 2464–2487. [Google Scholar] [CrossRef]
  22. Wang, T.; Wang, X.; Meng, Q. Joint Berth Allocation and Quay Crane Assignment under Different Carbon Taxation Policies. Transp. Res. Part B Methodol. 2018, 117, 18–36. [Google Scholar] [CrossRef]
  23. Cheimanoff, N.; Fontane, F.; Kitri, M.N.; Tchernev, N. Exact and Heuristic Methods for the Integrated Berth Allocation and Specific Time-Invariant Quay Crane Assignment Problems. Comput. Oper. Res. 2022, 141, 105695. [Google Scholar] [CrossRef]
  24. Iris, Ç.; Pacino, D.; Ropke, S. Improved Formulations and an Adaptive Large Neighborhood Search Heuristic for the Integrated Berth Allocation and Quay Crane Assignment Problem. Transp. Res. Part E Logist. Transp. Rev. 2017, 105, 123–147. [Google Scholar] [CrossRef]
  25. Xu, Q.; Mao, J.; Jin, Z. Simulated Annealing-Based Ant Colony Algorithm for Tugboat Scheduling Optimization. Math. Probl. Eng. 2012, 2012, 246978. [Google Scholar] [CrossRef]
  26. Wei, X.; Meng, Q.; Lim, A.; Jia, S. Dynamic Tugboat Scheduling for Container Ports. Marit. Policy Manag. 2023, 50, 492–514. [Google Scholar] [CrossRef]
  27. Sun, C.; Li, M.; Chen, L.; Chen, P. Dynamic Tugboat Scheduling for Large Seaports with Multiple Terminals. J. Mar. Sci. Eng. 2024, 12, 170. [Google Scholar] [CrossRef]
  28. Li, J.; Duan, X.; Xiong, Z.; Yao, P. Tugboat Scheduling Method Based on the NRPER-DDPG Algorithm: An Integrated DDPG Algorithm with Prioritized Experience Replay and Noise Reduction. Sustainability 2024, 16, 3379. [Google Scholar] [CrossRef]
  29. Zhong, H.; Zhang, Y.; Gu, Y. A Bi-Objective Green Tugboat Scheduling Problem with the Tidal Port Time Windows. Transp. Res. Part D Transp. Environ. 2022, 110, 103409. [Google Scholar] [CrossRef]
  30. Kang, L.; Meng, Q.; Tan, K.C. Tugboat Scheduling under Ship Arrival and Tugging Process Time Uncertainty. Transp. Res. Part E-Logist. Transp. Rev. 2020, 144, 102125. [Google Scholar] [CrossRef]
  31. Wei, X.; Jia, S.; Meng, Q.; Tan, K.C. Tugboat Scheduling for Container Ports. Transp. Res. Part E-Logist. Transp. Rev. 2020, 142, 102071. [Google Scholar] [CrossRef]
  32. Jia, S.; Li, S.; Lin, X.; Chen, X. Scheduling Tugboats in a Seaport. Transp. Sci. 2021, 55, 1370–1391. [Google Scholar] [CrossRef]
  33. Wang, X.; Liang, Y.; Wei, X.; Chew, E.P. An Adaptive Large Neighborhood Search Algorithm for the Tugboat Scheduling Problem. Comput. Ind. Eng. 2023, 177, 109039. [Google Scholar] [CrossRef]
  34. Yang, Z.-Y.; Cao, X.; Xu, R.-Z.; Hong, W.-C.; Sun, S.-L. Applications of Chaotic Quantum Adaptive Satin Bower Bird Optimizer Algorithm in Berth-Tugboat-Quay Crane Allocation Optimization. Expert Syst. Appl. 2024, 237, 121471. [Google Scholar] [CrossRef]
Figure 1. Schematic of integrated scheduling for tugboats, berths, and quay cranes.
Figure 1. Schematic of integrated scheduling for tugboats, berths, and quay cranes.
Jmse 12 01541 g001
Figure 2. Solution encoding format.
Figure 2. Solution encoding format.
Jmse 12 01541 g002
Figure 3. Initial solution generation process.
Figure 3. Initial solution generation process.
Jmse 12 01541 g003
Figure 4. Gantt chart of integrated scheduling for tugboats, berths, and quay cranes.
Figure 4. Gantt chart of integrated scheduling for tugboats, berths, and quay cranes.
Jmse 12 01541 g004
Figure 5. Comparison of algorithm parameters for different scale examples.
Figure 5. Comparison of algorithm parameters for different scale examples.
Jmse 12 01541 g005
Figure 6. Sensitivity analysis 3D plot.
Figure 6. Sensitivity analysis 3D plot.
Jmse 12 01541 g006
Figure 7. The relationship between stability and GAP.
Figure 7. The relationship between stability and GAP.
Jmse 12 01541 g007
Table 1. Comparison of the study details.
Table 1. Comparison of the study details.
ReferenceConditionsBerth Quay Crane ModelMethod
Yang et al. (2024) [34]Vessel berthing preferenceDiscreteQuay cranes cannot cross each other; quay crane quantity is time-invariantNonlinear mixed-integer programming modelAdaptive satin bower bird optimizer algorithm
This paperVessel berthing preference; vessel arrival time uncertaintyContinuousQuay cranes cannot cross each other; Quay crane quantity is time-variantLinear mixed-integer programming modelCPLEX,
CPLEX, adaptive large neighborhood search algorithm
Table 2. Model parameters and variables.
Table 2. Model parameters and variables.
NotationSets
v, v’Vessel index
tTime period index
qQuay index
sTugboat index
VSet of incoming vessels, V = {1,2, …|V|}
BSet of berths on the quay, B = {1,2, …|B|}
QSet of quay cranes, Q = {1,2, …|Q|}
TSet of time periods, T = {1, 2, …|T|}
SSet of tugboats, S = {1, 2, …|S|}
C1Parameter: vessel in-port fuel consumption rate
C2Parameter: vessel delay departure and waiting cost per time unit
C3Parameter: penalty cost per unit deviation from preferred berth
C4Parameter: tugboat fuel consumption rate
C5Parameter: quay crane usage rate
MParameter: a sufficiently large number
GAPParameter: the parameter represents the minimum buffer time interval between successive vessel berths in the same berth space.
SNvParameter: minimum number of tugboats required for vessel entry and exit
ATv Parameter: arrival time of vessel v
Lv Parameter: length of vessel v
PBv Parameter: preferred berth of vessel v
LDTv Parameter: latest departure time of vessel v
QUv Parameter: maximum number of quay cranes that can be allocated to vessel v
QLv Parameter: minimum number of quay cranes that can be allocated to vessel v
FTv Parameter: quay crane service hours required for vessel v
S1vVariable: time when vessel v starts berthing from the anchorage
ETvsVariable: time required for vessel s to enter the port when v tugboats are allocated
LTvsVariable: time required for vessel v to exit the port when s tugboats are allocated
S2vVariable: time when vessel v starts unberthing from the berth
BTv Variable: berthing time of vessel v
DTv Variable: departure time of vessel v
HTv Variable: loading and unloading time of vessel v at the quay
LQvVariable: index of the leftmost quay crane allocated to vessel v
RQvVariable: index of the rightmost quay crane allocated to vessel v
EHTVariable: time when vessel v finishes loading and unloading
Pv Variable: berthing position of vessel v
DSvVariable: deviation distance of vessel v from the preferred berth
RvVariable: delayed departure time of vessel v
Xvv′Variable: Xvv′ = 1 if vessel v’ berths later than vessel v depart, otherwise Xvv′ = 0
Yvv′Variable: Yvv′ = 1 if vessel v is completely berthed to the left of vessel v′, otherwise Yvv′ = 0
Zvv′Variable: Zvv′ = 1 if vessel v′ berths later than vessel v finishes loading and unloading, otherwise Zvv′ = 1
QvqVariable: AQvq = 1 if q quay cranes are allocated to vessel v, otherwise AQvq = 1
SAvsVariable: SAvs = 1 if s tugboats are allocated to incoming vessel v, otherwise SAvs = 0
SEvsVariable: SEvs = 1 if s tugboats are allocated to outgoing vessel v, otherwise SEvs = 0
RTvtVariable: RTvt if tugboats are allocated to incoming vessel v in period t, otherwise RTvt = 0
CTvtVariable: CTvt = 1 if tugboats are allocated to outgoing vessel v in period t, otherwise CTvt = 0
AU1vqVariable: auxiliary variable introduced to linearize the model
AU2vstVariable: auxiliary variable introduced to linearize the model
AU3vstVariable: auxiliary variable introduced to linearize the model
Table 3. Parameters of different types of vessels.
Table 3. Parameters of different types of vessels.
Vessel TypeLength (m)Quay Crane Handling Time (h)Minimum Required TugboatsMinimum Required Quay CranesMaximum Required Quay CranesProportion (%)
Small[100, 200][10, 50]1[2, 2][2, 3]30
Medium[200, 300][20, 80]2[2, 3][4, 5]30
Large [300, 400][80, 120]3[3, 4][6, 7]40
Table 4. Relationship between vessel entry/exit time and number of tugboats.
Table 4. Relationship between vessel entry/exit time and number of tugboats.
Tugboats12345
Type
Small32211
Medium-3221
Large--321
“-” indicates that the number of tugboats does not meet the vessel’s minimum number requirement.
Table 5. Parameters of small-scale example.
Table 5. Parameters of small-scale example.
VesselArrival TimeLengthPreferred BerthAgreed Departure TimeMinimum Quay CranesMaximum Quay CranesQuay Crane Hours Vessel Type
1144232330Medium
21517272466Medium
397153536109Large
413810353784Large
517519282422Medium
Table 6. Comparison of different quay crane transfer strategies.
Table 6. Comparison of different quay crane transfer strategies.
Size Time-Invariant Partially Time-Variant
OBJCTOBJCT
52357115.982310221.23
2125105.722124206.47
18535.76184810.12
158716.63158717.46
19255.5219078.85
Table 7. Comparison of ALNS and CPLEX solutions of small-scale instance.
Table 7. Comparison of ALNS and CPLEX solutions of small-scale instance.
SizeALNSCPLEXGap
/%
OBJF1F2F3F4F5CTOBJCT
523272109012240177563.822310221.230.73
21491948025260159075.772124206.471.17
184814404235146512.86184810.120
1588172101195121024.47158717.460.06
190714809235151535.4719078.850
6254324210016275191059.742523722.870.79
252624615020265184566.972493213.041.32
264524418036275191057.682607303.161.45
2141202024285163051.56213038.420.51
22012006026285163056.052194331.750.31
7271025012020320200067.862678794.591.19
280727016027325202564.0927821186.870.89
344429627048355247580.303410478.750.99
29112441302731022056.122892110.850.65
300324020028330220558.372970382.051.11
8293126014021370214059.66288410,295.841.62
303829416039385216066.54301113,966.440.92
429634847053405302085.37428412,785.590.28
357727616031370274073.893569135.970.22
370228823039410273583.0936749541.640.76
9368430819031430272564.70---
370033416041430273566.08---
4623382590564453150101.53---
403232216045440306571.7839994229.630.82
414432627043435307081.58---
‘-’ represents that the optimal solution could not be found within 5 h. Gap = [OBJ(ALNS) − OBJ(CPLEX)]/OBJ(CPLEX) × 100%.
Table 8. Comparison of solutions from different algorithms.
Table 8. Comparison of solutions from different algorithms.
SizeALNSLNSVNSGA
OBJCTOBJCTOBJCTOBJCT
5232763.82240076.334238990.412426172.10
214975.77232164.952155114.652210156.78
184812.86186010.87184912.751904143.21
158824.47159013.22159018.221654105.92
190735.47190919.19191825.291977123.45
10421161.144292102.08434882.034674412.34
4318102.28447470.134618154.685008347.56
4851109.82494161.665192221.525530468.78
433759.19456898.174591151.775074385.21
466364.45493467.95473761.415520492.90
20818684.07845397.258456161.759659987.65
8176122.59819990.91828094.329892865.78
8557114.738855127.548797116.9010,182934.21
8462121.39884472.90864774.3495891078.90
8785134.50881189.068828132.9810,366823.45
4018,158119.1618,601212.9618,610181.1521,4261956.34
17,597178.7117,737108.9017,846190.1321,2922078.92
18,598190.2019,821146.1320,259172.18623,4332345.67
18,956195.7820,414164.0419,318135.5524,0742123.89
18,591213.5720,941248.4018,643186.7922,4952290.10
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chu, L.; Zhang, J.; Chen, X.; Yu, Q. Optimization of Integrated Tugboat–Berth–Quay Crane Scheduling in Container Ports Considering Uncertainty in Vessel Arrival Times and Berthing Preferences. J. Mar. Sci. Eng. 2024, 12, 1541. https://doi.org/10.3390/jmse12091541

AMA Style

Chu L, Zhang J, Chen X, Yu Q. Optimization of Integrated Tugboat–Berth–Quay Crane Scheduling in Container Ports Considering Uncertainty in Vessel Arrival Times and Berthing Preferences. Journal of Marine Science and Engineering. 2024; 12(9):1541. https://doi.org/10.3390/jmse12091541

Chicago/Turabian Style

Chu, Liangyong, Jiawen Zhang, Xiuqian Chen, and Qing Yu. 2024. "Optimization of Integrated Tugboat–Berth–Quay Crane Scheduling in Container Ports Considering Uncertainty in Vessel Arrival Times and Berthing Preferences" Journal of Marine Science and Engineering 12, no. 9: 1541. https://doi.org/10.3390/jmse12091541

APA Style

Chu, L., Zhang, J., Chen, X., & Yu, Q. (2024). Optimization of Integrated Tugboat–Berth–Quay Crane Scheduling in Container Ports Considering Uncertainty in Vessel Arrival Times and Berthing Preferences. Journal of Marine Science and Engineering, 12(9), 1541. https://doi.org/10.3390/jmse12091541

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