Next Article in Journal
The Symbiotic Relationship between the Antarctic Limpet, Nacella concinna, and Epibiont Coralline Algae
Previous Article in Journal
Controls of Submarine Canyons Connected to Shore during the LGM Sea-Level Rise: Examples from Taiwan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discretization-Strategy-Based Solution for Berth Allocation and Quay Crane Assignment Problem

1
School of Traffic and Transportation Engineering, Central South University, Changsha 410075, China
2
School of Engineering, Deakin University, Waurn Ponds 3216, Australia
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2022, 10(4), 495; https://doi.org/10.3390/jmse10040495
Submission received: 8 March 2022 / Revised: 28 March 2022 / Accepted: 30 March 2022 / Published: 2 April 2022
(This article belongs to the Section Coastal Engineering)

Abstract

:
The continuous berth allocation and quay crane assignment problem considers the size of berths and ships, the number of quay cranes, the dynamic ships and non-crossing constraints of quay cranes. In this work, a mixed-integer linear programming model of this problem is established, aiming at minimizing the total stay time and delay penalty of ships. To solve the model, the continuous berth is separated into discrete segments via a proposed discretization strategy. Thereafter, a large neighborhood search algorithm composed of the random removal operator and relaxed sorting-based insertion operator and a backtracking comparison-based constraint repair strategy are proposed. The effectiveness of the model and algorithm presented is verified via real-life instances with different characteristics, and the performances of different combinations of removal operators and insertion operators in the large neighborhood search algorithmic framework are analyzed. Numerical results show that the large neighborhood search algorithm can optimally solve the small-scale instances in a reasonable time. Meanwhile, the results of large-scale instances show that the large neighborhood search algorithm incorporating the discretization strategy is more efficient than other genetic algorithms based on continuous optimization. With the proposed approach, high-quality berth and quay crane allocation results can be obtained efficiently.

1. Introduction

Maritime container transportation is one of the most important modes of transport in international trade, and about 80 percent of the global trade volume is carried out by sea and ports around the world [1]. With the increasing internationalization of container transportation and continuous growth of port trade volume, ports have become more and more busy and the competition between ports has become increasingly fierce. Berth and quay crane (QC) are the core operation resources of the forefront of container ports and making appropriate scheduling plans is of great significance to the effectiveness of container ports. In order to improve the efficiency of berth and QC operations at container terminals, this paper studies the berth allocation and quay crane assignment specific problem (BACASP).
Over the years, scholars have paid great attention to resource scheduling in ports and have obtained rich research results. Among them, the berth allocation problem is one of the most vital problems, which aims to allocate vessels to berths restricted by certain physical and technical constraints to maximize the efficiency of berth operations [2]. Bierwirth and Meisel [3] conducted a literature review on existing research results of berth allocation problem and classified the problem into continuous, discrete and mixed berth allocation problem according to the berth layouts. In the BAPD a quay is divided into many sections, called berths, where one vessel can be in service at a time. In the BAPC, there is no partition, and a vessel can be docked anywhere as long as there is enough space. Nishimura et al. [4] established a BAPD model by considering the waiting time and handling time of vessels. Subsequently, they further considered the vessels’ berthing priority and investigated the corresponding diversification of BAPD. The BAPD is inefficient in terms of berth utilization. With the rapid growth of transportation demand, discrete berth cannot meet the actual operational requirements. Imai et al. [5] studied the continuous berth allocation problem (BAPC), and obtained the upper and lower bounds of the problem by solving the corresponding simple discrete berth allocation problem (BAPD). Cordeau et al. [6] mapped the BAPD to a multi-depot vehicle routing problem with time windows (MDVRPTW) and presented a flow-based model for the BAPD. In addition, a strategy for discretizing the continuous berth space is proposed, and the BAPC and BAPD are unified into one solution framework. A tabu search algorithm is proposed to solve the BAPC and BAPD. The BAPC has been proved to be a NP-hard problem. Therefore, the research considerations for BAPC is dominated by heuristic algorithms, such as Genetic Algorithm (GA) [7,8,9], Memetic algorithm [10], ant colony optimization [11], particle swarm optimization algorithm [12], tubu search algorithm [13], differential evolution algorithm [14], etc.
In practice, the handling time of a vessel at the berth is related to the QCs assigned to the vessel, and incorporating berth allocation with quay assignment is essential for maximizing the efficiency of berth and quay operation. Park and Kim [15] first proposed the berth allocation and crane assignment problem (BACAP) and proposed a two-stage approach to solve the problem, wherein the first stage used sub-gradient optimization to determine the optimal berthing position and start time of berthing and number of QCs, while the second stage utilized a dynamic programming method to assign specific QCs. Meisel et al. [16] further considered the interference between QCs and the impacts of vessels’ berthing positions on its handling time, and proposed a mixed-integer linear model for the BACAP. Rodriguez-Molins et al. [17] proposed an innovative GRASP-based approach to solve the BACAP. Researchers in [18] proposed a mixed integer programming model for the BACAP and a multi-objective algorithm to achieve the balance between delayed time of vessels and handling energy consumption of QCs. Iris et al. [19] presented a BACAP model by considering the vessel speed-up cost, delayed completion penalty delayed departure penalty and QC cost. Thereafter, an algorithm based on large neighborhood search is proposed to solve the problem, which was confirmed capable of solving instances with up to 80 vessels. In their research, the number of QCs assigned to a vessel was assumed to be able to change dynamically. However, the costs of operations due to frequently moving and resetting of QCs were not considered. He et al. [20] proposed a mixed-integer programming model for BACAP and factors such as the work efficiency of QC operators, the operation cost of QCs and operation efficiency and performance-related monetary compensation between day and night, are considered.
The research mentioned above simplifies the BACAP by ignoring the restrictions that the QCs cannot cross each other along the track. Zhang et al. [21] and Türkoğulları et al. [22] argued that frequent change of QC assignment would result in considerable setup costs, and thus the QC reassignment was inhibited once the operation starts. Regarding this, Türkoğulları et al. [22] considered the non-crossing constraints of QCs, and proposed a mathematical model with a branch-cut-plane algorithm for solving small-scale berth allocation and quay crane assignment specific problems. Agra et al. [23] proposed a mixed-integer model for the BACASP, where heterogeneous QCs with different loading and unloading speeds and the non-crossing constraints of QCs are enforced. The model was solved via a branch-and-cut algorithm, which was proven to be able to effectively solve small-scale BACASPs to optimality. Thanos et al. [24] considered the movement range of QCs and established a BACASP model aiming at minimizing the transshipment distance of containers. Then a local search algorithm is proposed to solve this model. Bouzekri et al. [25] considered different water depths and multiple quays, and proposed a model integrating berth allocation and QC assignment problem.
Due to the high complexity, heuristic algorithms are usually employed for solving the BACAP [26], especially for dealing with large-scale problems. Apart from the methods mentioned above, other heuristics, especially the genetic algorithm [27], have also been applied to solve BACAP. In GA, the constraint handling strategy has a great influence on the quality of the obtained near-optimal solution. Among the studies on the BACAP, superiority of feasible solutions [28], repair method [29] and penalty function [30,31,32] are widely used. The superiority of feasible solutions strategy is easy implementation, but it more suitable for multi-objective models. The repair method refers to adjusting an infeasible solution to its nearest feasible solution. This strategy has high search efficiency. Moreover, the reason why the penalty function strategy has been widely used is its high flexibility and its easy implementation.
In this paper, a mixed-integer nonlinear programming model is proposed for BACASP to minimize the total port stay time of vessels and the penalties of delay with the spatial and temporal constraints such as the lengths of continuous berth, number of QCs and non-crossing of QCs, etc. In the traditional solution method, the berthing position and handling time of vessels are regarded as continuous decision variables, and the number of QCs is regarded as discrete decision variables. BACASP is regarded as a hybrid optimization problem and is solved by evolutionary algorithms such as genetic algorithm (GA) in [33]. This method contains continuous and discrete variables, making it difficult to design the coding information interaction method and to guarantee the convergence performance and stability of the algorithm [34]. To solve the problems associated with the existing methods, a strategy for discretizing the continuous berth space is proposed to solve the BAP [6,35]. The BAP is a subproblem of the BACASP to be addressed in our study.
In this work, an MINLP model is presented for the BACASP. The BACASP is converted to a discrete combinatorial optimization problem via a discretization strategy, and a large neighborhood search (LNS) incorporating a specific constraint-handling technique is proposed to efficiently solve it. The main contributions of this work are three-folds and summarized as follows: (1) A discretization-strategy-based method to solve the BACASP based on its combinatorial optimization property. (2) A specific LNS algorithm is designed to solve the BACASP. Under the framework of LNS algorithm, some destroy and repair operators are designed according to the problem’s characteristics to update the solution and a backtracking comparison-based (BCB) constraint-handling strategy is proposed to enhance the performance of the LNS for achieving feasible solutions of the BACASP, so as to achieve an efficient solution to the BACASP. (3) A large number of experiments are conducted to demonstrate the effectiveness of the proposed algorithm and constraint-handling strategy in solving the BACASP.

2. Problem Description

The problem studied in this paper is described below. Section 2.1 describes the notations and Section 2.2 introduces the problem model.
To facilitate the analysis, generally the BACASP can be described within a space-time two-dimensional coordinate, which is described in Figure 1. The horizontal axis represents the time, the vertical axis represents the berth, and the rectangles represent the vessels, the lengths of which in the vertical and horizontal directions ( l i and h i ) represent the length of the vessel and the handling time, respectively. The vessel can be berthed at any position in the continuous berth. The handling time of a vessel is affected by the number of assigned QCs and its berthing position when considering the interference between the QCs. As a result, it is assumed that each vessel corresponds to an optimal berthing position in which the interference between the QCs is minimal. The QC track is distributed along the coastline, where the QC can move to carry out loading and unloading operations on any vessel. However, due to the limitation of the track structure, the QCs cannot pass each other. Since different vessels cannot be served in the same berth at the same time, any two rectangles in the figure cannot overlap.
The objectives of this paper are to obtain the optimal berth allocation and QC assignment theme, so as to minimize the port stay times and delayed departure time of vessels.
For BACASP, the following assumptions are made in this paper.
  • There is no limitation of vessel draft and water depth.
  • Since the handling time of the QC is much greater than its movement time, the QC’s movement time is negligible.
  • The QC assigned to a vessel cannot serve other vessels during the berthing period of the vessel.
  • The safety distance between vessels is included in the length of the vessels.

2.1. Notation

Table 1 lists the variables used in our model for the BACASP along with their types and descriptions.

2.2. Model of BACASP

This section presents a model for the BACASP which is extended from Türkoğulları et al. (2014). In paper [22], the model assumes that the handling time of vessels is only related to the number of QCs. To be more realistic, we also consider the impacts of berthing preference and interference between QCs on the handling time of vessels. The objective and constraints of the model are as follows:
Min i N ( e i a i ) + ( e i d i ) +
s.t.
c i α h i ( 1 + Δ b i β ) m i , i N
e i = s i + h i , i N
Δ b i p i b i , i N
Δ b i b i p i , i N
e i s j + M ( 1 x i j ) , i , j N , i j
p i + l i p j + M ( 1 y i j ) , i , j N , i j
x i j + x j i + y i j + y j i 1 , i , j N , i j
L c i R c i , i N
c i = R c i L c i + 1 , i N
R c i L c j + M ( 1 y i j ) + M x i j + M x j i , i , j N , i j
s i a i , i N
p i + l i L , i N
c i q i m a x , i N
c i q i m i n , i N
R c i g , i N
x i j , y i j 0 , 1 , i , j N
s i , e i , p i , c i , R c i , L c i 0 , 1 , , , i N
The objective function (1) minimizes the port stay time and the delayed departure time, whose component for each vessel are: (i) the waiting time and handling time of vessels, (ii) the delayed departure time when the vessel is delayed to leave the port.
Constraint set (2) ensures that every vessel receives required QC hours considering productivity losses by QC interference and the chosen berthing position. It is worth noting that the constraint (2) in this paper is different from the constraint in [22]. They assume that the handling time of vessels is only related to the number of QCs. Considering the realistic scenario, the model also considers the impacts of berthing preference and productivity lost by QC interference. Constraint set (3) describes that the departure time of vessel is determined by the berthing start time and the handling time of vessel. Constraints (4) and (5) determine the deviations of the berthing position from the desired berthing position. Constraints (6)–(8) prevent any two vessels from overlapping in the space-time-diagram. Constraints (9) and (10) indicate the relationship of the specific number of QCs assigned to the same vessel. Constraint (11) means that for any two ships being served at the same time, the assigned QCs cannot cross each other. Constraints (12) indicate the vessels berthing after their arrival time. Constraints (13)–(18) are types and boundaries of each variable.

3. Methodology

This section describes the solution method for the BACASP. The framework of solution method is presented in Section 3.1. The main techniques of the method are described in Section 3.2.

3.1. Framework of Solution Method

LNS is a neighborhood search algorithm with strong global search ability, which is widely used in solving discrete optimization problems, and has achieved excellent performance in solving discrete problems [36,37]. In this study, a discretization strategy [6] is applied to divide the continuous berth into discrete segments, which transforms the BACASP into a discrete combinatorial optimization problem, and is further solved by the LNS algorithm. The frame of the LNS is shown in Algorithm 1.
The BACASP optimal solution needs to determine the order of services at the berths and the QCs allocated to the vessels. Let K be the set of berths after discretization and | K | denote the total number of berths (line 1 in Algorithm 1). For a particular berth k, let P k s and P k f be the left and right boundaries of berth k, respectively. We define a solution structure as follows in function (19):
X = S 1 ( B 1 ) S 1 ( B 2 ) S 1 ( B n 1 ) S 2 ( B 1 ) S 2 ( B 2 ) S 2 ( B n 2 ) S | K | ( B 1 ) S | K | ( B 2 ) S | K | ( B n | K | )
The columns of the solution structure represent the berth, and the rows represent the service order of vessels. X 0 , X n o w , X a c c e p t and X b e s t are the initial solution, the current solution, the accepted solution and the best solution, respectively. Set term S k ( B u ) denotes the u-th served vessel on berth k and the number of QCs assigned to the vessel. Hitherto, we prove that the optimal berthing position, berthing start time and number of cranes assigned to each vessel can be obtained as long as the solution structure S k ( B u ) is determined.
The berthing start time s B 1 of the first vessel on the berth is set as the arrival time of the vessel, and s B u of subsequent vessels must be later than the departure time of the previous one. The berthing start time of the vessel on berth k is:
s B u = a B u , u = 1 m a x ( a B u , a B u 1 + h B u 1 ) , 1 < u n k
Destroy operators remove τ vessels from the current solution structure (line 5 in Algorithm 1). τ is a random value in [ 1 , N · Φ ] where φ controls the number of vessels to be removed (%). A new solution structure can be obtained by implementing the destroy ( Z ) and repair ( R + ) operations (line 6 in Algorithm 1). There is conflict among vessels on different berths due to the continuity of the berths. The BCB strategy is used to adjust the infeasible solutions under the berthing sequence to form new solutions (line 7 in Algorithm 1). Afterwards, the QC allocation algorithm proposed by [22] is adopted to achieve optimal QC assignment results (line 8 in Algorithm 1). The basic idea is to first determine the number of QCs allocated to the vessels, and then the specific QCs are allocated in the order from left to right of the berthing position of the vessels. These processes are repeated until the stop criteria is met. The details of destroy and repair operators in Algorithm 1 are further discussed in the following subsections.
Algorithm 1: Framework of the LNS algorithm.
Input:Parameters of LNS
Output: X b e s t
1Determine the number of berth segments | K | , and P k s ,
  P k f                                    
2Generate an initial solution X 0
3 X b e s t X 0 , X a c c e p t X 0 , X n o w X 0
4While the stop criteria are not met do
5  Determine τ = R a n d o m [ 1 , N · Φ ]
6 X n o w = R + ( Z ( X n o w ) )
7  Adjust X n o w using the strategy BCB
8  Cranes assignment algorithm
9  Acceptance criteria for solution
10End while

3.2. Main Technique

3.2.1. Discretization Strategy

The discretization strategy was originally proposed by [6] to solve the BAP. This work extends the discretization strategy solution of BAP to address and solve the BACASP. The discretization strategy we proposed is designed to divide the continuous berth into discrete segments, and convert the BACASP to a discrete combinatorial optimization problem, which is efficiently solved by the proposed LNS. Figure 2a is a schematic diagram of the discretized berths. The short solid line divides the berths and the dotted line represents the berth boundary, indicating such as k, k 1 , and k + 1 . Each berth shares part of the adjacent berths.
Take berth k as an example. Berth k occupies the right half of the left berth and the left half of the right berth. The discretization strategy divides the berth segment, and the continuous docking of vessels has no effect. Figure 2b is the berthing sequence of berth k.
The discretization strategy has an impact on the final solution. There are some limitations in discretization strategy. The relationship between the length of discrete berths and vessels affects the scheduling solutions. If the discrete berth is too long, the occupancy rate of the berth will be reduced. If the discrete berth is too short, the discretization strategy will caused a fragmentation of the decision space. All berth segments should allow all types of vessels to dock during the course of the algorithm. Since the berths share a part of the adjacent berths, vessels allocated to different berths may violate the space constraints. Therefore, a constraint-handling strategy is proposed in this paper to solve this problem.

3.2.2. Destroy Operators

(1)
Related destroy operator
The related destroy operator was first proposed in [38] for vehicle routing problem to remove similar nodes. In this study, the degree of similarity between vessel i and j is computed through a relatedness function R ( i , j ) , where a higher value corresponds to less similar vessels. Since the arrival time and berthing position of the vessel have a greater impact on the objective value, the following Equation (21) is constructed.
R ( i , j ) = δ · | a i a j | + η · | p i p j |
where δ and η are the correlation coefficients of the arrival time and berthing position of the vessel, respectively. The operator randomly destroys a vessel i, calculates the R ( i , i ) ( i current solution X), and destroys vessel j with the maximum value of the relatedness function. This procedure is repeated until Φ vessels are removed.
(2)
Random destroy operator
The random destroy operator removes τ vessels randomly from the current solution. This operator expands the range of neighborhood search and avoids the LNS sitting at the local optimum.

3.2.3. Repair Operators

(1)
Deep greedy repair operator
The deep greedy repair operator was proposed by [39] for solving vehicle routing problems. In this study, we propose an operator suitable for this problem based on the idea of deep greedy repair operator. The deep greedy repair operator intends to repair the vessels to the positions with minimal objective value increment among all feasible positions. Let R p be the set of vessels that have been destroyed and X + i be the solution with vessel i reinserted. Let O b ( X ) be the objective value of current solution X. The implementation of this operator with respect to the BACASP is described in the following three steps. Firstly, examine whether set R p is empty and where the vessel (that belongs to R p ) can be repaired. For each vessel i R p , the deep greedy repair operator calculates and records its best insertion position p i , 1 with the lowest cost O b 1 ( X + i , p i , 1 ) (line 2 in Algorithm 2). Afterwards, vessel j = a r g m i n j R p O b 1 ( X + j , p j , 1 ) is inserted into its best insertion position p j , 1 of X. This procedure is repeated until R p is empty.
Algorithm 2: Deep greedy repair operator.
Input:The set of vessels removed R p
Output:The current solution X
1While R p do
2  Calculate O b 1 ( X + i , p i , 1 ) , O b 2 ( X + i , p i , 2 ) , of structure X + i with vessel i inserted to position p i , l
3  Insert vessel j = a r g m i n j R p O b 1 ( X + j , p j , 1 ) to its best insertion position
   p j , 1 into X
4  Remove vessel j from R p
5End while
(2)
Slack sorted repair operator
Mauri et al. [35] proposed the sorted greedy operator, in which vessels berthing in the same berth can only be served in a certain order, which may lead the search trapped in local optimum. As such, we propose the slack sorted repair operator considering the QC capacity demand and the remaining time to the due time of vessel. This operator first sorts the vessels in set R p according to the slack value V i in ascending order. Then, given the current solution, the operator repairs the first vessel into its best position in every berth. We define V i as follows:
V i = m i / ( d i a i )
where m i is required QC capacity of vessel i, and d i a i is the difference between the latest departure time and the arrival time of vessel. The handling time of the vessel is related to m i , and a larger m i indicates that the vessel needs longer time to departure.

3.2.4. BCB Strategy

Due to the limitations of the discretization strategy, the vessel allocation to these positions may not always satisfy the spatial constraints when the distances between the berthing positions are too small. Therefore, the BCB strategy is proposed in this paper to solve this problem.
We use LNS to solve the model combining the constraint handling strategy, which include the BCB strategy. The BCB strategy can solve the infeasible solutions caused by discretization strategy and reduce the waste of berths. Let E k be the set of unadjusted vessels on berth k. Let k be the berth on the left of berth k. The specific steps of BCB strategy are as follows:
Step 1: Select the leftmost berth of the coastline and mark it as k.
Step 2: If E k , continue; Otherwise, go to Step 7.
Step 3: Find the first vessel i E k and determine the berthing position p i of vessel i.
Step 4: If vessel i overlaps with vessel j ( j E k or j E k , and j i ), continue; otherwise, go to Step 6.
Step 5: Eliminate the overlap of vessel i with j and return to Step 4.
Step 6: Remove vessel i from E k and return to Step 2.
Step 7: If berth k is the rightmost berth of the coastline, stop; Otherwise, Continue.
Step 8: k = k + 1 . Return to Step 2.
Figure 3 shows the setting method of berthing position. The berthing position is determined in Step 3. Assume that Vessel 1 is at the berth k and b 1 is the best berthing position of vessel 1. If b 1 < P k s , then p 1 = P k s , as shown in Figure 3a. If P k s b 1 P k f and b 1 + l 1 P k f , then p 1 = b 1 , as shown in Figure 3b. Otherwise, p 1 = P k f l 1 , as shown in Figure 3c.
The overlap can be eliminated by delaying the berthing start time of the vessel and pushing back the berthing position of the vessel, and the influence of the two methods on the value of the objective function is calculated. Finally, the BCB strategy compares the objective value increments and chooses a way to make the increments smaller to eliminate the overlap.

3.2.5. Acceptance Criteria

With respect to the acceptance criteria (line 9 in Algorithm 1), a commonly used acceptance criteria (line 9 in Algorithm 1) in SA is applied in this paper [35]. Specifically, for a new solution X , a neighbor solution X with smaller objective value ( O b ( X ) < O b ( X ) ) is accepted, and with probability e ( O b ( X ) O b ( X ) ) / T otherwise. The temperature T starts at T s t a r t and decreased by Δ T in each iteration and resets to T s t a r t when the freezing temperature T e n d is reached.

4. Numerical Experiments

4.1. Generation of Instances

Two sets of problem instances are employed in this study for computational experiments. The first set of instances is the real data of 21 vessels reported by [22] at the Tianjin Five Continents International Container Terminal in Tianjin during a 200-h planning period. The terminal has a berth 1200 m long with 12 QCs. This set consists of 7 instances with different scales, denoted as i01–i07, and the numbers of vessels are 3, 6, 9, 12, 15, 18, and 21, respectively.
At present, there is no public large-scale example set of the BACASP. The second set of instances include vessels following the generation rules proposed by Türkoğulları et al. [22]. We generate the number of vessels N from 24 to 60. The length l i and the QC capacity demand m i of each vessel i are generated from uniform distribution U ( 3 , 8 ) and U ( 10 , 120 ) . The desired berthing position of vessel is generated by uniform distributions U ( 1 , L l i + 1 ) . The length of the vessel directly determines q i m i n and q i m a x . Vessels with length of 3, 4, 5, 6, 7, and 8 correspond to their q i m i n and q i m a x in the range [ 2 , 2 ] , [ 2 , 3 ] , [ 2 , 4 ] , [ 3 , 5 ] , [ 3 , 6 ] , and [ 3 , 7 ] . The due time is generated according to s i + m i / q i m a x + U ( 5 , 15 ) where m i / q i m a x is rounded up to the nearest integer ( m i / q i m a x ). Table 2 shows the means and variances of the Parameters of the sample.
All the experimental study is coded in C language and carried out on a computer with a 3.00 GHz CPU and 8 GB RAM. For each instance, the analysis is run for 20 consecutive times.

4.2. Parameter Setting

All parameters are used within the appropriate range, and extensive experiments based on the orthogonal design have been performed on different parameter combinations. Table 3 shows the parameters of the LNS algorithm, which include α , β , Φ , N t , T s t a r t , T e n d and Δ T . When the parameter changes within a range, and the most appropriate quantity is selected according to the influence of the variation of the parameter combination on the solution. Afterwards, a certain parameter combination is selected based on the best average result of the testing instances. All experiments are run for 20 consecutive times with different randomness seeds with the same parameter combination, as shown in Table 3.

4.3. Computational Results on the First Set of Instances

4.3.1. Analysis of Operators

To illustrate the performance of slack sorted repair operator, we increase pairwise combinations of different removal and insertion operators and the results are shown in Figure 4 and Figure 5. Figure 4 is the average and the best objective values of the solution obtained by different combinations of destroy–repair operators. Figure 5 depicts the standard deviation and CPU time of the solution obtained by different combinations of destroy–repair operators. Different operators are represented by the related destroy operator (Re), random destroy operator (Ra), deep greedy repair operator (Gr) and slack sorted repair operator (Sl). Figure 4 shows that the Ra-Sl combination obtains the best average and best values in most cases. Figure 5 The standard deviation shows an increasing trend with the increase of vessels’ size. In general, the value of the standard deviation of the Re-Gr combination is the largest, and the value of the standard deviation of the Ra-Sl combination is the smallest. The CPU time in Figure 5 shows that the combinations using Sl operator have shorter running time than the combinations using Gr operator. For scenario i07, the CPU time of the Re-Gr operator is 79.6 s, and the CPU time of the Re-Sl operator is 31.9 s. Compared with the Gr operator, the Sl operator saves ( ( 79.6 31.9 ) / 79.6 ) × 100 % = 59.9 % of the computing time. The reason is that the Sl operator can obtain the insertion order of the vessels by calculating the slack variables of the removed vessels. The Gr operator calculates the insertion cost of the removed vessels, and only obtains the currently best inserted vessels, and then continues the above calculation process until the insertion sequence of all vessels is completed.
In Figure 4, the difference among the best solutions obtained by different operator combinations is very small. To further investigate the difference between operator combinations, we implement them, respectively, based on the observed value obtained by wilcoxon rank sum (WRS) test [40]. At the significance level of 5%, this method detects whether there is a significant difference between the two groups. If yes, it indicates the bilateral detection probability p < 0.05 , rejects the null hypothesis, and the observed value h returns 1; If not, 0 is returned. Table 4 shows the WRS test of different combinatorial operators, where WRS 1 , WRS 2 and WRS 3 are the rank sum statistics of operators Ra-Sl and Ra-Gr, Re-Sl and Re-Gr, respectively. The results of Ra-Sl and Ra-Gr are significantly different in more than half of the examples, the results of Ra-Sl and Re-Sl are significantly different in all examples, and the results of Ra-Sl and Re-Gr are significantly different in most examples.

4.3.2. Analysis of BCB Strategy

Figure 6 shows the comparison results obtained under the BCB strategy and time-adjustment (TA) strategy, where the BCB strategy performs better in mean and variance overall. The time-adjustment strategy is lack of flexibility and cannot get better results in the large or complex scale instances. The BCB strategy combines adjustment of berthing start time and berthing position to ensure the stability and efficiency of the algorithm.

4.3.3. Analysis of LNS Algorithm

In order to confirm the merit of the proposed LNS for solving the BACASP, the results obtained by the LNS are compared with those obtained in Türkoğulları et al. [22], where the necessary and sufficient condition is proposed for generating an optimal solution of the BACASP, which is obtained by commercial solver CPLEX 12.2. In case this condition is not satisfied, the cutting plane algorithm is used until the condition holds.
The handling time of vessels in Türkoğulları et al. [22] is obtained by dividing the QC capacity demand of each vessel according to the number of cranes assigned, which is rounded up to the nearest integer. The costs attracted by vessels’ deviations from their best berthing positions and the impacts of interference among QCs have not been included in the handling time. Instead, the offset of the berthing position is considered in the objective function. To compare with the quality of solutions presented in Türkoğulları et al. [22], the handling time and objective in this set of instances are set to be the same as those in Türkoğulları et al. [22]. Specifically, we set parameters α and β in constraint (2) as 1 and 0, respectively, and the objective function then written as:
min i = 0 n γ 1 | s i a i | + γ 2 | p i b i | + γ 3 max ( 0 , a i + h i 1 d i )
The objective function contains the cost of waiting time of berthing, deviation from the desired berthing position and delayed departures of vessels. The corresponding coefficients γ 1 , γ 2 , and γ 3 are 1000, 1000 and 2000, respectively.
Table 5 shows the comparison results of the LNS algorithm and CPLEX solver for solving the BACASP. The “Best”, “Average” and “RT” in Table 5 denote the best objective value, average objective value and the average CPU time for 20 consecutive runs obtained by the LNS. The table shows that the LNS can solve all instances to optimality. The results also infer that the proposed LNS has good stability. For the instances i06 and i07, the method proposed in Türkoğulları et al. [22] is not able to obtain a feasible solution of the BACASP within the specified time (76,110 s) by CPLEX solver. As a result, the BACASP is relaxed to a BACAP, which is optimally solved by CPLEX. Afterwards, Türkoğulları et al. [22] checked the optimal solutions of the relaxed BACAP in i06 and i07, and found the solutions in aforementioned condition. In contrast, the proposed LNS can find the optimal solutions for all instances without relaxation within 41 s.

4.4. Computational Results on the Second Set of Instances

Due to the complexity of BACASP, exact methods are difficult to find the optimal in the efficient time, especially for dealing with a large-scale problem. Traditionally, BACASP is regarded as a continuous and discrete hybrid optimization problem, which is solved by the evolutionary algorithms, such as GA [33]. The genetic algorithm encodes the solution, interacts with the encoded information of different solutions, and solves the problem by combining the constraint processing method. Among them, the repair method (GA-RM) and the penalty function (GA-PF) strategy are used widely. GA-RM delays the berthing start time of overlapping vessels to get a feasible solution. The GA-PF constructs a fitness function that balances the objective function and constraints violation to obtain a feasible solution.
In order to further verify the discretization strategy and the performance of LNS algorithm in solving large-scale BACASP, the GA-RM, GA-PF and LNS algorithms are tested wiht the second set of examples, respectively. The results of the two algorithms are summarized in Table 6, where “BE”, “ME”, “SD” and “FS” donate the best objective value, the mean, the standard deviation of the objective value and the times feasible solutions are obtained in the 20 runs, respectively. When the scale of vessels is less than 40, GA-PF can obtain better solutions than GA-RM. When the number of vessels reaches 40, the quality of the solution obtained by GA-PF decreases significantly. The “SD” of GA-PF is the largest among all instances. The GA-PF has the advantage of high flexibility but disadvantage of difficulties of the parameters of the objective and the penalty. The penalty function of the GA-PF cannot effectively solve large-scale cases. The GA-RM has high feasible solution search efficiency, but there is a large systematic deviation in the search process. The average solution, the best solution and the standard deviation of the objective function obtained by the LNS algorithm are better than those obtained by GA-RM and GA-PF. Although this paper considers a few factors such as QCs interference and berthing deviation in the handling time, the LNS algorithm satisfies all constraints and obtains a feasible solution in each operation. In comparison, GA-PF cannot often obtain feasible solutions while the proposed method can.
The BACASP allocation scheme obtained by the above method can reflect the stay times of vessels, delayed departure times of vessels and the berth occupancy rate. The objective value of the model includes two parts: the total stay time and the delayed departure time of vessels. Taking the best solution obtained by LNS and GA-RM as an example. Table 7 shows the statistics of the vessel’s stay time in port and delayed departure time, as well as the difference of stay time Δ s t and the difference of delayed departure time Δ d t . The LNS algorithm can effectively reduce the stay time of vessels in port for 6 out of 7 instances. Furthermore, the LNS algorithm reduces the delayed departure time of vessels in all cases. In instance i08, the total times at port and delayed times of the 24 vessels decreased by 117 h and 72 h, respectively.
Berth occupancy rate reflects the degree of occupancy of vessels at berths. The calculation of the berth occupancy rate in this paper is as follows:
B O R = i = 0 n ( O c l i × O c t i ) L × T
where B O R is berth occupancy rate. O c l i and O c t i are the lengths of berths occupied by vessel i and occupied time of vessel i. L is the length of berth and T is the planning period.
The berth occupancy rate can be derived from information such as berthing time, start time, vessel length and handling time of the solution. As shown in Figure 7, the LNS algorithm can effectively improve the berth occupancy rate. Due to the limited resources of continuous berths and QCs, there are fluctuations in berth occupancy. Figure 8 shows the solution obtained by the LNS for the instance i13. It can be seen that the resources of the port are effectively utilized.

5. Conclusions

From a practical point of view, the BACASP allocation solution solved in this paper realizes the full utilization of port berth and QC resources on the premise of ensuring the service level and idleness reduction, which is conducive to avoiding ship congestion, improving the competitiveness of the port and exerting the maximum economic benefits for ports.
In this paper, the BACASP considers the non-crossing constraints of QCs, the mutual interference between QCs and the impact of berthing position on handling time. The final solution meets all constraints, and the allocation scheme is close to real-time operations.
The effect of different operator combinations has significant differences in solving the BACASP within the LNS algorithm in this work. Compared with the deep greedy repair operator, slack sorted repair operator can obtain higher-quality solutions in a shorter time. From the CPU time with different operator combinations in Section 4.3.1. The slack sorted repair operator saves 59.9% of the computing time. In this paper, the BACASP is transformed into a discrete optimization problem by a new discretization strategy. The effective search of solution can be obtained by the LNS with the BCB strategy. Compared with the traditional continuous optimization combined with GA, the proposed method is more effective and can ensure the efficiency and stability of the solution. In addition, the berth and quay crane allocation scheme obtained by the LNS reduces the stay time of vessels and improves the occupancy rate of berth resources.
This study has opened up several directions for future research. First, this study mainly focuses on the BACASP in a certain environment, but there are some uncertain factors in port such as unscheduled vessels and breakdown of QCs, which deserves further investigations. Second, although the LNS has been confirmed to be able to achieve high-quality solutions for the BACASP, the gap between the obtained solution and the optimal solution is unknown. As such, it will be interesting to develop accurate and efficient exact methods, such as branch-and-price and branch-and-cut, to achieve optimal solutions for the BACASP.

Author Contributions

Conceptualization, M.T. and B.J.; methodology, M.T. and B.J.; software, M.T.; validation, M.T., B.J. and X.F.; formal analysis, B.J. and X.F.; investigation, M.T. and B.J.; resources, B.J.; data curation, M.T. and B.J.; writing—original draft preparation, M.T.; writing—review and editing, B.J., X.F. and S.S.Y.; visualization, M.T.; supervision, B.J. and S.S.Y.; project administration, B.J.; funding acquisition, B.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Natural Science Foundation of Hunan Province, China under Grant No. 2020JJ5780, and Central South University under Grant No. 202045007.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BACAPBerth Allocation and Crane Assignment Problem
BACASPBerth Allocation and Crane Assignment Specific Problem
BCBBacktracking comparison-based
LNSLarge neighborhood search
ReRelated destroy operator
RaRandom destroy operator
DeDeep greedy repair operator
SlSlack sorted repair operator
QCQuay crane
NSet of vessels, the number of vessels n = | N |
GSet of QCs, the number of quay cranes g = | G |
LLength of the quay given as a multitude of 1-m sections
a i Expected time of arrival of vessel i N
l i Length of vessel i N
d i Due time of vessel i N .
b i Desired berthing position of vessel i N
m i Quay crane hours demand of vessel i N
q i m i n Lower bound on the number of QCs that can serve vessel i N simultaneously
q i m a x Upper bound on the number of QCs that can serve vessel i N simultaneously
α Interference exponent of cranes.
β Berth deviation factor.
MA large positive number

References

  1. Global Shipping Has Been Hit by the Coronavirus. Now Goods Are Getting Stranded. Available online: https://edition.cnn.com/2020/02/05/business/shipping-coronavirus-impact/index.html (accessed on 26 March 2022).
  2. Lim, A. The berth planning problem. Oper. Res. Lett. 1998, 22, 105–110. [Google Scholar] [CrossRef]
  3. Bierwirth, C.; Meisel, F. A follow-up survey of berth allocation and quay crane scheduling problems in container terminals. Eur. J. Oper. Res. 2015, 244, 675–689. [Google Scholar] [CrossRef]
  4. Nishimura, E.; Imai, A.; Papadimitriou, S. Berth allocation planning in the public berth system by genetic algorithms. Eur. J. Oper. Res. 2001, 131, 282–292. [Google Scholar] [CrossRef]
  5. Imai, A.; Sun, X.; Nishimura, E.; Papadimitriou, S. Berth allocation in a container port: Using a continuous location space approach. Transp. Res. Part B Methodol. 2005, 39, 199–221. [Google Scholar] [CrossRef] [Green Version]
  6. Cordeau, J.F.; Laporte, G.; Legato, P.; Moccia, L. Models and tabu search heuristics for the berth-allocation problem. Transp. Sci. 2005, 39, 526–538. [Google Scholar] [CrossRef] [Green Version]
  7. Tsai, A.H.; Lee, C.N.; Wu, J.S.; Chang, F.S. Novel wharf-based genetic algorithm for berth allocation planning. Soft. Comput. 2017, 21, 2897–2910. [Google Scholar] [CrossRef]
  8. Seyedalizadeh Ganji, S.R.; Babazadeh, A.; Arabshahi, N. Analysis of the continuous berth allocation problem in container ports using a genetic algorithm. J. Mar. Sci. Technol. 2010, 15, 408–416. [Google Scholar] [CrossRef]
  9. Frojan, P.; Correcher, J.F.; Alvarez-Valdes, R.; Koulouris, G.; Tamarit, J.M. The continuous Berth Allocation Problem in a container terminal with multiple quays. Expert Syst. Appl. 2015, 42, 7356–7366. [Google Scholar] [CrossRef]
  10. Mauri, G.R.; de Andrade, L.N.; Lorena, L.A.N. A Memetic Algorithm for a Continuous Case of the Berth Allocation Problem. IJCCI 2011, 25, 105–113. [Google Scholar]
  11. Cheong, C.Y.; Tan, K.C. A multi-objective multi-colony ant algorithm for solving the berth allocation problem. In Advances of Computational Intelligence in Industrial Systems; Springer: Berlin/Heidelberg, Germany, 2008; pp. 333–350. [Google Scholar]
  12. Ting, C.J.; Wu, K.C.; Chou, H. Particle swarm optimization algorithm for the berth allocation problem. Expert Syst. Appl. 2014, 41, 1543–1550. [Google Scholar] [CrossRef]
  13. Tong, J.; Nachtmann, H. A tabu search approach to the cargo prioritisation and terminal allocation problem. OR Spectr. 2020, 12, 147–173. [Google Scholar]
  14. Şahin, C.; Kuvvetli, Y. Differential evolution based meta-heuristic algorithm for dynamic continuous berth allocation problem. Appl. Math. Model. 2016, 40, 10679–10688. [Google Scholar] [CrossRef]
  15. Park, Y.M.; Kim, K.H. A scheduling method for Berth and Quay cranes. OR Spectr. 2003, 25, 1–23. [Google Scholar] [CrossRef]
  16. Meisel, F.; Bierwirth, C. Heuristics for the integration of crane productivity in the berth allocation problem. Transp. Res. Part E Logist. Transp. Rev. 2009, 45, 196–209. [Google Scholar] [CrossRef]
  17. Rodriguez-Molins, M.; Salido, M.A.; Barber, F. A GRASP-based metaheuristic for the Berth Allocation Problem and the Quay Crane Assignment Problem by managing vessel cargo holds. Appl. Intell. 2003, 40, 273–290. [Google Scholar] [CrossRef]
  18. He, J. Berth allocation and quay crane assignment in a container terminal for the trade-off between time-saving and energy-saving. Adv. Eng. Inform. 2016, 30, 390–405. [Google Scholar] [CrossRef]
  19. 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]
  20. He, J.; Wang, Y.; Tan, C.; Yu, H. Modeling berth allocation and quay crane assignment considering QC driver cost and operating efficiency. Adv. Eng. Inform. 2021, 47, 101252. [Google Scholar] [CrossRef]
  21. Zhang, C.; Zheng, L.; Zhang, Z.; Shi, L.; Armstrong, A.J. The allocation of berths and quay cranes by using a sub-gradient optimization technique. Comput. Ind. Eng. 2010, 58, 40–50. [Google Scholar] [CrossRef]
  22. Türkoğulları, Y.B.; Taşkın, Z.C.; Aras, N.; Altınel, İ.K. Optimal berth allocation and time-invariant quay crane assignment in container terminals. Eur. J. Oper. Res. 2014, 235, 88–101. [Google Scholar] [CrossRef]
  23. Agra, A.; Oliveira, M. MIP approaches for the integrated berth allocation and quay crane assignment and scheduling problem. Eur. J. Oper. Res. 2018, 264, 138–148. [Google Scholar] [CrossRef] [Green Version]
  24. Thanos, E.; Toffolo, T.; Santos, H.G.; Vancroonenburg, W.; Berghe, G.V. The tactical berth allocation problem with time-variant specific quay crane assignments. Comput. Ind. Eng. 2021, 155, 107168. [Google Scholar] [CrossRef]
  25. Bouzekri, H.; Alpan, G.; Giard, V. Integrated Laycan and Berth Allocation and time-invariant Quay Crane Assignment Problem in tidal ports with multiple quays. Eur. J. Oper. Res. 2021, 293, 892–909. [Google Scholar] [CrossRef]
  26. Ayvaz, D.; Topcuoglu, H.R.; Gurgen, F. Performance evaluation of evolutionary heuristics in dynamic environments. Eur. J. Oper. Res. 2012, 37, 130–144. [Google Scholar] [CrossRef]
  27. Duan, J.; Liu, Y.; Zhang, Q.; Qin, J. Combined Configuration of Container Terminal Berth and Quay Crane considering Carbon Cost. Math. Probl. Eng. 2021, 2021, 6043846. [Google Scholar] [CrossRef]
  28. Ji, B.; Yuan, X.; Yuan, Y. Modified NSGA-II for solving continuous berth allocation problem: Using multiobjective constraint-handling strategy. IEEE Trans. Cybern. 2017, 47, 2885–2895. [Google Scholar] [CrossRef]
  29. Han, X.L.; Lu, Z.Q.; Xi, L.F. A proactive approach for simultaneous berth and quay crane scheduling problem with stochastic arrival and handling time. Eur. J. Oper. Res. 2010, 207, 1327–1340. [Google Scholar] [CrossRef]
  30. Zhang, X.; Sun, B.; Sun, J.; Gou, Z. The berth and quay cranes integrated scheduling based on redundancy policy. In Proceedings of the 2014 33rd Chinese Control Conference (CCC), Nanjing, China, 28–30 July 2014; Volume 2014, pp. 7595–7600. [Google Scholar]
  31. Rodriguez-Molins, M.; Ingolotti, L.; Barber, F.; Salido, M.A.; Sierra, M.R.; Puente, J. A genetic algorithm for robust berth allocation and quay crane assignment. Prog. Artif. Intell. 2014, 2, 177–192. [Google Scholar] [CrossRef] [Green Version]
  32. Shang, X.T.; Cao, J.X.; Ren, J. A robust optimization approach to the integrated berth allocation and quay crane assignment problem. Transp. Res. Part E Logist. Transp. Rev. 2016, 94, 44–65. [Google Scholar] [CrossRef]
  33. Correcher, J.F.; Alvarez-Valdes, R. A biased random-key genetic algorithm for the time-invariant berth allocation and quay crane assignment problem. Expert Syst. Appl. 2017, 89, 112–128. [Google Scholar] [CrossRef]
  34. Chen, L.; Shen, J.; Qin, L.; Chen, H. An improved ant colony algorithm in continuous optimization. J. Syst. Sci. Syst. Eng. 2003, 12, 224–235. [Google Scholar] [CrossRef]
  35. Mauri, G.R.; Ribeiro, G.M.; Lorena, L.A.N.; Laporte, G. An adaptive large neighborhood search for the discrete and continuous berth allocation problem. Comput. Oper. Res. 2016, 70, 140–154. [Google Scholar] [CrossRef]
  36. Pacino, D.; Van Hentenryck, P. Large neighborhood search and adaptive randomized decompositions for flexible jobshop scheduling. In Proceedings of the Twenty-second International Joint Conference on Artificial Intelligence, Barcelona, Spain, 16–22 July 2011; Volume 102, pp. 1997–2002. [Google Scholar]
  37. Grangier, P.; Gendreau, M.; Lehuédé, F.; Rousseau, L.M. A matheuristic based on large neighborhood search for the vehicle routing problem with cross-docking. Comput. Oper. Res. 2017, 84, 116–126. [Google Scholar] [CrossRef]
  38. Shaw, P. A New Local Search Algorithm Providing High Quality Solutions to Vehicle Routing Problems. 1997, pp. 1–12. Available online: http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.51.1273&rep=rep1&type=pdf (accessed on 5 March 2022).
  39. Ropke, S.; Pisinger, D. An adaptive large neighborhood search heuristic for the pickup and delivery problem with time windows. Transp. Sci. 2006, 40, 455–472. [Google Scholar] [CrossRef] [Green Version]
  40. Chakraborti, S.; Van der Laan, P.; Bakir, S.T. Nonparametric control charts: An overview and some results. J. Qual. Technol. 2001, 33, 304–315. [Google Scholar] [CrossRef]
Figure 1. Diagram of the BACASP.
Figure 1. Diagram of the BACASP.
Jmse 10 00495 g001
Figure 2. Diagram of discrete berths and berthing order. (a) Description of the discretized berths. (b) Description of the berthing sequence of berth k.
Figure 2. Diagram of discrete berths and berthing order. (a) Description of the discretized berths. (b) Description of the berthing sequence of berth k.
Jmse 10 00495 g002
Figure 3. The setting method of berthing position.
Figure 3. The setting method of berthing position.
Jmse 10 00495 g003
Figure 4. Average and the best objective values obtained by different combinations of destroy– repair operators.
Figure 4. Average and the best objective values obtained by different combinations of destroy– repair operators.
Jmse 10 00495 g004
Figure 5. Standard deviation and CPU times obtained by different combinations of destroy– repair operators.
Figure 5. Standard deviation and CPU times obtained by different combinations of destroy– repair operators.
Jmse 10 00495 g005
Figure 6. The results of BCB and time adjustment strategy.
Figure 6. The results of BCB and time adjustment strategy.
Jmse 10 00495 g006
Figure 7. Comparison of the berth occupancy rate of the LNS and GA-RM.
Figure 7. Comparison of the berth occupancy rate of the LNS and GA-RM.
Jmse 10 00495 g007
Figure 8. Best solution obtained by the LNS for instance i13.
Figure 8. Best solution obtained by the LNS for instance i13.
Jmse 10 00495 g008
Table 1. Variables used in model. - is a dimensionless value.
Table 1. Variables used in model. - is a dimensionless value.
VariablesTypeDefinitionUnit
s i IntegerBerthing start time of vessel i N h
p i IntegerBerthing position of vessel i N m
c i IntegerThe number of QCs serving vessel i N -
L c i IntegerThe leftmost QC of the service vessel i N -
R c i IntegerThe rightmost QC of the service vessel i N -
h i IntegerThe handling time of vessel i N h
e i IntegerDeparture time of vessel i N h
x i j BinarySet to 1 if the berthing start time of vessel j N is later than
the departure time of vessel i N , 0 otherwise
-
y i j BinarySet to 1 if vessel i N is fully berthed on the left side of vessel
j N , 0 otherwise
-
Δ b i IntegerDeviation between the desired and the berthing position of
vessel i N
m
Table 2. The means and variances of the Parameter of the sample.
Table 2. The means and variances of the Parameter of the sample.
Parameter l i m i q i min q i max b i a i d i
Mean5.466.42.44.410.6294.7320.6
Variance1.91119.60.21.930.326,827.826,934.1
Table 3. Parameters of the LNS algorithm.
Table 3. Parameters of the LNS algorithm.
ParameterDescriptionRangeTuned Value
α Interference exponent for the cranes(0.9, 0.8)0.9
β Berth deviation factor(0.01, 0.02)0.01
Φ The number of vessels to be removed(0.3, 0.35, 0.4)0.3
N t Number of iterations(2000, 5000, 10,000)5000
T s t a r t Start temperature10001000
T e n d Freezing temperature0.010.01
Δ T Cooling rate0.9750.975
Table 4. The WRS test of different combinatorial operators.
Table 4. The WRS test of different combinatorial operators.
InstancesWRS 1 WRS 2 WRS 3
phphph
i013.11 × 10 4 11.08  × 10 2 15.88  × 10 2 0
i027.39  × 10 1 09.90  × 10 5 12.54  × 10 4 1
i035.78  × 10 2 01.83  × 10 2 11.32  × 10 1 0
i045.18  × 10 2 06.24  × 10 4 18.70  × 10 5 1
i054.77  × 10 3 19.96  × 10 3 18.70  × 10 5 1
i066.38  × 10 3 18.60  × 10 5 18.80  × 10 5 1
i073.48  × 10 3 11.31  × 10 4 11.03  × 10 4 1
Table 5. Comparison results of the LNS algorithm and CPLEX solver for solving the BACASP.
Table 5. Comparison results of the LNS algorithm and CPLEX solver for solving the BACASP.
InstanceTürkoğulları et al. [22]LNS
BestCPLEX RT(s)Cutting Plane Algorithm RT(s)BestAverageRT(s)
i01200037211.3200020000.55
i0221,00053,68935.321,00021,0002.05
i0321,00021,38340.321,00021,0004.85
i0421,00057,81844.921,00021,0008.40
i0535,00076,11083.535,00040,75015.35
i0643,000105.343,00044,50025.60
i0743,00089.343,00045,20040.25
Table 6. The results of LNS algorithm and GA.
Table 6. The results of LNS algorithm and GA.
InstanceNLNSGA-RMGA-PF
BEMESDFSBEMESDFSBEMESDFS
i0824375379.03.820547581.430.920504692.5120.616
i0928493497.84.120739790.733.420573822.4133.613
i1032563570.74.120855888.736.220676897.6133.516
i1136618626.55.1209221008.680.220790936.5137.09
i1240898913.26.3209801060.263.02017742278.4308.310
i135012391278.817.92014481501.242.02019812673.3550.119
i146014821535.824.72018981993.271.92020733249.6735.310
Table 7. The statistical results of stay time and delayed departure time of vessels.
Table 7. The statistical results of stay time and delayed departure time of vessels.
InstanceNLNSGA-RM Δ st Δ dt
Stay TimeDelayed
Departure Time
RT(s)Stay TimeDelayed
Departure Time
RT(s)
i08243621351.554798582.311772
i09284801384.1626113127.1146100
i103255013134.25722133165.6172120
i113660513213.9783139232.1178126
i124088315295.8844136293.7−39121
i1350116772736.151191257555.224185
i14601406761500.81535363876.9129287
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tang, M.; Ji, B.; Fang, X.; Yu, S.S. Discretization-Strategy-Based Solution for Berth Allocation and Quay Crane Assignment Problem. J. Mar. Sci. Eng. 2022, 10, 495. https://doi.org/10.3390/jmse10040495

AMA Style

Tang M, Ji B, Fang X, Yu SS. Discretization-Strategy-Based Solution for Berth Allocation and Quay Crane Assignment Problem. Journal of Marine Science and Engineering. 2022; 10(4):495. https://doi.org/10.3390/jmse10040495

Chicago/Turabian Style

Tang, Min, Bin Ji, Xiaoping Fang, and Samson S. Yu. 2022. "Discretization-Strategy-Based Solution for Berth Allocation and Quay Crane Assignment Problem" Journal of Marine Science and Engineering 10, no. 4: 495. https://doi.org/10.3390/jmse10040495

APA Style

Tang, M., Ji, B., Fang, X., & Yu, S. S. (2022). Discretization-Strategy-Based Solution for Berth Allocation and Quay Crane Assignment Problem. Journal of Marine Science and Engineering, 10(4), 495. https://doi.org/10.3390/jmse10040495

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