Probabilistic Distance, Capacity Clustering Location Model of a Semi-Obnoxious Facility, a Real Case of Tafo, Kumasi, Ghana

Affiliation(s)

^{1}
Department of Mathematics & Statistics, University of Energy & Natural Resources, Sunyani, Ghana.

^{2}
Department of Mathematics & Statistics, Kwame Nkrumah University of Science & Technology, Kumasi, Ghana.

^{3}
Department of Mathematics & Statistics, Catholic University College of Ghana, Sunyani, Ghana.

ABSTRACT

The collection of solid waste from third class communities in most devel-oping countries is by skip containers, however, the location of these facilities has been done arbitrary without any mathematical considerations as to the number of customers the facility is serving, the distance one has to travel to access it and thereby making some of these residences to dump their refuse in gutters, streams and even burn them. In this paper we proposed an improved probabilistic distance, capacity clustering location model which takes into consideration the weight of solid waste from a customer and the capacity of the skip container to locate the skip container to serve a required number of customers based on the capacity constraint of the container. The model was applied on a real world situation and compared with the existing practice in terms of average distance customers had to travel to access the facility. Our results gave a well shorter average travel distance by customers, gave a number of skip containers needed in an area based on their waste generation per capita.

The collection of solid waste from third class communities in most devel-oping countries is by skip containers, however, the location of these facilities has been done arbitrary without any mathematical considerations as to the number of customers the facility is serving, the distance one has to travel to access it and thereby making some of these residences to dump their refuse in gutters, streams and even burn them. In this paper we proposed an improved probabilistic distance, capacity clustering location model which takes into consideration the weight of solid waste from a customer and the capacity of the skip container to locate the skip container to serve a required number of customers based on the capacity constraint of the container. The model was applied on a real world situation and compared with the existing practice in terms of average distance customers had to travel to access the facility. Our results gave a well shorter average travel distance by customers, gave a number of skip containers needed in an area based on their waste generation per capita.

1. Introduction

A Facility Location Problem (FLP) is “defined as in the positioning of a set of points (facilities in our case) within a given location space on the basis of the distribution of demand points (users) to be allocated to the facilities” [1] . Practically applied either in private or in public sector, these problems deal with strategic and long-term decisions involving huge investment costs. In general, once a facility is positioned, it produces effect, positive or negative, on the users (actual or potential), whose intensity is often thought-about depending on the mutual distance. Of course if the effects are positive (desirable facilities), effective positions for the facilities are expected as close as possible to the demand points. It is the case of public utility sites such as schools, hospitals, shops, banks, metro stations and so on. On the contrary, in case of negative effects (undesirable or obnoxious facilities) users wish that facilities are as far as possible. Examples of this kind concern power or nuclear plants, rubbish dumps. However, also in these cases, when facilities are too far from the demand points, additional costs have to be paid in terms of logistic costs. For this reason, it is necessary to find compromised solutions able to balance the different aspects [1] . Obnoxious facilities are facilities in which users no longer consider the facility desirable and try to have it close as possible but avoid the facility and stay away from it. In this paper, we would be looking at waste bins and their locations. Location of a semi-obnoxious facility such as solid waste containers (skip containers) has received scant attention especially in developing countries where customers would have to carry their waste to dump into these containers. Arbitrary location of such facilities in the communities has come with its own problems to society, some of which are the long distances, one has to cover to access the facility forcing people to dispose of solid waste into drains in their neighborhood, containers getting full and overflow in no time and thereby compelling people to throw their waste around the bin thereby creating unhealthy feeding grounds for domestic animals and possibly an outbreak of a disease (like reported here [2] ) and unfriendly nature of the skip containers making it too difficult for children to dump solid waste into the bin. Each one of these problems increases the volume of filth in our communities, towns and cities rather than decreasing filth. It is rather therefore only reasonable that waste bins and skip containers can no longer be sited in arbitrary sites but must have behind them a concrete suggestion from a cluster model whose formulations and solution algorithms address the issue varying widely in terms of fundamental assumptions, mathematical complexity and computational performance [3] . In this paper we proposed an improved probabilistic distance, capacity clustering location model which takes into consideration the weight of solid waste from a customer and the capacity of the skip container to locate the skip container to serve a required number of customers based on the capacity constraint of the container taking Kumasi city as the case study. The results gave a shorter distance average travel distance for users or customers and also proposed the number of skip containers needed in the various zones considering the amount of waste generated in these zones.

2. Study Area

The research was conducted in Tafo, one of the nine sub-metropolitan assemblies in Kumasi, the second populated city in Ghana, West Africa. The nine sub-metropolitan areas are of Kumasi are: Suame, Manhyia, Bantama, Kwadaso, Nkyiaeso, Asokwa, Oforikrom, Subin and Tafo. The study area is the smallest of the nine sub metropolitan areas in terms of land area but it is the second highest generator of solid waste after Subin. The study area has a population of 157,226 with eleven communities within its domain. Four of these communities are categorized under class two (areas where we have fairly good road network but difficult for compactor trucks to collect waste) whiles the remaining seven are under class three (areas where we have bad road network and therefore impossible for proper collection of waste by heavy vehicles). The study area shares boundary with Manhyia to the east, Suame to the west and Subin to the south. The area generates about eighty-eight tones of solid waste a day. Our study considered five of the seven third class zones namely Old Tafo, Pankrono Dome, Pankrono West, Tafo Adompom and Ahenbronum constituting about 52% of the area population. The area has the second largest market in the Kumasi metropolis, has a number of government assisted schools from basic to secondary levels and one government hospital. Figure 1 shows the map of Kumasi Metropolis where the study has been done.

3. Literature Review

Location of facilities be it semi or obnoxious has had its fair share in terms of research in literature, that notwithstanding there are still many areas that needs attention to improve the difficulties service customers go through to access a facility. This section takes a look into various research works that has gone into

Figure 1. Map of Kumasi Metropolis.

facility location of obnoxious and semi-obnoxious alike. A semi-obnoxious facility is a facility which is useful to the end users but has some level of unwelcome effect that produces environmental concerns to the same end user. [4] defined an obnoxious facility as one that generates a disservice to the people nearby while producing an intended product or service. [5] first introduced the rectilinear maximin problem for locating an obnoxious facility. They developed a solution procedure by dividing the feasible region into rectilinear sub-regions and solved using a linear programming problem for each of these sub regions. [6] presented a multi-objective mixed integer linear model for undesirable facility location. The objectives considered are cost minimization and equity maximization. [7] [8] presented a multi-objective mixed integer programming model for the location of hazardous material facilities including the choice of variable technology with cost, risk and equity as the main objective functions. [9] proposed a dynamic multi-period multi-objective capacitated mixed integer programming model for the location of sanitary landfills. [10] studied the location of semi obnoxious facilities as a discrete location problem on a network. Several bi-criteria models are presented considering two conflicting objectives, the minimization of the obnoxious effect and the maximization of the accessibility of the communities to the closest open facility. These objectives tried to optimize the average value over the entire communities’ worst value is considered in two different ways, trying to optimize its average value over all the communities or trying to optimize its worst value. They used Euclidean distance to evaluate the obnoxious effect and the shortest distance path to evaluate the accessibility.

[11] presented a model for the location of semi obnoxious facilities and simultaneously routing the undesirable materials between facilities and communities alike. [12] developed two bi-criteria models for single allocation hub location problems. In both models they used total time as the first criteria to be minimized and used time minimization to process flow entering the hubs. They considered, total time as the second criteria in the first model and, in the second model, the maximum service time for the hubs is minimized. [13] studied a simplified model where clients will be assigned to (unreliable) facilities and reliable backup facilities if needed. [14] presented both scenario-based stochastic programming and a nonlinear mixed integer programming model and show that they are generally equivalent. Also, a constant-ratio approximation algorithm for the case where all failure rates are identical is proposed. [15] considered problems with a fortification budget in the first stage so that unreliable facilities can be fortified by hardening operations. Recently, this line of research is extended to investigate more general reliable network design problems. [16] considered a reliable multiple-echelon logistics network design problem where disruptions can happen in multiple echelons. [17] studied reliable hub-and-spoke network design problems in which hubs could be disrupted and affected flows will be rerouted through survived operational hubs. In their paper, [18] considered a multi-facility location problem in the presence of a line barrier with the starting point of the barrier uniformly distributed. Their objective was to locate n new facilities among m already existing facilities minimizing the summation of the weighted expected rectilinear barrier distances of the locations of new facilities and new and existing facilities. In their scholarly presentation, he sought to solve the problem of location where there are restrictions so that the restrictions prohibit establishment of new facilities in some specified regions. They designed a mixed-integer nonlinear programming model, which was conveniently transformed into a mixed-integer quadratic programming model. They also proposed two meta-heuristic algorithms, namely the genetic algorithm and the imperialist competitive algorithm for optimization for large problems. In their paper [9] , an integer programming model was proposed to help decision makers in choosing the sites to locate the unsorted waste collection bins in a residential town, as well as the capacities of the bins to be located at each collection site. Their model helped in assessing tactical decisions through constraints that force each collection area to be capacitated enough to fit the expected waste to be directed to that area, while taking into account Quality of Service constraints from the citizens’ point of view. Additionally, they proposed an effective constructive heuristic approach whose aim is to provide a good solution quality in an extremely reduced computational time. Computational results on data related to the city of Nardò, in the south of Italy, show that both exact and heuristic approaches provide consistently better solutions than that currently implemented, resulting in a lower number of activated collection sites, and a lower number of bins to be used [9] .

3.1. Maximum Insertion Probabilistic Model

A location problem is to locate a facility, or facilities, to serve optimally a given set of customers. The customers are given by their coordinates and demands. The coordinates are points in R^{n} (usually n = 2), and the demands are positive numbers q_{i}. Assuming N customers, the data of the problem is a set of points (coordinates)
$X=\left\{{X}_{1},{X}_{2},\cdots ,{X}_{N}\right\}$ in
${R}^{n}$ and a corresponding set of positive weights (demands)
$\left\{{q}_{1},{q}_{2},\cdots ,{q}_{N}\right\}$. Using the Euclidean distance
$d\left(X,Y\right)=\Vert X-Y\Vert $ between two points X, Y in R^{n}. If the customers are served by K facilities for a given K, we denote
${X}_{K}$ to be the set of customers assigned to the K^{th}-facility. The weighted sum of distance travelled by these customers is
$\underset{{X}_{i}\in {\chi}_{K}}{\sum}\frac{1}{{q}_{i}}\Vert {X}_{i}-{c}_{K}\Vert$ where c_{K} is the location of the K^{th}-facility. The customers
$X=\left\{{X}_{1},{X}_{2},\cdots ,{X}_{N}\right\}$, their demands
$\left\{{q}_{1},{q}_{2},\cdots ,{q}_{N}\right\}$ is to be located to facilities with centres
$\left\{{c}_{1},{c}_{2},\cdots ,{c}_{K}\right\}$ so as to minimize the weighted sum of distances travelled by all the customers to access the bin.
$\underset{{c}_{1},\cdots ,{c}_{K}}{\mathrm{min}}\underset{{X}_{1},\cdots ,{X}_{K}}{\mathrm{min}}{\displaystyle \underset{k=1}{\overset{K}{\sum}}{\displaystyle \underset{{X}_{1}\in {\chi}_{K}}{\sum}\frac{1}{{q}_{i}}\Vert {X}_{i}-{c}_{K}\Vert}}$ and each container has a capacity Q, such that the sum of demands assigned to the K^{th}—facility cannot exceed it, i.e.

$\underset{{X}_{i}\in {\chi}_{K}}{\sum}{q}_{i}}\le {Q}_{K$ (3.1)

Let a data set
$D\subset {R}^{p}$ be partitioned into K clusters
$\left\{{c}_{k}:k=1,2,\cdots ,K\right\}$,
$D={\displaystyle \underset{k=1}{\overset{K}{\cup}}{C}_{k}}$ and let C_{k} be the centre of the cluster C_{k} with capacity Q_{K}.

With each data point $X\in D$ and a cluster centre ${C}_{k}$, we denote:

· a distance ${d}_{k}\left(X,{c}_{k}\right)$ by ${d}_{k}(X)$

· a probability of membership in
${c}_{k}$ by
${P}_{k}\left(X\right)$ for each
$X\in D$ and each cluster C_{k},

$\frac{{P}_{k}\left(X\right){d}_{k}\left(X\right)}{{Q}_{k}}=\text{Constant,}$ (3.2)

depending on X and independent of the cluster k.

Given the cluster centers $\left\{{C}_{1},{C}_{2},\cdots ,{C}_{k}\right\}$, let X be a data point and let $\left\{{d}_{k}\left(X\right):k=1,2,\cdots ,k\right\}$ be its distances from the given centers.

From Equation (3.2) using i and k ${P}_{i}\left(X\right)\frac{{d}_{i}\left(X\right)}{{Q}_{i}}={P}_{k}\left(X\right)\frac{{d}_{k}\left(X\right)}{{Q}_{k}}$

$\text{Since}{\displaystyle \underset{i=1}{\overset{K}{\sum}}{P}_{i}}\left(X\right)=1,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{P}_{k}\left(X\right){\displaystyle \underset{i=1}{\overset{K}{\sum}}\left(\frac{\frac{{d}_{k}\left(X\right)}{{Q}_{k}}}{\frac{{d}_{i}\left(X\right)}{{Q}_{i}}}\right)}=1$

${P}_{k}\left(X\right)=\frac{{\displaystyle \underset{j\ne k}{\prod}\frac{{d}_{j}\left(X\right)}{{Q}_{j}}}}{{\displaystyle \underset{i=1}{\overset{K}{\sum}}{\displaystyle \underset{j\ne i}{\prod}\frac{{d}_{j}\left(X\right)}{{Q}_{j}}}}},\text{\hspace{0.17em}}k=1,2,\cdots ,K$

For two clusters, i.e. K = 2

${P}_{1}=\frac{\frac{{d}_{2}\left(X\right)}{{Q}_{2}}}{\frac{{d}_{1}\left(X\right)}{{Q}_{1}}+\frac{{d}_{2}\left(X\right)}{{Q}_{2}}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{P}_{2}=\frac{\frac{{d}_{1}\left(X\right)}{{Q}_{1}}}{\frac{{d}_{1}\left(X\right)}{{Q}_{1}}+\frac{{d}_{2}\left(X\right)}{{Q}_{2}}}$ (3.3)

For three clusters i.e. K = 3

${P}_{1}=\frac{\frac{{d}_{2}\left(X\right){d}_{3}\left(X\right)}{{Q}_{2}\times {Q}_{3}}}{\frac{{d}_{1}\left(X\right){d}_{2}\left(X\right)}{{Q}_{1}\times {Q}_{2}}+\frac{{d}_{1}\left(X\right){d}_{3}\left(X\right)}{{Q}_{1}\times {Q}_{3}}+\frac{{d}_{2}\left(X\right){d}_{3}\left(X\right)}{{Q}_{2}\times {Q}_{3}}}$ (3.4)

Denoting the constant term in Equation (3.2) by D(X), function in X, $\frac{{P}_{k}\left(X\right){d}_{k}\left(X\right)}{{Q}_{k}}=D(X)$

$\begin{array}{l}{P}_{k}\left(X\right)=\frac{D\left(X\right)}{\frac{{d}_{k}\left(X\right)}{{Q}_{k}}}\text{,}\text{\hspace{0.17em}}k=1,2,3,\cdots ,K\\ \text{Since}{\displaystyle \underset{i=1}{\overset{K}{\sum}}{P}_{k}}\left(X\right)=1\end{array}$

$D\left(X\right)=\frac{{\displaystyle \underset{j=1}{\overset{K}{\prod}}\frac{{d}_{j}\left(X\right)}{{Q}_{j}}}}{{\displaystyle \underset{i=1}{\overset{K}{\sum}}{\displaystyle \underset{j\ne i}{\prod}\frac{{d}_{j}\left(X\right)}{{Q}_{j}}}}}$ (3.5)

For two clusters, i.e. K = 2,

$D\left(X\right)=\frac{\frac{{d}_{1}\left(X\right){d}_{2}\left(X\right)}{{Q}_{1}{Q}_{2}}}{\frac{{d}_{1}\left(X\right)}{{Q}_{1}}+\frac{{d}_{2}\left(X\right)}{{Q}_{2}}}$ (3.6)

For three clusters, i.e. K = 3;

$D\left(X\right)=\frac{\frac{{d}_{1}\left(X\right){d}_{2}\left(X\right){d}_{3}\left(X\right)}{{Q}_{1}{Q}_{2}{Q}_{3}}}{\frac{{d}_{1}\left(X\right)}{{Q}_{1}}+\frac{{d}_{2}\left(X\right)}{{Q}_{2}}+\frac{{d}_{3}\left(X\right)}{{Q}_{3}}}$ (3.7)

Therefore for the entire data set D is the sum of (3.6) over all points, and is a function of the K cluster centers,

$F\left({c}_{1},{c}_{2},\cdots ,{c}_{k}\right)={\displaystyle \underset{i=1}{\overset{N}{\sum}}\frac{{\displaystyle \underset{k=1}{\overset{K}{\prod}}{d}_{k}\left({X}_{i},{c}_{k}\right)}}{{\displaystyle \underset{t=1}{\overset{K}{\sum}}{\displaystyle \underset{j\ne t}{\prod}{d}_{j}\left({X}_{i},{c}_{j}\right)}}}}$ (3.8)

3.2. An Extremal Principle

Consider a case of two clusters, using X to be the given data point with distances ${d}_{1}\left(X\right),\text{\hspace{0.17em}}{d}_{2}\left(X\right)$ to the cluster centers and a known cluster sizes ${Q}_{1}$ and ${Q}_{2}$ known. The probabilities in (3) are the optimal solution ${P}_{1},{P}_{2}$ of the extremal problem.

$\begin{array}{l}\text{Minimize}\left\{\frac{{d}_{1}\left(x\right){P}_{1}^{2}}{{Q}_{1}}+\frac{{d}_{2}\left(x\right){P}_{2}^{2}}{{Q}_{2}}\right\}\\ \text{Subjectto}\text{\hspace{0.17em}}{P}_{1}+{P}_{2}=1\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{P}_{1},{P}_{2}\ge 0\end{array}$ (3.9)

The Lagrangian of the problem is:

$L\left({P}_{1},{P}_{2},\lambda \right)=\frac{{d}_{1}\left(x\right){P}_{1}^{2}}{{Q}_{1}}+\frac{{d}_{2}\left(x\right){P}_{2}^{2}}{{Q}_{2}}-\lambda \left({P}_{2}+{P}_{1}-1\right)$ (3.10)

Setting the partial derivatives with respect to ${P}_{1},{P}_{2}$ to zero, we have ${P}_{1}{d}_{1}\left(x\right)={P}_{2}{d}_{2}\left(x\right)$.

Substitute the probabilities in (3.4) in (3.10) we have ${L}^{\ast}\left({P}_{1}\left(x\right),{P}_{2}\left(x\right),\lambda \right)=\frac{\frac{{d}_{1}\left(x\right){d}_{2}\left(x\right)}{{Q}_{1}{Q}_{2}}}{\frac{{d}_{1}\left(x\right)}{{Q}_{1}}+\frac{{d}_{2}\left(x\right)}{{Q}_{2}}}$ which is the same as the joint distance function obtained in Equation (3.9).

The extremal problem for the entire data set $D=\left\{{X}_{1},{X}_{2},\cdots ,{X}_{N}\right\}\subset {R}^{n}$

Minimize

$\underset{i=1}{\overset{N}{\sum}}\left(\frac{{d}_{1}\left({x}_{i}\right){p}_{1}{\left({x}_{i}\right)}^{2}}{{Q}_{1}}+\frac{{d}_{2}\left({x}_{i}\right){p}_{2}{\left({x}_{i}\right)}^{2}}{{Q}_{2}}\right)$ (3.11)

$\begin{array}{l}\text{Subjectto:}{p}_{1}\left({x}_{i}\right)+{p}_{2}\left({x}_{i}\right)=1\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{p}_{1}\left({x}_{i}\right),\text{\hspace{0.17em}}{p}_{2}\left({x}_{i}\right)\ge 0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}i=1,2,3,\cdots ,N\end{array}$

where ${p}_{1}\left({x}_{i}\right),{p}_{2}\left({x}_{i}\right)$ are the cluster probabilities at ${x}_{i}$ and ${d}_{1}\left({x}_{i}\right),{d}_{2}\left({x}_{i}\right)$ are the corresponding distances. The problem separates into N as in (3.11) and its optimal value is $\underset{i=1}{\overset{N}{\sum}}\frac{\frac{{d}_{1}\left({x}_{i}\right){d}_{2}\left({x}_{i}\right)}{{Q}_{1}{Q}_{2}}}{\frac{{d}_{1}\left({x}_{i}\right)}{{Q}_{1}}+\frac{{d}_{2}\left({x}_{i}\right)}{{Q}_{2}}}$ the sum of the joint distance function of all points.

3.3. Cluster Centers

Writing Equation (3.11) as a function of the cluster centers ${C}_{1},{C}_{2}$

$f\left({C}_{1},{C}_{2}\right)={\displaystyle \underset{i=1}{\overset{N}{\sum}}\left(\frac{{d}_{1}\left({x}_{i},{C}_{1}\right){p}_{1}{\left({x}_{i}\right)}^{2}}{{Q}_{1}}+\frac{{d}_{2}\left({x}_{i},{C}_{2}\right){p}_{2}{\left({x}_{i}\right)}^{2}}{{Q}_{2}}\right)}$.

For Euclidean distance ${d}_{k}\left(x,{c}_{k}\right)=\Vert x-{c}_{k}\Vert ,k=1,2$ so that

$f\left({C}_{1},{C}_{2}\right)={\displaystyle \underset{i=1}{\overset{N}{\sum}}\left(\frac{\Vert {x}_{i}-{C}_{1}\Vert {p}_{1}{\left({x}_{i}\right)}^{2}}{{Q}_{1}}+\frac{\Vert {x}_{i}-{C}_{2}\Vert {p}_{2}{\left({x}_{i}\right)}^{2}}{{Q}_{2}}\right)}$ (3.12)

and look for centers ${C}_{1},{C}_{2}$ that minimizes f and the probabilities ${p}_{1}\left({x}_{i}\right),{p}_{2}\left({x}_{i}\right)$ are given for $i=1,2,3,\cdots ,N$ and with an assumption that the cluster centers ${C}_{1},{C}_{2}$ do not coincide with any of the points $i=1,2,3,\cdots ,N$

From the assumption above, the gradient of Equation (3.12) with respect to ${c}_{k}$ is

$\begin{array}{c}{\nabla}_{c}f\left({c}_{1},{c}_{2}\right)=-{\displaystyle \underset{i=1}{\overset{N}{\sum}}\frac{{x}_{i}-{c}_{k}}{\Vert {x}_{i}-{c}_{k}\Vert}}\times {p}_{k}{\left({x}_{i}\right)}^{2}\\ =-{\displaystyle \underset{i=1}{\overset{N}{\sum}}\frac{{x}_{i}-{c}_{k}}{{d}_{k}\left({x}_{i},{c}_{k}\right)}}\times {p}_{k}{\left({x}_{i}\right)}^{2},\text{\hspace{0.17em}}k=1,2\end{array}$

Setting the gradient to zero, and grouping like terms, we have

$\underset{i=1}{\overset{N}{\sum}}\left(\frac{{p}_{k}{\left({x}_{i}\right)}^{2}}{{d}_{k}\left({x}_{i},{c}_{k}\right)}\right){x}_{i}}=\left({\displaystyle \underset{i=1}{\overset{N}{\sum}}\frac{{p}_{k}{\left({x}_{i}\right)}^{2}}{{d}_{k}\left({x}_{i},{c}_{k}\right)}}\right){c}_{k$

${C}_{k}={\displaystyle \underset{i=1}{\overset{N}{\sum}}\left(\frac{{u}_{k}\left({x}_{i}\right)}{{\displaystyle \underset{j=1}{\overset{N}{\sum}}{u}_{k}\left({x}_{j}\right)}}\right){x}_{i}}$ (3.13)

where ${u}_{k}=\frac{{p}_{k}{\left({x}_{i}\right)}^{2}}{{d}_{k}\left({x}_{i},{c}_{k}\right)}$ for $k=1,2$

Giving the minimizers as

${C}_{1}={\displaystyle \underset{i=1}{\overset{N}{\sum}}\left(\frac{{u}_{1}\left({x}_{i}\right)}{{\displaystyle \underset{j=1}{\overset{N}{\sum}}{u}_{1}\left({x}_{j}\right)}}\right){x}_{i}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{C}_{2}={\displaystyle \underset{i=1}{\overset{N}{\sum}}\left(\frac{{u}_{2}\left({x}_{i}\right)}{{\displaystyle \underset{j=1}{\overset{N}{\sum}}{u}_{2}\left({x}_{j}\right)}}\right){x}_{i}}$ (3.14)

where

${u}_{1}=\frac{{p}_{1}{\left({x}_{i}\right)}^{2}}{{d}_{1}\left({x}_{i},{c}_{1}\right)},\text{\hspace{0.17em}}\text{}{u}_{2}=\frac{{p}_{2}{\left({x}_{i}\right)}^{2}}{{d}_{2}\left({x}_{i},{c}_{2}\right)}$ (3.15)

For a function of K cluster centers $f\left({C}_{1},{C}_{2},\cdots ,{C}_{K}\right)={\displaystyle \underset{k=1}{\overset{K}{\sum}}{\displaystyle \underset{i=1}{\overset{N}{\sum}}\left(\frac{{d}_{k}\left({x}_{i},{C}_{k}\right){p}_{k}{\left({x}_{i}\right)}^{2}}{{Q}_{k}}\right)}}$ an analog of (3.13) then by the results of Equation (3.2) the minimizers of f is

${C}_{k}={\displaystyle \underset{i=1}{\overset{N}{\sum}}\left(\frac{{u}_{k}\left({x}_{i}\right)}{{\displaystyle \underset{j=1}{\overset{N}{\sum}}{u}_{k}\left({x}_{j}\right)}}\right){x}_{i}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{with}\text{\hspace{0.17em}}{u}_{k}=\frac{{p}_{k}{\left({x}_{i}\right)}^{2}}{{d}_{k}\left({x}_{i},{c}_{k}\right)}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}k=1,2,\cdots ,K$ (3.16)

4. Algorithm for Probabilistic Distance Location Model

The major steps involved in the formation of the algorithm are described below.

1) Calculation of number of skip containers

The number of skip containers (k) needed in each zone is dependent on demands (q_{i}) of the customers and the capacity of the (Q) of the container;
$k=\frac{{\displaystyle \underset{i=1}{\overset{N}{\sum}}{q}_{i}}}{Q}.$

2) Selection of initial cluster centers

· The x-coordinates of all the customers are arranged in ascending order. The highest x-coordinates is selected with its demand recorded, the next highest x-coordinates value is recorded with its demand until the capacity constraints of the container is satisfied.

· The centroid of these points is found by taken the average of the extreme x-values and that of the y-values.

· The process is continued until all the initial centres are assigned.

3) Updating of the cluster centers

i) Compute distances from C_{i} for all
$X\in D.$

ii) Update the centre ${{C}^{\prime}}_{i}$ using Equation (16)

iii) If the difference between $\underset{i=1}{\overset{K}{\sum}}\left|{C}_{i}^{n-1}-{C}_{i}^{n-2}\right|$ and $\underset{i=1}{\overset{K}{\sum}}\left|{C}_{i}^{n}-{C}_{i}^{n-1}\right|$ is less than $\epsilon $ stop else return to ii.

4) Data Analysis & Result Using the Proposed Model

The data for the study was obtained jointly from Waste management division of Kumasi Metropolitan Assembly (KMA), Physical Planning unit of KMA, Ghana Statistical service (regional office) and a six (6) month continuous field work. The study area was divided into four zones due to physical boundaries such as streams, huge bridges which cannot be used by vehicles. Zone one 1) has seven hundred and forty-nine (749) houses with and has 1100 bins, zone two 2) with five hundred and sixty-one (561) houses has 828 bins, zone three 3) has five hundred and forty-two (542) houses with 792 bins and zone four had six hundred and twenty-three (623) houses with 789 bins of 140 liters respectively. The model was coded in C++ and ran on Intel core 470 Gb computer.

Cluster Centre from the Zones

A 14 m^{3} skip containers were used for the location centers, with a compaction factor of 0.4, a container can hold 166 of 140 litre bins. With the number of bins in zone one, the number of bins needed is calculated as
$\frac{{\displaystyle \underset{i=1}{\overset{n}{\sum}}{q}_{i}}}{{v}_{j}}=\frac{1100}{166}=6.62\cong 7$.

Coordinates of all the houses obtained by a Geographic Information System were taken and implemented using the probabilistic model. The final cluster (location centers) and their coordinates for each of the seven location centers by the model after 500 iterations is given in the tables below.

The maps mainly indicate the clusters obtained by the maximum insertion with their cluster centers where the waste collected should be sent to. The red color spots among the clusters of points represent the cluster centers for each cluster (which is differentiated by colors).

From Figure 2, there are seven sub-clusters in the seven bounded areas shown in the map below. Each sub-cluster serves users or customers in its bounded area. All the wastes collected in the area are sent to the associated sub-cluster center. The details of Figure 2 are stored in Table 1.

From Figure 3, there are five sub-clusters in for the five bounded areas shown in the map below. Each sub-cluster serves users or customers in its bounded area. All the wastes collected in the area is sent to the associated

Figure 2. Final cluster centers in zone 1.

Table 1. Summary of Zone 1 for Figure 2.

Figure 3. Final cluster centers in zone 2.

sub-cluster center. The detail of Figure 3 is stored in Table 2.

From Figure 4, there are five sub-clusters in for the five bounded areas shown in the map below. Each sub-cluster serves users or customers in its bounded area. All the wastes collected in the area are sent to the associated sub-cluster center. The detail of Figure 4 is stored in Table 3.

From Figure 5, there are five sub-clusters in for the five bounded areas shown in the map below. Each sub-cluster serves users or customers in its bounded area. All the wastes collected in the area are sent to the associated sub-cluster center. The detail of Figure 5 is stored in Table 4.

This table contains the summary of sub-cluster centers, number of customers in a sub-cluster and the number of bins used by customers in each sub-cluster in Zone 1. It is associated to Figure 2 and also contains the average distance from travelled from any location in the bounded area to its associated sub-cluster center. The cluster center points recorded in this table represents the geographical locations of each sub-cluster given by the maximum insertion method in and it is labeled with a red point in Figure 2. The cluster centers serve all

Table 2. Summary of Zone 2 for Figure 3.

Figure 4. Final cluster centers in zone 3.

Table 3. Summary of Zone 3 for Figure 4.

users or customers in the bounded area. Thus, wastes collected in each bounded area in Zone 1 are sent to the sub-cluster center. There are seven sub-clusters in Zone 1 and labeled 1 - 7 in the table and in Figure 2.

This table contains the summary of sub-cluster centers, number of customers in a sub-cluster and the number of bins used by customers in each sub-cluster in Zone 2. It is associated to Figure 3 and also contains the average distance from travelled from any location in the bounded area to its associated sub-cluster center. The cluster center points recorded in this table represents the geographical locations of each sub-cluster given by the maximum insertion method in and it is labeled with a red point in Figure 3. The cluster centers serve all

Figure 5. Final cluster centers in zone 4.

Table 4. Summary of Zone 4 for Figure 5.

users or customers in the bounded area. Thus, wastes collected in each bounded area in Zone 2 are sent to the sub-cluster center. There are five sub-clusters in Zone 2 and labeled 1 - 5 in the table and in Figure 3.

This table contains the summary of sub-cluster centers, number of customers in a sub-cluster and the number of bins used by customers in each sub-cluster in Zone 3. It is associated to Figure 4 and also contains the average distance from travelled from any location in the bounded area to its associated sub-cluster center. The cluster center points recorded in this table represents the geographical locations of each sub-cluster given by the maximum insertion method in and it is labeled with a red point in Figure 4. The cluster centers serve all users or customers in the bounded area. Thus, wastes collected in each bounded area in Zone 3 are sent to the sub-cluster center. There are five sub-clusters in Zone 3 and labeled 1 - 5 in the table and in Figure 4.

This table contains the summary of sub-cluster centers, number of customers in a sub-cluster and the number of bins used by customers in each sub-cluster in Zone 4. It is associated to Figure 5 and also contains the average distance from travelled from any location in the bounded area to its associated sub-cluster center. The cluster center points recorded in this table represents the geographical locations of each sub-cluster given by the maximum insertion method in and it is labeled with a red point in Figure 5. The cluster centers serve all users or customers in the bounded area. Thus, wastes collected in each bounded area in Zone 4 are sent to the sub-cluster center. There are five sub-clusters in Zone 4 and labeled 1 - 5 in the table and in Figure 5.

5. Conclusion

The adoption of the improved probabilistic distance location model to the study area of 2475 houses and about 3509 of 140 litre bins has shown that the area needs twenty-two (22) skip containers, seven in zone one (1), five skip containers each in zones two (2), three (3), and four (4) respectively as compared to the current 15 skip containers in the area. The model has also clearly identified the required number of customers to be assigned to each skip container based on the capacity of each customer and the capacity of the skip container. Our results have given an average walking distance of 92.86 m for a customer to access the skip container as compared to the existing average of 234.71 m to access a skip container.

Cite this paper

Otoo, D. , Sebil, C. , Kessie, J. and Larbi, E. (2019) Probabilistic Distance, Capacity Clustering Location Model of a Semi-Obnoxious Facility, a Real Case of Tafo, Kumasi, Ghana.*American Journal of Operations Research*, **9**, 146-160. doi: 10.4236/ajor.2019.93009.

Otoo, D. , Sebil, C. , Kessie, J. and Larbi, E. (2019) Probabilistic Distance, Capacity Clustering Location Model of a Semi-Obnoxious Facility, a Real Case of Tafo, Kumasi, Ghana.

References

[1] Barbati, M. (2013) Models and Algorithms for Facility Location Problems with Equity Considera-tions.

[2] http://78.47.45.183/agency/ultimate-fm/20180215/72633487/filth-in-tafo-ahenbronom-as-residents-fear-outbreak-of-diseases

[3] Shiripour, S., Mahdavi, I., Amiri-Aref, M., Mohammadnia-Otaghsara, M. and Mahdavi-Amiri, N. (2012) Multi-Facility Location Problems in the Presence of a Probabilistic Line Barrier: A Mixed Integer Quadratic Programming Model. International Journal of Production Research, 50, 3988-4008.

https://doi.org/10.1080/00207543.2011.579639

[4] Erkut, E. and Neuman, S. (1989) Analytical Models for Locating Undesirable Facilities. European Journal of Operational Research, 40, 275-291.

https://doi.org/10.1016/0377-2217(89)90420-7

[5] Erkut, E. and Neuman, S. (1992) A Multiobjective Model for Locating Undesirable Facilities. Annals of Operations Research, 40, 209-227.

https://doi.org/10.1007/BF02060478

[6] Drezner and Wesolowsky, G.O. (1983) The Location of an Obnoxious Facility with Rectangular Distance. Journal of Regional Science, 23, 241-248.

https://doi.org/10.1111/j.1467-9787.1983.tb00800.x

[7] Wyman, M. and Kuby, M. (1993) A Multiobjective Loca-tion-Allocation Model for Assessing Toxic Waste Processing Technologies. Studies in Locational Analysis, 4, 193-196.

[8] Wyman, M. and Kuby, M. (1995) Proactive Optimization of Toxic Waste Transportation, Location and Tech-nology. Location Science, 3, 167-185.

https://doi.org/10.1016/0966-8349(95)00014-3

[9] Ghiani, G., Laganà, D., Manni, E. and Triki, C. (2012) Capacitated Location of Collection Sites in an Urban Waste Management System. Waste Management, 32, 1291-1296.

https://doi.org/10.1016/j.wasman.2012.02.009

[10] Costa, M.G., Captivo, M.E. and Clímaco, J. (2008) Capacitated Single Allocation Hub Location Problem—A Bicriteria Approach. Computers & Operations Research, 35, 3671-3695.

https://doi.org/10.1016/j.cor.2007.04.005

[11] Cappanera, P., Gallo, G. and Maffioli, F. (2004) Discrete Facility Location and Routing of Obnoxious Facilities. Discrete Applied Mathematics, 133, 3-28.

https://doi.org/10.1016/S0166-218X(03)00431-1

[12] Fonseca, M.C. and Captivo, M.E. (2007) Analysis of Bicriteria Semi Obnoxious Location Models. Paper Presented at Tercer Encuentro de la Red Iberoamericana de Evaluación Decisión Multicriterio, R.E.D.-M. 2007, Culiacán.

[13] Shen, Z., Zhan, R. and Zhang, J. (2011) The Reliable Facility Location Problem: Formulations, Heuristics, and Approximation Algorithms. INFORMS Journal on Computing, 23, 470-482.

https://doi.org/10.1287/ijoc.1100.0414

[14] Li, X., Tomasgard, A. and Barton, P.I. (2011) Nonconvex Generalized Benders Decomposition for Stochastic Separable Mixed-Integer Nonlinear Programs. Journal of Optimization Theory and Applications, 151, 425.

https://doi.org/10.1007/s10957-011-9888-1

[15] Li, Q., Zeng, B. and Savachkin, A. (2012) Reliable Facility Location Design under Disruptions. Computers & Operations Research, 40, 901-909.

https://doi.org/10.1016/j.cor.2012.11.012

[16] An, Y., Zeng, B., Zhang, Y. and Zhao, L. (2014) Reliable p-Median Facility Location Problem: Two-Stage Robust Models and Algorithms. Transportation Research Part B: Methodological, 64, 54-72.

https://doi.org/10.1016/j.trb.2014.02.005

[17] An, Y., Zhang, Y. and Zeng, B. (2015) The Reliable Hub-and-Spoke Design Problem: Models and Algorithms. Transportation Research Part B: Methodological, 77, 103-122.

https://doi.org/10.1016/j.trb.2015.02.006

[18] Farahani, R.Z., SteadieSeifi, M. and Asgari, N. (2010) Multiple Criteria Facility Location Problems: A Survey. Applied Mathematical Modelling, 34, 1689-1709.

https://doi.org/10.1016/j.apm.2009.10.005

[1] Barbati, M. (2013) Models and Algorithms for Facility Location Problems with Equity Considera-tions.

[2] http://78.47.45.183/agency/ultimate-fm/20180215/72633487/filth-in-tafo-ahenbronom-as-residents-fear-outbreak-of-diseases

[3] Shiripour, S., Mahdavi, I., Amiri-Aref, M., Mohammadnia-Otaghsara, M. and Mahdavi-Amiri, N. (2012) Multi-Facility Location Problems in the Presence of a Probabilistic Line Barrier: A Mixed Integer Quadratic Programming Model. International Journal of Production Research, 50, 3988-4008.

https://doi.org/10.1080/00207543.2011.579639

[4] Erkut, E. and Neuman, S. (1989) Analytical Models for Locating Undesirable Facilities. European Journal of Operational Research, 40, 275-291.

https://doi.org/10.1016/0377-2217(89)90420-7

[5] Erkut, E. and Neuman, S. (1992) A Multiobjective Model for Locating Undesirable Facilities. Annals of Operations Research, 40, 209-227.

https://doi.org/10.1007/BF02060478

[6] Drezner and Wesolowsky, G.O. (1983) The Location of an Obnoxious Facility with Rectangular Distance. Journal of Regional Science, 23, 241-248.

https://doi.org/10.1111/j.1467-9787.1983.tb00800.x

[7] Wyman, M. and Kuby, M. (1993) A Multiobjective Loca-tion-Allocation Model for Assessing Toxic Waste Processing Technologies. Studies in Locational Analysis, 4, 193-196.

[8] Wyman, M. and Kuby, M. (1995) Proactive Optimization of Toxic Waste Transportation, Location and Tech-nology. Location Science, 3, 167-185.

https://doi.org/10.1016/0966-8349(95)00014-3

[9] Ghiani, G., Laganà, D., Manni, E. and Triki, C. (2012) Capacitated Location of Collection Sites in an Urban Waste Management System. Waste Management, 32, 1291-1296.

https://doi.org/10.1016/j.wasman.2012.02.009

[10] Costa, M.G., Captivo, M.E. and Clímaco, J. (2008) Capacitated Single Allocation Hub Location Problem—A Bicriteria Approach. Computers & Operations Research, 35, 3671-3695.

https://doi.org/10.1016/j.cor.2007.04.005

[11] Cappanera, P., Gallo, G. and Maffioli, F. (2004) Discrete Facility Location and Routing of Obnoxious Facilities. Discrete Applied Mathematics, 133, 3-28.

https://doi.org/10.1016/S0166-218X(03)00431-1

[12] Fonseca, M.C. and Captivo, M.E. (2007) Analysis of Bicriteria Semi Obnoxious Location Models. Paper Presented at Tercer Encuentro de la Red Iberoamericana de Evaluación Decisión Multicriterio, R.E.D.-M. 2007, Culiacán.

[13] Shen, Z., Zhan, R. and Zhang, J. (2011) The Reliable Facility Location Problem: Formulations, Heuristics, and Approximation Algorithms. INFORMS Journal on Computing, 23, 470-482.

https://doi.org/10.1287/ijoc.1100.0414

[14] Li, X., Tomasgard, A. and Barton, P.I. (2011) Nonconvex Generalized Benders Decomposition for Stochastic Separable Mixed-Integer Nonlinear Programs. Journal of Optimization Theory and Applications, 151, 425.

https://doi.org/10.1007/s10957-011-9888-1

[15] Li, Q., Zeng, B. and Savachkin, A. (2012) Reliable Facility Location Design under Disruptions. Computers & Operations Research, 40, 901-909.

https://doi.org/10.1016/j.cor.2012.11.012

[16] An, Y., Zeng, B., Zhang, Y. and Zhao, L. (2014) Reliable p-Median Facility Location Problem: Two-Stage Robust Models and Algorithms. Transportation Research Part B: Methodological, 64, 54-72.

https://doi.org/10.1016/j.trb.2014.02.005

[17] An, Y., Zhang, Y. and Zeng, B. (2015) The Reliable Hub-and-Spoke Design Problem: Models and Algorithms. Transportation Research Part B: Methodological, 77, 103-122.

https://doi.org/10.1016/j.trb.2015.02.006

[18] Farahani, R.Z., SteadieSeifi, M. and Asgari, N. (2010) Multiple Criteria Facility Location Problems: A Survey. Applied Mathematical Modelling, 34, 1689-1709.

https://doi.org/10.1016/j.apm.2009.10.005