Mining Method Selection by Integrated Ahp and Promethee Method

Selecting the best mining method among many alternatives is a multicriteria decision making problem. The aim of this paper is to demonstrate the implementation of an integrated approach that employs AHP and PROMETHEE together for selecting the most suitable mining method for the " Coka Marin " underground mine in Serbia. The related problem includes five possible mining methods and eleven criteria to evaluate them. Criteria are accurately chosen in order to cover the most important parameters that impact on the mining method selection, such as geological and geotechnical properties, economic parameters and geographical factors. The AHP is used to analyze the structure of the mining method selection problem and to determine weights of the criteria, and PROMETHEE method is used to obtain the final ranking and to make a sensitivity analysis by changing the weights. The results have shown that the proposed integrated method can be successfully used in solving mining engineering problems.


INTRODUCTION
defines underground mining as a technique of recovering minerals from deposits below the earth's surface.
Important indicators of underground mining, such as net efficiency production, excavation cost, ore losses, ore dilution and the final financial effects depend on selected and applied mining method.Certainly the most important goal of the applied mining method is the achievement of a lower excavation cost and consequently a greater financial profit.However, the mining method selection cannot be solely based on these criteria.
Safe working conditions, ore excavation efficien cy and ore dilution are also important characteristics of a selected mining method, which significantly affect the financial effects of the applied mining method.
Numerous considerations indicate that the undergorund mining method selection depends on the large number of relevant factors.These factors can be classified into three main groups, i.e.: mining-geological factors, such as: ground conditions, hangingwall and footwall, ore thickness, general shape, dip, plunge, depth below the surface, grade distribution, quality of resource, etc. mining-technical factors, such as: annual productivity, applied equipment, environmental considerations, mine recovery, flexibility of methods, machinery and mining rate, and economic factors, such as: capital cost, operating cost, mineable ore tons, orebody grades and mineral value.DEJAN BOGDANOVIC, DJORDJE NIKOLIC and IVANA ILIC In practice, there are cases in which the mininggeological factors allow the application of a certain mining method, but its application is not justified from the aspect of financial effects.There are also cases in which a certain mining method considers the application of certain types of machinery, but it is not justified from the aspect of mining-technical factors.
Anyway, mining method selection is a timeconsuming and difficult process, requiring great advanced knowledge and experience.The process can be a hard task for mining engineers and managers.For a proper and effective evaluation, the decision maker may need to analyze a large amount of data and to consider many factors.There are several methodologies that were developed in the past for mining method selection.Bajkonurov (1969), Budko (1971), Imenitov (1970) and Popov (1970) have suggested a procedure for mining method selection that comprises two phases (preliminary mining method selection and the selection of the most suitable mining method from a group of applicable methods).The selection of the most suitable method from a group of applicable methods is made by the techno-economic model procedure.This model is based on the estimation of expected financial effects, which can be aquired by the implementation of each method from a group of applicable methods.The mining method with the best financial effects is selected.
Later, Boshkov andWright (1973), Morrison (1976), Laubscher (1981) and Hartman (1987) have suggested a series of approaches for mining method selection.However, these approaches were inadequate for the development of a methodology that would automatically choose a mining method.
The first numerical approach for mining method selection was suggested by Nicholas (1981Nicholas ( , 1992)).This methodology rates different mining methods based on the rankings of specific input parameters.The mining method with the highest sum result is selected.Nicholas has suggested some modifications that involve the weighting of various categories, such as that of ore geometry, ore zone, hangingwall and footwall.Miller et al. (1995) have developed the UBC methodology.The UBC mining method selection is a modification of the Nicolas methodology that better represents typical Canadian mining design practices.The main weaknesses of these approaches are in the fact that the importance of each selection criteria has not been considered.
An up-to-date approach considers mining method selection as a multi-criteria decisionmaking (MCDM) problem with a finite number of alternatives that have to be ranked considering many different and conflicting criteria.The advantage of these methods is that they can account for both financial and non-financial impacts.Among these methods, the most popular ones are scoring models, analytic hierarchy process -AHP, analytic network process -ANP, axiomatic design -AD, utility models, TOPSIS, ELECTRE and PROMETHEE.It is essential to develop all the elements related to the situation of MCDM in detail before selecting an appropriate MCDM method to solve the problem under consideration (Bufardi et al. 2004, Mergias et al. 2007).
Multi-criteria decision making (MCDM) methods, such as AHP and Fuzzy AHP, which are used for mining method selection problems in the literature, make the evaluations using the same evaluation scale and preference functions on the criteria basis.According to this, Ataei et al. (2008) have used the AHP approach for mining method selection.Otherwise, Bitarafan and Ataei (2004) have used different fuzzy methods as an innovative tool for criteria aggregation in mining decision problems.Also, Alpay and Yavuz (2009) have suggested a combination of AHP and fuzzy logic methods for underground mining method selection.The PROMETHEE method has some strength in the comparison with existing methods thanks to its simplicity and capacity to AHP AND PROMETHEE METHOD FOR MINING METHOD EVALUATION approximate the way that human mind expresses and synthesizes preferences when facing multiple contradictory decision perspectives.Decision Lab is a software that supports this method, and it also makes a sensitivity analysis for the results.This method provides a visual and powerful tool called Geometrical Analytic for Interactive Aid (GAIA) plane to identify conflicts among criteria and to group the alternatives (Albadvi et al. 2007).
In this paper, an AHP-PROMETHEE integrated approach for the selection of the most suitable mining method will be introduced, and the implementation process will be explained with a real world example.We shall use the AHP method to analyze the structure of the mining method selection problem and determine the weights of criteria.After that, we shall use the PROMETHEE method for final ranking.In this process the criteria, which have the greatest effect on the mining method selection, are determined by a sensitivity analysis.
This paper is divided into five sections.In Sec tion "Introduction", the studied problem is discussed.Section "AHP and PROMETHEE methods" briefly describes the two proposed approaches.In Section "AHPPROMETHEE integrated methodology", a proposed AHP-PRO-METHEE integrated approach for mining method selection is presented in detail.How the proposed approach is used on a real example is explained in Section "A numerical example of proposed methodology".In last Section "Conclusion", the concluding remarks are is discussed.

AHP METHOD
The Analytic Hierarchy Process (AHP) is a method developed by Saaty (1980) to support multicriteria decision making.This method is used in the decisionmaking process to help people setting priorities among the alternatives and making better decisions by taking into account qualitative and quantitative aspects of the decision.Lee et al. (2001) defines the AHP as a quantitative technique that facilitates the structuring of a complex multi-attribute problem and provides an objective methodology that is applied to a wide variety of decisions in the human judgment process.
The AHP involves decomposing a complex MCDM problem into a multi-level hierarchical structure of objectives, criteria and alternatives.The decomposition into a hierarchy is based on previous studies, research and empirical experiences.Once the hierarchy has been developed, one moves to assess the relative importance of decision criteria, then compare the decision alternatives with respect to each criterion, and finally determine the overall priority for each decision alternative and the overall ranking of the decision alternatives.
The assessment of the relative importance of decision criteria and the comparison of decision alternatives with respect to each criterion is done by a pair-wise comparison, which involves the following three tasks: developing a comparison matrix at each level of the hierarchy, starting from the second level and going down; computing the relative weights for each element of the hierarchy; and estimating the consistency ratio to check the consistency of the judgment.
Let {A 1 , A 2 , . . ., A n } be n alternatives, and {w 1 , w 2 , . . .,w n } be their current weights.The pairwise comparison is conducted by asking a decisionmaker or an expert questions, such as which criterion is more important with regard to the decision goal and by what scale (1-9), as shown in Table I.The answers to these questions are forming a pairwise comparison matrix that can be defined as follows:

Weight calculation
The matrix of the pairwise comparision A = [a i j ] represents the value of the expert's preference among individual pairs of alternatives (A i versus A j for all i, j = 1, 2, . . ., n).
After this, the decision-maker compares pairs of alternatives for all the possible pairs, and the comparison matrix A is obtained, where the element a i j shows the preference weight of A i obtained by comparison with A j .
The a i j elements estimate the ratios w i /w j , where w is the vector of current weights of the alternative (which is our goal).
The matrix has reciprocal properties, which are a ji = 1/a i j .
After all pair-wise comparison matrices are formed, the vector of weights w = [w 1 ,w 2 , . . .,w n ] is computed on the basis of Satty's eigenvector procedure in two steps.First, the pair-wise comparison matrix, A = [a i j ] nxn , is normalized, and then the weights are computed.Normalization (3) for all j = 1, 2, . . ., n. (4) For a valid comparison, we need to check the consistency of the pair-wise matrix (CI ).
where λ max is an important validating parameter in AHP and is used as a reference index to screen information by calculating the Consistency Ratio (CR) of the estimated vector.CR can be calculated using the following equation: where RI is the random consistency index obtained from a randomly generated pair-wise comparison matrix.
If CR < 0.1, then the comparisons are acceptable.If CR ≥ 0.1, the values of the ratio are indicative of inconsistent judgments.In this case, the original values in the pair-wise comparison matrix A should be reconsidered and revised.
Determining the overall priority for each decision alternative and the overall ranking of decision alternatives is done by synthesizing the results over all levels.The weighted priorities of the decision alternatives are added components in order to obtain an overall weight (wAi ) or priority of each alternative over the entire hierarchy.

PROMETHEE (Preference Ranking Organization
Method for Enrichment Evaluation) is an outranking method for a finite set of alternatives (Brans et al. 1984).These options include the choice of an appropriate preference function and the weighting given to each variable.The preference function defines how one object is to be ranked relative to another, and translates the deviation between the AHP AND PROMETHEE METHOD FOR MINING METHOD EVALUATION evaluations of two samples on a single parameter into a preference degree.The preference degree represents an increasing function of the deviation; hence, smaller deviations will contribute to weaker degrees of preference and larger ones to stronger degrees of preference.Six preference functions represented by specific shapes are available in the PROMETHEE method.Each shape is dependent on two thresholds, Q and P. Q is an indifference threshold representing the largest deviation that is considered negligible, and the preference threshold P represents the smallest deviation that is considered as decisive.P cannot be smaller than Q.The Gaussian threshold S is a middle value that is only used with the Gaussian preference function (Brans 1982, Brans andVincke 1985).
The PROMETHEE method is based on the calculation of positive flow (Ф + ) and negative flow (Ф − ) for each alternative according to the given weight for each criterion.The positive outranking flow expresses how much each alternative is outranking all the others.The higher the positive flow (Ф + → 1), the better the alternative.The negative outranking flow expresses how much each alternative is outranked by all the others.The smaller the negative flow (Ф − → 0), the better the alternative.The PROMETHEE II complete ranking is based on a calculation of net outranking flow value (Ф) that represents the balance between the positive and negative outranking flows.The higher the net flow, the better the alternative (Brans andMareschal 1994, Anand andKodali 2008).The stepwise procedure for achieving this outranking method is presented below: Step 1. Establishment of an impact matrix/double entry table.An impact matrix for the selected criteria ( j = 1 . . .n) and alternatives (i = 1 . ..m) can be established by using cardinal (quantitative) and ordinal (qualitative) data.
Step 2. Application of the preference function P(a, b).For each criterion, the selected preference function P(a, b) is applied to decide how much the outcome a is preferred to b.
Step 3. Calculation of an overall or global preference index Pi (a, b) that represents the intensity of preference of a over b.
PROMETHEE I provides a partial ranking of the alternatives and more realistic information about incomparability, while PROMETHEE II provides a complete ranking of the alternatives by calculating the net flow: Both partial ranking (PROMETHEE I) and complete ranking (PROMETHEE II) of the set of considered alternatives can be proposed to the stakeholders in order to solve the decisionmaking problem.
Finally, the relative position of the alternatives in terms of contributions to the various criteria are given by the geometrical analysis for interactive aid (GAIA).

AHP-PROMETHEE INTEGRATED METHODOLOGY
Both PROMETHEE and AHP methods have strengths and weaknesses, and can be singly used for mining method selection.In this paper, our idea is to lessen the weakness and enhance the strengths of these two methods by the process of integration and combination of their proceedings.Macharis et al. (2004) have analyzed the strengths and weaknesses of both PROMETHEE and AHP methods.They have made the comparative analysis of the following elements in both methods: the underlying value judgments, the structuring of the problem, the treatment of inconsistencies, the determination of weights, the evaluation elicitation, the management of the rank reversal problem, the support of group decisions, the availability of software packages and the possibility to visualize the problem.Based on this comparative analysis, we have concluded that a number of favorable characteristics of the AHP method could enhance PROMETHEE, namely at the level of structuring the decision problem and determining weights.The criteria weights, obtained by AHP, have a higher level of coherence, correlation, consistency and accuracy than weights determined on the basis of intuition or a domain specialist's knowledge, which is mostly used in the PROMETHEE method.
In the Data gathering stage, alternative mining methods and the criteria that will be used in their evaluation are determined, and the decision hierarchy is formed.
In the second stage, AHP computations are used for forming pairwise comparison matrices in order to determine the criteria weights.The individual evaluations are made by using the scale provided in Table I to determine the values of the elements of pairwise comparison matrices.Criterium Decision Plus software is used for computations in this process.
In the third stage -PROMETHEE computations, mining method priorities are found.Firstly, preference functions and parameters are determined by the authors.After this, the partial ranking with PROMETHEE I and the complete ranking with PROMETHEE II and GAIA plane are determined.Decision Lab software is used in this process.
In the last stage -Decision making, the best mining method is selected according to the rankings and GAIA plane obtained by PROMETHEE I and II.
The schematic representation of the proposed approach is presented in Figure 1.
The total ore reserves of "Coka Marin" are calculated as 1 160 000 tons, with an average grade of 1.0 'per cent ' Cu,2.6 'per cent' Pb,6.4 'per cent' Zn,2.71 gram per ton Au and 24.22 gram per ton Ag.Table II.shows physical and mechanical characteristics of "Coka Marin" deposit.

DATA GATHERING
On the basis of the physical and mechanical characteristics of this deposit (Table II), mining-technical and economic factors, the following mining methods are considered: room and pillar (A 1 ), room and pillar with fill (A 2 ), shrinkage stoping (A 3 ), cut and fill (A 4 ) and sublevel caving (A 5 ).
In Table III the criteria that have impact on the mining method selection are given.These criteria comprise most of all relevant factors for the undergorund mining method selection.
There are three levels in the decision hierarchy.The overall goal of the decision process -the selection of the best mining method -is on the first level of the hierarchy.The criteria are on the second level, and alternatives are on the third level of the hierarchy.

AHP COMPUTATIONS
On the basis of the decision hierarchy for mining method selection problem, the criteria to be used in the evaluation process are assigned weights by using AHP method.
Let A = {A 1 , A 2 , A 3 , A 4 , A 5 } be the set of possible mining alternatives, and C = {C 1 ,C 2 , . . ., C 11 } the set of selection criteria.A decision-maker is asked to form individual pairwise comparison matrix by using the scale given in Table I.Table IV shows the 11 × 11 pair-wise comparison matrix constructed to express the decision-makers' empirical estimate of the level of importance for each individual criterion.The maximum eigenvector was obtained from this matrix using the Criterium Decision Plus software.
Table V presents the results obtained from the computations based on the pairwise comparison matrix.
The excavation costs (C 9 ), thickness (C 1 ), ore excavation efficiency (C 7 ) and work safety (C 10 ) are determined as the most important criteria in the mining method selection process by AHP.The Consistency Ratio of the pairwise comparison matrix is calculated as 0.071 < 0.1.So, the weights are shown to be consistent, and they can be used in the decision making process.

PROMETHEE COMPUTATIONS
On the basis of the evaluation criteria, mining method alternatives are evaluated and the evaluation matrix is formed.In this process, some criteria have a quantitative structure, while others have a qualitative structure or an uncertain structure that cannot be accurately measured.Accordingly, some criteria (C 6 , C 7 and C 8 ) are expressed on quantitative scales, while others are qualitative.The usage of both quantitative and qualitative scales ensures that all criteria are properly processed on the best way.The qualitative evaluation has been done by an expert on a 5-point scale -Table VI.The worst category is very poor (numerical value 1), and the best category is very high (numerical value 5).
The evaluations of these five alternatives, according to the previously stated criteria, i.e., evaluation matrix, are displayed in Table VII.
After the evaluation matrix is determined, mining method alternatives are evaluated by the Firstly, the partial ranking is determined via PROMETHEE I (Fig. 3) on the basis of the flow values in Table VIII.PROMETHEE I uses positive and negative flow values to find the partial ranking.
Mining method A 2 is determined as the worst alternative according to the PROMETHEE I partial ranking.A 5 , A 1 , A 3 and A 4 alternatives are preferred to A 2 alternative.A 5 alternative is preferred to A 1 and A 4 alternatives.Also, A 3 alternative is preferred to A 1 and A 4 alternatives.On the other hand, A 5 and A 3 (also A 1 and A 4 ) are incomparable alternatives.It is obvious that PROMETHEE I did not provide information about the best alternative.
The best alternative is identified by PROMETHEE II complete ranking (Fig. 4).Net flow values given in the last column of Table VIII are used in this process.
A 5 alternative (sublevel caving) is selected as the best alternative based on the information provided by PROMETHEE II, and the other    alternatives are ranked in the order of A 3 (shrinkage stoping), A 1 (room and pillar), A 4 (cut and fill) and A 2 (room and pillar with fill).
The advantage of Decision Lab software package reflects in the use of GAIA (Geometrical Analysis for Interactive Assistance).Considering that 1 value is satisfactory (1 = 87.92%), the validity of using this tool can be discussed, where, 1 presents the measure of the amount of information being preserved by the defined model.In the real world, the applications of the value of 1 has always been larger than 60%, and in most cases larger than 80% (Brans and Mareschal 1994).
Based on the GAIA plane, it is possible to AHP AND PROMETHEE METHOD FOR MINING METHOD EVALUATION    easily determine the discriminative strength of each criterion, as well as the aspects of consistency and inconsistency as the quality of each alternative by every criterion.The alternatives are shown by triangles, and the criteria are presented as axes with square endings.The eccentricity of the position of square criteria is representing the strength of influence of that criterion, while the similarity in preference among certain criteria is defined with almost the same direction of axes of these criteria.
For the ranking, it is possible to determine the agreement among criteria C 1 , C 2 , C 3 , C 4 , C 5 , C 6 , C 9 and C 10 , while the criteria C 7 , C 8 and C 11 are not in compliance with other criteria.Also, the alternative position (triangles) determines the strength or weakness in relation to the alternative criteria.The closer orientation axis of the criteria to an alter-  native shows the better position of the alternative according to that criteria.The alternative A 5 (Cluster A) in Figure 5 can be determined as the best option because it is the closest alternative regarding the axis directions of criteria with the greatest impact (C 1 , C 3 , C 9 and C 10 ), and is directed to the nearest position of the decision stick pi, which defines the compromise solution in accordance with the weighted criteria.Contrary to it, the alternative A 2 is the worst one because it is good only by one weak criterion C 8 (Cluster B), and it is opposes to the direction of the decision stick pi which results have also been obtained through the PROMETHEE ranking.

DECISION-MAKING
According to the integrated AHP and PROME-THEE computations, it is decided that the sublevel caving (A 5 ) is the most suitable mining method for the "Coka Marin" underground mine.It was also analyzed how the variation in the criteria weights after the decision will affect the ranking.The sensitivity analysis has been performed, and the resulting "stability intervals" values are given in Table IX.Table IX gives for each criterion the limits within the weights' values, which can vary without changing the PROMETHEE II complete ranking.From the result of the sensitivity analysis, it is clear that Ore excavation efficiency (C 7 ) have the greatest impacts on the complete ranking.

CONCLUSION
In this paper, a decision approach is provided for mining method selection problem.Mining method selection is the most important decision made in a mine project.The selection of a suitable mining method for an ore deposit requires the consideration of a numerous criteria, such as mining-geological factors, economic factors and miningtechnical factors.This selection problem is based on the comparisons of mining method alternatives according to the identified criteria.
An integrated AHP and PROMETHEE decision making method has been used in the proposed approach for excavating the "Coka Marin" polymetallic ore deposit.The proposed approach differs from the present mining method selection in the literature.Here, AHP is used to assign weights to the criteria for mining method selection, while PROMETHEE is used for the complete ranking of the alternatives.The weights obtained from AHP are used in PROMETHEE computations, and the alternative priorities are determined based on these weights.In this study it was shown that the calculation of the criteria weights is important in the PROMETHEE method, and they could change the ranking.
The proposed integrated method can help decisionmakers to easily choose and analyze factors and attributes.The strengths of this approach over the existing methods are as follows: the PROMETHEE method takes into account the preference function of each criterion, which is determined by the decision-makers; each criterion is evaluated on a different basis, and it is possible to make better decisions.PROMETHEE I provides a partial ranking to identify the alternatives that cannot be compared and the alternatives that are indifferent, while PROMETHEE II provides a complete ranking for alternatives.The GAIA plane is a powerfull tool that suggests a differentiation power to the criteria, similar criteria, independent criteria and opposite criteria.By utilizing the PROMETHEE method to make a sensitivity analysis of the results, the most effective criterion in decision making is determined.These opportunities are not available in current methods, such as AHP, fuzzy AHP, ELECTRE, TOPSIS, etc.
Therefore, on the basis of the obtained results by the proposed integrated AHP and PROMETHEE method, the most suitable mining method for extracting the "Coka Marin" polymetallic ore is the sublevel caving method (alternative A 5 ).
The proposed model has only been implemented on a mining method selection problem in the RTB DEJAN BOGDANOVIC, DJORDJE NIKOLIC and IVANA ILIC Bor company, but the company management has found the proposed model satisfactory and implementable.Palavras-chave: AHP, mineração, processo de decisão baseado em múltiplos critérios, PROMETHEE.

TABLE I Pair-wise Comparison Scale for AHP preference.
--DEJAN BOGDANOVIC, DJORDJE NIKOLIC and IVANA ILIC

TABLE IX Stability intervals .
AHP AND PROMETHEE METHOD FOR MINING METHOD EVALUATION