Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 507953, 11 pages

http://dx.doi.org/10.1155/2013/507953

Research Article

A Capacitated Location-Allocation Model for Flood Disaster Service Operations with Border Crossing Passages and Probabilistic Demand Locations

Seyed Ali Mirzapour,¹ Kuan Yew Wong,¹ and Kannan Govindan²

¹Department of Manufacturing and Industrial Engineering, Faculty of Mechanical Engineering, Universiti Teknologi Malaysia, 81310 UTM Skudai, Malaysia

²Department of Business and Economics, University of Southern Denmark, 5230 Odense, Denmark

Received 16 August 2013; Revised 14 October 2013; Accepted 14 October 2013

Academic Editor: Tsan-Ming Choi

Copyright © 2013 Seyed Ali Mirzapour et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Potential consequences of flood disasters, including severe loss of life and property, induce emergency managers to find the appropriate locations of relief rooms to evacuate people from the origin points to a safe place in order to lessen the possible impact of flood disasters. In this research, a p-center location problem is considered in order to determine the locations of some relief rooms in a city and their corresponding allocation clusters. This study presents a mixed integer nonlinear programming model of a capacitated facility location-allocation problem which simultaneously considers the probabilistic distribution of demand locations and a fixed line barrier in a region. The proposed model aims at minimizing the maximum expected weighted distance from the relief rooms to all the demand regions in order to decrease the evacuation time of people from the affected areas before flood occurrence. A real-world case study has been carried out to examine the effectiveness and applicability of the proposed model.

1. Introduction

Due to various occurrences of disasters such as floods, hurricanes, and earthquakes which lead to enormous property damages and human injuries, disaster management has recently become a crucial issue. Within disaster management tasks, finding the appropriate locations of relief rooms can help emergency managers to ensure public safety and well-being in tragic situations. Thus, transportation specialists play a prominent role in the emergency management process which focuses on people evacuation from disaster events in order to mitigate injury and loss of life. In this regard, there has been an increasing interest among practitioners and scholars in the field of emergency facility location problems. Compared to all disasters, flood is more possible to be predicted and prevented. Therefore, it is vital for flood disaster managers to find the appropriate locations of relief rooms or emergency medical services and allocate them to populations in order to provide efficient services.

A restricted region is a limitation in an area in which the geographic features obstruct the construction of relief rooms. Generally, a restricted planar area is divided into three categories: forbidden regions, congested regions, and regions where neither placement nor travelling is permitted through them. Lakes, mountains, highways, and military zones are some practical examples of these restricted regions. Hamacher and Nickel [1] have performed an extensive review of facility location problems with restrictions. These limitations were incorporated by Klamroth [2] into a model in which a fixed line barrier in a region divides it into two half-planes and two passages located along the line barrier provide communication between these subregions. The model proposed by Klamroth [2] is a specific formulation, and it may not be appropriate for cases that have regional probability distributions of customers or different objective functions. In the real world, when we are dealing with establishing permanent emergency service facilities for future disasters over a long term planning horizon, it is essential to consider the regional probability distribution of customers due to the uncertainties of events and partial information and data. Therefore, a more applicable model which provides better solutions would be beneficial for overcoming this weakness.

To the best of our knowledge, few research activities which consider barriers in emergency facility location problems have been accomplished, due to the computational complexities associated with these problems. Moreover, far too little attention has been paid to developing a practical mathematical model of emergency facility location that considers probabilistic customer regions and a line barrier with border crossings simultaneously. In this research, in order to address the abovementioned problems, we develop a mixed integer nonlinear programming (MINLP) model of a capacitated facility location-allocation problem (CFLAP) in an area in which a fixed line barrier and probabilistic customer regions are taken into account.

The rest of this paper is organized as follows: Section 2 provides a comprehensive literature review on related topics. Section 3 gives an explanation of the problem, and Section 4 presents the proposed mathematical model. In order to show the applicability of the proposed model, a suitable case study is illustrated in Section 5, and the model is validated in Section 6. Finally, Section 7 presents some concluding remarks and outlines a number of future research directions.

2. Literature Review

2.1. Emergency Facility Location Problems

The emergency facility location models are categorized into three different types: P-median, P-center, and maximal covering location problems. Since the objective of this research is to design an effective response strategy for disaster management organizations in order to reduce casualties, a planar center location problem (PCLP) would best suit our purpose. Actually, whenever the average cost of servicing the customers is less important than ensuring that no customer receives poor quality of services, the minimax location problem is typically proposed [3]. The procedure of the PCLP model, which is in fact a minimax problem, looks for the center in order to minimize the maximum service time, cost, or loss. Research in this area started with Sylvester [4] which was the first paper that considered the Weber problem with the PCLP model. Elzinga and Hearn [5] then introduced the Euclidean center problem with equal weights and offered an efficient solution technique for their proposed model. Francis et al. [6] made it possible to stress the importance of different customers by assigning unequal weights to them. Mehrez et al. [7] applied a single objective optimization to determine the location of a new hospital. The solutions of this problem were generated by a qualitative approach based on the evaluation of experts’ judgments. After that, Daskin et al. [8] proposed an error-bound driven demand point aggregation for the minimax problem with rectilinear distance. Hurtado et al. [9] came out with some constraints on the minimax problem and used a time algorithm to solve it.

In the last two decades, the field of disaster management has gained increasing attention. Many authors have conducted several studies on predisaster preparedness and postdisaster responsiveness, including Altay and Green III [10], Simpson and Hancock [11], Ichoua [12], and Shishebori and Jabalameli [13]. Caunhye et al. [14] and Galindo and Batta [15] performed an extensive review of optimization models in emergency logistics for disaster operations management. Sherali et al. [16] proposed a capacitated facility location-allocation model in order to minimize the congested evacuation time by finding the optimal locations of evacuation shelters under hurricane or flood conditions. Chang et al. [17] studied the problem of locating relief rooms and allocating relief resources for flood emergency preparation. They proposed two stochastic models under possible flood scenarios in which the demand locations are uncertain and depend on the level of flood in different scenarios. In addition, Sheu [18] presented an earthquake relief model that forecasts the demand of affected regions and coordinates the provision of relief supplies. Mete and Zabinsky [19] dealt with an optimization approach for locating and distributing medical supplies under demand and transportation uncertainties for disaster situations. Some other recent relevant studies on this issue include Jiang et al. [20], Rawls and Turnquist [21], Yi et al. [22], Halper and Raghavan [23], Wu et al. [24], Xu et al. [25], and Duran et al. [26].

2.2. Facility Location Problems under Uncertainty

When addressing the strategic decision of establishing permanent facilities to prepare for future disasters over a long term planning horizon, information is uncertain due to the time lag; so considering the location of demands as a known point with certainty leads to poor model performance [27]. The first research in the probabilistic Weber problem was done by Cooper [28] who focused on a single facility location problem. After that, Katz and Cooper [29] approximated the expected distance of transportation via Euclidean distance in which the customer locations were assumed to have a normal distribution. They proposed an algorithm to minimize the total expected distance between the new facility and existing facilities. Later on, for the case of rectilinear distance, Wesolowsky [30] considered three different probability distributions of bivariate normal, bivariate symmetric exponential, and bivariate uniform for customer locations in solving the probabilistic Weber problem.

Although some studies have recently been done for the case of multifacility minimax location problems, much less attention has been given to the probabilistic formulation of minimax problems. Stochastic demand locations were first considered by Carbone and Mehrez [31] for the minimax location problem where the locations of demand points were normally distributed with predetermined means and variances. An emergency location problem with a weighted minimax objective function has been presented by Berman et al. [32] in which the weights associated with each demand point were not given and each of them was assumed to have a uniform distribution. They proved that the objective function is convex for parameters of the uniform distribution, and, thus, the model can be solved by using standard optimization techniques. The minimax regret location problem has been investigated by Averbakh and Bereg [33] in the case of uncertain weights and coordinates of demand points. They proposed a model with two objective functions of the median and center problems by using rectilinear distance. Foul [34] applied a mathematical model with the minimax objective function to find the unique location where the demand points follow a bivariate uniform distribution. Then, Pelegrín et al. [35] proposed a comprehensive framework for the 1-center problem where the weights of demands have an arbitrary probability distribution. They have studied a variety of problems using a number of objective functions.

After that, Durmaz et al. [36] presented a model on the probabilistic capacitated multifacility Weber problem in a two-dimensional region by assuming that the locations of customers have some multivariate probability distributions. Later on, Canbolat and Von Massow [37] shared their work on locating emergency facilities with random demands. They proposed an approach that explicitly addresses the uncertainty with respect to where an emergency will occur by minimizing the maximum risk through the minimization of the maximum expected distance. Then, Hosseinijou and Bashiri [38] formulated an appropriate model to represent the minimax single facility location problem in which demand areas are weighted and their coordinates have a bivariate uniform distribution. They proposed a model which finds the optimal location of the transfer point such that the maximum expected weighted distance from the facility point to all demand points through the transfer point is minimized.

2.3. Facility Location Problems with Barriers

In many cases of facility location problems, difficulties in solution methodologies occur when barriers limit the establishment of new facilities in a specified region. The Weber problem with a circular barrier was first presented by Katz and Cooper [39]. Since then, researchers have studied a variety of location problems in the presence of barriers using different types of objective function. They showed that the objective functions of these problems are nonconvex and offered some heuristic methods for solving them. Klamroth and Wiecek [40] focused on solving a multiobjective median problem (MOMP) with the minisum objective function by assuming a line barrier in an area and showed that their proposed model is nonconvex. Many researchers have applied different heuristic algorithms for the Weber problem that considers the convex or nonconvex polyhedral barriers in a region (e.g., [41–43]). The location problem has also been extended by Klamroth [44] who suggested that the single circular barrier can be subdivided into some feasible convex regions. By using this method, the nonconvex objective function of the problem has been transformed into some convex objective functions in each region. When the number of subdivided convex regions increases, this method is not suitable, since dividing the convex regions becomes awkward. To overcome this complexity, a genetic algorithm has been proposed by Bischoff and Klamroth [45] based on the Weiszfeld technique to solve the model with a large number of demand points. Bischoff et al. [46] presented two alternative location-allocation heuristics for the multifacility Weber problem with barriers.

The locations of optimum points for a given number of passages in two cases of capacitated and incapacitated facility location problems have been discussed by Huang et al. [48]; their model tried to minimize the transportation cost. Increasingly, more authors have addressed the planar center location problem in the presence of barriers in which the objective function is to find a facility that minimizes the maximum distance (e.g., [49–51]). A different perspective of using a wave front approach was proposed by Frieß et al. [52] for the minimax location problem in the presence of barriers. Canbolat and Wesolowsky [53] developed a model by considering a probabilistic line barrier that occurs randomly on the horizontal route. Their objective was to locate a single new facility to minimize the sum of the expected distances from all customer locations. After that, Canbolat and Wesolowsky [54] presented a new experimental approach for the Weber problem by using the Varignon frame. Their proposed method tried to minimize the sum of the weighted distances from the facility to customer points by considering barriers in a region. Then, Canbolat and Wesolowsky [47] presented a planar single facility location problem in the presence of border crossing. They developed the models for minisum and minimax problems in the presence of a fixed line barrier which divides a region into two subregions.

In short, although previous studies have proposed different approaches to address the facility location-allocation problem in the presence of a barrier, they have somehow neglected the probabilistic nature of customer locations. Moreover, very limited research activities have been carried out on multifacility location problems; albeit in many real-world applications, more than a single facility is required to be established. Therefore, an appropriate mathematical model is proposed in this paper in order to address these issues.

3. Problem Definition and Formulation

3.1. Facility Location Problem in the Presence of a Line Barrier

It is important to introduce some preliminary definitions before formally defining the problem. Based on Klamroth [2], the facility location problem in the presence of a line barrier is described as follows.

Let be a line and be a set of points on . Then, is called a fixed line barrier with passages. The feasible region is defined as the union of two closed half-planes and on both sides of . is the coordinate of a border crossing passage. Since the connections are located on a fixed horizontal line barrier, we set , . A set of existing demand regions with locations of , and positive weights of is given in . Then, is the -norm barrier distance between and in the presence of a fixed line barrier where communication between two subplanes is allowed only through the passages. The basic minimax model with barriers can be stated as follows:

Any demand region that is in the same half-plane will not be affected by the presence of the barrier. The other demand regions that are not in the same half-plane with the new facility will connect to the new facility through just one of the passages along the line barrier. For this problem, the rectilinear distance, , is considered to compute the distance between the new emergency facilities and all demand regions.

3.2. Stochastic Weighted Emergency Facility Location Problem in the Presence of a Line Barrier

Assume that each demand region has a coordinate of such that , are independent random variables and each of them has a bivariate uniform distribution. This assumption is realistic in real-world cases wherein the geographical distribution of population in different regions is uniform in square or rectangular areas [38]. Furthermore, a uniform distribution provides a “building block” to deal with a situation which is not precisely approximated [34].

The problem is to find the optimum location of the new facility such that the maximum expected weighted distance from the facility to all demand regions, , is minimized:

where represents the expected value. The above mentioned problem is illustrated in Figure 1. The weights of emergency needs in each demand region are defined as a function of two factors. The first factor is related to the qualitative criterion of demand region which can be the quality of road construction from a demand region to the new facility point whereas the next one is concerned with the quantitative criterion of demand region which can be a function of population, potential of emergency situation occurrence, and other issues. The weight of the qualitative factor has an inverse relation with which means that will reduce whenever the quality of road construction increases and this causes the new facility to be located farther from demand region . The weight of the quantitative factor has a direct relation with ; for instance, whenever the population of demand region increases, will increase too and the new emergency facility will tend to be located nearer to this customer area. In order to use both quantitative and qualitative factors in the total weight of demand region , they must be defined in the same range to be comparable. One method of doing so is normalizing both types of factors to be in the domain of zero to one. Based on the experience of emergency managers, both quantitative and qualitative factors have the same significance.

Figure 1: Demand locations with a fixed line barrier and two passages on the plane.

In this paper, the case of uniformly distributed demand regions in several rectangles is considered. It uses the rectilinear distance concept for measuring distances. Imagine a number of demand regions in which , are independent random variables in and with probability density functions and , respectively. Each area is a demand region, and it is desirable to find the optimum locations of emergency services among all these regions. The objective is to find the optimal locations of emergency facilities to provide efficient services for people. The mathematical expression for the expected distance between new facility and demand region is presented as

Let and be denoted by and , respectively. Then, we have

The closed-form expression for the expected distance was developed by Foul [34], which can be written as

can be defined like by substituting , , with , , , respectively. With this expression, the problem is formulated as follows:

Equation (7) can be transformed into an equivalent constrained nonlinear program which is presented here:

4. Proposed Mathematical Model

In this section, a MINLP model is proposed for CFLAP. The proposed MINLP model aims at finding the suitable locations of some emergency facilities in the presence of a fixed line barrier such that the expected weighted rectilinear distance between the emergency facilities and the farthest demand regions is minimized. The following assumptions and notations are used in formulating the MINLP model for CFLAP.

Assumptions. (i) The demand regions’ coordinates have a bivariate uniform distribution because customer populations (residential areas) are normally dispersed or structured in square or rectangular areas.

(ii) There is a line barrier in the area which is fixed.

(iii) Two emergency service facilities are required to be established. This will ensure a more balanced allocation of the demand regions to the emergency service facilities since the entire region is divided into two half-planes by a line barrier.

Notations

:Number of new facilities:Number of demand regions or customer regions:Number of passages:Positive weight of each existing demand region :Binary variable: 1, if demand region is assigned to new facility , 0, otherwise:Binary variable: 1, if demand region is served by new facility through passage , 0, otherwise:Binary variable: 1, if demand region and new facility are located in different half-planes, 0, otherwise:Binary variable: 1, if new facility is located in the upper half-plane, 0, otherwise:Binary parameter: 1, if demand region is located in the upper half-plane, 0, otherwise:Maximum capacity of new facility :Binary variable: 1, if , 0, otherwise:Binary variable: 1, if , 0, otherwise:Binary variable: 1, if , 0, otherwise:Binary variable: 1, if , 0, otherwise:Binary variable: 1, if , 0, otherwise:Binary variable: 1, if , 0, otherwise:Binary variable: 1, if , 0, otherwise:Binary variable: 1, if , 0, otherwise:Binary variable: 1, if , 0, otherwise:Binary variable: 1, if , 0, otherwise:Binary variable: 1, if , 0, otherwise:Binary variable: 1, if , 0, otherwise.

Based on the above assumptions and notations, the objective function and constraints will be described.

4.1. Objective Function

The objective function tries to minimize while satisfying constraint (9). Since needs to be minimized, the shortest possible path through one of the border crossings along the barrier will be selected. The objective function can be formed as follows:

Note that the expected distance from the passage coordinates along the line barrier to the demand regions is denoted by and which are equal to and and they can be defined like and by substituting with and with , respectively.

4.2. Constraints

The constraints of this problem are formulated as follows.

4.2.1. Expected Distance

Based on the obtained closed-form expression in (6), the expected distance is computed in different ways depending on the location of the new facility and the interval of demand regions. Since this model will be solved by the GAMS software (full academic version 22.1), it is required to convert the formulas of into a few constraints by introducing new binary variables of , , and . This makes it easier for the software to perform the computation process. Likewise, new binary variables of , , and for are introduced. and can be manipulated like and by replacing , , with , , and , , with , , , respectively. These constraints can be written as follows:

The first three constraints (constraints (10), (11), and (12)) are related to . Whenever the variable is less than (), the binary variable in constraint (10) will be equal to one; otherwise, it makes the binary variable to be equal to zero. Constraint (11) ensures that the expected distances between facilities and demand regions are calculated by the third formula of (6) whenever the variable is more than (); otherwise will be equal to zero. If neither nor is equal to one, then the binary variable must be equal to one due to constraint (12). Constraint (12) also guarantees that just one of these cases can occur for the same new facility and customer region; assignment of two or three binary variables to one will lead to infeasibility and is not allowed. The rest of the constraints (constraints (13) to (21)) follow the same interpretation for , , and , respectively.

4.2.2. The Position of New Facilities and Customer Regions

Constraint (22) indicates whether the new facility is located in the upper half-plane or not:

Constraint (23) shows whether the customer region and the new facility are both located in the same half-plane or not:

Constraints (22) and (23) help the objective function to compute the expected distance between the new facilities and the customer regions in different cases.

4.2.3. Demand Allocation

According to constraint (24), each demand region must be served just by one of the new facilities:

The next constraint assures that the shortest path through just one of the passages will be selected to serve the allocated customer region in a different half-plane:

4.2.4. Capacity

Constraint (26) guarantees that each new facility cannot serve more than its corresponding capacity:

The modified objective function with respect to constraints (10) to (21) looks as follows:

5. A Case Illustration

To examine the practicality and effectiveness of the proposed model for CFLAP, a case is illustrated for the emergency facility location problem. One of the flood disaster management activities in the Amol city, Mazandaran province, northern part of Iran, was selected as a case study. In recent years, severe floods have occurred in the northern part of Iran. In 1999, floods in this area left 34 dead and 15 missing after more than a week and presumed dead. More than 100 people were injured and various areas were damaged with estimated losses of 15 million dollars [55]. Recently, the frequency of flood occurrence in the Mazandaran province, Iran, has increased and an appropriate model is required to find the suitable locations of relief rooms in order to decrease the evacuation time of people from the affected areas to a safe place before flood occurrence to reduce the possible impact of flood disasters.

The Amol city is located at 36°23′N latitude and 52°20′E longitude in northern Iran. This city is situated on the Haraz river bank, and it is less than 20 kilometers south of the Caspian sea and 10 kilometers north of the Alborz mountain. Based on the 2010 census, its population was around 224,000 with 61,085 families. As shown in Figure 2, the Haraz river divides the city into two subplanes, and two bridges are located along the river to connect the two subregions with each other. With respect to the geographic information system (GIS) map of the city, it has been divided into ten demand areas. As we mentioned before, the weights of emergency needs in each demand region are a function of both qualitative and quantitative criteria. The population of each area, as a quantitative criterion, was obtained from the statistics center of Iran and is shown in Table 1; the data for quality of road construction were gained from the GIS map of Amol which was prepared by the Department of Housing and Urban Development of Mazandaran Province. Table 2 presents the data for the qualitative factor.

Table 1: Population data in different regions.

Population

1
13464
0.057132

2
26139
0.110916

3
13984
0.059339

4
25432
0.107916

5
23916
0.101483

6
8888
0.037715

7
12972
0.055044

8
35072
0.148822

9
38392
0.162910

10
37405
0.158722

Table 2: Data for quality of road construction.

HQ (%)
GQ (%)
AQ (%)
MNR (%)
MJR (%)
Score (%) of each area

1
0.38
40.91
37.5
20.19
0
275.46
0.089512

2
1.65
39.76
35.76
11.53
9.41
281.62
0.091514

3
0
11.76
28.43
25.49
34.31
382.32
0.124236

4
1.37
6.65
23.29
39.33
29.35
388.61
0.126280

5
2.96
20
28.89
24.44
23.7
345.89
0.112398

6
0.72
11.15
40.66
24.75
22.7
357.5
0.116171

7
2.78
59.45
30.87
6.45
0
240.09
0.078018

8
5.97
68.1
20.71
3.36
1.37
224.59
0.072981

9
1.54
70.51
26.15
1.79
0
228.17
0.074145

10
15.79
14.04
14.04
10.53
45
353.11
0.114744

Figure 2: The illustration of the demand regions and passages.

For the sake of normalizing the quantitative weight of each demand region, the population number of each area is divided by the total population. For the qualitative factor, the quality of road construction is categorized into 5 groups which are high quality (HQ), good quality (GQ), acceptable quality (AQ), minor repair required (MNR), and major repair required (MJR). Since this factor is inversely related with , the best quality construction receives the significance of 1 and the worst quality construction gets the maximum significance score which is equal to 5. The weight of the qualitative factor is computed similarly as the quantitative criterion. Since the significance of both criteria is the same, the final weight of each area is obtained by adding together half of the weights of the qualitative and quantitative factors.

The coordinates of ten demand regions with their corresponding final weights are tabulated in Table 3, while Table 4 gives the coordinates of the passages. It is noted that the establishment of two relief rooms with a maximum capacity percentage of 0.545 and 0.46 is required to cover all the demand regions. The proposed MINLP model contains 292 integer variables and 4 noninteger variables in terms of its complexity. Hence, it has been implemented in the GAMS 22.1 software using the BARON solver in order to identify the optimal values of and . BARON was created at the University of Illinois which implements global optimization algorithms of a branch and bound type and solves nonconvex problems to obtain global solutions [56]. Computations were performed using a system with a Core 2 duo 2.1 GHz CPU and 2 GB RAM and the computation time for solving the model was about 11 minutes. The time complexity grows exponentially with the increasing number of passages or new facilities.

Table 3: Parameter values for the demand regions.

1

0.073322

2

0.101215

3

0.0917875

4

0.117098

5

0.1069405

6

0.076943

7

0.066531

8

0.1109015

9

0.1185275

10

0.136733

Table 4: The coordinates of passages.

1
6.9
8.5

2
10.1
8.5

The optimal locations of the emergency facilities and the value of the objective function are reported in Table 5. In addition, Table 6 shows the optimal corresponding allocation clusters of the two new emergency facilities. For a better understanding, Figure 3 illustrates the optimum locations of the new facilities and their corresponding allocation. Based on the results, we recommend that demand regions 1, 2, 8, 9, and 10 are served by the first new facility and the rest of the demand areas are served by the second new facility. The results will be useful for emergency managers to decrease the evacuation time in affected areas before flood occurrence.

Table 5: Optimal solutions for the locations of relief rooms.

9.8
13.312
8.6
6.139
0.682

Table 6: Optimal allocation of demand areas to the two new facilities.

1

0

1

0

0

1

0

1

0

1

0

1

0

1

1

0

1

0

1

0

Figure 3: Illustration of the computation results.

6. Validation of the Model

Model validation is essential for evaluating the proposed model to make sure that the chosen new facility locations truly minimize the maximum value of distance from the stochastic demand regions to the chosen new facility points within the defined area. Based on the validation procedure of Hosseinijou and Bashiri [38], ten demand points were randomly generated within each demand region. Then, the objective functions of the recommended emergency facility locations were compared with the objective function of the proposed model.

Some efforts of validating the model are presented here. Ten scenarios were created from five randomly selected areas in Table 3. A random data toolbox in the Minitab software was used for generating random point data in each area. The weight of each point in the scenarios is equal to the weight of its allocated area. Table 7 presents the randomly selected areas for validation of the proposed model. In each of the scenarios, the actual value of the objective function is computed since the customer location is a point. The maximum capacity percentage of the two new facilities is set as 0.425 and 0.58, and the optimum expected value of the objective function for the five selected areas is equal to 1.11. The validation results of the proposed model are tabulated in Table 8. The actual value of the objective function for each scenario is listed in column 7; column 8 shows the error of our model, which is the difference between the actual value in column 7 and the expected value of the objective function for the five selected regions.

Table 7: Coordinates of randomly selected areas.

1

0.159

3

0.198

5

0.234

7

0.145

9

0.264

Table 8: Validation results of the proposed model.

Scenario

× [11.4, 17.1]

× [5.9, 8]

× [4.7, 8]

× [8.9, 11.4]

× [13.7, 18.2]

1

1.160
0.05

2

1.147
0.037

3

1.092
0.018

4

1.189
0.079

5

1.049
0.061

6

1.425
0.315

7

1.062
0.048

8

1.179
0.069

9

1.144
0.034

10

1.087
0.023

The final analysis of the model validation is reported in Table 9. The first row is related to the average of the objective function values for the ten scenarios which are extracted from column 7 of Table 8. The second row indicates the average error between the expected objective function value of the proposed model, and the obtained average amount in the first row. The average error percentage of the model is presented in the last row. As can be seen in Table 9, the average error percentage is equal to 3.87% which shows that the proposed CFLAP model is a good estimator of the actual objective function value and its solution is near to optimal. Hence, the proposed stochastic model is valid.

Table 9: Analysis of validation results for the proposed model.

1.153

0.043

3.87%

In order to highlight the advantages of the proposed model as compared to the most relevant existing literature, a comparison has been made with the model developed by Canbolat and Wesolowsky [47]. The results are summarized in Table 10.

Table 10: Comparison results.

Canbolat and Wesolowsky [47]
This study

Deterministic coordinates of demand locations
Stochastic coordinates of demand locations

Single facility location
Multifacility location

Numerical example
Real-world application

As shown, the work of Canbolat and Wesolowsky [47] was dedicated to a single facility location problem where deterministic demand coordinates were utilized. In addition, their model was tested using a numerical example. In contrast, this study has developed a stochastic model that considers probabilistic demand locations for a multifacility location-allocation problem. The applicability of the proposed model has been evaluated using an actual case study, thus justifying its usefulness in the real-world context.

7. Conclusions

In this study, we have considered the planar center location-allocation problem for locating some emergency facilities in the presence of a fixed line barrier in a region. Although some studies have been done on facility location problems with barriers, little attention has been devoted to the probabilistic formulation of these problems where demand coordinates are distributed according to a bivariate probability distribution. In this paper, we have proposed a MINLP model for CFLAP by considering probabilistic demand regions and a fixed line barrier simultaneously, thus making the model more realistic. This study aims at finding the optimum locations of two relief rooms and their corresponding allocation which leads to minimizing the maximum expected weighted distance from the relief rooms to all the demand regions. A real-world case study has been carried out in order to show the effectiveness of the proposed model and the problem has been solved by the GAMS 22.1 software in order to find the optimum solution.

As an approximation, uniform distributions have been used to represent demand region locations. Thus, one possible area for future research is to formulate the problem with the exact probability distributions based on goodness-of-fit tests. In addition, considering the minimax regret objective function for hazardous location problems is another possible extension that can be studied.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

H. W. Hamacher and S. Nickel, “Restricted planar location problems and applications,” Naval Research Logistics, vol. 42, no. 6, pp. 967–992, 1995.

K. Klamroth, “Planar weber location problems with line barriers,” Optimization, vol. 49, no. 5-6, pp. 517–527, 2001.

A. Krishna, G. S. Dangayach, and R. Jain, “Service recovery: literature review and research issues,” Journal of Service Science Research, vol. 3, no. 1, pp. 71–121, 2011.

J. J. Sylvester, “A question in the geometry of situation,” Quarterly Journal of Pure and Applied Mathematics, vol. 1, no. 1, pp. 79–80, 1857.

J. Elzinga and D. W. Hearn, “Geometrical solutions for some minimax location problems,” Transportation Science, vol. 6, no. 4, pp. 379–394, 1972.

R. L. Francis, L. F. McGinnis Jr., and J. A. White, Facility Layout and Location: An Analytical Approach, Prentice-Hall, Englewood Cliffs, NJ, USA, 2nd edition, 1992.

A. Mehrez, Z. Sinuany-Stern, T. Arad-Geva, and S. Binyamin, “On the implementation of quantitative facility location models: the case of a hospital in a rural region,” Journal of the Operational Research Society, vol. 47, no. 5, pp. 612–625, 1996.

M. S. Daskin, S. M. Hesse, and C. S. Revelle, “α-Reliable p-minimax regret: a new model for strategic facility location modeling,” Location Science, vol. 5, no. 4, pp. 227–246, 1997.

F. Hurtado, V. Sacristán, and G. Toussaint, “Some constrained minimax and maximin location problems,” Studies in Locational Analysis, no. 15, pp. 17–35, 2000.

N. Altay and W. G. Green III, “OR/MS research in disaster operations management,” European Journal of Operational Research, vol. 175, no. 1, pp. 475–493, 2006.

N. C. Simpson and P. G. Hancock, “Fifty years of operational research and emergency response,” Journal of the Operational Research Society, vol. 60, no. 1, pp. S126–S139, 2009.

S. Ichoua, “Relief distribution networks: design and operations,” in Supply Chain Optimization, Design, and Management: Advances and Intelligent Methods, I. Minis, V. Zeimpekis, G. Dounias, and N. Ampazis, Eds., pp. 171–186, IGI Global, Hershey, Pa, USA, 2011.

D. Shishebori and M. S. Jabalameli, “Improving the efficiency of medical services systems: a new integrated mathematical modeling approach,” Mathematical Problems in Engineering, vol. 2013, Article ID 649397, 13 pages, 2013.

A. M. Caunhye, X. Nie, and S. Pokharel, “Optimization models in emergency logistics: a literature review,” Socio-Economic Planning Sciences, vol. 46, no. 1, pp. 4–13, 2012.

G. Galindo and R. Batta, “Review of recent developments in OR/MS research in disaster operations management,” European Journal of Operational Research, vol. 230, no. 2, pp. 201–211, 2013.

H. D. Sherali, T. B. Carter, and A. G. Hobeika, “A location-allocation model and algorithm for evacuation planning under hurricane/flood conditions,” Transportation Research Part B, vol. 25, no. 6, pp. 439–452, 1991.

M.-S. Chang, Y.-L. Tseng, and J.-W. Chen, “A scenario planning approach for the flood emergency logistics preparation problem under uncertainty,” Transportation Research Part E, vol. 43, no. 6, pp. 737–754, 2007.

J.-B. Sheu, “An emergency logistics distribution approach for quick response to urgent relief demand in disasters,” Transportation Research Part E, vol. 43, no. 6, pp. 687–709, 2007.

H. O. Mete and Z. B. Zabinsky, “Stochastic optimization of medical supply location and distribution in disaster management,” International Journal of Production Economics, vol. 126, no. 1, pp. 76–84, 2010.

W. Jiang, L. Deng, L. Chen, J. Wu, and J. Li, “Risk assessment and validation of flood disaster based on fuzzy mathematics,” Progress in Natural Science, vol. 19, no. 10, pp. 1419–1425, 2009.

C. G. Rawls and M. A. Turnquist, “Pre-positioning of emergency supplies for disaster response,” Transportation Research Part B, vol. 44, no. 4, pp. 521–534, 2010.

P. Yi, S. K. George, J. A. Paul, and L. Lin, “Hospital capacity planning for disaster emergency management,” Socio-Economic Planning Sciences, vol. 44, no. 3, pp. 151–160, 2010.

R. Halper and S. Raghavan, “The mobile facility routing problem,” Transportation Science, vol. 45, no. 3, pp. 413–434, 2011.

L. Wu, Q. Tan, and Y. Zhang, “Delivery time reliability model of logistic network,” Mathematical Problems in Engineering, vol. 2013, Article ID 879472, 5 pages, 2013.

J. Xu, J. Gang, and X. Lei, “Hazmats transportation network design model with emergency response under complex fuzzy environment,” Mathematical Problems in Engineering, vol. 2013, Article ID 517372, 22 pages, 2013.

S. Duran, M. A. Gutierrez, and P. Keskinocak, “Pre-positioning of emergency items for CARE international,” Interfaces, vol. 41, no. 3, pp. 223–237, 2011.

W. Klibi, S. Ichoua, and A. Martel, “Prepositioning emergency supplies to support disaster relief: a stochastic programming approach,” Centre interuniversitaire de recherche sur les réseaux d'entreprise, la logistique et le transport (CIRRELT), Canada, 2013

L. Cooper, “A random locational equilibrium problem,” Journal of Regional Sciences, vol. 14, no. 1, pp. 47–54, 1974.

I. N. Katz and L. Cooper, “An always-convergent numerical scheme for a random locational equilibrium problem,” SIAM Journal of Numerical Analysis, vol. 11, no. 4, pp. 683–692, 1974.

G. O. Wesolowsky, “The weber problem with rectangular distances and randomly distributed destinations,” Journal of Regional Sciences, vol. 17, no. 1, pp. 53–60, 1977.

R. Carbone and A. Mehrez, “The single facility minimax distance problem under stochastic location of demand,” Management Science, vol. 26, no. 1, pp. 113–115, 1980.

O. Berman, J. Wang, Z. Drezner, and G. O. Wesolowsky, “A probabilistic minimax location problem on the plane,” Annals of Operations Research, vol. 122, no. 1–4, pp. 59–70, 2003.

I. Averbakh and S. Bereg, “Facility location problems with uncertainty on the plane,” Discrete Optimization, vol. 2, no. 1, pp. 3–34, 2005.

A. Foul, “A 1-center problem on the plane with uniformly distributed demand points,” Operations Research Letters, vol. 34, no. 3, pp. 264–268, 2006.

B. Pelegrín, J. Fernández, and B. Tóth, “The 1-center problem in the plane with independent random weights,” Computers & Operations Research, vol. 35, no. 3, pp. 737–749, 2008.

E. Durmaz, N. Aras, and İ. K. Altınel, “Discrete approximation heuristics for the capacitated continuous location—allocation problem with probabilistic customer locations,” Computers & Operations Research, vol. 36, no. 7, pp. 2139–2148, 2009.

M. S. Canbolat and M. Von Massow, “Locating emergency facilities with random demand for risk minimization,” Expert Systems with Applications, vol. 38, no. 8, pp. 10099–10106, 2011.

S. A. Hosseinijou and M. Bashiri, “Stochastic models for transfer point location problem,” International Journal of Advanced Manufacturing Technology, vol. 58, no. 1–4, pp. 211–225, 2012.

I. N. Katz and L. Cooper, “Facility location in the presence of forbidden regions, I: formulation and the case of Euclidean distance with one forbidden circle,” European Journal of Operational Research, vol. 6, no. 2, pp. 166–173, 1981.

K. Klamroth and M. M. Wiecek, “A bi-objective median location problem with a line barrier,” Operations Research, vol. 50, no. 4, pp. 670–679, 2002.

Y. P. Aneja and M. Parlar, “Algorithms for weber facility location in the presence of forbidden regions and/or barriers to travel,” Location Science, vol. 3, no. 2, pp. 134–134, 1995.

S. E. Butt and T. M. Cavalier, “An efficient algorithm for facility location in the presence of forbidden regions,” European Journal of Operational Research, vol. 90, no. 1, pp. 56–70, 1996.

R. G. McGarvey and T. M. Cavalier, “A global optimal approach to facility location in the presence of forbidden regions,” Computers and Industrial Engineering, vol. 45, no. 1, pp. 1–15, 2003.

K. Klamroth, “Algebraic properties of location problems with one circular barrier,” European Journal of Operational Research, vol. 154, no. 1, pp. 20–35, 2004.

M. Bischoff and K. Klamroth, “An efficient solution method for weber problems with barriers based on genetic algorithms,” European Journal of Operational Research, vol. 177, no. 1, pp. 22–41, 2007.

M. Bischoff, T. Fleischmann, and K. Klamroth, “The multi-facility location-allocation problem with polyhedral barriers,” Computers & Operations Research, vol. 36, no. 5, pp. 1376–1392, 2009.

M. S. Canbolat and G. O. Wesolowsky, “A planar single facility location and border crossing problem,” Computers & Operations Research, vol. 39, no. 12, pp. 3156–3165, 2012.

S. Huang, R. Batta, K. Klamroth, and R. Nagi, “The K-connection location problem in a plane,” Annals of Operations Research, vol. 136, no. 1, pp. 193–209, 2005.

P. Nandikonda, R. Batta, and R. Nagi, “Locating a 1-center on a Manhattan plane with “arbitrarily” shaped barriers,” Annals of Operations Research, vol. 123, no. 1–4, pp. 157–172, 2003.

P. M. Dearing, K. Klamroth, and R. Segars Jr., “Planar location problems with block distance and barriers,” Annals of Operations Research, vol. 136, no. 1, pp. 117–143, 2005.

A. Sarkar, R. Batta, and R. Nagi, “Placing a finite size facility with a center objective on a rectangular plane with barriers,” European Journal of Operational Research, vol. 179, no. 3, pp. 1160–1176, 2007.

L. Friebß, K. Klamroth, and M. Sprau, “A wavefront approach to center location problems with barriers,” Annals of Operations Research, vol. 136, no. 1, pp. 35–48, 2005.

M. S. Canbolat and G. O. Wesolowsky, “The rectilinear distance weber problem in the presence of a probabilistic line barrier,” European Journal of Operational Research, vol. 202, no. 1, pp. 114–121, 2010.

M. S. Canbolat and G. O. Wesolowsky, “On the use of the Varignon frame for single facility weber problems in the presence of convex barriers,” European Journal of Operational Research, vol. 217, no. 2, pp. 241–247, 2012.

A. R. Gharagozlou, A. R. A. Ardalan, and M. H. Diarjan, “Environmental planning for disaster management by using GIS, (a case study about flood in Mazandaran),” The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, vol. 37, B8, pp. 323–326, 2008.

N. Sahinidis and M. Tawarmalani, “Baron,” in GAMS—The Solver Manuals, pp. 19–32, GAMS Development Corporation, Washington, DC, USA, 2008.