Container Loading Optimization in Rail – Truck Intermodal Terminals Considering Energy Consumption

Rail–truck intermodal terminals are an important type of dry port and play a vital role in inland freight transport. This paper addresses the container loading problem in rail–truck intermodal terminals considering energy consumption under the sustainability concept. We analyze the effect factors of energy efficiency for container loading operations and develop an optimization model to minimize the total handling time and container reshuffling. A genetic algorithm is designed to obtain the optimal container loading sequence. Computational experiments on a specific Chinese rail–truck intermodal terminal were conducted to evaluate the performance of our approach. Results show our approach has a good performance for different sizes, and the total handing time, reshuffling times and energy consumption of the handling task are prominently decreased.


Introduction
Intermodal freight transport involves the transportation of freight in an intermodal container or vehicle using multiple modes of transportation (e.g., rail, ship, truck) without any handling of the freight itself when changing modes.The method reduces cargo handling, and so improves security, reduces damage and loss, and allows freight to be transported faster [1].In intermodal freight transport systems, all containers are transshipped by rail, ship and truck in intermodal terminals.The handling efficiency of the terminal significantly impacts the performance and service quality of intermodal freight transport.Therefore, most terminals focus on optimizing loading, unloading and stockpiling operations to improve efficiency.
With rapidly increasing concern about environment pollution, carbon dioxide emissions from transportation becomes a significant environmental threat.Sustainable transportation development is an important research field in sustainability, and green transportation is proposed by many countries to sustain the development of economic globalization.As key nodes of the freight transportation network, intermodal terminals have numerous heavy-duty pieces of equipment.These equipment operations can cause vast quantities of energy consumption and pollutant emission, so intermodal terminals play an important role in sustainable transportation development.Terminals have a responsibility for energy saving and emission reduction.Therefore, terminals need to consider not only handling efficiency but also energy consumption when they optimize loading, unloading and stockpiling operations.
Recently, dry ports received more concern from researchers and practitioners around the world and have a vital role in inland freight transport.The rail-truck intermodal terminal is an important kind of dry port that is configured with advanced equipment.In rail-truck intermodal transportation systems, container trains are used to transport massive quantities of containers for long distances, trucks are responsible for short distance pick-up and delivery activities, and containers are rapidly transshipped in terminals.
With the rapid development of worldwide shipping and inland freight transportation corridors, inland container transportation volume has experienced a sharp increase, which generates a higher demand for handling efficiency of rail-truck intermodal terminals.However, the current handling strategy of terminal is inefficient, which means that the equipment efficiency cannot be fully developed.Therefore, rail-truck intermodal terminals must optimize equipment utilization to improve handling efficiency, as well as considering sustainable development to reduce energy consumption.For now, containers in rail-truck intermodal terminals are mainly outbound containers, and the main handling operation is container loading.Therefore, it is necessary for rail-truck intermodal terminals to optimize container loading.
In this paper, the container loading optimization problem in rail-truck intermodal terminal is considered.We analyze the energy consumption in outbound container loading operations and determine the key factors for improving the energy efficiency.An outbound container loading optimization model is proposed by considering energy consumption and a genetic algorithm is developed to obtain an approximate optimal loading sequence.The rest of paper is organized as follows.The relevant literature is reviewed in next section.The container loading optimization problem is described in the third section and formulated in the fourth section.A genetic algorithm is developed in the fifth section.Computational results are reported in the sixth section and the final section covers the conclusion.

Literature Review
Intermodal transportation has been largely studied in recent literature, mainly focusing on intermodal terminal operations, intermodal transportation network design, intermodal transportation routes optimization and synchronization of operations [2][3][4].As the core node of the intermodal transportation system, intermodal terminals provide equipment and facilities for container transfer between ship, rail and truck.Most intermodal terminal operations research has specifically focused on container ports.Literature on container port operations is multi-faceted, having addressed issues such as berth scheduling [5][6][7][8], quay crane scheduling [9][10][11][12], stowage planning and sequencing [13,14], storage activities in the yard [15,16], and allocation and dispatching of yard cranes and transporters [17][18][19].Many corresponding models and algorithms were developed for operational planning and scheduling in container ports.As the research further develops, problem formulation became more complex as more uncertain factors are considered, and the research focus develops from one operation optimization to integrated optimization of multiple operations [20][21][22].
In contrast with container ports, specific literature on operation optimization of rail-truck intermodal terminal is relatively scarce.Although rail-truck terminals and container ports have similar transfer equipment, the specific operation procedures and rules have significant differences between these two types of intermodal terminals.To compare operations in container ports and rail-truck intermodal terminals, the most significant distinction is that the ship handling area and container yard of the container port are compressed into one handling area in rail-truck intermodal terminals.The rail-mounted gantry crane of rail-truck intermodal terminals is simultaneously responsible for loading-unloading operations and storage activities, which are separately performed by quay cranes and yard cranes in container ports.Thus, the operations related to rail-mounted gantry cranes are more complex, especially the outbound container loading operations, which need to consider container reshuffling in the loading process.Relevant research achievements in container ports cannot be directly applied in rail-truck intermodal terminals.
Recently, operation optimization of rail-truck intermodal terminals has gained more attention, and the existing literature mainly focuses on the storage space allocation problem (SSAP) and the rail-mounted gantry crane scheduling problem (RGCSP).The SSAP of inbound containers was formulated as a two-stage optimization model: first to balance the workload of inbound containers, then reducing the overlapping amounts [23].The container assignment problem of rail-truck transshipment terminals was formulated as a two-stage optimization model for minimizing overlapping amounts and operation distance [24].An exact solution procedure was developed to determine disjunct yard areas of varying size for multiple gantry cranes in polynomial runtime, and to ensure the workload for a given pulse of trains is equally distributed among cranes [25].A dynamic programming approach was proposed to determine yard areas for gantry cranes to accelerate train processing speed [26].The RGCSP was formulated as an optimization model whose objective is to determine an optimization handling sequence in order to minimize rail-mounted gantry crane idle load time in handling tasks [27].An optimization model was proposed for the RMGC scheduling problem based on a dual cycle mode and a genetic algorithm was designed to obtain the optimization handling sequence [28].
For the study of sustainable intermodal transportation development.a multi-period mixed integer nonlinear single objective optimization problem was proposed to minimize transportation, hub location, rerouting, environmental and social costs with near optimal shipment quantities and hub allocations as the prime decisions [29].An extensive survey for environmental sustainability in freight transportation is developed to help fill some of the gaps in the theory and to enhance practice [30].A geospatial intermodal freight transport model was used to examine the environmental, economic, and time-of-delivery tradeoffs associated with freight transportation in the Great Lakes region and examine opportunities for marine vessels to replace a portion of heavy-duty trucks for containerized freight transport [31].A mixed-integer mathematical programming model was presented for a multi-objective, multi-mode and multi-period sustainable load planning problem by considering import/export load flows to satisfy the transport demands of customers and many other related issues [32].For the study on energy consumption in intermodal terminals, the yard crane (YC) scheduling problem was formulated as a mixed integer programming model whose two objectives minimize the total completion delay of all task groups and the total energy consumption of all YCs [17].An optimal model was built with consideration of key factors such as the crane moving distance, turning distance and the practical operation rules, which are directly related to the total energy consumption [33].The problem of integrated quay crane (QC) scheduling, internal truck (IT) scheduling and YC scheduling was formulated as a mixed integer programming model where the objective is to minimize the total departure delay of all vessels and the total transportation energy consumption of all tasks [21].
Based on the literature review above, we can draw three conclusions: (i) most existing studies focused on operations optimization in container ports, and scarce literature has focused on rail-truck intermodal terminals.Because of differences between two types of intermodal terminals in operation procedures and rules, the existing studies are hard directly apply to rail-truck intermodal terminals.(ii) The energy consumption of intermodal terminals has been the subject of much more attention in recent years, but current studies on energy consumption all focus on container ports; specific literature considering energy consumption in rail-truck intermodal terminal is scarce.(iii) In the scarce literature related to rail-truck intermodal terminals, studies only focused on improving handling efficiency and did not consider energy consumption while optimizing operations.

Problem Description
All outbound containers are handled in the main operation area of rail-truck intermodal terminals, which is configured with arrival-departure lines, truck operation lanes, inbound container yard, outbound container yard and rail-mounted gantry cranes.Figure 1 gives a schematic representation of the main operation area in a Chinese rail-truck intermodal terminal.container yard, outbound container yard and rail-mounted gantry cranes.Figure 1 gives a schematic representation of the main operation area in a Chinese rail-truck intermodal terminal.According to the length of arrival-departure lines and the number of rail-mounted gantry cranes, the main operation area is equally divided, and each rail-mounted gantry crane is responsible for a fixed area.Thus, the scope of our study is limited to one fixed handling area.
Outbound containers in rail-truck intermodal terminals can be classified into two types based on their status at different handling stages.The first type is the outbound containers brought into the terminal by trucks, abbreviated as OCT.OCT1 are allocated in the container yard and OCT2 are directly loaded onto rail vehicles.The other is outbound containers that are already in the container yard waiting to being loaded onto rail vehicles, abbreviated as OCC.Because the loading operations are only for OCT2 and OCC, the object of our study is these two types of outbound containers.
Outbound container loading operations can be described as moving from one initial loading position (container yard or truck) to a matching end position (rail vehicle), and an empty move back to the next outbound container initial loading position.An outbound container handling task is finished when all outbound containers loaded onto rail vehicles.A sample of a container loading task is shown in Figure 2. Before the handling task begins, OCC have already been allocated in the container yard.During the handling task, OCT2 are brought into the terminal by trucks and park in the assigned space.All outbound containers have been assigned specific initial and end loading positions before they are loaded.A sample of the loading operation is shown in Figure 3.According to the length of arrival-departure lines and the number of rail-mounted gantry cranes, the main operation area is equally divided, and each rail-mounted gantry crane is responsible for a fixed area.Thus, the scope of our study is limited to one fixed handling area.
Outbound containers in rail-truck intermodal terminals can be classified into two types based on their status at different handling stages.The first type is the outbound containers brought into the terminal by trucks, abbreviated as OCT.OCT 1 are allocated in the container yard and OCT 2 are directly loaded onto rail vehicles.The other is outbound containers that are already in the container yard waiting to being loaded onto rail vehicles, abbreviated as OCC.Because the loading operations are only for OCT 2 and OCC, the object of our study is these two types of outbound containers.
Outbound container loading operations can be described as moving from one initial loading position (container yard or truck) to a matching end position (rail vehicle), and an empty move back to the next outbound container initial loading position.An outbound container handling task is finished when all outbound containers loaded onto rail vehicles.A sample of a container loading task is shown in Figure 2.  According to the length of arrival-departure lines and the number of rail-mounted gantry cranes, the main operation area is equally divided, and each rail-mounted gantry crane is responsible for a fixed area.Thus, the scope of our study is limited to one fixed handling area.
Outbound containers in rail-truck intermodal terminals can be classified into two types based on their status at different handling stages.The first type is the outbound containers brought into the terminal by trucks, abbreviated as OCT.OCT1 are allocated in the container yard and OCT2 are directly loaded onto rail vehicles.The other is outbound containers that are already in the container yard waiting to being loaded onto rail vehicles, abbreviated as OCC.Because the loading operations are only for OCT2 and OCC, the object of our study is these two types of outbound containers.
Outbound container loading operations can be described as moving from one initial loading position (container yard or truck) to a matching end position (rail vehicle), and an empty move back to the next outbound container initial loading position.An outbound container handling task is finished when all outbound containers loaded onto rail vehicles.A sample of a container loading task is shown in Figure 2. Before the handling task begins, OCC have already been allocated in the container yard.During the handling task, OCT2 are brought into the terminal by trucks and park in the assigned space.All outbound containers have been assigned specific initial and end loading positions before they are loaded.A sample of the loading operation is shown in Figure 3. Before the handling task begins, OCC have already been allocated in the container yard.During the handling task, OCT 2 are brought into the terminal by trucks and park in the assigned space.All outbound containers have been assigned specific initial and end loading positions before they are loaded.A sample of the loading operation is shown in Figure 3.According to the loading operations of outbound containers, the energy consumption of container loading in the terminal can be divided into loading crane energy consumption, empty crane energy consumption and reshuffling crane energy consumption.The details of energy consumption are shown in Table 1.

Loading crane energy consumption
This aspect is generated by the loading crane moving from the container yard or trucks to rail vehicles.Because all outbound containers have specific handling positions, the loading crane energy consumption is a fixed value in one handling task.In general, this kind of energy consumption is determined and never changed after the handling task is generated.

Empty crane energy consumption
This aspect is generated by the empty crane moving between two outbound container loading operations.This energy consumption is determined by the outbound container handling sequence.Different handling sequences have a great influence on empty crane energy consumption, so optimizing the handling sequence can directly reduce this kind of energy consumption and improve energy efficiency of outbound container loading operations.

Reshuffling crane energy consumption
This aspect is generated by container reshuffling operations in the loading process.According to the different status of the upper container, reshuffling crane energy consumption can be divided into two types.In the first type, the upper container does not belong to the same handling task, so it must be reshuffled to ensure the handling task can go on smoothly.This type of energy consumption is unavoidable.The other type is caused by the overlapping of containers in the same handling task, and this type of energy consumption can be reduced or avoided by optimizing the handling sequence.
As observed in Table 1, the key issue for improving the energy efficiency of outbound container loading in rail-truck intermodal terminals is to reduce the empty move distance of rail-mounted gantry cranes and container reshuffling operations in the handling task.Therefore, the optimization objective of our study is to optimize the outbound container loading sequence to minimize the handling move distance and container reshuffling.Decreasing the energy consumption from empty According to the loading operations of outbound containers, the energy consumption of container loading in the terminal can be divided into loading crane energy consumption, empty crane energy consumption and reshuffling crane energy consumption.The details of energy consumption are shown in Table 1.

Loading crane energy consumption
This aspect is generated by the loading crane moving from the container yard or trucks to rail vehicles.Because all outbound containers have specific handling positions, the loading crane energy consumption is a fixed value in one handling task.In general, this kind of energy consumption is determined and never changed after the handling task is generated.

Empty crane energy consumption
This aspect is generated by the empty crane moving between two outbound container loading operations.This energy consumption is determined by the outbound container handling sequence.Different handling sequences have a great influence on empty crane energy consumption, so optimizing the handling sequence can directly reduce this kind of energy consumption and improve energy efficiency of outbound container loading operations.

Reshuffling crane energy consumption
This aspect is generated by container reshuffling operations in the loading process.According to the different status of the upper container, reshuffling crane energy consumption can be divided into two types.In the first type, the upper container does not belong to the same handling task, so it must be reshuffled to ensure the handling task can go on smoothly.This type of energy consumption is unavoidable.The other type is caused by the overlapping of containers in the same handling task, and this type of energy consumption can be reduced or avoided by optimizing the handling sequence.
As observed in Table 1, the key issue for improving the energy efficiency of outbound container loading in rail-truck intermodal terminals is to reduce the empty move distance of rail-mounted gantry cranes and container reshuffling operations in the handling task.Therefore, the optimization objective of our study is to optimize the outbound container loading sequence to minimize the handling move distance and container reshuffling.Decreasing the energy consumption from empty crane movements and reshuffling crane movements can be effective in improving the energy efficiency of outbound container loading in rail-truck intermodal terminals.

Coordinate Frame Transforming for Handling Area
In order to expediently formulate the problem, the handling area is transformed into three-dimensional coordinates.In the X direction, the departure line, rows of the container yard and the truck operation lane are orderly encoded from 1 to 2+n; n is the number of rows in the container yard.In the Y direction, the bays of the handling area are orderly encoded from 1 to m; m is the number of bays in the container yard.Figure 4 gives an encoding example of the X/Y direction.
Sustainability 2019, 11, x FOR PEER REVIEW 6 of 16 crane movements and reshuffling crane movements can be effective in improving the energy efficiency of outbound container loading in rail-truck intermodal terminals.

Coordinate Frame Transforming for Handling Area
In order to expediently formulate the problem, the handling area is transformed into three-dimensional coordinates.In the X direction, the departure line, rows of the container yard and the truck operation lane are orderly encoded from 1 to 2+n; n is the number of rows in the container yard.In the Y direction, the bays of the handling area are orderly encoded from 1 to m; m is the number of bays in the container yard.Figure 4 gives an encoding example of the X/Y direction.In the Z direction, the layers are orderly encoded from low to high.Because all handling positions in the departure line and truck operation lane have only one layer, the z-value of these positions are all 1.A sample of coordinate frame transforming is shown in Figure 5.After coordinate frame transforming, the handling positions in the handling task are transformed into coordinate points in a three-dimensional coordinate system, and the loading operations can be depicted as movements among different points.For describing the sequential relationship of operations in the container loading task, we use moving time to replace moving distance, and the corresponding objective is to minimize whole handling time and container reshuffling in the container loading task.crane movements and reshuffling crane movements can be effective in improving the energy efficiency of outbound container loading in rail-truck intermodal terminals.

Coordinate Frame Transforming for Handling Area
In order to expediently formulate the problem, the handling area is transformed into three-dimensional coordinates.In the X direction, the departure line, rows of the container yard and the truck operation lane are orderly encoded from 1 to 2+n; n is the number of rows in the container yard.In the Y direction, the bays of the handling area are orderly encoded from 1 to m; m is the number of bays in the container yard.Figure 4 gives an encoding example of the X/Y direction.In the Z direction, the layers are orderly encoded from low to high.Because all handling positions in the departure line and truck operation lane have only one layer, the z-value of these positions are all 1.A sample of coordinate frame transforming is shown in Figure 5.After coordinate frame transforming, the handling positions in the handling task are transformed into coordinate points in a three-dimensional coordinate system, and the loading operations can be depicted as movements among different points.For describing the sequential relationship of operations in the container loading task, we use moving time to replace moving distance, and the corresponding objective is to minimize whole handling time and container reshuffling in the container loading task.After coordinate frame transforming, the handling positions in the handling task are transformed into coordinate points in a three-dimensional coordinate system, and the loading operations can be depicted as movements among different points.For describing the sequential relationship of operations in the container loading task, we use moving time to replace moving distance, and the corresponding objective is to minimize whole handling time and container reshuffling in the container loading task.

Assumptions
The following assumptions are introduced for the problem formulation: (i) Each loading operation only involves one container once.(ii) Operation positions are assumed to be known and fixed after the loading task is generated.(iii) All loading operations are non-preemptive; that is, once a loading operation starts, it must be completed without any pause or shift.
(iv) The containers in the model are assumed to be of the same size.

Notations and Variables
The parameters and variables used for the mathematical formulations are defined in Table 2. Start time of task i from (a, e, k) to (b, f , Finish time of task i from (a, e, k) to (b, f , l) X i j 1, if task j immediately begins after task i has finished; 0, otherwise.1. Tasks S and T are considered to be the initial and finish status, Thus, when task i is the first task, X Si = 1.When task i is the last task, X iT = 1 r i (a,e,k+m) Reshuffling times of task i caused by its m th upper container B i j (a,e,k),(c,g,l) 0, if task i begins earlier than task j; 1, otherwise.(a, e, k) is the starting position of task i, and (c, g, l) is the starting position of task j

Objective Function
The objective function of the outbound container loading optimization model is written as follows: The objective function is to minimize total handling time and container reshuffling times in the container loading task.

Constraints
The constraints of the outbound container loading optimization model are introduced as follows to ensure the practical feasibility of the solution: (a,e,k),(a,e,k+m) , ∀i, j ∈ T, ∀(a, e, k), (a, e, k + m) ∈ P B ij (a,e,k),(a,e,k+m) ≤ 1 − X ij , ∀i, j ∈ T, ∀(a, e, k), (a, e, k + m) ∈ P (6) r i (a,e,k+m) ≥ 0, ∀i ∈ T, ∀(a, e, k), (a, e, k + m) ∈ P (9) Constraint ( 2) and Constraint (3) are preorder and subsequent operation constrains.Constraint (2) represents that there is at most one subsequent operation for each container loading operation.Constraint (3) represents that there is at most one preorder operation for each container loading operation.Constraint (4) indicates the continuous time relationship between two consecutive loading operations.Constraint (5) defines the reshuffling times of a task container caused by its upper container.Constraint ( 6) defines the relationship of two continuous handling tasks.Constraint (7) and Constraint (8) are the beginning and finishing operation constrains.They represent that the whole container loading task can only have one beginning operation and one finishing operation.

A Genetic Algorithm for the Problem
The outbound container loading optimization ties into the crane scheduling problem, which has been proved to be non-deterministic polynomial [34], and the formulations presented above cannot be solved exactly in a reasonable time.Therefore, we developed a genetic algorithm (GA) to obtain the approximate optimal outbound container loading sequence.
In genetic algorithm implementation, the chromosome representation and genetic operator design are two vital steps for quickly accessing the feasible space and an effective movement toward the optimal solution.The proposed genetic algorithm is illustrated in Figure 6 and the main implementation steps are elaborated as follows.9 be solved exactly in a reasonable time.Therefore, we developed a genetic algorithm (GA) to obtain the approximate optimal outbound container loading sequence.
In genetic algorithm implementation, the chromosome representation and genetic operator design are two vital steps for quickly accessing the feasible space and an effective movement toward the optimal solution.The proposed genetic algorithm is illustrated in Figure 6 and the main implementation steps are elaborated as follows.

Chromosome Representation
The chromosome is designed by permutation encoding.Each chromosome represents a possible loading sequence for the outbound containers and consists of T  number of genes.Each gene represents a loading operation, and the value of the gene is the indices of this loading operation.The gene sequence is the loading sequence implemented from left to right.
The operation indices correlate with operation positions; one operation index has a fixed start position and finish position, and the two positions cannot be changed in GA implementation.The genetic operators are only used for the operation index.Figure 7 shows an example of chromosome representation.The chromosome represents a handling task with five loading operations, which are handled in the sequence of

Chromosome Representation
The chromosome is designed by permutation encoding.Each chromosome represents a possible loading sequence for the outbound containers and consists of T number of genes.Each gene represents a loading operation, and the value of the gene is the indices of this loading operation.The gene sequence is the loading sequence implemented from left to right.
The operation indices correlate with operation positions; one operation index has a fixed start position and finish position, and the two positions cannot be changed in GA implementation.The genetic operators are only used for the operation index.Figure 7 shows an example of chromosome representation.The chromosome represents a handling task with five loading operations, which are handled in the sequence of 4 → 1 → 3 → 2 → 5 .

9
the approximate optimal outbound container loading sequence.
In genetic algorithm implementation, the chromosome representation and genetic operator design are two vital steps for quickly accessing the feasible space and an effective movement toward the optimal solution.The proposed genetic algorithm is illustrated in Figure 6 and the main implementation steps are elaborated as follows.

Chromosome Representation
The chromosome is designed by permutation encoding.Each chromosome represents a possible loading sequence for the outbound containers and consists of T  number of genes.Each gene represents a loading operation, and the value of the gene is the indices of this loading operation.The gene sequence is the loading sequence implemented from left to right.
The operation indices correlate with operation positions; one operation index has a fixed start position and finish position, and the two positions cannot be changed in GA implementation.The genetic operators are only used for the operation index.Figure 7 shows an example of chromosome representation.The chromosome represents a handling task with five loading operations, which are handled in the sequence of

Evaluation of Fitness Value
According to the objective function mentioned above, the fitness function can be designed as follows.The fitness value of the chromosome could be calculated based on the fitness function.

Genetic Operators Design
For the selection operator, we chose roulette wheel selection, which selects potentially useful chromosomes for recombination based on the selection probability of each individual chromosome.

Crossover operator
The genes in the chromosome represent indices of all loading operations in the handling task, so the gene values cannot be lost and repeated in the offspring.To avoid infeasibility of the offspring generated, an order crossover operator is adopted.The crossover process is shown as follows.
(i) Firstly, we randomly select a segment in the parents, and place the segment from one parent to the front of the other parent.(ii) In order to avoid duplication, we scan the offspring from left to right and delete the repeated gene values in the substring.(iii) The crossover operation is finished when the offspring generated by the order crossover operator are both feasible.
A sample of the crossover operation is shown in Figure 8.

10
According to the objective function mentioned above, the fitness function can be designed as follows.The fitness value of the chromosome could be calculated based on the function.

Genetic Operators Design
For the selection operator, we chose roulette wheel selection, which selects potentially useful chromosomes for recombination based on the selection probability of each individual chromosome.

Crossover operator
The genes in the chromosome represent indices of all loading operations in the handling task, so the gene values cannot be lost and repeated in the offspring.To avoid infeasibility of the offspring generated, an order crossover operator is adopted.The crossover process is shown as follows.
(i) Firstly, we randomly select a segment in the parents, and place the segment from one parent to the front of the other parent.
(ii) In order to avoid duplication, we scan the offspring from left to right and delete the repeated gene values in the substring.
(iii) The crossover operation is finished when the offspring generated by the order crossover operator are both feasible.
A sample of the crossover operation is shown in Figure 8.

Mutation operator
As mentioned above, for ensuring the integrity of the handling tasks, missing and duplicating gene values must be forbidden in the offspring.Thus, an inversion mutation operator is used.The operator randomly picks two mutation points in the chromosome and then inverts the gene values of the mutation points.This mutation operator can guarantee the generation of a feasible chromosome.A sample of this mutation operation is shown in Figure 9.

Mutation operator
As mentioned above, for ensuring the integrity of the handling tasks, missing and duplicating gene values must be forbidden in the offspring.Thus, an inversion mutation operator is used.The operator randomly picks two mutation points in the chromosome and then inverts the gene values of the mutation points.This mutation operator can guarantee the generation of a feasible chromosome.A sample of this mutation operation is shown in Figure 9.

Stopping Criterion
The algorithm iteration will stop when the iterative times meet the specified value.Then, the best individual from the current population will be reported as the final solution.

Stopping Criterion
The algorithm iteration will stop when the iterative times meet the specified value.Then, the best individual from the current population will be reported as the final solution.

Computational Experiments
To verify the proposed approach for container loading optimization, several experiments were performed.Computational experiments are applied to the data from a specific Chinese rail-truck intermodal terminal.For evaluating the performance of our approach, a comparison is made between our approach and the current approach applied in rail-truck intermodal terminals.The current approach is straightforward and does not involve any optimization techniques.The outbound containers are orderly loaded onto rail vehicles from left to right in the departure line.Furthermore, to further evaluate the performance of our approach, experiments on different sample sizes are conducted.
The details of the specific rail-truck intermodal terminal are depicted as follows.One departure line can stop 120 rail vehicles for twenty-foot equivalent The outbound container yard has three rows and 120 bays, and the maximum layer number is two.There are four rail-mounted gantry cranes in the main operation area, and each crane is responsible for 30 bays.The maximum operations A  For the computational sample above, experiments are conducted for 100 independent runs.Based on the preliminary tests, the algorithm parameters are set as follows: population size is 50, crossover rate is 0.95, mutation rate is 0.15, and the maximum number of generations is 200.Then, we chose empty move time and reshuffling times as evaluation indicators, and made a comparison between our approach and the current approach.The comparison results are shown in Table 4.There are five GAP shown in Table 4. GAP a shows the optimization of the total handling time between our approach and the current approach.GAP b shows the optimization of the reshuffling times between our approach and the current approach.GAP c shows the optimization of the total energy consumption between our approach and the current approach.GAP d presents the percentage of empty move energy consumption in total energy consumption.GAP e presents the percentage of reshuffling energy consumption in total energy consumption.These GAP are calculated as follows.
GAP b = (T r obtained f rom CA − T r obtained f rom OA)/T r obtained f rom CA (13) GAP e = E r obtained f rom OA/E W obtained f rom OA In Table 4, the total energy consumption E W , energy consumption of an empty move E e and energy consumption of reshuffling E r , are calculated as follows: where C p is the power consumption per hour of crane usage and C r is the average power consumption per reshuffling.In this paper, C p is set at 300 kW/h and C r is set at 15 kW.As observed in Table 4, the GAP of total handling time is 7.1% and the GAP of reshuffling times is 60%.The reshuffling times obtained by our approach result from the overlapping of upper containers, which do not belong to the same handling task.The reshuffling operations caused by the overlapping between containers in the same handling task are optimized by our approach.Our approach observably reduces the total handling time and container reshuffling times in the outbound container loading process.There is a 22% decrease in total energy consumption after loading optimization by our approach.
In order to further evaluate the performance of our approach, experiments on different sample sizes are conducted.In general, there are 2-4 rail-mounted gantry cranes in the main operation area of this Chinese rail-truck intermodal terminal.Therefore, the quantitative range of loading operations in one handling task is from 30 to 60.In this range, computational experiments are carried out to verify the effectiveness of our approach for different sample sizes.For each task size, experiments are implemented for 100 independent runs.The results are shown in Table 5.As observed in Table 5, handling tasks with different sample sizes can be solved in a short time.The total computational time is about 1-2 minutes to generate a loading sequence for a container train combined with 120 rail vehicles, which can meet the actual scheduling requirements in rail-truck intermodal terminals.In Figure 10, based on optimizing the outbound container loading sequence, the total handling time, reshuffling times and energy consumption of the handling task are decreased.
13 is 60%.The reshuffling times obtained by our approach result from the overlapping of upper containers, which do not belong to the same handling task.The reshuffling operations caused by the overlapping between containers in the same handling task are optimized by our approach.Our approach observably reduces the total handling time and container reshuffling times in the outbound container loading process.There is a 22% decrease in total energy consumption after loading optimization by our approach.
In order to further evaluate the performance of our approach, experiments on different sample sizes are conducted.In general, there are 2-4 rail-mounted gantry cranes in the main operation area of this Chinese rail-truck intermodal terminal.Therefore, the quantitative range of loading operations in one handling task is from 30 to 60.In this range, computational experiments are carried out to verify the effectiveness of our approach for different sample sizes.For each task size, experiments are implemented for 100 independent runs.The results are shown in Table 5.As observed in Table 5, handling tasks with different sample sizes can be solved in a short time.The total computational time is about 1-2 minutes to generate a loading sequence for a container Based on the experimental results mentioned above, our approach is different from previous studies from a sustainability perspective, as we consider not only improving handling efficiency, but also reducing energy consumption in container loading operation optimization.The high energy efficiency of operations in rail-truck intermodal terminals is significant for sustainable development of freight transportation and contribute to energy conservation and emission reduction.

Conclusions
The container loading optimization problem in rail-truck intermodal terminals was considered in this paper.By analyzing the energy consumption of loading operations in the terminal, the key issue to improve the energy efficiency of outbound container loading in rail-truck intermodal terminals was to optimize the loading sequence for reducing the energy consumption from empty crane movements and reshuffling crane movements.An optimization model was presented to minimize total handling time and container reshuffling times in container loading tasks.Based on a GA, computational experiments on data from a specific rail-truck intermodal terminal in China were implemented.The results showed that optimizing the loading sequence can significantly reduce the total handling time, reshuffling times and energy consumption, and had a good performance for different sizes.The proposed approach can be useful and practical to help the operators in rail-truck intermodal terminals make an outbound container loading plan.

Figure 1 .
Figure 1.Schematic representation of the main operation area.

Figure 2 .
Figure 2. A sample of a container loading task.

Figure 1 .
Figure 1.Schematic representation of the main operation area.

Sustainability 2019 , 16 4
11,  x FOR PEER REVIEW 4 of container yard, outbound container yard and rail-mounted gantry cranes.Figure1gives a schematic representation of the main operation area in a Chinese rail-truck intermodal terminal.

Figure 1 .
Figure 1.Schematic representation of the main operation area.

Figure 2 .
Figure 2. A sample of a container loading task.

Figure 2 .
Figure 2. A sample of a container loading task.

Figure 3 .
Figure 3.A sample of loading operation.

Figure 3 .
Figure 3.A sample of loading operation.

Figure 4 .
Figure 4.An encoding example of the X/Y direction.

Figure 5 .
Figure 5.A sample of coordinate frame transforming.

Figure 4 .
Figure 4.An encoding example of the X/Y direction.In the Z direction, the layers are orderly encoded from low to high.Because all handling positions in the departure line and truck operation lane have only one layer, the z-value of these positions are all 1.A sample of coordinate frame transforming is shown in Figure5.

Figure 4 .
Figure 4.An encoding example of the X/Y direction.

Figure 5 .
Figure 5.A sample of coordinate frame transforming.

Figure 5 .
Figure 5.A sample of coordinate frame transforming.

Figure 7 .
Figure 7.A sample of chromosome representation.

Figure 7 .
Figure 7.A sample of chromosome representation.Figure 7. A sample of chromosome representation.

Figure 7 .
Figure 7.A sample of chromosome representation.Figure 7. A sample of chromosome representation.

Figure 8 .
Figure 8.A sample of the crossover operation.

Figure 9 .
Figure 9.A sample of the mutation operation.

Figure 9 .
Figure 9.A sample of the mutation operation.

Figure 10 .
Figure 10.GAP for different sample sizes.

Figure 10 .
Figure 10.GAP for different sample sizes.

Table 1 .
Details of energy consumption.

Table 3 .
task with a sample size of 30.

Table 4 .
Comparison results between our approach and the current approach.Notes: T W is total handling time, T r is reshuffling times, E W is total energy consumption, E e is energy consumption of an empty move, E r is energy consumption of reshuffling.a, b, c, d, e are the indices of GAP.

Table 5 .
Performance for different sample sizes.

Table 5 .
Performance for different sample sizes.