SciELO - Scientific Electronic Library Online

 
 número20Inventariação do património geomorfológico cárstico do Parque Nacional da Furna Feia (RN, Brasil) como proposta para uso geoturísticoImpacto da pandemia da Covid-19 nos hábitos de utilização dos parques verdes urbanos: o caso de Condeixa-a-Nova índice de autoresíndice de assuntosPesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

  • Não possue artigos similaresSimilares em SciELO

Compartilhar


GOT, Revista de Geografia e Ordenamento do Território

versão On-line ISSN 2182-1267

GOT  no.20 Porto jun. 2020  Epub 30-Dez-2020

https://doi.org/10.17127/got/2020.20.007 

Artigos Originais

Reducing Costs of Preventive Lightning Systems by Locational Optimization

Reduzindo custos de sistemas preventivos de descargas atmosféricas por otimização locacional

Matheus, Libórioi 

Willian, Maiaii 

Martins, Carlosiii 

Ekel, Petriii 

Laudares, Sandroiv 

Patrícia Bernardesi 

iPontifícia Universidade Católica de Minas Gerais, Departamento de Administração Av. Itaú, 525 - Dom Cabral, 30535-012, Belo Horizonte, Minas Gerais, Brasil (m4th32s@gmail.com)(patriciabernardes@pucminas.br)

iiUniversidade Federal de Minas Gerais, Programa de Pós-graduação em Engenharia Elétrica 31270-901, Belo Horizonte, Brasil (willianmaia1@yahoo.com.br)

iiiPontifícia Universidade Católica de Minas Gerais, Departamento de Engenharia Elétrica 30535-012, Belo Horizonte, Brasil (capsm@pucminas.br) (petr.ekel2709@gmail.com)

ivPontifícia Universidade Católica de Minas Gerais, Departamento de Geografia 30535-012, Belo Horizonte, Brasil (laudares@pucminas.br)


ABSTRACT

The installation of lightning detection sensor systems involves a large volume of resources. However, the absence of these systems can lead to irreparable damage, such as the death of a person. How to reduce the costs of implementing these systems to protect workers and citizens or even public and private equipment? Usually, this problem is not solved on a formal level, but based on the intuition of the involved specialists, without applying optimization techniques. Taking this into account, this research explores the theoretical structure of locational optimization methods to propose a simplified heuristic method that does not require any programming skills for its application. Specifically, its use directly resulted in the cost reduction of preventive lightning systems in an area containing 41 mines. All mines must be covered by the sensor system which has a coverage radius of 10 km. The experimental results demonstrate that the application of this method generates more effective results than those obtained using known methods of locational optimization. The method proposed in this research contributes directly to reduce costs of preventive lightning systems and indirectly to protect workers in mining areas as well as to protect populations in urban areas where the frequency and incidence of lightning is especially high.

Palavras-chave: lightning discharge; locational problems; preventive lightning detection systems; spatial optimization.

RESUMO

A instalação de sistemas de sensores de detecção de descargas atmosféricas envolve um grande volume de recursos. No entanto, a ausência desses sistemas pode levar a danos irreparáveis, como a morte de uma pessoa. Como reduzir os custos de implantação desses sistemas para se proteger trabalhadores e cidadãos ou mesmo equipamentos públicos e privados? Normalmente, esse problema não é resolvido em um nível formal e sim com base na intuição dos especialistas envolvidos, sem aplicar técnicas de otimização. Levando isso em consideração, essa pesquisa explora a estrutura teórica dos métodos de otimização locacional para propor um método heurístico simplificado que não exige qualquer habilidade de programação para sua aplicação. Especificamente, seu uso foi demonstrado diretamente pela redução de custos de sistemas preventivos de raios em uma área contendo 41 minas. Todas as minas devem ser cobertas pelo sistema de sensores, que tem um raio de cobertura de 10 km. Os resultados experimentais demonstram que a aplicação desse método gera resultados mais efetivos que do que os obtidos utilizando métodos conhecidos de otimização locacional. O método proposto nesta pesquisa contribui diretamente para reduzir custos do sistema preventivo de raio e indiretamente para proteger trabalhadores em áreas de mineração ou mesmo para proteger populações em áreas urbanas onde a frequência e a incidência de descargas atmosféricas são especialmente elevadas.

Keywords: descargas atmosféricas; problemas de localização; sistemas preventivos de detecção de raios; otimização espacial.

1. Introduction

If, on the one hand, atmospheric discharges detection systems are quite expensive (Gomes et al., 2012). On the other hand, cargo and equipment exposed to lightning in mining areas are also quite expensive. In addition, the safety of workers in the mining industry must be considered (Berger, 2007).

The geographic analysis of the elements that make up the context of this work, through the exchange of theories and concepts between mathematical and geographic sciences, can contribute to more rational decisions (Laudares, Laudares & Libório, 2016). For example, it can serve for minimizing the costs of deploying atmospheric discharge detection systems, protecting equipment, and ensuring the safety of workers in mines. Figure 1 illustrates the problem of the atmospheric discharge detection system allocation for 41 mines of a mining company. The question that arises with this problem is: where to install the atmospheric discharge detection systems to protect all mines with the fewest equipment?

Google maps

Figure 1 Mines exposed to atmospheric discharge.  

This research reviews several methods of locational optimization (Church & ReVelle, 1974; Farahani, Asgari, Heidari, Hosseininia & Goh, 2012; Bergantiños, Gómez-Rúa, Llorca, Pulido & Sánchez-Soriano, 2020). Based on this review, an optimization method is developed and applied. The results obtained with its utilization are compared with the solutions generated by the use of other two optimization methods.

The solution contributes to the reduction of the costs of implantation of atmospheric discharge detection systems. This contribution may be especially relevant for the expansion of integrated networks of lightning detection sensors, such as, for example, the sensor network of the Metropolitan Region of Paraíba Vale (Magina et al., 2016) which is responsible for protecting the population and public facilities.

2. Thunderstorm Warning Systems

Atmospheric discharge is a natural, seasonal, and destructive phenomenon, occurring because of the accumulation of electrical charges in the atmosphere, usually when storms occurs (Visacro, 2005; Pinto Jr, 2005). The atmospheric discharge can be defined as an intense electric discharge reaching kilometers of extension, with peaks of the current intensity, being able to surpass 1 kA (Visacro ,2005). Atmospheric discharges occur when the electric field, produced by electric charges, exceeds the insulation capacity of the air. At this moment a rapid movement of electrons from a region with negative charges begins to reach a region with positive charges, producing the atmospheric discharge (Magina, Pinto Jr & Naccarato, 2016; Lopez, Perez, Herrera, Aranguren & Porras, 2012). As a result, atmospheric discharges from the cloud may in some cases reach the earth's surface (Pinto Jr, 2005).

Brazil is cited by Pinto Jr (2005) and Magina et al. (2016) as the country with the highest incidence of atmospheric discharges in the world, being reached by about 50 million lights per year each year. In particular, the lightning discharges are frequent, reaching, mainly, people who work outdoors (Pinto Jr, 2005). These accidents cause deaths in 20% of cases. Between 2001 and 2014, it was registered 1,792 deaths, representing 2% of the deaths due to lightning in the world (Berger, 2007).

The study of lightning discharge in urban areas has been gaining space in the specialized literature as urban heating and local physiographic specificities have been related to changes in the frequency and quantity of lightning discharge in urban areas (Stallins & Rose, 2008). Studies show that the increased lifetime of storms caused by the effect of urban heat islands and the concentration of pollution explain the increase in the number of lightning discharge in urban areas (Farias, Pinto Jr, Naccarato & Pinto, 2009). As a result, urban areas have the highest average of lightning discharge per square kilometer, with the frequency of lightning discharges being positively correlated with heat fluctuations (Sokol & Rohli, 2018).

Thunderstorm warning systems (TWSs) fit precisely in this context which the large amount of lightning discharges observed in Brazil exposes people to risks and public and private facilities to potential damage.

In mining industry, protection structures protect the equipment allocated to outdoor works, for example, construction works, drilling, environmental assessments, maintenance and operation of electricity or telecommunications networks, among others (Cooray, 2009; Leite, Modena & Sueta, 2007). The safety of workers in these operations can also be achieved through the use of structures such as detection and warning systems for atmospheric discharges (Lopez et al., 2012; Magina et al., 2016).

TWSs can serve as such system (Lundquist & Scuka, 1971). Unlike passive systems, which protect equipment through fixed bars, cages, lightning rods or grounding techniques (Cooray, 2009), TWSs detect and alert in advance possible atmospheric discharges (Lopez et al., 2012; Magina et al., 2016). Each sensor of TWS covers 10 km radius (Lundquist & Scuka, 1971).

Detection and warning systems monitor the electromagnetic field and the electrostatic field of the storm, detecting possible atmospheric discharges (Magina et al., 2016). Monitoring is performed by a system composed of a set of sensors interconnected to a central processing unit (Lopez et al., 2012). The sensors collect information from the storms that are sent and analyzed by the central processing unit. This analysis indicates the probability of an atmospheric discharge occurring. Risk situations trigger warning signals (Lundquist & Scuka, 1971).

The warning signals, once triggered, create a window, for places like this one, of 32 minutes Lopez et al., 2012). This allows performing worker safety procedures (Berger, 2007; Gomes, Ab Kadir & Cooper, 2012) and protecting equipment (Cooray, 2009; Leite et al., 2007).

3. Optimal location problems

The problem of this research can be formulated as a problem of an optimal location identification. The first studies and theoretical results related to locational optimization problems can be found in Weber (1929). This work, dedicated to the theory of location in the industry, considers geographic factors in the search of the location that minimizes the costs of installation of factories. Later, methods derived from gravitational models emerged, such as the Gravity Point (Reilly, 1931), Inverse Distance Weighting (Converse, 1949), and the Area of Influence (Applebaum, 1966). Although widely used in locational problems in Geography, these methods do not perform any type of optimization.

In Ortiz-Astorquiza, Contreras, and Laporte (2018), several locational optimization methods are used to solve problems associated with the location of facilities. Optimal location methods have been applied to solve different problems in different fields. For example, to indicate the ideal locations of a sensor network for monitoring spatial phenomena (Nguyen, Kodagoda, Ranasinghe & Dissanayake, 2015) or to position sensors for surveillance and monitoring of environmental areas (Mathew & Surana, 2012). The results of (HassanzadehFard & Jalilian, 2018) are associated with determining of places of the installation and sizes of reactive power sources in electrical distribution systems. The authors of (Kuschel et al., 2019) consider a problem of defining the smallest number of radars, providing the maximal area aiming to maximize coverage of the monitored area, observing the positions allowed for their installation (Kuschel et al., 2019). The results of (Tong & Murray, 2017) are directed at configuring the locations of a hospital network to minimize distances between service units and locations, noting demand level at each location (Tong & Murray, 2017).

Another example is to identify optimal locations to install warehouses, observing that the installation costs are composed of fixed part of investment and variable part of investment that varies according to the number of customers to be serviced by the warehouse (Ivanov et al., 2019).

Church and ReVelle (1974) conceptualize the locational optimization problem as the maximization of facilities coverage, having as restrictions a critical distance and a specific number of installations. However, the most common locational optimization problem in actual life is demonstrated in Farahani et al. (2012) and Bergantiños et al. (2020). In these works, the presented coverage problem considers that all facilities must be serviced by at least one type of equipment. This equipment has, in turn, a maximum coverage radius. For example, the costs of deploying emergency facilities can be minimized, satisfying a maximum coverage radius, and the entire set of cities or locations must be met. In Owen and Daskin (1998), the best places to trade, should cover the entire consumer market, and its long-term growth, with the lowest operating costs.

The combinatorial optimization problems associated with localization have their common basis to the application of approximate methods (Blesa & Blum, 2018), in particular, the approximate methods of discrete mathematical programming (Cooper, 1963; Ivanov et al., 2019). The corresponding objective functions have a large number of local minimums, and the problem global solution is to be reached (Gill, Murray & Wright, 2019). The difference between these solutions is shown by Blesa and Blum (2018) as the optimization gap. This gap can be reduced with reallocations or substitutions (Ivanov et al., 2019).

The heuristic methods can be structured and constructed according to the boundary conditions of each problem (Salhi, 2017). Although stochastic methods (Owen & Daskin, 1998), linear programming (Ivanov et al., 2019) and genetic algorithms (Kuschel et al., 2019) can be used in optimal location problems, a heuristic designed specifically for the problem of minimizing the costs of implementing lightning systems can present more satisfactory results.

For example, the greedy methods presented by Church and ReVelle (1974), Owen and Daskin (1998) and Expósito-Márquez et al. (2019). The applications of these methods allows one, for example, to minimize the number of installations, covering the demand, without exceeding the maximum allowed distance using the Greedy Adding (GA) (Church & ReVelle, 1974) or Greedy Adding with Substitution (GAS) algorithm (Expósito-Márquez et al., 2019). The GAS algorithm is similar to the concepts of the traveling salesman and warehouse location, where optimization is also local and global optimization is not guaranteed (Church & ReVelle, 1974).

4. Materials and Methods

4.1. Allocation and Reallocation (AR) method

The Allocation and Relocation (AR) method developed in this research can be performed in three phases: data preparation, processing, and reprocessing.

In this case study, we had to minimize the costs involved in implementing a system for the lightning discharges detection in 41 mines from a private mining company. We highlight that the locations that can install the system have interconnection with the company's automation system, access infrastructure and availability of electricity.

The geographic coordinates of the possible installation points of the sensor system for the protection of the mines were collected. With these coordinates (longitude and latitude), the distances between the mines were calculated (Table 1).

Table 1 Distances between mines (Km) 

𝑑𝑖𝑗 n1 n2 n3 n4 n5
n1 0 12 7 16 17
n2 12 0 8 9 9
n3 7 8 0 12 13
n4 16 9 12 0 9
n5 17 9 13 9 0

In order to calculate the distance 𝑑𝑖𝑗 Formula 1 was applied, where 𝑥𝑖 is the mine 𝑖 longitude, and 𝑥𝑗 and is mine j longitude and 𝑐 is Earth’s curvature adjustment (see Field, 2018).

(1) 

It can be attested that a mine is protected when the maximum system coverage radius (10 km) is less or equal the distance between two mines (Formula 2).

(2)  

Through Formula 3 it is possible to calculate the system coverage potential for each mine, that is βi:

(3)  

The use of Formulas 2 and 3 permits one to represent the system coverage potential for each mine, reflected in table 2 creation.

Table 2 System coverage potential per mine  

α𝑖𝑗 n1 n2 n3 n4 n5 βi
n1 1 0 1 0 0 2
n2 0 1 1 1 1 4
n3 1 1 1 0 0 3
n4 0 1 0 1 1 3
n5 0 1 0 1 1 3

In the AR method the sensors are initially allocated to mines with the lowest βi value. Therefore, when the mine coverage is equal to 1, we allocate a sensor for this mine and remove it from the resource available.

But, as seen, all βi are greater than 1. At the same time, the mine with the smallest β in Table 3 is n1. Next, it was seen that n3 is potentially reached by a sensor installed in n1. Therefore, the protection system allocation was evaluated if n1 is ADEQUATE. To do this it was checked if β3 ≤ β1.

Table 3 System's potential coverage per mine 

α𝑖𝑗 n1 n2 n3 n4 n5 βi
n1 1 0 1 0 0 2
n2 0 1 1 1 1 4
n3 1 1 1 0 0 3
n4 0 1 0 1 1 3
n5 0 1 0 1 1 3

Since the condition β3 ≤ β1 is FALSE, it is possible to confirm that regarding AR method the protection system allocation in n1 is NOT ADEQUATE. Therefore, protection system should be allocated from n1 to n1. Thus, the mines covered by the system installed in n3 (n1, n2 and n3) are taken from Table 4 and then the process is repeated.

Table 4 System's potential coverage per mine 

α𝑖𝑗 n1 n2 n3 n4 n5 βi
n1 1 0 1 0 0 2
n2 0 1 1 1 1 4
n3 1 1 1 0 0 3
n4 0 1 0 1 1 2
n5 0 1 0 1 1 2

The procedure is restarted identifying in Table 4 the mine with the lowest βi value. In this case, the sensor can be allocated in the mines n4 or n5. Therefore, we evaluate whether the allocation of the protection system in n4 or in n5 is ADEQUATE. For this, we check if β4 ≤ β5 and/or if β5 ≤ β4. As these relationships are TRUE, we can state that the allocation of the system in n4 or on n5 is ADEQUATE. The procedure repeats until all values are taken off the table.

4.2. Greedy Adding methods

Running the greedy adding method, we identify and allocate a system of sensors in the mine with greater value. Then we remove from the table all the mines served by the allocated system. The second system will always be allocated to the mine with the highest potential until all values are removed from the table.

Running GAS we identify and allocate the system in the mine with greater value. We repeat the procedure to find the location of the second and third sensor. After that, we removed the sensor from the second location and defined a new location for this sensor, looking for the largest mine. Then, we compare the two solutions generated, choosing the one with the largest number of covered mines. In summary, with the GAS method, we change the default locations with each new entry. Iterations are associated with replacing a chosen facility when its change improves maximum coverage.

5. Results and Discussions

In this section the results of the application of the methods: AR, GA and GAS are presented and compared, and then evaluated from the perspective of the assumptions adopted, resulting in a smaller number of systems, with the greatest possible redundancy and the ease of replication. Figure 2 shows the results of applying the GA method, which requires 11 sensor system to cover all the 41 mines.

The results of using GAS method, shown in Figure 3, indicate that 10 sensor system are required to fully cover the mines.

The results from using AR method presented in Figure 4, indicate that 10 sensor system are required to cover the total of mines.

To evaluate the system redundancy premise, we calculate the percentage of mines covered by more than one sensor system or redundancy rate. In order to do this, we identified, in the mine coverage potential table, the mines where the systems were allocated. Next, we select the lines corresponding to these mines, identifying the mines covered by more than one sensor system. From these cases we calculate the redundancy rates for each method. Therefore, GA, GAS and AR methods have respectively 32%, 22% and 32% of the mines covered by more than one sensor system.

Finally, to evaluate the premise of the ease of method replication we calculate the number of iterations required to process them and execute them. The number of iterations can be understood as a relation between effort (time) and complexity (risk of error). This relation is verified in Figure 5, in which we see the evolution of the number of interactions that each method demands to find the solution of the problem, where is the number of sensor system and the number of interactions.

Google maps.

Figure 2 Representing GA method applied.  

Google maps

Figure 3 Representing GAS method applied 

Google maps

Figure 4 Representing AR method applied.  

Figure 5 Iterations and optimum solution per method. 

Considering that this iteration requires the approximate effort of 3 minutes, we estimate that the time taken to solve the problem through the GA, GAS and AR methods is 33, 165 and 57 minutes, respectively.

The results from the application of each method, which cover all the mines regarding the premises of minimizing cost, greater redundancy, and greater ease of replication, are summarized in table 5.

Table 5 Methods comparison analysis 

Premise GA GAS AR
Solution cost 11 10 10
Redundancy rate 32 22 32
Replication difficulty 11 55 19

Following Salhi's (2017) understanding, it was demonstrated that methods designed to deal with specific boundary conditions of the problem of interest produce better results.

The evolution of the number of iterations, which reflects the effort or complexity required for the execution of the AR method is linear with respect to the number of alternatives, allowing its replication for problems with a greater number of mines. Among the three methods, AR is the one with the lowest standard deviation between the systems installation sites. This is a consequence of a better geographic distribution of the systems. Therefore, the AR was obtained, with fewer systems in relation to GA, and a higher rate of redundancy. Finally, the AR was also efficient in the minimum equipment allocation. We found by AR, with greater ease and redundancy, a solution that minimizes the costs of the system implantation equals to that found in the GAS method.

5.1. Contributions, limitations, and future research

This research brings general contributions that will be described in this section. First, it reveals that heuristic methods can be constructed according to the boundary conditions of each problem (Salhi, 2017), but that this construction requires some type of programming (e.g. Blesa & Blum, 2018; Ivanov et al., 2019; Kuschel et al., 2019) that limit a broader replication of these methods in locational problems. Second, it offers a heuristic method that: 1) can be easily replicated without any programming; 2) produces an optimal solution to an area coverage problem with minimal resources; 3) can be adapted to similar locational problems, such as, for example, hospital location (e.g. Salgado et al., 2019), location of children’s schools (e.g. Libório, et al., 2020), location of supermarkets (e.g. Baños, Wandosell & Parra, 2016), micro power plants location optimization (e.g. Borges et al., 2019) among others; and 4) stands as an alternative to methods traditionally used in the area of Geography e.g. Reilly’s (1931) Gravity Point, Converse’s (1949) Inverse Distance Weighting and, Applebaum’s (1966) Area of Influence.

Two specific contributions can be cited. First, the application of the AR method reduces the mining company’s costs with TWS to protect workers and equipment exposed to lightning discharges. Second, cheaper TWS are of interest to governments concerned with protecting populations and public facilities in urban areas where there is a greater frequency and quantity of lightning discharges due to the presence of heat islands and pollution (Stallins & Rose, 2008; Farias et al., 2009; Sokol & Rohli, 2018).

This research has three main limitations. First, the method did not consider the relevance of upward electrical discharges in triggering the system alerts. Second, there is no redundancy constraints. Third, although the number of iterations to obtain the optimal solution by the AR method does not grow exponentially as by the GAS method, processing a large number of alternatives can be quite laborious. Exploring solutions for these limitations would be an interesting suggestion for future works. Researchers are also encouraged to implement the AR method for other problems such as environmental, engineering, logistics, marketing or urban problems.

6. Conclusions

The objective of this research was to develop a method for reducing costs of preventive lightning systems in an area containing 41 mines, considering that: the system's sensors have a range of 10 km; all mines must be covered by at least one sensor; the total number of sensors to be installed must be as small as possible; the method must be easy to replicate and should not require any programming; the redundancy of sensors per mine is not mandatory, but it is a desirable feature.

The achieved results show that the AR method developed in this research fulfilled the proposed goals. First, because all 41 mines were covered by at least one sensor. Second, because the solution obtained by the AR presented the best arrangement of sensors number, redundancy rate and replication difficulty. Third, because AR can be easily replicated without the need for programming, using spreadsheets such as Microsoft Excel. The number of sensors that determine the cost of the lightning detection system decreases from 11 to 10 sensors when using the AR and GAS methods. The sensor redundancy rate per mine increases from 22 to 32 when using the AR method instead of the GAS method. The difficulty of replication, the relation between time effort and complexity or risk of error, falls from 55 iterations using the GAS method to 19 iterations using the AR method.

Acknowledges

This work was supported by the Coordination for the Improvement of Higher Education Personnel (CAPES) under Grant 0001; Brazilian National Council for Scientific and Technological under Grant 311032/2016-8; and Companhia Vale - Brazil (Vale), Research in Planning and Decision Making Support under Grant P&D 2019-2022.

References

Applebaum, W. (1966). Methods for determining store trade areas, market penetration, and potential sales. Journal of Marketing Research, 127-141. [ Links ]

Baños, R., Wandosell, G., & Parra, M. C. (2016). Web GIS to enhance relational capital: the case of general merchandise retailers. Journal of Knowledge Management, 20(3), 578-593. [ Links ]

Bergantiños, G., Gómez-Rúa, M., Llorca, N., Pulido, M., & Sánchez-Soriano, J. (2020). Allocating costs in set covering problems. European Journal of Operational Research, 284(3), 1074-1087. [ Links ]

Berger, G. (2007). Lightning-caused accidents and injuries to humans. In Proceeding from IX International Symposium on Lightning Protection. Foz do Iguacu. [ Links ]

Blesa, M. J., & Blum, C. (2018). On Solving the Maximum Disjoint Paths Problem with Ant Colony Optimization. Handbook of Approximation Algorithms and Metaheuristics: Contemporary and Emerging Applications, 2, 1. [ Links ]

Borges, V. V., Libório, M. P., Silva, I. S., Laudares, S., de Oliveira, B. A. G., & Silva Martins, C. A. P. (2019). Ferramenta de geovisualização web para a implantação de centrais geradoras hidrelétricas. Espaço Energia, 30, 1-8. [ Links ]

Church, R., & ReVelle, C. (1974, December). The maximal covering location problem. In Papers of the Regional Science Association (Vol. 32, No. 1, pp. 101-118). Springer-Verlag. [ Links ]

Converse, P. D. (1949). New laws of retail gravitation. Journal of marketing, 14(3), 379-384. [ Links ]

Cooper, L. (1963). Location-allocation problems. Operations research, 11(3), 331-343. [ Links ]

Cooray, V. (2009). Lightning protection. The Institution of Engineering and Technology. [ Links ]

Expósito-Márquez, A., Expósito-Izquierdo, C., Brito-Santana, J., & Moreno-Pérez, J. A. (2019). Greedy randomized adaptive search procedure to design waste collection routes in La Palma. Computers & Industrial Engineering, 137, 106047. [ Links ]

Farahani, R. Z., Asgari, N., Heidari, N., Hosseininia, M., & Goh, M. (2012). Covering problems in facility location: A review. Computers & Industrial Engineering, 62(1), 368-407. [ Links ]

Farias, W. R. G., Pinto Jr, O., Naccarato, K. P., & Pinto, I. R. C. A. (2009). Anomalous lightning activity over the Metropolitan Region of São Paulo due to urban effects. Atmospheric Research, 91(2-4), 485-490. [ Links ]

Field, K. (2018). Cartography. Esri Press. [ Links ]

Gill, P. E., Murray, W., & Wright, M. H. (2019). Practical optimization. Society for Industrial and Applied Mathematics. [ Links ]

Gomes, C., Ab Kadir, M. Z. A., & Cooper, M. A. (2012, September). Lightning safety scheme for sheltering structures in low-income societies and problematic environments. In 2012 International Conference on Lightning Protection (ICLP) (pp. 1-11). IEEE. [ Links ]

HassanzadehFard, H., & Jalilian, A. (2018). Optimal sizing and location of renewable energy based DG units in distribution systems considering load growth. International Journal of Electrical Power & Energy Systems, 101, 356-370. [ Links ]

Ivanov, D., Tsipoulanidis, A., & Schönberger, J. (2019). Facility location planning and network design. In Global Supply Chain and Operations Management (pp. 155-202). Springer, Cham. [ Links ]

Kuschel, H., Cristallini, D., & Olsen, K. E. (2019). Tutorial: Passive radar tutorial. IEEE Aerospace and Electronic Systems Magazine, 34(2), 2-19. [ Links ]

Laudares, S., Laudares, J. B., & Libório, M. P. (2016). Geographic Information Systems Used as a Practical Way of Teaching Mathematics. Journal of Geographic Information System, 8(5), 608-617. [ Links ]

Leite, D. M., Modena, J., & Sueta, H. E. (2007). Present situation of structures lightning protection standardization and comparisons between the valid Brazilian standard and the 62305 IEC series. In IX International Symposium on Lightning Protection (pp. 1-5). [ Links ]

Libório, M. P., Bernardes, P., Ekel, P. I., Ramalho, F. D., & Santos, A. C. G. (2020). Geomarketing e o problema da questão locacional nos estudos de marketing. Revista Brasileira de Marketing, 19(2), 448-469. [ Links ]

Lopez, J., Perez, E., Herrera, J., Aranguren, D., & Porras, L. (2012, September). Thunderstorm warning alarms methodology using electric field mills and lightning location networks in mountainous regions. In 2012 International Conference on Lightning Protection (ICLP) (pp. 1-6). IEEE. [ Links ]

Lundquist, S. A. S., & Scuka, V. (1971). U.S. Patent No. 3,611,365. Washington, DC: U.S. Patent and Trademark Office. [ Links ]

Magina, F. C., Pinto Jr, O., & Naccarato, K. P. (2016). Atmospheric Electric Field Sensors Network Integration in Brazil. IEEE Latin America Transactions, 14(7), 3056-3064. [ Links ]

Mathew, G., & Surana, A. (2012, December). A static coverage algorithm for locational optimization. In 2012 IEEE 51st IEEE Conference on Decision and Control (CDC) (pp. 806-811). IEEE. [ Links ]

Nguyen, L. V., Kodagoda, S., Ranasinghe, R., & Dissanayake, G. (2015). Information-driven adaptive sampling strategy for mobile robotic wireless sensor network. IEEE Transactions on Control Systems Technology, 24(1), 372-379. [ Links ]

Ortiz-Astorquiza, C., Contreras, I., & Laporte, G. (2018). Multi-level facility location problems. European Journal of Operational Research, 267(3), 791-805. [ Links ]

Owen, S. H., & Daskin, M. S. (1998). Strategic facility location: A review. European journal of operational research, 111(3), 423-447. [ Links ]

Pinto Jr, O. (2005). A arte da guerra contra os raios. Oficina de Textos. [ Links ]

Reilly, W. J. (1931). The law of retail gravitation. WJ Reilly. [ Links ]

Salgado, R. R., Libório, M. P., Silva Martin, C. A. P., & Bernardes, P. (2019). WebGIS de geomarketing: solução de dois estágios para um problema de localização do varejo. Raega-O Espaço Geográfico em Análise, 46(1), 75-89. [ Links ]

Salhi, S. (2017). Heuristic search: The emerging science of problem solving. Springer. [ Links ]

Sokol, N. J., & Rohli, R. V. (2018). Land cover, lightning frequency, and turbulent fluxes over Southern Louisiana. Applied geography, 90, 1-8. [ Links ]

Stallins, J. A., & Rose, L. S. (2008). Urban lightning: Current research, methods, and the geographical perspective. Geography Compass, 2(3), 620-639. [ Links ]

Tong, D., & Murray, A. T. (2017). Location analysis: Developments on the horizon. In Regional Research Frontiers -Vol. 2 (pp. 193-208). Springer, Cham. [ Links ]

Visacro F. S. (2005). Descargas atmosféricas: uma abordagem de engenharia. Artliber. [ Links ]

Weber, A. (1929). Theory of the Location of Industries. University of Chicago Press. [ Links ]

Received: July 30, 2020; Accepted: November 19, 2020

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License