A heuristic approach for multi-product capacitated single-allocation hub location problems

Tese de mestrado, Estatística e Investigação Operacional, Universidade de Lisboa, Faculdade de Ciências, 2015Em redes onde o fluxo entre nodos é muito elevado (como pode ser o caso do transporte de pessoas e mercadorias ou até mesmo fluxo de dados numa rede), torna-se menos dispendioso criar pontos onde se concentram os fluxos provenientes das diferentes origens para depois serem consolidados e redistribuídos até aos destinos. A esses pontos dá-se o nome de hubs. O problema de localização de hubs consiste na localização de hubs numa rede e na alocação de todos os nodos da rede a esses hubs, de modo a que se possa encaminhar os fluxos entre os pares origem-destino a menos que sejam hubs. A rede constituída pelos hubs é normalmente definida como completa e não se permitem ligações diretas entre os pares origem-destino. Para além disso, assume-se que existe um factor de desconto para o fluxo que circula entre hubs. Neste tipo de redes (hub-and-spoke networks) podem aparecer duas variantes, no que diz respeito à alocação dos nodos aos hubs: single-allocation e multiple-allocation. No primeiro caso, permite-se apenas uma ligação de cada nodo não hub a um hub de modo a que todo o fluxo com origem e destino a cada nodo saia e chegue a esse nodo através de apenas um hub. No caso em que se tem multiple-allocation, cada nodo poderá ser afecto a mais do que um hub e o fluxo que chega e sai desse nodo poderá usar mais do que um hub. Algumas variantes que se poderão considerar para este problema incluem restrições de capacidade nos hubs (restrições que limitam a capacidade de um hub processar uma certa quantidade de fluxo de origem, limitações na capacidade total, limitações no processamento de fluxo que sai do hub, etc.), restrições de capacidade nos arcos, problemas multi-periódicos, presença de incerteza, o número de hubs ser fixo, o tipo de objectivo (minimizar custos, minimizar distâncias entre hubs, etc.) entre outras. A necessidade de aproximar este tipo de problemas aos casos que se observam no mundo real leva à inclusão de cada vez mais restrições dando origem a mais variantes do problema. Neste trabalho, será abordado o problema de localização de hubs na variante single-Allocation, com restrições de capacidade em relação ao fluxo que cada hub é capaz de processar. Para além disso, considera-se fluxos relativos a mais do que um tipo de produto. Este problema é designado por Problema Multi-produto de Localização de Hubs com Capacidade1. Cada hub poderá ser dedicado a processar apenas um tipo de produto, poderá processar mais do que um, ou mesmo todos. A rede de hubs é completa para cada produto mas, no entanto, se se considerar a rede de hubs para todos os produtos, esta poderá não ser completa. Como constatado em Correia et al. [17], no caso em que cada hub processa todos os tipos de produto, resolver o problema multi-produto ao invés de se resolver vários problemas, um para cada produto em separado, dá origem a melhores resultados. A complexidade inerente a este tipo de problemas leva a que sejam classificados como problemas NP-Hard pois não existem algoritmos que sejam capazes de os resolver em tempo polinomial. Por esta razão faz sentido desenvolver algoritmos heurísticos de modo a se conseguir obter, em tempo útil, soluções para instâncias maiores do problema . Como referido em Meyer et al. [51], em problemas de localização de hubs, duas soluções com valores objectivo muito semelhantes poderão ser estruturalmente muito diferentes, e portanto, através um mecanismo de pesquisa local poderá ser muito difícil a passagem de uma boa solução para outra melhor. Por esta razão, neste trabalho opta-se por uma heurística que se baseia num método em que se constroem soluções repetidamente. Para a construção das soluções, considerando que um processo de construção do tipo Greedy poderia dar origem a um número limitado de soluções e que as componentes da solução que são escolhidas por último são as piores, optou-se pelo desenvolvimento de um algoritmo de Ant Colony Optimization (ACO). Esta meta-heurística baseia-se no comportamento apresentado pelas formigas quando estas procuram alimento. Quando uma formiga deixa a colónia em busca de alimento, no seu trajeto, deposita um químico (feromona) que pode ser detectado por outras formigas. Quanto maior a concentração de feromona, maior a atração de cada formiga por esse trajeto e, portanto, os trajetos com maiores concentrações de feromonas serão percorridos por mais formigas. Por outro lado, se o caminho de ida e volta até ao alimento for mais curto, mais vezes será percorrido e maior será a concentração de feromona nesse caminho. O resultado destes dois tipos de reforço positivo nas concentrações de feromona nos trajetos percorridos pelas formigas, aliados ao facto de que existe evaporação do químico (a concentração de feromona diminui nos caminhos menos percorridos ao longo do tempo) dá origem aos \carreirinhos" de formigas que se podem observar na natureza e que normalmente representam o caminho mais curto entre o alimento e a Colónia de formigas. Considere-se o problema em questão em que se tem n nodos e p produtos. Para a representação das soluções, em vez de se considerar uma matriz binária n χ n χ p, onde o valor 1 representa uma afetação, considerou-se uma matriz n χ p, em que cada entrada representa, para cada produto, o hub ao qual o nodo foi afecto. O caso em que um nodo é afecto a si mesmo indica que esse nodo é hub para o produto correspondente. Este tipo de representação permite reduzir o tamanho da matriz e diminuir o uso da memória computacional. Antes da construção de uma solução, é aplicado um pré-processamento que vai evitar, com base nas restrições do problema, que certas componentes da solução sejam consideradas durante o processo de construção da solução. Deste modo, reduz-se o espaço de procura de soluções e algum esforço computacional. Para a construção de uma solução, escolhe-se o tamanho da colonia (o número de formigas que pertencem à colónia) e cada formiga vai escolhendo, sucessivamente, componentes da solução através de uma regra pseudo-aleatória onde algumas componentes da solução são escolhidas de um modo greedy e outras são escolhidas através de roulette wheel selection. A cada componente da solução é atribuído um valor inicial de feromona e, à medida que cada formiga vai adicionando componentes à solução, o valor da feromona associado à componente adicionada vai decrescendo, o que resulta na diminuição da probabilidade de que essa componente seja escolhida pela próxima formiga, dando origem à diversificação do conjunto de soluções construído por cada colónia. No fim, depois de todas as formigas terem construído uma solução, escolhe-se a melhor solução e reforça-se a concentração de feromona na melhor solução construída pela colónia. Se, por acaso, uma formiga der origem a uma solução não admissível, a solução construída por essa formiga não é considerada. Para mais detalhe em relação a este processo consultar Dorigo et al. [20]. Este tipo de algoritmo permite a inclusão de métodos de pesquisa local de modo a que a solução obtida por cada colónia seja melhorada. Com o objectivo de obter um algoritmo mais eficiente, escolheu-se incluir esta possibilidade e procedeu-se ao reforço da concentração de feromona após feita uma pesquisa local. Na pesquisa local efectuada, usaram-se três tipos de vizinhança. Um deles fecha os hubs dedicados que só servem a si próprios e realoca-os a outros já abertos para esse mesmo produto. Outro, escolhe aleatoriamente um nodo alocado a um hub dedicado para um dado produto e realoca-o a outro hub dedicado ao mesmo produto. Um terceiro, escolhe um hub aleatoriamente e transforma-o num nodo, realocando-o a outro hub dedicado ao mesmo tipo de produto. De modo a obter soluções iniciais melhores, explora-se a possibilidade de atribuir valores iniciais de feromona mais altos às componentes de solução pertencentes à solução da relaxação linear, na proporção do valor correspondente no caso das variáveis 0-1. Uma outra variação explorada consiste em fazer o reforço do valor de feromona às componentes da solução, apenas quando esta é a melhor de todas encontrada até ao momento, permitindo que haja evaporação de certas componentes de solução que poderão estar a ser escolhidas consecutivamente e permitindo que se escape mais facilmente de óptimos locais. Após implementação do algoritmo procede-se à fase dos testes computacionais em instâncias do problema com 10, 20, 25 e 40 nodos, 1, 2 e 3 produtos e hubs que processam 1, 2 e 3 produtos. As instâncias usadas nos testes computacionais pertencem ao Australian Post data set e foram adaptados por Correia et al. [17] de modo a que se tivesse dados para mais do que um tipo de produto.In this thesis, an heuristic procedure is proposed for the the multi-product capacitated single-allocation hub location problem. When addressing a problem in which it is necessary to determine the transportation of large commodity flows between many origin-destination (O-D) pairs, instead of using direct links, it becomes more efficient to design the networks in such a way that some of the nodes become consolidation centers or hubs. The Multi-Product Capacitated Single-Allocation Hub Location Problem (MP-CSAHLP according to Correia et al. [17]), is a NP-Hard problem in which several types of ow are considered, making it possible to consider the case when multiple types of products are to be shipped between each O-D pair. It can be seen as an extension of the classical Capacitated Single-Allocation Hub Location Problem. In the problem investigated in this work, no more than one hub can be located in each node and the hubs can be either dedicated (each hub can only handle one type of product) or non-dedicated (one hub can handle more than one type product). The hubs have capacity limitations regarding the incoming flow. Furthermore, the hub network is complete for each product but, when considering the hub network as a whole, it does not necessarily have to be complete. The goal is to locate the hubs in the network, allocate the non-hub nodes to the opened hubs and route the flow between each O-D pair. The objective is to minimize the total ow routing cost plus the setup costs of the hubs and costs of preparing the hubs to handle the different types of products. In order to obtain feasible solutions to the above problem, an Ant Colony Optimization procedure is proposed, which is a constructive, population-based meta-heuristic based in the foraging behavior of ants. Indirect communication between the ants through pheromones reflects the colony search experience. High-quality solutions are found as an outcome of the global cooperation among all the ants of the colony. A preprocessing procedure is also proposed in which some solution components are forbidden based on the problems restrictions. Such preprocessing reduces the search space and thus may reduce the computational effort. The proposed heuristic uses a single ant colony, which simultaneously chooses the hubs and allocates the nodes to the hubs. Once these solutions are found, the routing of the flow is computed in a short amount of time, using the optimization models for the MP-CSAHLP in which some variables (location and allocation) are fixed. The results show that the proposed heuristic has the potential to find good quality solutions for the MP-CSAHLP and that its performance can be improved with finer parameter tuning, longer runs and more intense local search

Scalable Facility Location for Massive Graphs on Pregel-like Systems

We propose a new scalable algorithm for facility location. Facility location is a classic problem, where the goal is to select a subset of facilities to open, from a set of candidate facilities F , in order to serve a set of clients C. The objective is to minimize the total cost of opening facilities plus the cost of serving each client from the facility it is assigned to. In this work, we are interested in the graph setting, where the cost of serving a client from a facility is represented by the shortest-path distance on the graph. This setting allows to model natural problems arising in the Web and in social media applications. It also allows to leverage the inherent sparsity of such graphs, as the input is much smaller than the full pairwise distances between all vertices. To obtain truly scalable performance, we design a parallel algorithm that operates on clusters of shared-nothing machines. In particular, we target modern Pregel-like architectures, and we implement our algorithm on Apache Giraph. Our solution makes use of a recent result to build sketches for massive graphs, and of a fast parallel algorithm to find maximal independent sets, as building blocks. In so doing, we show how these problems can be solved on a Pregel-like architecture, and we investigate the properties of these algorithms. Extensive experimental results show that our algorithm scales gracefully to graphs with billions of edges, while obtaining values of the objective function that are competitive with a state-of-the-art sequential algorithm

Locating and Protecting Facilities Subject to Random Disruptions and Attacks

Recent events such as the 2011 Tohoku earthquake and tsunami in Japan have revealed the vulnerability of networks such as supply chains to disruptive events. In particular, it has become apparent that the failure of a few elements of an infrastructure system can cause a system-wide disruption. Thus, it is important to learn more about which elements of infrastructure systems are most critical and how to protect an infrastructure system from the effects of a disruption. This dissertation seeks to enhance the understanding of how to design and protect networked infrastructure systems from disruptions by developing new mathematical models and solution techniques and using them to help decision-makers by discovering new decision-making insights. Several gaps exist in the body of knowledge concerning how to design and protect networks that are subject to disruptions. First, there is a lack of insights on how to make equitable decisions related to designing networks subject to disruptions. This is important in public-sector decision-making where it is important to generate solutions that are equitable across multiple stakeholders. Second, there is a lack of models that integrate system design and system protection decisions. These models are needed so that we can understand the benefit of integrating design and protection decisions. Finally, most of the literature makes several key assumptions: 1) protection of infrastructure elements is perfect, 2) an element is either fully protected or fully unprotected, and 3) after a disruption facilities are either completely operational or completely failed. While these may be reasonable assumptions in some contexts, there may exist contexts in which these assumptions are limiting. There are several difficulties with filling these gaps in the literature. This dissertation describes the discovery of mathematical formulations needed to fill these gaps as well as the identification of appropriate solution strategies

Dual-Based Local Search for Deterministic, Stochastic and Robust Variants of the Connected Facility Location Problem

In this dissertation, we propose the study of a family of network design problems that arise in a wide range of practical settings ranging from telecommunications to data management. We investigate the use of heuristic search procedures coupled with lower bounding mechanisms to obtain high quality solutions for deterministic, stochastic and robust variants of these problems. We extend the use of well-known methods such as the sample average approximation for stochastic optimization and the Bertsimas and Sim approach for robust optimization with heuristics and lower bounding mechanisms. This is particular important for NP-complete problems where even deterministic and small instances are difficult to solve to optimality. Our extensions provide a novel way of applying these techniques while using heuristics; which from a practical perspective increases their usefulness

Allocation Strategies in Hub Networks

Cataloged from PDF version of article.In this paper, we study allocation strategies and their effects on total routing costs in hub networks. Given a set of nodes with pairwise traffic demands, the p-hub median problem is the problem of choosing p nodes as hub locations and routing traffic through these hubs at minimum cost. This problem has two versions; in single allocation problems, each node can send and receive traffic through a single hub, whereas in multiple allocation problems, there is no such restriction and a node may send and receive its traffic through all p hubs. This results in high fixed costs and complicated networks. In this study, we introduce the r-allocation p-hub median problem, where each node can be connected to at most r hubs. This new problem generalizes the two versions of the p-hub median problem. We derive mixed-integer programming formulations for this problem and perform a computational study using well-known datasets. For these datasets, we conclude that single allocation solutions are considerably more expensive than multiple allocation solutions, but significant savings can be achieved by allowing nodes to be allocated to two or three hubs rather than one. We also present models for variations of this problem with service quality considerations, flow thresholds, and non-stop service. (C) 2011 Elsevier B.V. All rights reserved

Bilevel facility location problems: theory and applications.

In this doctoral thesis we focus on studying facility location problems considering customer preferences. In these problems, there is a set of customers or users who demand a service or product that must be supplied by one or more facilities. By facilities it is understood some object or structure that offers some service to customers. One of the most important assumptions is that customers have established their own preferences over the facilities and should be taken into account in the customer-facility assignment. In real life, customers choose facilities based on costs, preferences, a predetermined contract, or a loyalty coefficient, among others. That is, they are free to choose the facilities that will serve them. The situation described above is commonly modeled by bilevel programming, where the upper level corresponds to location decisions to optimize a predefined criteria, such as, minimize location and distribution costs or maximize the demand covered by the facilities; and the lower level is associated to -customer allocation- to optimize customer preferences. The hierarchy among both levels is justified because the decision taken in the upper level directly affects the decision’s space in the lower level

Model and solution methods for some hub location problems

In this thesis we study some hub location problems in the context of transportation networks. These are combinatorial optimization problems appearing in situations where there is a need of transporting some traffic, like items, people, and information, from many origins to many destinations. Instead of sending these flows using a direct shipment between all pairs of nodes in the network, a subset of these nodes is selected to use as hubs, with the aim of consolidating and distribute the flows. Thus, hubs induce a subnetwork that sends the traffic more efficiently and at a cheaper cost, allowing economies of scale when large amounts of traffic between nodes on this subnet are transported. We study different variants of hub location problems that try to model several real world situations and characteristics. In all of them, we aim to minimize the cost of sending traffic through the transportation network.In this thesis we study some hub location problems in the context of transportation networks. These are combinatorial optimization problems appearing in situations where there is a need of transporting some traffic, like items, people, and information, from many origins to many destinations. Instead of sending these flows using a direct shipment between all pairs of nodes in the network, a subset of these nodes is selected to use as hubs, with the aim of consolidating and distribute the flows. Thus, hubs induce a subnetwork that sends the traffic more efficiently and at a cheaper cost, allowing economies of scale when large amounts of traffic between nodes on this subnet are transported. We study different variants of hub location problems that try to model several real world situations and characteristics. In all of them, we aim to minimize the cost of sending traffic through the transportation network

Generation and optimisation of real-world static and dynamic location-allocation problems with application to the telecommunications industry.

The location-allocation (LA) problem concerns the location of facilities and the allocation of demand, to minimise or maximise a particular function such as cost, profit or a measure of distance. Many formulations of LA problems have been presented in the literature to capture and study the unique aspects of real-world problems. However, some real-world aspects, such as resilience, are still lacking in the literature. Resilience ensures uninterrupted supply of demand and enhances the quality of service. Due to changes in population shift, market size, and the economic and labour markets - which often cause demand to be stochastic - a reasonable LA problem formulation should consider some aspect of future uncertainties. Almost all LA problem formulations in the literature that capture some aspect of future uncertainties fall in the domain of dynamic optimisation problems, where new facilities are located every time the environment changes. However, considering the substantial cost associated with locating a new facility, it becomes infeasible to locate facilities each time the environment changes. In this study, we propose and investigate variations of LA problem formulations. Firstly, we develop and study new LA formulations, which extend the location of facilities and the allocation of demand to add a layer of resilience. We apply the population-based incremental learning algorithm for the first time in the literature to solve the new novel LA formulations. Secondly, we propose and study a new dynamic formulation of the LA problem where facilities are opened once at the start of a defined period and are expected to be satisfactory in servicing customers' demands irrespective of changes in customer distribution. The problem is based on the idea that customers will change locations over a defined period and that these changes have to be taken into account when establishing facilities to service changing customers' distributions. Thirdly, we employ a simulation-based optimisation approach to tackle the new dynamic formulation. Owing to the high computational costs associated with simulation-based optimisation, we investigate the concept of Racing, an approach used in model selection, to reduce the high computational cost by employing the minimum number of simulations for solution selection

Optimal Size and Location of Warehouses under Risk of Failure

Nowadays, every company faces challenges that seem to be loaded with a contradiction: how to reduce operations and transportation costs while increasing customer satisfaction levels. Designing a supply chain network is an effective solution to such an issue. Supply chain network design involves making decisions about the number, sizes, and locations of the facilities in a supply chain. The focus of this study is how to choose appropriate warehouse locations and sizes in supply chain network design. The study is divided into two parts. In the second part, the risk of warehouse failure is considered while in the first part, it is not. Three sets of mathematical optimization models for warehouse location and branch assignment were developed. The first set of mathematical optimization models covered the case of warehouse location without risk. Two sets of decision variables were introduced to determine the locations for new warehouses and assign warehouses to branches. The second set of mathematical optimization models covered the warehouse location problem under the risk of warehouse failure. Again, two sets of decision variables were introduced. The first set of decision variables helped in determining the locations for new warehouses, and the second set helped in assigning a primary and a backup warehouse to each branch. The backup warehouse to be used in case of failure of the primary warehouse. The third set of mathematical optimization models covered the case in which some warehouses can be fortified to become totally risk-free. Each branch was either assigned to a primary fortified warehouse only or to a primary warehouse that was not fortified and a secondary fortified warehouse. Fortification model required an additional variables indicating which warehouses to be fortified. Warehouses with multiple capacity levels and multiple part category types were considered, which is a contribution to the topic of warehouse disruption risk. Specialized warehouses were also considered in this dissertation, which is another contribution of this dissertation. Some linearization and relaxation methods were used to help in solving the three models. Further, a solution methodology was presented based on the solution to scenario subproblems that are more easily, i.e., more quickly, solved. This requires an algorithm to determine the scenarios. Each scenario represents the number and sizes of warehouses needed to be built. The scenarios are novel in that they do not specify a subset of warehouses to be opened, but rather they specify the number of warehouses of each size to be opened. The results showed the effectiveness of the proposed solution methodology by application to an example based on a case study of a Canadian company; and a created example based on European cities

Hub Network Design and Discrete Location: Economies of Scale, Reliability and Service Level Considerations

In this thesis, we study three related decision problems in location theory. The first part of the dissertation presents solution algorithms for the cycle hub location problem (CHLP), which seeks to locate p-hub facilities that are connected by means of a cycle, and to assign non-hub nodes to hubs so as to minimize the total cost of routing flows through the network. This problem is useful in modeling applications in transportation and telecommunications systems, where large setup costs on the links and reliability requirements make cycle topologies a prominent network architecture. We present a branch and-cut algorithm that uses a flow-based formulation and two families of mixed-dicut inequalities as a lower bounding procedure at nodes of the enumeration tree. We also introduce a greedy randomized adaptive search algorithm that is used to obtain initial upper bounds for the exact algorithm and to obtain feasible solutions for large-scale instances of the CHLP. Numerical results on a set of benchmark instances with up to 100 nodes confirm the efficiency of the proposed solution algorithms. In the second part of this dissertation, we study the modular hub location problem, which explicitly models the flow-dependent transportation costs using modular arc costs. It neither assumes a full interconnection between hub nodes nor a particular topological structure, instead it considers link activation decisions as part of the design. We propose a branch-and-bound algorithm that uses a Lagrangean relaxation to obtain lower and upper bounds at the nodes of the enumeration tree. Numerical results are reported for benchmark instances with up to 75 nodes. In the last part of this dissertation we study the dynamic facility location problem with service level constraints (DFLPSL). The DFLPSL seeks to locate a set of facilities with sufficient capacities over a planning horizon to serve customers at minimum cost while a service level requirement is met. This problem captures two important sources of stochasticity in facility location by considering known probability distribution functions associated with processing and routing times. We present a nonlinear mixed integer programming formulation and provide feasible solutions using two heuristic approaches. We present the results of computational experiments to analyze the impact and potential benefits of explicitly considering service level constraints when designing distribution systems