Storage Space Allocation of Inbound Container in Railway Container Terminal

Efficient storage strategy of railway container terminals is important in balancing resource utilization, reducing waiting time, and improving handling efficiency. In this paper, we consider the formulation and solution algorithm for storage space allocation problem of inbound containers in railway container terminal. The problem is formulated as two-stage optimization models, whose objectives are balancing the workload of inbound containers and reducing the overlapping amounts. An algorithm implement process based on rolling horizon approach is designed to solve the proposed models. Computational experiments on an actual railway container terminal show that the proposed approach is effective to solve space allocation problem of inbound container and is significant for the operation and organization of railway container terminals.


Introduction
With rapid development of passenger-dedicated railway in China, more transportation capacities of railway freight corridor are released, which puts forward huge opportunities and challenges for railway container transportation.At present, operation and organization of railway container terminals cannot meet developing demands of container transportation.To enable container to rapidly transfer between rail and truck, modern transshipment technologies are required to improve the resource utilization in railway container terminals.
As key resources of railway container terminals, storage spaces are responsible for temporarily storing containers which are moved over long distances by container trains and short distances by trucks.Storage space allocation is important in balancing space utilization, improving the efficiency of container handling, and reducing the turnaround time of containers.The storage space allocation is a vital basis and constraint for other resources utilization in railway container terminal.
Inbound and outbound container operations in railway container terminal are different.Inbound containers arrive predictably in large quantities, are temporarily stored in container yard, and depart one by one in an unpredictable sequence.On the contrary, outbound containers arrive in a random sequence and depart predictably.This paper focuses on the inbound containers in railway container terminals.
The rest of this paper is organized as follows: the relevant literature is reviewed in the next section.The storage space allocation problem of inbound containers in railway container terminals is described in Section 3 and formulated in Section 4.An algorithm implement process based on rolling horizon approach is developed in Section 5. Computational results are reported in Section 6 and finally Section 7 covers the conclusion.

Literature Review
The storage space allocation problem of inbound containers in railway container terminal belongs to the storage space allocation problem (SSAP) which is defined as the temporary allocation of the inbound/outbound containers to the storage blocks at each time period with aim of balancing the workload between blocks in order to minimize the storage/retrieval times of containers [1].
The SSAP was firstly formulated for a container terminal in Hong Kong [2].The problem was solved by a rolling horizon approach.In each planning horizon, the problem is decomposed into two levels and formulated each level as a mathematical programming model.In order to determine the storage location of arriving export containers by considering its weight, a dynamic programming model is formulated to minimize the number of relocation movements expected for the loading operation [3].The process of determining the storage locations for outbound containers was divided into two stages: space allocation stage and stage of locating individual containers [4].A storage location assignment problem for outbound containers of maritime terminal was decomposed into two stages.The problem in the first stage is solved by a mixed integer programming model, while a hybrid sequence stacking algorithm is applied to solve the problem in the second stage [5].
An approach for allocating storage space to groups of outbound containers in port container terminals was proposed with considering the impacts of various space-reservation strategies on the productivity of the loading operation [6].For finding the best allocation of containers in a yard bay in order to minimize the number of reshuffles, a domaindependent heuristically guided planner for obtaining the optimized reshuffling plan was proposed with a stacking state and a container demand known [7].In order to improve the operations efficiency of retrieving inbound containers in a modern automatic container terminal, inbound container space allocation models were proposed to optimally allocate the arrival inbound containers so as to minimize the expected container retrieval time [8].A novel approach using antbased control was proposed for allocating containers to storage blocks in a marine container terminal with the competing objectives of balancing the workload among yard blocks and minimizing the distance traveled of trucks between yard blocks and berths [9].
In order to solve the storage space problem of outbound containers for utilizing space efficiently and make loading operations more efficient, two heuristic algorithms are suggested based on the duration-of-stay of containers and the subgradient optimization technique, respectively [4].An extended version of SSAP was proposed and an efficient genetic algorithm was developed to solve the extended problem in a container terminal [1].Construction algorithms and a tabu search heuristic are presented for dynamic space allocation problem to optimize the space/resource assignments during the implementation of project activities [10].Three tabu search heuristics are presented for dynamic space allocation problem.The first heuristic is a simple basic tabu search heuristic.The second heuristic adds diversification and intensification strategies to the first and the third heuristic is a probabilistic tabu search heuristic [11].
A hybrid insertion algorithm is designed for solving the problem which integrates the yard truck scheduling and the storage allocation to minimize the weighted sum of total delay of requests and the total travel time of yard trucks [12].For the problem of determining the stacking positions for incoming containers in automated container terminals, an online search algorithm was proposed to dynamically adjust and optimize a stacking policy by continuously generating and evaluating variants of stacking policies [13].A decision support system was present to manage container stacking problem, berth allocation problem, and the quay crane assignment problem in a coordinated way.A domainoriented heuristic planner for calculating the number of reshuffles needed to allocate the containers in the appropriate place [14].A construction and a hybrid algorithm (HGT) based on the GRASP and tabu search metaheuristics were proposed to solve the dynamic space allocation problem, where project duration is divided into a number of consecutive periods, each of them associated with a number of activities [15].
According to the literature review above, most studies focused on the storage space allocation problem in maritime container terminals.Specific literature on railway container terminal is scarce.Because storage strategy of inbound and outbound containers in railway container terminals is obviously different from maritime container terminals, existing studies are hardly applied in railway container terminals.In this paper, we consider the storage space allocation problem of inbound container in railway container terminals.Given that container arrival-departure time and operation sequences are known, two-stage storage space allocation models are proposed, whose objectives are to balance workloads among container blocks and assign containers to the optimum positions.

Problem Description
The Chinese railway container terminals have advanced arrival-departure lines, storage spaces, and handling equipment.Figure 1 gives a schematic representation of a typical railway container terminal in China.Our study is based on the configuration and layout of the representation in Figure 1.
As observed in Figure 1, container yard of railway container terminal is composed by inbound container yard, outbound container yard, and auxiliary container yard.Since most of inbound container allocation operations occur in inbound container yard, we set inbound container yard as the study scope of this paper.
According to the inbound container status in different handling stages, containers to be handled in inbound container yard can be classified into the following three types.
(i) Inbound containers on container train wait for unloading and allocating to the inbound container yard, abbreviated as ICT.
(ii) ICT containers temporarily stored in inbound container yard wait for loading to trucks for customers, abbreviated as ICTY.
(iii) ICT containers are unloaded and directly loaded to trucks for customers, abbreviated as ICTT.
According to optimization objectives of SSAP, most literatures decomposed the storage space allocation problem into two stages.The first stage is to balance containers workloads among blocks and evenly allocate containers to each block.The second stage is the slot allocation for containers which are allocated to blocks based on the first stage optimization results.In this paper, our study decomposed the storage space allocation problem of inbound container in railway container terminal into two stages based on the optimization objectives of SSAP.
(i) The first stage is workload balance similar to the exciting study, whose objective is to balance workloads of inbound containers among blocks and evenly allocate inbound containers to each block.
(ii) The second stage is container slot allocation based on the first stage optimization results, whose objective is to minimize overlapping amounts of ICT.

Problem Formulation
In this section, according to the problem description above, the storage space allocation problem of inbound container in railway container terminal is formulated as two-stage optimization models based on rolling horizon approach.At each planning epoch, we plan for a fixed horizon in immediate future and execute the plan accordingly up to the next planning epoch; then we formulate a new plan based on the latest information; this pattern goes on continually [2].The rolling of planning horizon is shown in Figure 2. The workload balance is implemented in each planning epoch, and the container slot allocation is implemented in each period of planning epoch.(1) There is enough resource, that is, gantry crane and container yard space, to handle the allocation workload at the considered block.
(2) The arriving and departure time of containers are known in advance and there is no time delay in scheduling period.
(3) Handling sequence of containers is assumed to be known.
(4) The containers in the model are assumed to be of one size.ICT  : number of ICT that are unloaded from trains in period  and to be picked up in period +, 1 ≤  ≤ , 0 ≤  ≤  − ;

Workload
ICT  : number of ICT in block  that are unloaded from trains in period , 1 ≤  ≤ , 1 ≤  ≤ ; ICT  : number of ICT in block  that are unloaded from trains in period  and to be picked up in period  + , 1 ≤  ≤ , 1 ≤  ≤ , 0 ≤  ≤  − ; ICTT  : number of ICTT that are unloaded and directly loaded to trucks in period , 1 ≤  ≤ .

Objective Function.
By attention to the problem description in Section 3, the objective function of workload balance model is written as follows: Objective function is to balance workloads of inbound containers among blocks and evenly allocate ICT to blocks in each planning epoch.Constraint (2) ensures that the total number of ICT that are unloaded from trucks in period  and to be loaded onto rail vehicles in period  +  is the sum of these containers assigned to all the blocks.Constraint (3) ensures that the total number of ICT in block  that are unloaded from trucks in period  is the sum of these containers loaded onto rail vehicles in period  +  in this block.
Constraint (4) indicates that the number of ICTY in block  during period  is the sum of initial number of ICTY in block  and the number of ICTY transferred from the ICT that unloaded in the block in planning epoch.

Objective Functions.
By attention to the problem description in Section 3, the objective functions of container slot allocation model are written as follows: The second stage objective function ( 9) minimizes overlapping amount which is caused by ICT allocated in considered block at same planning period.

Constraints.
The constraints of container slot allocation model are introduced as follows to ensure the practical feasibility of the solution.
(1) Overlapping amounts constrains of ICT: Constraint (11) ensures that each container cannot be allocated upon the empty container slot.
Constraint (12) indicates the allocation preferences constraint of ICT.It ensures that the allocation positions of ICT are close to the truck operation lane in order to reduce loading time of ICTY.

Solution Algorithm
In order to solve the two-stage models presented above, an algorithm implement process based on rolling horizon approach is proposed in this section.The workload balance model is converted to a linear integer programming model and a heuristic algorithm is designed to solve the container slot allocation model.The implement process is shown in Figure 3.

HA Implementation for Container Slot Allocation Model.
The notations of HA are shown in Table 1 and the procedure of HA is shown as follows.
Step 1.According to the initial block information at the beginning of period, get a feasible allocated set of container slots  by removing infeasible container slots.Parameter initialization: let V = 1,  = ,  = 1, = {},  = {},  ICT V = , and go to Step 2.
Step 2. Allocate the V ICT.Search empty container slots in block row from  to  − ( − 1).If empty container slot exists, go to Step 3. If empty container slot does not exist, go to Step 4.
Step 3.   Step 7. Allocation of the V ICT has finished.Let  = ,  = 1, and V = V + 1; then if V ≤ , go to Step 2; otherwise go to Step 8.
Step 8. Calculate the overlapping amounts based on  and output the set of container slot allocation .Procedure terminates.

Computational Experiment
To illustrate the proposed model and algorithm for space allocation problem of inbound container in railway container terminal, computational experiments are performed by using the actual data from a specific railway container terminal  in China [16].In order to evaluate the improvement of our approach, a comparison is made between our approach and random allocation algorithm which is currently used in railway container terminals.Furthermore, to evaluate the effectiveness and practicability of our approach, numerical experiments of 7 days and 30 days are carried out.To implement the proposed algorithm, the parameters related to the specific railway container terminal are needed.There are four blocks in the inbound container yard.Each block is composed by 30 bays, 6 rows, and 2 layers.The numerical experiments are performed based on a personal computer with Intel Core i5-2450M @ 2.50 GHz processors and 4 GB RAM.
Because most of ICTY are picked up no more than two days after they allocated to blocks, we choose 3 days as a planning horizon, 1 day as a planning epoch, and 6   Based on the computational example above, experiment is conducted, and a comparison between our approach (OA) and random allocation algorithm (RAA) is made to evaluate the performance of OA on space allocation problem of inbound container in railway container terminal.The computational results are shown in Table 4 and Figure 4.
As observed in Table 4 and Figure 4, the imbalance amounts and overlapping amounts obtained by OA are both fewer than the amounts obtained by CA, the average GAP of imbalance amounts is 42.7%, and the average GAP of overlapping amounts is 48.8%.In order to evaluate the effectiveness and practicability of our approach, numerical experiments of 7 days and 30 days are conducted.The computational results of 7 days are shown in Table 5 and Figure 5 and of 30 days are shown in Figure 6.
As observed in Table 5 and Figures 5 and 6, the performance of our approach is satisfactory in solving different size instances.The results of computational experiments indicate that our approach is efficient to solve space allocation problem of inbound container in railway container terminals.

Conclusion
In this paper, we considered the space allocation problem of inbound container in railway container terminals with the arrival-departure time and operation sequence of containers known.Two-stage optimization models were proposed; the first stage is workload balance model, whose objective is to balance workloads of inbound containers among blocks and evenly allocate inbound containers to each block.The second stage is container slot allocation model, whose objective is to minimize overlapping amounts of ICT.An algorithm implement process based on rolling horizon approach is designed to solve the proposed models.Computational experiments on an actual railway container terminal show that the models and algorithm proposed in this paper are effective to balance workloads of inbound containers and reduce the overlapping amounts.
In future, proposing a stochastic programming model for container slot allocating problem by considering random factors is a possibility for further research.

Figure 1 :
Figure 1: Schematic representation of a railway container terminal.

4. 1 .
Assumptions.The following four assumptions are introduced for the formulation of the problem.
Balance Model 4.2.1.Notations and Variables.The notations and variables of workload balance model are defined as follows: layer identifier of container slot, 1 ≤  ≤ ; (, , ): container slot of  row,  bay, and  layer;   : departure time of container in container slot of  row,  bay, and  layer;  ICT V : overlapping amounts of the V ICT; : an infinitesimal number;   : if (, , ) has container,   = 1.Otherwise,   = 0;  V : if the V container is allocated to (, , ),  V = 1.Otherwise,  V = 0;  ICT ,(−) : overlapping of (, ,  − ) after ICT was allocated to (, , ).If   <  (−) ,  ICT ,(−) = 1.Otherwise,  ICT : total number of bays in the considered block; : total number of rows in considered block; : maximum layer number of stack; : bay identifier of container slot, 1 ≤  ≤ ; : row identifier of container slot, 1 ≤  ≤ ; : [2] workload balance model proposed in Section 4 is a nonlinear model, because of the objective function of model.In order to obtain a solution, a conversion method should be used to convert the model to a linear model.According to the conversion mentioned in[2], we define   = max {} (ICT  + ICTY  ) and   = min {} (ICT  + ICTY  ).Then the workload balance model can be converted as the linear integer programming model as follows: go to Step 2. Otherwise, let  ICT Step 4. If (, , ) ∉ , go to Step 6. Otherwise allocate container to (, , ) and judge the feasibility of solution.If the solution is not feasible, go to Step 8.If the solution is feasible, go to Step 5.

Table 2 :
Inbound container information of 4 planning periods in 1 day (TEU).

Table 4 :
Comparison between OA and RAA in 1 day.: GAP 1 = (imbalance amounts obtained from RAA − imbalance amounts obtained from OA) * 100/imbalance amounts obtained from RAA. GAP 2 = (overlapping amounts obtained from RAA − overlapping amounts obtained from OA) * 100/overlapping amounts obtained from RAA. Notes

Table 5 :
Comparison between OA and RAA in 7 days.: GAP 1 = (imbalance amounts obtained from RAA − imbalance amounts obtained from OA) * 100/imbalance amounts obtained from RAA. GAP 2 = (overlapping amounts obtained from RAA − overlapping amounts obtained from OA) * 100/overlapping amounts obtained from RAA. Notes