Large Neighborhood Search for LNG Inventory Routing
|
|
- Mariah Carter
- 5 years ago
- Views:
Transcription
1 Large Neighborhood Search for LNG Inventory Routing Vikas Goel ExxonMobil Upstream Research Company Kevin C. Furman ExxonMobil Upstream Research Company Jin-Hwa Song formerly at ExxonMobil Corporate Strategic Research Amr S. El-Bakry ExxonMobil Production Company Abstract Liquefied Natural Gas (LNG) is steadily becoming a common mode for commercializing natural gas. Due to the capital intensive nature of LNG projects, the optimal design of LNG supply chains is extremely important from a profitability perspective. Motivated by the need for a model that can assist in the design analysis of LNG supply chains, we address an LNG inventory routing problem where optimized ship schedules have to be developed for an LNG project. In this paper, we present an arc-flow formulation based on the MIP model of Song and Furman [17]. We also present a set of construction and improvement heuristics to solve this model efficiently. The heuristics are evaluated based on a set of realistic test instances that are very large relative to the problem instances seen in recent literature related to this problem. Extensive computational results indicate that the proposed methods are computationally efficient in finding optimal or near optimal solutions and are substantially faster than state-of-the-art commercial optimization software. 1 Introduction Natural gas is expected to be the world s fastest growing fossil fuel with consumption increasing at an average rate of 1.6% per year from 2008 to 2035 [1]. In many cases, large reserves of natural gas are located in areas of little or no local demand. As a result, natural gas must be transported over long distances, either via pipelines or shipped as liquefied natural gas (LNG) in specially designed ships. To be transported as LNG, natural gas produced from a reservoir is first processed to remove impurities. The purified natural gas is liquefied by Corresponding author; vikas.goel@exxonmobil.com 1
2 cooling it to 163 C at a liquefaction terminal (also referred to as a production terminal). The LNG is stored at dedicated storage facilities at the production terminal before being loaded on LNG ships that transport the LNG to regasification terminals (also referred to as regas terminals). At the regas terminal, LNG is stored before being regasified (i.e., vaporized) and injected into the local natural gas pipeline grid. For a more thorough introduction to the LNG business, the reader is referred to Tusiana and Shearer [19]. LNG accounts for a growing share of world natural gas trade with world natural gas liquefaction capacity expected to double between 2009 and 2035 [1]. However, new LNG production projects are complex and extremely capital intensive. As a result, the optimal design of the supply chain of an LNG project is extremely important from a profitability perspective. From an operations perspective, managing an LNG project involves negotiating a delivery schedule for each customer on an annual basis. This Annual Delivery Program (ADP) is negotiated to best accommodate the expected requirements of each party in a given year. Due to the negotiated nature of the schedule, the ADP can impose inefficiencies in the utilization of the shipping fleet and the overall LNG supply chain. In this paper we present an optimization model that optimizes the LNG shipping schedule and inventory management for an LNG project. The main motivation behind this model is to provide an objective approach for the analysis of key LNG supply chain design decisions such as LNG fleet size and composition, and terminal infrastructure during the concept and detailed design phases for new LNG projects. The proposed model can also be used to develop ADPs in operational LNG projects and to quantify the inefficiency in a schedule caused by the negotiation process. The model can also be used to conduct what-if analysis to test how robust a schedule is to disruption events. LNG inventory routing can be considered as a special case of maritime inventory routing problems (MIRP). MIRP combines inventory management and ship routing, which are typically treated separately in industrial practice. The reader is referred to the papers by Christiansen et al. [8, 7] and Andersson et al. [3] for a thorough review of maritime inventory routing literature. Christiansen and Fagerholt [5] define a basic maritime inventory routing problem as the transportation of a single product that is produced at loading ports and consumed at unloading ports where each port has a given inventory storage capacity and a production or consumption rate. The number of visits to a port and the quantity of product to be loaded or unloaded are not predetermined. As discussed by Andersson et al. [3], nearly every paper concerning combined inventory management and routing addresses a new version of the problem. This is in contrast to the work on classical routing problems such as the vehicle routing problem (VRP) that have widely accepted definitions and assumptions. The LNG inventory routing problem (LNG IRP) studied in this article is based on a real-world application and shares the fundamental properties of a single product MIRP. However, the LNG IRP addressed here includes several variations including variable production and consumption rates, LNG specific contractual obligations, and berth constraints. Most importantly, the LNG IRP seeks to generate schedules where each ship makes several voyages over a time horizon with both the number of voyages and the time horizon being considerably larger than those considered by a typical MIRP. Andersson et al. [4] provide an excellent overview of the business cases and common characteristics for LNG inventory routing. Grønhaug and Christiansen [12] develop discrete time arc-flow and path-flow models for 2
3 operational LNG inventory routing and test these models on a number of operational planning cases with time horizons up to 60 days. Grønhaug et al. [13] develop a branch-and-price method for the discrete time path-flow formulation of the LNG IRP and apply it to an extended set of test cases with time horizons up to 75 days with improved computational performance over the previous effort [12]. While the above work focuses on optimizing inventory management at production terminals only together with the ship schedules, Fodstad et al. [9] develop a discrete-time arc-flow model for LNG IRP that also considers optimizing decisions related to sales downstream of the regas terminals. Although atypical of the LNG industry, they also include the capability for split loads in their model. The model presented by Fodstat et al. is implemented in the LNGScheduler software [10] and shows improved run times compared to those reported by Grønhaug et al. [13] over the same test set. Fodstad et al. [10] apply this model to a set of instances with a time horizon of 181 days. There has recently been some development of software applications for combined inventory management and ship scheduling for the LNG supply chain. van de Broecke and Adams [20] and Stremersch et al. [18] describe the design of such decision support tools. Fodstad et al. [9, 10] discuss the model, implementation and some computational testing for LNGScheduler. The paper by Rakke et al. [15] appears to be a first attempt to address problems of developing ADPs for large LNG projects. In this paper an arc-flow model and rolling horizon heuristic method are applied to problems with up to 46 ships and a one year planning horizon. Halvorsen-Weare and Fagerholt [14] develop a decomposition-based heuristic method to solve an arc-flow model of the the LNG IRP, and address instances with a moderate number of ships and terminals with a time horizon ranging from 30 to 360 days. While these authors have addressed larger problems than previously reported in literature, their work limits optimization of ship schedules together with inventory management at the production terminals only. In this paper, we address large scale LNG IRPs where ship schedules are optimized together with inventory management on both the production and regas terminals. This approach provides the ability to analyze supply chain designs and develop ship schedules from a broader general-interest perspective. Christiansen et al. [6] develop a greedy construction heuristic and a solution improving genetic algorithm for an MIRP from the cement industry. They note that efficient neighborhood operators for improvement heuristics remain a challenge. The methods described in this paper are inspired by the fundamental modeling and algorithmic concepts developed by Song and Furman [17] for the MIRP problem defined by Furman et al. [11]. The two ship neighborhood operator proposed by Song and Furman [17] is inspired by the classical 2-opt neighborhood search where the shipping schedules of a pair of ships are optimized at each iteration of the search. In this article, we study the application of the two-ship neighborhood operator and a new neighborhood operator to extremely large-scale instances of the LNG inventory routing problem. We also propose approaches to improve the efficiency of the two-ship search operator. Although the solution techniques presented here have been developed and tested in the context of this specific LNG IRP, these can be generally applied to a wide variety of combined inventory management and routing problems. For a review of very large neighborhood search (VLNS) methods, the reader is referred to the excellent survey by Ahuja et al. [2]. The rest of the paper is organized as follows. The formal problem description is introduced 3
4 in Section 2. Section 3 presents a mixed-integer programming model for the problem. Section 4 discusses various solution techniques used in the large scale neighborhood search procedure. In Section 5, computational results are presented for various algorithms described in this paper. Section 6 illustrates how various practical analyses are enabled by such a model. Finally, in Section 7 the conclusions and future research directions are discussed. 2 Problem Description We consider an LNG IRP from the perspective of an vertically integrated oil and gas company. The goal of the problem is to find an optimal schedule for a heterogeneous pool of ships that deliver LNG from a set of production terminals to a set of regasification terminals such that constraints related to inventory storage, port operations and contractual obligations are satisfied. While the problem aims at finding an optimal schedule for a given supply chain design that specifies fleet composition and size, terminal storage and berth facilities, the main motivation for solving this problem is based on a need to develop a model that can assist in analyzing alternate supply chain designs. Formally, the problem considers a set of production and regas terminals with fixed storage capacities. Each production terminal has a given production profile. On the other hand, operations at the regas terminals are more flexible such that regasification rates can be adjusted within specified bounds on a daily basis. Profiles for the minimum and maximum regas rates over the planning horizon are specified for each regas terminal. In addition, total LNG demand over the planning horizon at each regas terminal is specified. Finally, each terminal has a limited number of berths to load and unload cargoes. Each loading or unloading activity spans over one time period and occupies a berth for that duration. The supply chain includes a set of heterogeneous ships. We restrict our attention to full load, full discharge problems. In other words, we consider the case where each ship will fully load at a single production terminal and fully discharge at a single regas terminal (i.e., there is no partial unloading a.k.a. split delivery). Several previous papers have considered split delivery of LNG cargoes [9, 10, 12, 13, 4]. From a planning and design analysis perspective, this assumption is reasonable since this is typically considered the preferred mode of operation within the LNG industry. Due to its low boiling point, the LNG loaded on a ship gradually boils-off (vaporizes) while the ship travels or while it waits at a terminal. The boil-off gas is commonly burnt as fuel and causes the volume delivered to be less than the volume loaded on the ship. We refer to the volumes loaded and discharged by a ship as the load and discharge capacities, respectively, for the ship. The difference between these capacities represents the average boil-off loss based on a typical voyage between two terminals. While actual boil-off amounts experienced during operations may differ from the above approximation, we consider our approximation reasonable since boil-off amounts represent a small fraction of the total cargo size and minor deviations from this approximation are not expected to significantly impact long term supply chain design analysis. We consider the terminal storage and contractual demand constraints to be soft constraints. In other words, production in excess of storage capacity at a production terminal is considered as lost production. Similarly, regasification in excess of existing inventory lev- 4
5 els at a regas terminal is considered as a stockout where the excess has to be purchased from the spot market. The objective is to minimize the total penalties associated with lost production, stockouts, and unmet contractual demands at regas terminals. This problem has several complicating characteristics including: heterogeneous fleets, seasonal travel times, limited berth availability at terminals, seasonality in production rates, and variable regas rates. As mentioned above, a solution to the above problem can support key capital decisions regarding the sizing of terminal infrastructure and fleets. A solution to the above problem can also be used to evaluate the effects of pooling shipping fleets, to quantify key operations metrics including fleet utilization, and evaluate the ability of a supply chain design to meet contractual obligations. In defining the mathematical model in the next section, all of the characteristics of the problem description are fully developed. 3 Model The proposed mathematical model is formulated as a mixed integer programming (MIP) problem. 3.1 Sets and Parameters We first define a number of sets and parameters that are used in the mathematical model. L R : Production terminals. L = {1, 2,..., L }. : Regas terminals. R = {1, 2,..., R }. J : All terminals. J = L R. V : Vessels (Ships). v = {1, 2,..., V }. J v : Set of terminals that ship v can visit. T : Planning horizon. T = {1, 2,..., T } c j : Storage capacity at terminal j. c v,j : Volume loaded (discharged) by ship v at production (regas) terminal j. D r : Demand for LNG over planning horizon at regas terminal r. b j : Number of berths at terminal j. p l,t : Production rate during time period t at production terminal l. I 0 j : Initial inventory at terminal j. j 0 v : Initial destination terminal for ship v. t 0 v : Arrival date at initial destination terminal for v. d L r,t : Minimum regas rate of regas terminal r. d U r,t : Maximum regas rate of regas terminal r. τ j,j,t : travel time (including berth time) from terminal j to terminal j for voyage starting at time t. w j,t : Weight for penalizing lost production (stockout) at production (regas) terminal j. 5
6 w D r : Weight for penalizing unmet demand at regas terminal r. 3.2 Network This section describes the time-space network formulation for the LNG IRP described above. This network model is closely related to the model of Song and Furman [17]. As presented by Savelsbergh and Song [16], a time-space network formulation can be viewed as an integer multi-commodity network flow formulation where a ship is a commodity and a node represents a possible visit to a terminal at a particular time. Figure 1: Example of Time-Space Network Structure The time-space network in Figure 1 includes source (SRC) and sink (SNK) nodes to represent initial and final locations for the ships. The network also includes nodes (referred to as regular nodes) that represent each terminal at a given time period. For each ship v, there are five types of arcs; an arc from the SRC to SNK node represents that a ship is not utilized; an arc from SRC to a regular node represents the arrival of a ship to its initial destination; an arc from a regular node to SNK represents the final departure of the ship; waiting arcs allow a ship to wait at a terminal without occupying a berth. These arcs connect regular nodes corresponding to the terminal at successive time periods. Finally, a travel arc from a regular node n 1 to a regular node n 2 represents the loading (or unloading) activity at the terminal corresponding to node n 1 immediately followed by travel to the terminal corresponding to node n 2. If the schedule for a ship has enough slack, the ship can either wait at a production terminal after loading a cargo, or travel to the destination regas terminal and wait before it unloads a cargo, or split the waiting time between the two terminals. From a ship scheduling perspective, these are symmetric solutions. Using a travel arc to represent a loading or unloading activity immediately followed by travel eliminates this symmetry. The formal node and arc definitions are presented below. N : Set of nodes. N = {(j, t) j J v } {SRC} {SNK} 6
7 N : Set of regular nodes. N = N\{SRC, SNK}. A φ v : Arc from SRC node to SNK node for ship v. = {v, (SRC), (SNK)} A I v : Set of arcs from source to initial destination for ship v. = {v, (SRC), (j 0 v, t 0 v)} A D v : Set of waiting arcs for ship v. = {v, (j, t), (j, t + 1) j J v } A T v : Set of travel arcs for ship v. = {v, (l, t), (r, t + τ l,r,t ) l J v L, r J v R} {v, (r, t), (l, t + τ r,l,t ) l J v L, r J v R} A F v : Set of arcs from regular nodes to SNK node for ship v. = {v, (j, t), (SNK) j J v } A v : Set of arcs for ship v. = A φ v A I v A D v A T v A F v A : Set of arcs. A = v A v δ n + : Set of outgoing arcs from node n. δ n : Set of incoming arcs to node n. Note that to ensure consistency, we assume that the travel times satisfy the following condition τ j,j,t τ j,j,t t T, j, j J, j j This condition ensures that the network does not allow a ship to arrive at its destination terminal earlier by delaying its departure Decision Variables A feasible schedule for this LNG IRP specifies whether a ship travels on any particular arc within the network model. The solution also specifies the regas rates at each regas terminal during each time period. All of the other variables can be considered state variables as their values will be a direct result of the binary variables on the arcs and the regas rates. The following variables have been defined in the mixed integer programming formulation presented below. I j,t : Inventory level at terminal j at the end of time period t. d r,t : Regas rate during time period t at regas terminal r. o l,t : Lost production at production terminal l during time period t. s r,t : Stockout at regas terminal r during time period t. δ D r : Unmet demand at regas terminal r during planning horizon; x a : Binary variable for arc a. 7
8 3.3 MIP Formulation The mixed integer programming formulation which we will refer to as model P follows. min w l,t o l,t + w r,t s r,t + wr D δr D (1) (l,t) N l L (r,t) N r R r R s.t. x a x a = 0, v V, n N, (2) a A v δ n + a A v δ + SRC a A v δ n x a = 1, v V, (3) x a = 1, v V, (4) a A v δ SNK I l,t = I l,t 1 + p l,t v I r,t = I r,t 1 d r,t + v a (A T v AF v ) δ+ (l,t) c v,l x a o l,t l L, t T (5) c v,r x a + s r,t r R, t T, (6) v δ D r a (A T v A F v ) δ n + D r t a (A T v AF v ) δ+ (r,t) x a b j, j J, t T, n = (j, t) N (7) v a (A T v AF v ) δ+ (r,t) c v,r x a, r R (8) d L r d r,t d U r, r R, t {1, 2,..., T }, (9) 0 I j,t c j, j J, t T (10) o l,t 0 l L (11) s r,t 0 r R (12) δ D r 0 r R (13) x a {0, 1} a A (14) The objective is to minimize the sum of weighted lost production, stockout, and unmet demands. Constraints (2)-(4) are network-flow conservation constraints for each ship. Constraints (5) and (6) model the inventory balance at the production and regas terminals, respectively. Lost production and stockout variables provide slack for these constraints. Constraint (7) ensures that the number of ships loading or unloading a cargo at a terminal within a time period does not exceed the number of berths. Note that a ship can wait at a terminal without occupying a berth. Constraint (8) enforces a lower bound on the unmet demand variable based on total deliveries scheduled for a regas terminal. Since the objective function seeks to minimize the total unmet demand, this constraint will be tight in every optimal solution that has a positive unmet demand. The final set of constraints (9)-(14) ensures that all the variables satisfy their specific bounds. It should be noted that for the above model to calculate the lost production and stockout amounts correctly, the penalty parameters for lost production and stockouts should be monotonically decreasing in time, i.e., w j,t > w j,t+1. This ensures that a solution will not involve lost production (stockout) until the inventory reaches capacity (falls to zero). In other words, without this condition the model could report solutions where lost production 8
9 or stockouts are reported in advance of the actual event. For example, consider a solution where production terminal l experiences a loss of ɛ units in time period t. If w l,t 1 = w l,t, an equivalent solution can be constructed where the loss is distributed over periods t and t 1. This would be unrealistic if the terminal inventory at t 1 is lower than the storage capacity. A monotonically decreasing penalty weight ensures that the lost production and stockout calculations are correct. In its present form, the model is generic with respect to the length of the planning horizon and the time discretization. Based on planning conventions in the LNG industry and typical loading and discharging durations, a year-long planning horizon with one day time discretization is expected to be a reasonable scale for the problem. An important characteristic of this formulation is that it is easy to accommodate various practical constraints of the problem. Many of these constraints can be handled by adding, removing, and/or fixing nodes and arcs. For example, planned maintenance at a terminal that disallows any loading or unloading activities can be accounted for by removing all travel arcs from that terminal for the duration of the activity. Due to the complex nature of the above problem, commercially available Mixed Integer Programming solvers cannot solve our real world problems within a reasonable amount of time. In the following section, we discuss some of the solution methods developed in order to overcome this obstacle. 4 Solution Methods We develop local neighborhood search heuristics with the goal of finding good solutions to the proposed model in a short amount of time. This approach can also be used to provide a warm-start for an exact solution method. Algorithm 1 proposes a three-step heuristic method for the above model. The first step of the heuristic uses a construction heuristic to build an initial feasible solution to the model. The construction heuristic builds a feasible solution by iteratively scheduling ship departures for the entire fleet on a daily basis. Ship departures on any given day are scheduled using an urgency-based greedy heuristic. Two improvement heuristics that employ different neighborhood structures are then used in sequence to improve this solution. The Time-Window Improvement Heuristic seeks to improve the solution obtained from the construction heuristic by allowing for departure dates for each voyage in the current solution to change within a small time-window defined around the current departure dates. The Two-Ship Improvement Heuristic seeks to improve the existing solution by allowing a subset of the ships to change schedules over the entire horizon. The subproblems generated by both improvement heuristics can be solved as Mixedinteger programs (MIPs) and solved using a commercial solver. Algorithm 1 Solution Heuristic Build initial solution using Construction Heuristic (CH) Improve solution using Time window improvement heuristic (TWIH) Improve solution using Two-ship improvement heuristic (2SIH) 9
10 4.1 Construction Heuristic (CH) The construction heuristic is based on a greedy approach that seeks to deliver as much LNG as possible and as soon as possible from the production terminals to the regas terminals. In addition, the construction heuristic seeks to minimize lost production and stockouts by prioritizing voyages to terminals that have the most urgent demand for shipping capacity during the next few time periods. The heuristic uses this approach to schedule ship departures at all terminals on a daily basis. A feasible solution is generated by fixing daily ship departure decisions iteratively over the entire planning horizon. The construction heuristic (shown in Algorithm 2) begins by initializing the terminal inventory levels at the beginning of the planning horizon, and setting the initial arrival dates and locations for each ship. The heuristic builds a solution where the daily regas rate at any regas terminal is equal to the mean regas rate for the terminal on the respective day. The construction heuristic then iterates over the time periods to fix ship departures from each terminal in any given time period. The heuristic uses a similar approach for scheduling departures from production and regas terminals. We present the explanation for scheduling departures from production terminals. Variations for scheduling departures from regas terminals are presented parenthetically. In order to schedule departures at terminal k during time period t, the heuristic first calculates the closing inventory level at terminal k for time period t using equation (5) (equation (6) for regas terminals) if no ship departures were to take place from terminal k during time period t. The closing inventory is used to identify the set of ships V k that can be dispatched from terminal k during time period t. The set V k includes ships that are located at production (regas) terminal k at period t and that can be loaded (discharged) based on their load (discharge) capacities and the closing inventory at terminal k in period t. The heuristic schedules departures of ships in set V k from production (regas) terminal k as long as a berth is available at terminal k, and the terminal inventory permits loading (discharging) some ship in V k. In scheduling the next departure, the heuristic first identifies the destination terminal with the most urgent need for shipping capacity. The need for shipping capacity at a production (regas) terminal is quantified as the difference between the total planned production (regasification) at the terminal during a future time horizon, and the total shipping capacity already scheduled to reach that terminal during that time horizon. The terminal ˆk with the most urgent need for shipping capacity is selected as the destination for the next departure from terminal k. Note that not all ships may be allowed to travel to all terminals. Further, our model restricts ship to travel from a production (regas) terminal to a regas (production) terminal only. Therefore, the selection of the next destination is restricted to terminals such that at least one ship in V k can travel from terminal k to terminal ˆk. In order to maximize the delivery of LNG, the heuristic selects the ship ˆv with the largest loading (discharging) capacity among all ships in V k that can travel from production (regas) terminal k to terminal ˆk. The heuristic then schedules a departure for ship ˆv from terminal k to terminal ˆk in period t and the closing inventory level at terminal k is updated accordingly. This approach is repeated until no more departures can be scheduled from terminal k in period t. A complete schedule for the entire time horizon is built by applying this approach for all terminals over all time periods. 10
11 Finally, lost production (stockout) volumes for all production (regas) terminals can be calculated based on the production (regasification) profiles together with the shipping schedule. Algorithm 2 Construction Heuristic Initialize ship locations and terminal inventory levels Fix daily regas rates to dl r,t +du r,t 2 for all regas terminals r for t = 1, 2,..., T do for all terminals k do Evaluate closing inventory level for period t assuming no ship departures at terminal k during period t V k := Set of ships that are located at terminal k in period t and that can be dispatched based on ship capacity and terminal inventory while another berth is available and V k is not empty do K := Set of terminals k such that at least one ship in V k can be dispatched from terminal k to terminal k Select destination terminal ˆk K that has the most urgent demand for shipping capacity Select ship ˆv that has the largest capacity among ships that can be dispatched from k to ˆk Schedule departure of ship ˆv from terminal k to ˆk in period t Update closing inventory for period t based on scheduled departure end while Evaluate stockout or lost production at k at t end for end for 4.2 Time-Window Improvement Heuristic (TWIH) The time-window improvement heuristic is a large neighborhood search method. This heuristic searches in the neighborhood consisting of all solutions in which any voyage departure is delayed or advanced by at most m days, compared to the departure date for that voyage in the input solution. An optimal solution in this neighborhood can be obtained by solving a modified version of Model P that excludes travel arcs for all voyages that depart outside the time windows defined by the above neighborhood. This modified version of model P can be solved using a commercial MIP solver. 4.3 Two-Ship Improvement Heuristic (2SIH) The two-ship improvement heuristic (described in Algorithm 3) is based on the variable neighborhood search method. At each iteration, the heuristic selects a pair of ships and defines the search neighborhood to consist of the set of feasible schedules for these ships. 11
12 Schedules for all other ships are considered fixed. A subproblem is solved to optimize schedules within this neighborhood. The algorithm iteratively selects pairs of ships and solves the corresponding subproblems for the defined neighborhoods. The heuristic is terminated when a better solution is not found for ( ) V 2 successive iterations. Since the objective function for this problem is lower bounded by zero, the heuristic can also be terminated as soon as a solution with a zero objective value is obtained. From an implementation perspective, the optimal solution for a subproblem induced by a pair of ships can be obtained by solving a reduced version of Model P where variables for all other ships are fixed based on the current solution. We refer to the reduced model as a two-ship subproblem. Algorithm 3 Two Ship Improvement Heuristic Initialize solution; U B := objective value of solution count := 0 while count ( ) V 2 and UB > 0 do Select ships v and v Fix schedules for all ships other than v and v Solve two-ship subproblem if better solution is found then Update solution; U B := objective value of solution count := 0 else count := count + 1 end if end while Ship-pair Selection The efficiency of the two-ship improvement heuristic depends on the sequence of two-ship subproblems solved during the algorithm. By selecting ship-pairs that are most likely to lead to a better solution, the number of subproblems solved and the overall CPU time can be minimized. We present three schemes for selecting two-ship subproblems during the two-ship improvement heuristic. Lexicographic selection For a problem where the fleet of ships is represented by the ordered set { } v 1, v 2, v 3,..., v V, we consider the selection scheme to be lexicographic if it selects two-ship subproblems in the following order: (s 1, s 2 ), (s 1, s 3 ),..., (s 1, s V ), (s 2, s 3 ),..., (s V 1, s V ). The lexicographic ordering scheme may be effective in a practical setting where the relative importance of ships is inherent in the order in which they are included in the fleet. Metrics based selection In this method, we use the current solution together with problem specific insights to identify the two-ship subproblem that is likely to achieve the largest improvement in objective function. Specifically, we define two metrics based on the current 12
13 solution for each ship-pair for which a subproblem can be solved next. These metrics represent estimates for the maximum potential reduction and for the likelihood of achieving that reduction in the objective function if the corresponding two-ship subproblem were optimized. The ship-pair that has the highest likelihood of achieving a reduction among the ship-pairs that have the largest potential for reduction is selected. Algorithm 4 provides a detailed specification for the selection method. Let πv 1 represent the improvement in the objective function if the schedule for ship v is re-optimized. We estimate πv 1 as the total production and stockouts that would be eliminated if ship v could be rescheduled such that the sequence of terminals visited by v remains unchanged but the ship arrives at all terminals as soon as possible; i.e., the re-optimized schedule does not involve any demurrage. The estimate for reduction in the objective function by solving the two-ship subproblem for ship-pair (v 1, v 2 ) (represented by Π 1 (v 1,v 2 )) is estimated as the greater of πv 1 1 and πv 1 2. We evaluate the likelihood of achieving a reduction in the objective function by solving the two-ship subproblem for ship-pair (v 1, v 2 ) assuming that the sequence of terminals visited by v 1 and v 2 in the re-optimized solution will remain unchanged compared to the current solution. However, departure dates for the two ships from a terminal might be swapped or re-ordered in order to reduce overall lost production and stockouts. In order to quantify this likelihood, we evaluate π(v 2 1,v 2 ) as the number of times in the current solution where the departure for ship v 1 is delayed due to the departure for v 2 from the same terminal. In evaluating π(v 2 1,v 2 ), we count cases where ship v 1 is forced to demur at a terminal while ship v 2 departs from that terminal. Cases where ship v 2 is the last ship to depart from a terminal before v 1 arrives at that terminal and is forced to demur are also included in the calculation of π(v 2 1,v 2 ). Finally, the likelihood of achieving a reduction in the objective function by solving the two-ship subproblem for ship-pair (v 1, v 2 ) (represented by Π 2 (v 1,v 2 )) is approximated as the number of times the ship that is estimated to lead to bigger reduction in objective function is delayed due to the other ship in the pair. Note that since Π 1 (v 1,v 2 ) and Π2 (v 1,v 2 ) are symmetric with respect to (v 1, v 2 ), an arbitrary ordering of ships is assumed such that v 1 < v 2 can be defined. The set P consists of ship-pairs (v 1, v 2 ) for which Π 1 (v 1,v 2 ) is within the top α percent of the range of Π 2 (, ). The ship-pair in P that has the largest estimated likelihood of achieving a reduction in the objective function is selected. The ship associated with the largest value for this metric is selected as v 2. Model based selection The restriction that at most two ships can change their schedules can be represented by constraints (15)-(17). Binary variable y v represents whether or not ship v can change schedule. Constraint 15 restricts at most two ships to change their schedule. Constraint (16) determines the value of y v based on whether the schedule for ship v in solution represented by variable x is different from that in a given solution x. y v 2 (15) v y v = 1 a δ v x a x a v (16) y v {0, 1} v (17) 13
14 Algorithm 4 Metrics based two-ship subproblem selection for all ships v do πv 1 := total lost production and stockouts that would be eliminated if all demurrage in current schedule for ship v were eliminated with the sequence of terminals visited by ship v remained unchanged end for for all ship pairs (v 1, v 2 ) such that v 1 < v 2 do Π 1 (v 1,v 2 ) := max { } πv 1 1, πv 1 2 end for for all ship pairs (v 1, v 2 ) do π(v 2 1,v 2 ) := number of times in current solution where departure for v 1 is delayed due to departure for v 2 from the same terminal end for for all ship pairs (v 1, v 2 ) such that v 1 < v 2 do Π 2 (v 1,v 2 ) := { π 2 (v1,v 2 ) : if π 1 v 1 π 1 v 2 π 2 (v 2,v 1 ) : if π 1 v 1 > π 1 v 2 end for Π 1,max := max (v 1,v 2 ) Π1 (v 1,v 2 ), Π1,min := min (v 1,v 2 ) Π1 (v 1,v 2 ) P := { (v 1, v 2 ) Π 1 (v 1,v 2 ) [ Π 1,max α 100 (Π1,max Π 1,min ), Π 1,max]} (v1, v2) := arg max (v 1,v 2 ) P Π2 (v 1,v 2 ) return (v1, v2) 14
15 Dimensions Cont. 0-1 Best known Problem ( L, R, V ) Variables Variables Constraints solution P1 (1,2,6) P2 (1,3,8) P3 (2,1,10) P4 (3,1,13) P5 (1,4,15) P6 (4,1,4) P7 (1,6,6) P8 (1,1,14) P9 (1,2,17) P10 (1,4,27) P11 (1,5,14) P12 (1,4,18) P13 (1,8,40) P14 (1,10,69) Table 1: Problem Instances We define a restricted version of model P (referred to as rp) by incorporating constraints (15)-(17). The solution to rp will identify the two-ship subproblem that will lead to the biggest improvement in the objective function relative to the existing solution x. The model based selection approach involves selecting the two ships that are associated with the largest values for y ( ) in the solution of the linear relaxation of model rp. For large problems, solving the linear relaxation for model rp can be expensive. It is easy to see that the constraint set for Model rp includes a network-substructure that is separable over the set of ships. However, the constraints representing terminal inventory and berth limits, together with constraints that limit the number of ships that can change schedule cause the model to be coupled over the set of ships. Decomposition methods, such as Dantzig-Wolfe or Lagrangean decomposition can be employed for solving the linear relaxation of rp efficiently when solution in the full space is inefficient. 5 Computational Results In this section, we present computational results to demonstrate the performance of the proposed heuristic. The results are based on an implementation of the model in AIMMS. All optimization models have been solved with CPLEX 11.1 on a Dell 7500 Windows PC with a 2.93 GHz dual-quad core processor and 24 GB RAM. 15
16 5.1 Problem Set Table 1 shows the sizes of the 14 problems on which the proposed heuristics were evaluated. The table reports the dimensions ( L, R, V ) for each problem where L, R and V represent the number of production terminals, regas terminals and ships, respectively. All the problems are defined over a one year planning horizon. These problems are realistic test cases allowing for a limited amount of pooling of the shipping fleet among multiple shipping lanes restricted via the allowed arcs in the network. The number of continuous and binary variables, the number of constraints and the best known solution for each problem are also reported. We classify the problems in to three categories based on the ease of finding good solutions using CPLEX. easy problems: includes P1, P2, P3, P4, P5. CPLEX generates the best known solutions for these problems in less than 1 CPU hour medium problems: includes P6, P7, P8, P9, P10. CPLEX generates a good (within 15% of best known) solution for these problems within 1-5 CPU hours hard problems: includes P11, P12, P13, P14. CPLEX cannot generate a good solution for these problems within 10 CPU hours 5.2 CPLEX vs. Proposed Heuristic We first compare the performances of CPLEX and the proposed heuristic. For this comparison, the results presented for the heuristic are based on randomly selecting the next two-ship subproblem to be solved at each iteration of the two-ship improvement heuristic (2SIH). The proposed heuristic is therefore referred to as H-Randm. To account for the effect of randomization, results for H-Randm are based on 100 runs for problems P1-P12 and 50 runs for problems P13-P14. Our implementation for the time-window improvement heuristic (TWIH) is based on defining a time-window that allows each ship departure to be delayed or advanced by at most four time-periods. Each two-ship subproblem is solved with a node-limit of 1000 nodes. The time-window improvement heuristic is solved with a node-limit of 5 nodes. Results for CPLEX are based on runs with CPU time limits of 5 hours, 10 hours and 20 hours for easy, medium and hard problems, respectively. Table 2 compares the solutions generated by CPLEX with those generated by H-Randm. For ease of exposition, we report the solution quality in terms of the percent gap of a solution with respect to the best known solution. Specifically, for a solution X and best known solution ˆX, X ˆX the gap is computed as ˆX 100. Since the best known solutions for problems P4 and P11 have a zero objective value, the solution quality for these cases is reported in terms of the actual objective values. Table 2 reports the percent gaps relative to the best known solution for the first feasible (FF) and the best feasible solutions (BF) found by CPLEX, together with the percent gap between BF and the best bound when CPLEX terminates. The table also reports percent gaps relative to the best known solution for the solutions found by H-Randm at the end of the three stages in the H-Randm; i.e., construction heuristic (CH), time-window 16
17 Problem CPLEX H-Randm FF BF CPLEX Gap CH CH + TWIH CH + TWIH + 2SIH Min Avg Max (%) (%) (%) (%) (%) (%) (%) (%) P P P P4* [ ] [0.00] 0.00 [399.17] [288.24] [0.00] [0.00] [0.00] P P P P P P P11* [ ] [ ] [504.34] [504.34] [0.00] [2.62] [130.87] P P P Table 2: Solution quality comparison: CPLEX vs H-Randm improvement heuristic (TWIH) and two ship improvement heuristic (2SIH). Since 2SIH involves randomization, the minimum, maximum and average percent gaps relative to the best known solution are reported at the end of 2SIH. Table 3 compares the CPU times required by CPLEX with those required by H-Randm. The CPU times reported in columns that are marked as scld are scaled by the average CPU time taken by H-Randm. Columns that report actual CPU times are marked as act. The table reports CPU times taken by CPLEX to generate the first feasible (TFF) and the best feasible (TBF) solutions, time required to generate a feasible solution better than the average solution for H-Randm (TRH), and the total run time for CPLEX. The table also reports the cumulative CPU time taken to complete the three stages of H-Randm, including the minimum, maximum and average CPU times for the overall heuristic. Table 2 shows that CPLEX can prove optimality for only 3 problems and generate the best known solution for 8 problems within the CPU time limits. We compare the performance of CPLEX with H-Randm for the three problem classes individually. For the easy problems, CPLEX is able to generate the best known solutions for all problems. For problems P1, P3 and P4, H-Randm generates solutions within 1.26% of the best known solution. However, CPLEX is much slower (2.55x to 70x) in generating its best solution than the average CPU time taken by H-Randm. Further, CPLEX is 1.78x-19.02x slower in generating a solution that is at least as good as that generated by H-Randm. For problems P2 and P5, H-Randm generates solutions with objective value 8.5% to 25.6% worse than the best known solution and CPLEX. However, CPLEX is much slower in gener- 17
18 Problem CPLEX H-Randm TFF TBF TRH Total CH CH + TWIH CH + TWIH + 2SIH Min Avg Max (scld) (scld) (scld) (act) (scld) (scld) (scld) (act) (scld) P P P P P P P P P P P P P P Table 3: CPU time comparison: CPLEX vs H-Randm ating its best solution (14.4x-25.3x) for these problems than H-Randm. Further, in the case of P2, CPLEX is 17.1x slower than H-Randm in generating a solution that is at least as good as that generated by H-Randm. Overall, it can be concluded that for the easy problems, CPLEX generates better solutions within the time limit provided but H-Randm generates good solutions much faster. Note that although the speed-up factors above are based on the average performance of H-Randm, similar conclusions can be drawn based on the worst case performance of H-Randm. For 3 out of 5 medium problems, both CPLEX and H-Randm are able to generate the best known solutions. However, CPLEX is 68.6x to 131.9x slower than H-Randm. In the case of P10, H-Randm always generates a better solution than CPLEX although CPLEX requires 3.3x more CPU time than the average run-time of H-Randm. Finally, the solutions reported by H-Randm for P7 demonstrate significant variance in quality (0.04%-26.27% gap) although 67% runs of H-Randm generates a better solution than CPLEX. For the hard problems, the best feasible solutions generated by CPLEX are clearly poor compared to the best known solutions (466.62% to % gap for P12, P13 and P14) and also much worse than the average (and worst) solutions generated by H-Randm. Further, CPLEX is slower (1.03x to 143.7x) in generating its best feasible solution than the average CPU time required the H-Randm. In fact, for 3 out of 4 problems, CPLEX is slower in generating its best feasible solution than the worst performance of H-Randm. Note that comparing TBF or TRH with the total CPU time for H-Randm gives a conservative estimate of the speedup delivered by H-Randm since the total CPU time for H-Randm 18
19 also includes the time required to prove convergence in addition to the time required to find the best solution. Table 2 also shows that CH and TWIH together (referred to as CH+TWIH) form a very effective heuristic for generating initial solutions. The CH+TWIH method provides better solutions than the first feasible solution generated by CPLEX for 11 out of 14 problems. For 3 problems, CH+TWIH generates the best known solution while for all hard problems that method generates better solutions than the best solution generated by CPLEX. Further, CH+TWIH is much faster than CPLEX in generating feasible solutions. Finally, it is worth noting the large variations in CPU time for H-Randm for cases where the heuristic returns the same solution across all runs. Specifically, the ratio of maximum and minimum CPU time taken by H-Randm ranges from 1.00 (problem P6) to 16.6 (problem P4). This variation is amplified for problems where the heuristic finds a solution with a zero objective value since the heuristic can stop immediately without having to complete a full cycle of solves of two ship subproblems to prove that no further improvements can be achieved. These large variations indicate the potential for reducing overall CPU time through an intelligent ordering of two ship subproblems during 2SIH. 5.3 Sub-problem Selection in 2SIH Next we evaluate the effectiveness of the three schemes presented in Section for selecting the next two-ship subproblem during 2SIH. We refer to the heuristic when combined with the lexicographic, metrics-based and model-based selection schemes as H-Lexic, H-Mtrcs, and H-Model, respectively. The performance of these two-ship selection schemes is compared against H-Randm. The same parameters and termination criteria for the construction heuristic, the TWIH, and the two-ship subproblems as those presented above are used for all versions of the heuristic. In addition, we partition the set of iterations (two-ship subproblems solved) during 2SIH in to a set of cycles, where each cycle refers to a set of ( ) V 2 successive iterations. To achieve search diversification during 2SIH, we do not allow for a subproblem for the same ship-pair to be solved more than once during any cycle. This restriction is applied for all versions of the heuristic. For H-Model, this constraint is included explicitly in model rp. In order to reduce the overall CPU time used for selecting two-ship subproblems in H-Model, we solve the relaxation of model rp only when a new solution is obtained. Similarly, in the case of H-Mtrcs, metrics for quantifying the improvement that a two-ship subproblem will yield are re-used until a better solution is obtained. We use α = 20% for identifying the candidate ship-pairs for selection during H-Mtrcs. Note that results for P6, P8 and P9 are not included in this analysis since the best known solutions for these problems are generated by (CH+TWIH) itself. Also, results reported for the three largest problems (P10, P14 and P13) with H-Model are based on the use of Dantzig- Wolfe decomposition to solve the linear relaxation of rp for selecting two-ship subproblems during 2SIH. Results reported for H-Model for all other problems are obtained by solving rp using CLPEX. Table 4 compares the quality of solutions obtained by the heuristic method with the four ordering schemes. Solution quality is reported in terms of the % gap of the objective value relative to the best known solution (except in the case of P4 and P11 since the best known 19
Inventory Routing Problem for the LNG Business
Inventory Routing Problem for the LNG Business Roar Grønhaug 1 Marielle Christiansen 1 Guy Desaulniers 2,3 Jacques Desrosiers 2,4 1 Norwegian University of Science and Technology, 2 GERAD, 3 Ecole Polytechnique
More informationSolving Operational and Tactical Problems in the LNG Value Chain
Solving Operational and Tactical Problems in the LNG Value Chain Truls Flatberg, Marte Fodstad, Oddvar Kloster, Arnt-Gunnar Lium, Vibeke Nørstebø, Frode Rømo, Kristin Tolstad Uggen and Adrian Werner INFORMS
More informationThe LNG Inventory Routing Problem with Pick-Up Contracts
The LNG Inventory Routing Problem with Pick-Up Contracts Marielle Christiansen Roar Grønhaug Agenda The Liquefied Natural Gas Supply Chain Problem characteristics Regasification terminals Ship movements
More informationNetwork Flows. 7. Multicommodity Flows Problems. Fall 2010 Instructor: Dr. Masoud Yaghini
In the name of God Network Flows 7. Multicommodity Flows Problems 7.1 Introduction Fall 2010 Instructor: Dr. Masoud Yaghini Introduction Introduction In many application contexts, several physical commodities,
More informationA Mathematical Model for Driver Balance in Truckload Relay Networks
Georgia Southern University Digital Commons@Georgia Southern 12th IMHRC Proceedings (Gardanne, France 2012) Progress in Material Handling Research 2012 A Mathematical Model for Driver Balance in Truckload
More informationA heuristic for rich maritime inventory routing problems
A heuristic for rich maritime inventory routing problems Oddvar Kloster, Truls Flatberg, Geir Hasle Seminar NA / UNSW, Sydney, Australia July 5 2011 1 Outline SINTEF Introduction Model Algorithms Test
More informationAn Optimization Algorithm for the Inventory Routing Problem with Continuous Moves
An Optimization Algorithm for the Inventory Routing Problem with Continuous Moves Martin Savelsbergh Jin-Hwa Song The Logistics Institute School of Industrial and Systems Engineering Georgia Institute
More informationStorage Allocation and Yard Trucks Scheduling in Container Terminals Using a Genetic Algorithm Approach
Storage Allocation and Yard Trucks Scheduling in Container Terminals Using a Genetic Algorithm Approach Z.X. Wang, Felix T.S. Chan, and S.H. Chung Abstract Storage allocation and yard trucks scheduling
More informationVehicle Routing with Cross Docks, Split Deliveries, and Multiple Use of Vehicles. Arun Kumar Ranganathan Jagannathan
Vehicle Routing with Cross Docks, Split Deliveries, and Multiple Use of Vehicles by Arun Kumar Ranganathan Jagannathan A thesis submitted to the Graduate Faculty of Auburn University in partial fulfillment
More informationA Vehicle Routing in Glass Waste Collection
A Tino Henke, M. Grazia Speranza, Gerhard Wäscher University of Magdeburg, Germany University of Brescia, Italy Clausthal-Zellerfeld, May Introduction and Motivation A Glass Waste Collection Problem Introduction
More informationModelling the mobile target covering problem using flying drones
Modelling the mobile target covering problem using flying drones Luigi Di Puglia Pugliese 1 Francesca Guerriero 1 Dimitrios Zorbas 2 Tahiry Razafindralambo 2 1 Department of Mechanics, Energy and Management
More informationFleet Sizing and Empty Freight Car Allocation
Fleet Sizing and Empty Freight Car Allocation Philipp Hungerländer, Sebastian Steininger 13th July 2018 Abstract Empty freight car allocation problems as well as fleet sizing problems depict highly important
More informationSolving a Supply-Delivery Scheduling Problem with Constraint Programming
Solving a Supply-Delivery Scheduling Problem with Constraint Programming Katherine Giles and Willem-Jan van Hoeve Tepper School of Business, Carnegie Mellon University 5000 Forbes Avenue, Pittsburgh, PA
More informationAN ADAPTIVE LARGE NEIGHBORHOOD SEARCH ALGORITHM FOR A SELECTIVE & PERIODIC INVENTORY ROUTING PROBLEM
AN ADAPTIVE LARGE NEIGHBORHOOD SEARCH ALGORITHM FOR A SELECTIVE & PERIODIC INVENTORY ROUTING PROBLEM Deniz Aksen 1, Onur Kaya 2, F. Sibel Salman 2, Özge Tüncel 2 1 College of Adm. Sci. and Econ., Koç University,
More informationHigh-performance local search for solving real-life inventory routing problems
High-performance local search for solving real-life inventory routing problems Thierry Benoist 1, Bertrand Estellon 2, Frédéric Gardi 1, Antoine Jeanjean 1 1 Bouygues e-lab, Paris, France 2 Laboratoire
More informationA TABU SEARCH METAHEURISTIC FOR ASSIGNMENT OF FLOATING CRANES
1 st Logistics International Conference Belgrade, Serbia 28 - November 13 A TABU SEARCH METAHEURISTIC FOR ASSIGNMENT OF FLOATING CRANES Dragana M. Drenovac * University of Belgrade, Faculty of Transport
More informationOptimization in offshore supply vessel planning
Optimization in offshore supply vessel planning Elin E. Halvorsen-Weare and Kjetil Fagerholt This work has been published in Optimization and Engineering. The final publication is available at Springer
More informationMetaheuristics for scheduling production in large-scale open-pit mines accounting for metal uncertainty - Tabu search as an example.
Metaheuristics for scheduling production in large-scale open-pit mines accounting for metal uncertainty - Tabu search as an example Amina Lamghari COSMO Stochastic Mine Planning Laboratory! Department
More informationTitle: A Column Generation Algorithm for the Log Truck Scheduling Problem.
Title: A Column Generation Algorithm for the Log Truck Scheduling Problem. Authors: Myrna Palmgren Department of Optimization Linkoping University S-58183 Linkoping, Sweden e-mail: mypal@mai.liu.se Mikael
More informationHybrid search method for integrated scheduling problem of container-handling systems
Hybrid search method for integrated scheduling problem of container-handling systems Feifei Cui School of Computer Science and Engineering, Southeast University, Nanjing, P. R. China Jatinder N. D. Gupta
More informationDynamic Vehicle Routing and Dispatching
Dynamic Vehicle Routing and Dispatching Jean-Yves Potvin Département d informatique et recherche opérationnelle and Centre interuniversitaire de recherche sur les réseaux d entreprise, la logistique et
More informationRouting in maritime logistics
Routing in maritime logistics Truls Flatberg and Oddvar Kloster SINTEF 1 Outline Maritime routing Pickup and delivery variations Free delivery location Predefined number of visits Inter arrival gap Generic
More informationCongestion Reduction through Efficient Empty Container Movement
Congestion Reduction through Efficient Empty Container Movement Final Report METRANS Project 15-28 August 15, 2017 Principal Investigator: Maged M. Dessouky Ph.D. Graduate Student: Santiago Carvajal Daniel
More informationOptimizing a Containership Stowage Plan. using a modified Differential Evolution algorithm
Optimizing a Containership Stowage Plan using a modified Differential Evolution algorithm Speaker: Dr. Yun Dong ydong@tli.neu.edu.cn Supervisor: Pro. Lixin Tang Lixintang@mail.neu.edu.com The Logistics
More informationXXXII. ROBUST TRUCKLOAD RELAY NETWORK DESIGN UNDER DEMAND UNCERTAINTY
XXXII. ROBUST TRUCKLOAD RELAY NETWORK DESIGN UNDER DEMAND UNCERTAINTY Hector A. Vergara Oregon State University Zahra Mokhtari Oregon State University Abstract This research addresses the issue of incorporating
More informationISE480 Sequencing and Scheduling
ISE480 Sequencing and Scheduling INTRODUCTION ISE480 Sequencing and Scheduling 2012 2013 Spring term What is Scheduling About? Planning (deciding what to do) and scheduling (setting an order and time for
More informationEnterprise-Wide Optimization for Operations of Crude-Oil Refineries: Closing the Procurement and Scheduling Gap
Enterprise-Wide Optimization for Operations of Crude-Oil Refineries: Closing the Procurement and Scheduling Gap Brenno C. Menezes a*, Ignacio E. Grossmann a, Jeffrey D. Kelly b a Carnegie Mellon University,
More informationRehandling Strategies for Container Retrieval
Rehandling Strategies for Container Retrieval Tonguç Ünlüyurt and Cenk Aydin Sabanci University, Faculty of Engineering and Natural Sciences e-mail: tonguc@sabanciuniv.edu 1 Introduction In this work,
More informationCollaborative Logistics
Collaborative Logistics Martin Savelsbergh Ozlem Ergun Gultekin Kuyzu The Logistics Institute Georgia Institute of Technology 35th Annual Conference of the Italian Operations Research Society Lecce, September
More informationProcedia - Social and Behavioral Sciences 189 ( 2015 ) XVIII Annual International Conference of the Society of Operations Management (SOM-14)
Available online at www.sciencedirect.com ScienceDirect Procedia - Social and ehavioral Sciences 189 ( 2015 ) 184 192 XVIII Annual International Conference of the Society of Operations Management (SOM-14)
More informationMetaheuristics. Approximate. Metaheuristics used for. Math programming LP, IP, NLP, DP. Heuristics
Metaheuristics Meta Greek word for upper level methods Heuristics Greek word heuriskein art of discovering new strategies to solve problems. Exact and Approximate methods Exact Math programming LP, IP,
More informationA Particle Swarm Optimization Algorithm for Multi-depot Vehicle Routing problem with Pickup and Delivery Requests
A Particle Swarm Optimization Algorithm for Multi-depot Vehicle Routing problem with Pickup and Delivery Requests Pandhapon Sombuntham and Voratas Kachitvichayanukul Abstract A particle swarm optimization
More informationINTEGER PROGRAMMING BASED SEARCH
INTEGER PROGRAMMING BASED SEARCH A Thesis Presented to The Academic Faculty by Michael R Hewitt In Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in the School of Industrial
More informationUniversity Question Paper Two Marks
University Question Paper Two Marks 1. List the application of Operations Research in functional areas of management. Answer: Finance, Budgeting and Investment Marketing Physical distribution Purchasing,
More informationAirline Disruptions: Aircraft Recovery with Maintenance Constraints
1 Airline Disruptions: Aircraft Recovery with Maintenance Constraints Niklaus Eggenberg Dr. Matteo Salani and Prof. Michel Bierlaire In collaboration with APM Technologies Funded by CTI Switzerland 2 Dr.
More informationCapacitated vehicle routing problem for multi-product crossdocking with split deliveries and pickups
Available online at www.sciencedirect.com Procedia - Social and Behavioral Sciences 62 ( 2012 ) 1360 1365 WC-BEM 2012 Capacitated vehicle routing problem for multi-product crossdocking with split deliveries
More informationMulti-Commodity Network Flow Based Approaches for the Railroad Crew Scheduling Problem
Multi-Commodity Network Flow Based Approaches for the Railroad Crew Scheduling Problem Balachandran Vaidyanathan Department of Industrial & Systems Engineering, University of Florida, Gainesville, FL 32611
More informationHybrid MIP method for a Pickup and Delivery Problem with Time Windows and Dock Service Constraints
International Journal of Operations Research International Journal of Operations Research Vol. 8, No. 1, 36 56 (2011) Hybrid MIP method for a Pickup and Delivery Problem with Time Windows and Dock Service
More informationCongestion Reduction Through Efficient Empty Container Movement
Congestion Reduction Through Efficient Empty Container Movement August 2017 A Research Report from the National Center for Sustainable Transportation Maged Dessouky, University of Southern California Santiago
More informationExtended Model Formulation of the Proportional Lot-Sizing and Scheduling Problem. Waldemar Kaczmarczyk
Decision Making in Manufacturing and Services Vol. 5 2011 No. 1 2 pp. 49 56 Extended Model Formulation of the Proportional Lot-Sizing and Scheduling Problem with Lost Demand Costs Waldemar Kaczmarczyk
More informationVehicle Routing Tank Sizing Optimization under Uncertainty: MINLP Model and Branch-and-Refine Algorithm
Vehicle Routing Tank Sizing Optimization under Uncertainty: MINLP Model and Branch-and-Refine Algorithm Fengqi You Ignacio E. Grossmann Jose M. Pinto EWO Meeting, Sep. 2009 Vehicle Routing Tank Sizing
More informationProcedia - Social and Behavioral Sciences 109 ( 2014 ) Selection and peer review under responsibility of Organizing Committee of BEM 2013.
Available online at www.sciencedirect.com ScienceDirect Procedia - Social and Behavioral Sciences 109 ( 2014 ) 779 783 2 nd World Conference On Business, Economics And Management-WCBEM 2013 A hybrid metaheuristic
More informationColumn Generation Methods for Disrupted Airline Schedules
1 Column Generation Methods for Disrupted Airline Schedules Niklaus Eggenberg Dr. Matteo Salani and Prof. Michel Bierlaire In collaboration with APM Technologies Funded by CTI Switzerland Index 2 Index
More informationLocomotive Fuelling Problem (LFP) in Railroad Operations. Bodhibrata Nag 1 & Katta G.Murty 2
1 Locomotive Fuelling Problem (LFP) in Railroad Operations Bodhibrata Nag 1 & Katta G.Murty 2 About 75% of the world s railroads operate with diesel fuel. Even though European railroads rely on electric
More informationVehicle Routing with Departure Smoothing
ERASMUS UNIVERSITY ROTTERDAM Erasmus School of Economics Vehicle Routing with Departure Smoothing Master s Thesis Econometrics and Management Science Author: N.M. van der Zon Student ID: 387441 Supervisors:
More informationDISPATCHING TRANSPORT VEHICLES IN MARITIME CONTAINER TERMINALS
DISPATCHING TRANSPORT VEHICLES IN MARITIME CONTAINER TERMINALS by Pyung-Hoi Koo Department of Systems Management and Engineering, Pukyong National University, Busan, Korea Yongsoro 45, Namgu, Busan, South
More informationHeuristic Techniques for Solving the Vehicle Routing Problem with Time Windows Manar Hosny
Heuristic Techniques for Solving the Vehicle Routing Problem with Time Windows Manar Hosny College of Computer and Information Sciences, King Saud University, Riyadh, Saudi Arabia mifawzi@ksu.edu.sa Keywords:
More informationThe vehicle routing problem with demand range
DOI 10.1007/s10479-006-0057-0 The vehicle routing problem with demand range Ann Melissa Campbell C Science + Business Media, LLC 2006 Abstract We propose and formulate the vehicle routing problem with
More informationStrategic Design of Robust Global Supply Chains: Two Case Studies from the Paper Industry
Strategic Design of Robust Global Supply Chains: Two Case Studies from the Paper Industry T. Santoso, M. Goetschalckx, S. Ahmed, A. Shapiro Abstract To remain competitive in today's competitive global
More informationGROUP elevator scheduling is important to transportation
352 IEEE TRANSACTIONS ON AUTOMATION SCIENCE AND ENGINEERING, VOL. 7, NO. 2, APRIL 2010 Optimization of Group Elevator Scheduling With Advance Information Jin Sun, Student Member, IEEE, Qian-Chuan Zhao,
More informationDiscrete and dynamic versus continuous and static loading policy for a multi-compartment vehicle
European Journal of Operational Research 174 (2006) 1329 1337 Short Communication Discrete and dynamic versus continuous and static loading policy for a multi-compartment vehicle Yossi Bukchin a, *, Subhash
More informationGenerational and steady state genetic algorithms for generator maintenance scheduling problems
Generational and steady state genetic algorithms for generator maintenance scheduling problems Item Type Conference paper Authors Dahal, Keshav P.; McDonald, J.R. Citation Dahal, K. P. and McDonald, J.
More informationMulti-Period Vehicle Routing
Multi-Period Vehicle Routing Bruce L. Golden Robert H. Smith School of Business University of Maryland Presented at TRANSLOG Workshop Reñaca, Chile December 2009 1 Colleagues and Co-Authors Damon Gulczynski,
More information^ Springer. The Logic of Logistics. Theory, Algorithms, and Applications. for Logistics Management. David Simchi-Levi Xin Chen Julien Bramel
David Simchi-Levi Xin Chen Julien Bramel The Logic of Logistics Theory, Algorithms, and Applications for Logistics Management Third Edition ^ Springer Contents 1 Introduction 1 1.1 What Is Logistics Management?
More informationReal-Time Control Strategies for Rail Transit
Real-Time Control Strategies for Rail Transit Outline: Problem Description and Motivation Model Formulation Model Application and Results Implementation Issues Conclusions 1 Problem Context High frequency
More informationTramp Ship Routing and Scheduling - Models, Methods and Opportunities
Downloaded from orbit.dtu.dk on: Oct 14, 2018 Tramp Ship Routing and Scheduling - Models, Methods and Opportunities Vilhelmsen, Charlotte; Larsen, Jesper; Lusby, Richard Martin Publication date: 2015 Document
More informationSolving the Empty Container Problem using Double-Container Trucks to Reduce Congestion
Solving the Empty Container Problem using Double-Container Trucks to Reduce Congestion Santiago Carvajal scarvaja@usc.edu (424) 333-4929 Daniel J. Epstein Department of Industrial and Systems Engineering
More informationA New Formulation for the Combined Maritime Fleet Deployment and Inventory Management Problem
A New Formulation for the Combined Maritime Fleet Deployment and Inventory Management Problem Bo Dong 1 ( ), Tolga Bektaş 2, Saurabh Chandra 3, Marielle Christiansen 4, and Kjetil Fagerholt 5,6 1 Norwegian
More informationDeveloping a hybrid algorithm for distribution network design problem
Proceedings of the 2011 International Conference on Industrial Engineering and Operations Management Kuala Lumpur, Malaysia, January 22 24, 2011 Developing a hybrid algorithm for distribution network design
More informationThe Time Window Assignment Vehicle Routing Problem
The Time Window Assignment Vehicle Routing Problem Remy Spliet, Adriana F. Gabor June 13, 2012 Problem Description Consider a distribution network of one depot and multiple customers: Problem Description
More informationModeling and solving the train load planning problem in seaport container terminals
11 IEEE International Conference on Automation Science and Engineering Trieste, Italy - August 24-27, 11 ThB2.2 Modeling and solving the train load planning problem in seaport container terminals Daniela
More informationOptimizing Package Express Operations in China
Optimizing Package Express Operations in China Baris Yildiz 1 and Martin Savelsbergh 2 1 Department of Industrial Engineering, Koc University, Istanbul, Turkey 2 H. Milton Stewart School of Industrial
More informationDYNAMIC FREIGHT SELECTION FOR REDUCING LONG-HAUL ROUND TRIP COSTS
DYNAMIC FREIGHT SELECTION FOR REDUCING LONG-HAUL ROUND TRIP COSTS Arturo E. Pérez Rivera & Martijn Mes Department of Industrial Engineering and Business Information Systems University of Twente, The Netherlands
More informationABSTRACT. Timetable, Urban bus network, Stochastic demand, Variable demand, Simulation ISSN:
International Journal of Industrial Engineering & Production Research (09) pp. 83-91 December 09, Volume, Number 3 International Journal of Industrial Engineering & Production Research ISSN: 08-4889 Journal
More informationA Branch-and-Price Algorithm for the Vehicle Routing Problem with Roaming Delivery Locations
A Branch-and-Price Algorithm for the Vehicle Routing Problem with Roaming Delivery Locations Gizem Ozbaygin 1,2, Oya Karasan 2, Martin Savelsbergh 1, and Hande Yaman 2 1 H. Milton Stewart School of Industrial
More informationEvolutionary Algorithms
Evolutionary Algorithms Evolutionary Algorithms What is Evolutionary Algorithms (EAs)? Evolutionary algorithms are iterative and stochastic search methods that mimic the natural biological evolution and/or
More informationScheduling Quay Crane and Yard Equipment A case study of terminal in Central region
Thousand TEU Scheduling Quay Crane and Yard Equipment A case study of terminal in Central region Phong Ho Thanh Department of Industrial and Systems Engineering International University Vietnam National
More informationA Note on the Truck and Trailer Routing Problem
A Note on the Truck and Trailer Routing Problem Shih-Wei Lin a, Vincent F. Yu b,*, Shuo-Yan Chou b a Department of Information Management, Chang Gung University No. 259, Wen-Hwa 1st Road, Kwei-Shan, Tao-Yuan
More informationDecomposition Approaches for Optimal Production Distribution Coordination of Industrial Gases Supply Chains
Decomposition Approaches for Optimal Production Distribution Coordination of Industrial Gases Supply Chains Miguel Zamarripa, Pablo A. Marchetti, Ignacio E. Grossmann Department of Chemical Engineering
More informationLONG-HAUL FREIGHT SELECTION FOR LAST-MILE COST REDUCTION
LONG-HAUL FREIGHT SELECTION FOR LAST-MILE COST REDUCTION Arturo E. Pérez Rivera & Martijn Mes Department of Industrial Engineering and Business Information Systems University of Twente, The Netherlands
More informationMulti-Period Vehicle Routing with Call-In Customers
Multi-Period Vehicle Routing with Call-In Customers Anirudh Subramanyam, Chrysanthos E. Gounaris Carnegie Mellon University Frank Mufalli, Jose M. Pinto Praxair Inc. EWO Meeting September 30 th October
More informationStrategic Planning Models and Approaches to Improve Distribution Planning in the Industrial Gas Industry
Strategic Planning Models and Approaches to Improve Distribution Planning in the Industrial Gas Industry Leily Farrokhvar Dissertation submitted to the Faculty of the Virginia Polytechnic Institute and
More informationSolving an integrated Job-Shop problem with human resource constraints
Solving an integrated Job-Shop problem with human resource constraints PMS 10 - Tours (France) O. Guyon 1.2, P. Lemaire 3, É. Pinson 2 and D. Rivreau 2 1 École des Mines de Saint-Étienne 2 LISA - Institut
More informationA Decomposition Approach for the Inventory Routing Problem
A Decomposition Approach for the Inventory Routing Problem Ann Melissa Campbell Martin W.P. Savelsbergh Tippie College of Business Management Sciences Department University of Iowa Iowa City, IA 52242
More informationStrategies for Coordinated Drayage Movements
Strategies for Coordinated Drayage Movements Christopher Neuman and Karen Smilowitz May 9, 2002 Abstract The movement of loaded and empty equipment (trailers and containers) between rail yards and shippers/consignees
More informationResearch on Multi-Objective Co-operation Optimization of Shipyard Berth Allocation System
Research on Multi-Objective Co-operation Optimization of Shipyard Berth Allocation System Tang Weiwei 1*, Yang Kunrong 2 1 School of Civil Engineering and Transportation South China University of Technology
More informationA Framework for Integrating Planning Activities in Container Terminals
A Framework for Integrating Planning Activities in Container Terminals August 30th, 2007 S. H. Won and K. H. Kim Dept. of Industrial Engineering, Pusan National University, South Korea Contents 1 Introduction
More informationThe Tactical Berth Allocation Problem with QC Assignment and Transshipment Costs
The Tactical Berth Allocation Problem with QC Assignment and Transshipment Costs Models and Heuristics Ilaria Vacca Transport and Mobility Laboratory, EPFL joint work with Giovanni Giallombardo, Luigi
More informationScheduling and Coordination of Distributed Design Projects
Scheduling and Coordination of Distributed Design Projects F. Liu, P.B. Luh University of Connecticut, Storrs, CT 06269-2157, USA B. Moser United Technologies Research Center, E. Hartford, CT 06108, USA
More informationSimulated Annealing Algorithm for Vehicle Routing Problem with Transshipment
Simulated Annealing Algorithm for Vehicle Routing Problem with Transshipment Sukanya Thonethong #1, Jirachai Buddhakulsomsiri #2 # Logistics and Supply Chain Systems Engineering Sirindhorn International
More informationIterative train scheduling in networks with tree topologies: a case study for the Hunter Valley Coal Chain
22nd International Congress on Modelling and Simulation, Hobart, Tasmania, Australia, 3 to 8 December 2017 mssanz.org.au/modsim2017 Iterative train scheduling in networks with tree topologies: a case study
More informationGuy Desaulniers Eric Prescott Gagnon Louis Martin Rousseau Ecole Polytechnique, Montreal
Guy Desaulniers Eric Prescott Gagnon Louis Martin Rousseau Ecole Polytechnique, Montreal 1 Introduction Vehicle routing problem with time windows Motivation Large neighborhood search Hybrid LNS and Column
More informationA Simulation Study of Bulk Cement Product Maritime Tranportation Considering Vessel s Maintenance
P P P P School [1] A Simulation Study of Bulk Cement Product Maritime Tranportation Considering Vessel s Maintenance 1 2 3 4 Nurhadi SiswantoP P, ElisabethP P, Effi LatiffiantiP P, Ahmad RusdiansyahP P,
More informationBalanced Billing Cycles and Vehicle Routing of Meter Readers
Balanced Billing Cycles and Vehicle Routing of Meter Readers by Chris Groër, Bruce Golden, Edward Wasil University of Maryland, College Park University of Maryland, College Park American University, Washington
More informationSimultaneous Sensor Selection and Routing of Unmanned Aerial Vehicles for Complex Mission Plans
Simultaneous Sensor Selection and Routing of Unmanned Aerial Vehicles for Complex Mission Plans Frank Mufalli, Rajan Batta, and Rakesh Nagi Department of Industrial & Systems Engineering, University at
More informationHours of service regulations in road freight transport: an optimization-based international assessment. Thibaut Vidal
Hours of service regulations in road freight transport: an optimization-based international assessment Thibaut Vidal Seminar, Universidade Federal Fluminense, March 15 th, 2013 Context of this research
More informationCROSS-DOCKING: SCHEDULING OF INCOMING AND OUTGOING SEMI TRAILERS
CROSS-DOCKING: SCHEDULING OF INCOMING AND OUTGOING SEMI TRAILERS 1 th International Conference on Production Research P.Baptiste, M.Y.Maknoon Département de mathématiques et génie industriel, Ecole polytechnique
More informationVehicle Routing Problems in Supply Chain Management
Lecture 12. Vehicle Routing Problems in Supply Chain Management With the emergence of international markets and the growth of globalization, the management of supply chains has gained increased attention.
More informationWe consider a distribution problem in which a set of products has to be shipped from
in an Inventory Routing Problem Luca Bertazzi Giuseppe Paletta M. Grazia Speranza Dip. di Metodi Quantitativi, Università di Brescia, Italy Dip. di Economia Politica, Università della Calabria, Italy Dip.
More informationRouting and Scheduling in Tramp Shipping - Integrating Bunker Optimization Technical report
Downloaded from orbit.dtu.dk on: Jan 19, 2018 Routing and Scheduling in Tramp Shipping - Integrating Bunker Optimization Technical report Vilhelmsen, Charlotte; Lusby, Richard Martin ; Larsen, Jesper Publication
More informationStrengths & Drawbacks of MILP, CP and Discrete-Event Simulation based Approaches for Large-Scale Scheduling
Strengths & Drawbacks of MILP, CP and Discrete-Event Simulation based Approaches for Large-Scale Scheduling Pedro M. Castro Assistant Researcher Laboratório Nacional de Energia e Geologia Lisboa, Portugal
More informationManagement Science Letters
Management Science Letters 2 (202) 7 80 Contents lists available at GrowingScience Management Science Letters homepage: www.growingscience.com/msl Investigating transportation system in container terminals
More informationIntegrated Location Routing and Scheduling Problems
Integrated Location Routing and Scheduling Problems Zeliha Akça zelihaakca@lehigh.edu Rosemary T. Berger rosemary.berger@verizon.net Theodore K. Ralphs tkr2@lehigh.edu Department of Industrial and Systems
More informationContainer packing problem for stochastic inventory and optimal ordering through integer programming
20th International Congress on Modelling and Simulation, Adelaide, Australia, 1 6 December 2013 www.mssanz.org.au/modsim2013 Container packing problem for stochastic inventory and optimal ordering through
More informationA Log-Truck Scheduling Model Applied to a Belgian Forest Company. Jean-Sébastien Tancrez 1
A Log-Truck Scheduling Model Applied to a Belgian Forest Company Jean-Sébastien Tancrez 1 Louvain School of Management, Université catholique de Louvain, Mons, Belgium {js.tancrez@uclouvain.be} Abstract.
More informationOptimal Communications Systems and Network Design for Cargo Monitoring
Optimal Communications Systems and Network Design for Cargo Monitoring Daniel T. Fokum Ph.D. dissertation oral defense Department of Electrical Engineering & Computer Science University of Kansas December
More informationCHAPTER 5 SUPPLIER SELECTION BY LEXICOGRAPHIC METHOD USING INTEGER LINEAR PROGRAMMING
93 CHAPTER 5 SUPPLIER SELECTION BY LEXICOGRAPHIC METHOD USING INTEGER LINEAR PROGRAMMING 5.1 INTRODUCTION The SCMS model is solved using Lexicographic method by using LINGO software. Here the objectives
More informationSolving the Resource Allocation Problem in a Multimodal Container Terminal as a Network Flow Problem
Author manuscript, published in "Computational Logistics, Jürgen W. Böse, Hao Hu, Carlos Jahn, Xiaoning Shi, Robert Stahlbock, Stefan Voß (Ed.) (2011) 341-353" DOI : 10.1007/978-3-642-24264-9_25 Solving
More informationSystems : An Analysis. TG Crainic, J Damay, M Gendreau, R Namboothiri June 15, 2009
Strategic Planning of National/Regional Freight Transportation Systems : An Analysis TG Crainic, J Damay, M Gendreau, R Namboothiri June 15, 2009 Talk Outline Problem Motivation Multi-product Multi-modal
More informationINTEGRATING VEHICLE ROUTING WITH CROSS DOCK IN SUPPLY CHAIN
INTEGRATING VEHICLE ROUTING WITH CROSS DOCK IN SUPPLY CHAIN Farshad Farshchi Department of Industrial Engineering, Parand Branch, Islamic Azad University, Parand, Iran Davood Jafari Department of Industrial
More information