*Article* **Application of Fuzzy Analytic Hierarchy Process to Underground Mining Method Selection**

### **Sanja Baji´c \*, Dragoljub Baji´c , Branko Glušˇcevi´c and Vesna Risti´c Vakanjac**

Faculty of Mining and Geology, University of Belgrade, 11000 Belgrade, Serbia; dragoljub.bajic@rgf.bg.ac.rs (D.B.); branko.gluscevic@rgf.bg.ac.rs (B.G.); vesna.ristic@rgf.bg.ac.rs (V.R.V.)

**\*** Correspondence: sanja.bajic@rgf.bg.ac.rs

Received: 30 December 2019; Accepted: 14 January 2020; Published: 23 January 2020

**Abstract:** The paper proposes a problem-solving approach in the area of underground mining, related to the evaluation and selection of the optimal mining method, employing fuzzy multiple-criteria optimization. The application of fuzzy logic to decision-making in multiple-criteria optimization is particularly useful in cases where not enough information is available about a given system, and where expert knowledge and experience are an important aspect. With a straightforward objective, multiple-criteria decision-making is used to rank various mining methods relative to a set of criteria and to select the optimal solution. The considered mining methods represent possible alternatives. In addition, various criteria and subcriteria that influence the selection of the best available solution are defined and analyzed. The final decision concerning the selection of the optimal mining method is made based on mathematical optimization calculations. The paper demonstrates the proposed approach as applied in a case study.

**Keywords:** multiple-criteria decision-making; underground mines; mining methods; expert knowledge

### **1. Introduction**

It is well-known that in most cases, a large number of criteria and subcriteria for decision-making matrices are uncertain and decision makers are unable to arrive at exact numerical values for comparing decisions. As such, mathematical methods are needed to effectively treat uncertainty, vagueness, and subjectivity. Viewed from that perspective, fuzzy logic is a scientifically based approach that relies on experience and intuition (or expert judgment). The fuzzy analytic hierarchy process (FAHP) enables the evaluation and analysis of criteria using fuzzified evaluation scales based on Saaty's scale [1].

In recent years, scientists worldwide have introduced a number of new theories and procedures for selecting underground mining methods, which generally involve gray correlation and multiple-criteria decision-making (AHP, FAHP, TOPSIS, PROMETHEE, ELECTRE, and VIKOR). Multiple-criteria decision-making (MCDM) methods have been demonstrated as useful problem-solving tools in various fields of engineering [2–4]. FAHP is widely applied. Guo et al. [5] used FAHP to determine evaluation index weights when they assessed the stability of a worksite above an abandoned coalmine, which threatened the safety of a high-speed railroad line. Pipatprapa et al. [6] used structural equation modeling (SEM) and FAHP to investigate factors suitable for assessing the environmental performance of the food industry. Lee et al. [7] proposed an FAHP-based decision-making model for selecting the best location for a frontal solar facility, given that the electric power demand, fossil fuel depletion, and environmental awareness necessitate power supply from renewable sources. Chatterjee & Stevi´c [8] used FAHP for supplier selection in supply chain management. Božani´c et al. [9] compare the FAHP method to another method that uses the fuzzy approach in MCDM for ranking the locations for deep wading as a technique of crossing the river by the army tank units. Stankovi´c et al. [10] used FAHP

for determining the importance of the traffic accessibility criteria. Mallick et al. [11] applied FAHP in groundwater management of a semi-arid region.

An in-depth literature review revealed that much research has used MCDM techniques to define optimal mining methods for different ores. Özfırat [12] applied FAHP to assess the use of certain machinery in the Amasra coalmine, in order to boost production, downsize the workforce and, consequently, reduce the number of accidents. Chander et al. [13] propose a decision-making technique for the selection of the optimal underground bauxite mining method. Based on AHP and VIKOR multiple-criteria optimization techniques, their results show that the optimal mining method, in that case, was cut-and-fill. Balusa & Gorai [14] compare mining methods using five MCDM models (TOPSIS, VIKOR, ELECTRE, PROMETHEE II, and WPM). They employed AHP to determine the weights of effective criteria for the Tummalapalle uranium mine in India. The results indicate that the selected mining methods were not equally efficient. Balusa & Gorai [15] used FAHP to select suitable underground mining methods. Bogdanovic et al. [16] employed a combination of AHP and PROMETHEE to select the most suitable mining method for the Coka Marin underground mine in ˇ Serbia: AHP to analyze the structure of the problem and determine criteria weights, and PROMETHEE for final ranking and sensitivity analysis. Alpay & Yavuz [17] developed a decision-making support system for the Karaburun underground chromite mine in Eskisehir, Turkey. They applied AHP to find acceptable alternatives. Yazdani-Chamzini et al. [18] proposed a selection model for the optimal mining method at the Angouran mine, one of the main producers of zinc in Iran. They developed the model based on FAHP and FTOPSIS. Then Asadi et al. [19] used a TOPSIS model to select the optimal mining method for the Tazareh coalmine in Iran. Javanshirgiv & Safari [20] applied fuzzy TOPSIS to select the optimal mining method for the Kamar Mahdi mine in Iran. Ataei et al. [21] also used TOPSIS to do the same for the Jajarm mine in Iran. For this mine, Naghadehi et al. [22] proposed a combination of FAHP and AHP: FAHP to determine criteria weights and AHP to rank the mining methods. On the other hand, some researchers have employed MCDM models to address mine dewatering, which is a parallel process in mining operations. Baji´c et al. [23] describe the selection of the optimal groundwater control system for the open cast-mine Buvaˇc (Bosnia and Herzegovina), using soft optimization and fuzzy optimization (VIKOR and FAHP) techniques. For the same case study, Polomˇci´c et al. [24] performed mathematical optimization calculations applying fuzzy dynamic TOPSIS.

The present paper describes and tests a decision-making algorithm for the selection of the optimal underground mining method. The algorithm is applied in a real case study to the Borska Reka copper mine (Serbia). First, the relevant alternatives are identified and then the selection criteria are analyzed. This if followed by MCDM, to select the optimal mining method. Finally, the best choice is the method that maximizes the output of useful components and minimizes tailings. In addition, the optimal solution involves the shortest mining time and the lowest consumption of energy and materials, along with full safety at work and no adverse effect on mine development.

### **2. Case Study**

The FAHP-based methodology for decision-making on the optimal underground mining method was applied in a real case study. The study area is the Borska Reka copper mine in eastern Serbia (Figure 1). In terms of regional metallogeny, the Bor ore field and Borska Reka copper mine belong the so-called Bor Zone, which coincides with the Timok igneous complex. In geologic terms, the sediments are composed of volcanites and volcanoclastic rocks, quartz-diorite porphyritic rocks, hydrothermally altered volcanic and volcanoclastic rocks, pelites with tuffs and tuffites, conglomerates, sandstones, Quaternary alluvial sediments, and technogenic deposits.

The mineral composition of the ore from Borska Reka includes chalcopyrite, covellite, chalcosine, rutile, hematite, magnetite, sphalerite, galenite, tetrahedrite, tennantite, digenite, cubanite, and native gold. The prevalent ore is pyrite, the dominant copper mineral is chalcopyrite, and there are covellite, chalcosine, and bornite to a lesser extent. On the other hand, enargite and molybdenite are very rare. However, this ratio of copper minerals is not uniform across the ore body. Certain parts have

elevated concentrations of covellite, chalcosine and bornite, but they are rarely dominant. There are also frequent occurrences of rutile, magnetite and hematite, as well as sphalerite and galenite. Tetrahedrite tennantite, digenite, cubanite and native gold are very rare and occur sporadically. dominant. There are also frequent occurrences of rutile, magnetite and hematite, as well as sphalerite and galenite. Tetrahedrite tennantite, digenite, cubanite and native gold are very rare and occur sporadically.

*Symmetry* **2020**, *12*, 192 3 of 23

Past exploration has revealed that the Borska Reka ore body is among very large deposits in the geometric sense, with elevated copper concentrations. The ore body is at an angle of 45◦–55◦ . Its maximum length is ~1.410 m and maximum width 635 m. The ore body is deep; the average ultimate depth is ~920 m from the ground surface. Past exploration has revealed that the Borska Reka ore body is among very large deposits in the geometric sense, with elevated copper concentrations. The ore body is at an angle of 45°–55°. Its maximum length is ~1.410 m and maximum width 635 m. The ore body is deep; the average ultimate depth is ~920 m from the ground surface.

**Figure 1.** Geographic position of the study area: Borska Reka copper mine. **Figure 1.** Geographic position of the study area: Borska Reka copper mine.

### **3. Methodology 3. Methodology**

methods.

The underground mining decision-making algorithm is shown in Figure 2. FAHP is the optimization technique. In general, one of the limiting factors of conventional methods applied to select the optimal mining technology is often a lack of data. Mines are complex geologic systems and mining operations are dynamic as the size and depth of the mine constantly increase (in plain view and elevation). As such, mining requires continual adaptation to new conditions. The contribution to science of decision-making methods based on fuzzy logic is the ability to focus on overcoming uncertainties inherent in mining method selection. The underground mining decision-making algorithm is shown in Figure 2. FAHP is the optimization technique. In general, one of the limiting factors of conventional methods applied to select the optimal mining technology is often a lack of data. Mines are complex geologic systems and mining operations are dynamic as the size and depth of the mine constantly increase (in plain view and elevation). As such, mining requires continual adaptation to new conditions. The contribution to science of decision-making methods based on fuzzy logic is the ability to focus on overcoming uncertainties inherent in mining method selection.

On the other hand, compared to other methods that include the fuzzy approach, FAHP offers certain specific advantages in optimal underground mining method selection. Because of the depth of the ore deposit and imprecise data typical of such a geologic system, which make it impossible to accurately define all the physical, mechanical and geologic conditions, the entire mining process requires constant "learning" and gradual, hierarchical problem-solving, to achieve the set objective. FAHP involves a continual "learning process", along with discussion among experts and prioritizing. On the other hand, compared to other methods that include the fuzzy approach, FAHP offers certain specific advantages in optimal underground mining method selection. Because of the depth of the ore deposit and imprecise data typical of such a geologic system, which make it impossible to accurately define all the physical, mechanical and geologic conditions, the entire mining process requires constant "learning" and gradual, hierarchical problem-solving, to achieve the set objective. FAHP involves a continual "learning process", along with discussion among experts and prioritizing.

Consequently, the use of FAHP highlighted the quality of this technique based on expert judgment or, in other words, reflected the decision-makers' knowledge and experience in evaluating information, to arrive at an optimal decision concerning multiple alternative underground mining Consequently, the use of FAHP highlighted the quality of this technique based on expert judgment or, in other words, reflected the decision-makers' knowledge and experience in evaluating information, to arrive at an optimal decision concerning multiple alternative underground mining methods.

*Symmetry* **2020**, *12*, 192 4 of 23

**Figure 2.** Multiple-criteria decision-making (MCDM) model for the selection of optimal underground mining method. **Figure 2.** Multiple-criteria decision-making (MCDM) model for the selection of optimal underground mining method. mining method.

FAHP is a combination of the conventional AHP method [1] and the fuzzy set theory [25]. It is

**Figure 2.** Multiple-criteria decision-making (MCDM) model for the selection of optimal underground

FAHP is a combination of the conventional AHP method [1] and the fuzzy set theory [25]. It is implemented using triangular fuzzy numbers [26]. *TFN* (Figure 3) in set *R* is a triangular fuzzy number if its membership function ெሺሻ: → ሾ0,1ሿ is defined as follows: FAHP is a combination of the conventional AHP method [1] and the fuzzy set theory [25]. It is implemented using triangular fuzzy numbers [26]. *TFN* (Figure 3) in set *R* is a triangular fuzzy number if its membership function µ*M*(*x*) : *R* → [0, 1] is defined as follows: implemented using triangular fuzzy numbers [26]. *TFN* (Figure 3) in set *R* is a triangular fuzzy number if its membership function ெሺሻ: → ሾ0,1ሿ is defined as follows: <sup>⎧</sup> − <sup>−</sup>

$$\mu\_{\rm{TFN}}(\mathbf{x}) \;= \begin{cases} \frac{\mathbf{x}}{\mathbf{s} - \mathbf{l}} - \frac{\mathbf{l}}{\mathbf{s} - \mathbf{l}'} & \mathbf{x} \in [\mathbf{l}, \mathbf{s}]\\ \frac{\mathbf{x}}{\mathbf{s} - \mathbf{d}} - \frac{\mathbf{d}}{\mathbf{s} - \mathbf{d}'} & \mathbf{x} \in [\mathbf{s}, \mathbf{d}]\\ 0, & \mathbf{x} \notin [\mathbf{l}, \mathbf{d}] \end{cases}$$

0, ∉ ሾ, ሿ

⎩

where ≤≤. where *l* ≤ *s* ≤ *d*. where ≤≤.

**Figure 3.** Membership function of triangular fuzzy number (TFN). **Figure 3.** Membership function of triangular fuzzy number (TFN).

**Figure 3.** Membership function of triangular fuzzy number (TFN). The modification of AHP into FAHP is in that the relative importance of the optimality criteria is described by linguistic variables [27], determined by the expert, and modeled by triangular fuzzy The modification of AHP into FAHP is in that the relative importance of the optimality criteria is described by linguistic variables [27], determined by the expert, and modeled by triangular fuzzy numbers (TFN). In other words, fuzzy numbers describe the pairwise comparison matrices of the The modification of AHP into FAHP is in that the relative importance of the optimality criteria is described by linguistic variables [27], determined by the expert, and modeled by triangular fuzzy numbers (TFN). In other words, fuzzy numbers describe the pairwise comparison matrices of the

numbers (TFN). In other words, fuzzy numbers describe the pairwise comparison matrices of the

optimality criteria. The fuzzified Saaty scale, proposed by many authors [26,28–30] is used. One of them is shown in Table 1.


**Table 1.** Fuzzified scale [26,30].

Chang [26] made the first development steps and Deng [30] modified the method. Also, Baji´c et al. [23] applied fuzzy optimization to mine hydrogeology. Based on the above, the FAHP analysis was implemented in the following steps:

*(a)* First the problem relating to the selection of the underground mining method was examined and the alternatives and criteria/subcriteria that influence the selection of the optimal alternative were identified. This involved the selection of a team of experts and the "exploitation" of their knowledge and experience.

*(b)* Pairs of criteria (Equation (1)), subcriteria (Equation (2)) and alternatives (Equation (3)) were evaluated and compared using the FAHP scale (Table 1):

$$A = \begin{bmatrix} a\_{11} & a\_{12} & \cdots & a\_{1m} \\ a\_{21} & a\_{22} & \cdots & a\_{2m} \\ \dots & \dots & \dots & \dots \\ a\_{m1} & a\_{m2} & \cdots & a\_{mm} \end{bmatrix} \tag{1}$$

where: *aij* = 1 for every *i* = *j*, (*i*, *j* = 1, 2, . . . , *m*) and *aij* = <sup>1</sup> *aji*

$$A\_{j} = \begin{bmatrix} a\_{11} & a\_{12} & \cdots & a\_{1k\_{j}} \\ a\_{21} & a\_{22} & \cdots & a\_{2k\_{j}} \\ \dots & \dots & \dots & \dots \\ a\_{k\_{j}1} & a\_{k\_{j}2} & \cdots & a\_{k\_{j}k\_{j}} \end{bmatrix} \tag{2}$$

where criterion *C<sup>j</sup>* is composed of *k<sup>j</sup>* subcriteria,

$$\begin{array}{c} \mathbf{Y}\_{k} = \begin{bmatrix} a\_{11} & a\_{12} & \cdots & a\_{1N} \\ a\_{21} & a\_{22} & \cdots & a\_{2N} \\ \dots & \dots & \dots & \dots \\ a\_{N1} & a\_{N2} & \cdots & a\_{NN} \end{bmatrix} \tag{3} \\ \tag{3} \\ \tag{4}$$

where *N* is the number of alternatives relative to each of the *K* subcriteria; *k* = 1, 2, . . . ,*K*.

*(c)* The weights of all three matrices from step **b** are determined gradually, using fuzzy extent analysis [26] or fuzzy arithmetic and the extension principle [31]. All the resulting weights are normalized:

$$w\_i = \sum\_{j=1}^{m} a\_{lj} \otimes \left[ \sum\_{k=1}^{m} \sum\_{l=1}^{m} a\_{kl} \right]^{-1} \tag{4}$$

where *i* = 1, 2, . . . , *m*

$$w\_j' = \left(\sum\_{l=1}^{k\_j} a\_{il} \otimes \left[\sum\_{i=1}^{k\_j} \sum\_{l=1}^{k\_j} a\_{il}\right]^{-1}\right) \otimes w\_j \tag{5}$$

where *j* = 1, 2, . . . , *m*; *p* = 1, 2, . . . , *k<sup>j</sup>*

$$\mathcal{W} = \left( w\_1^1, w\_1^2, \dots, w\_1^{k\_1}; w\_2^1, w\_2^2, \dots, w\_2^{k\_2}; \dots; w\_j^1, w\_j^2, \dots, w\_j^{k\_j}; \dots; w\_m^1, w\_m^2, \dots, w\_m^{k\_m} \right) \tag{6}$$

where *W* are subcriteria weights, whose total "length" is *K*

$$\mathcal{W} = (\mathcal{W}\_1, \mathcal{W}\_2, \dots, \mathcal{W}\_K) \tag{7}$$

*(d)* The next step is the application of the aggregation principle, to reduce two hierarchy tiers (criteria and subcriteria) to a single tier:

$$K = \sum\_{j=1}^{m} k\_j \tag{8}$$

where *C*1,*C*2, . . . ,*C<sup>m</sup>* is a set of *m* criteria, each with its subcriteria; *kj*—number of subcriteria of the *j*-th criterion.

*(e)* The fuzzy decision matrix and fuzzy performance matrix are now calculated. The fuzzy decision matrix results from calculations of the fuzzy extent analysis from step *c* for the alternatives:

$$\begin{array}{rcl} X &=& \begin{bmatrix} \mathfrak{x}\_{11} & \mathfrak{x}\_{12} & \cdots & \mathfrak{x}\_{1K} \\ \mathfrak{x}\_{21} & \mathfrak{x}\_{22} & \cdots & \mathfrak{x}\_{2K} \\ \cdots & \cdots & \cdots & \cdots \\ \mathfrak{x}\_{N1} & \mathfrak{x}\_{N2} & \cdots & \mathfrak{x}\_{NK} \end{bmatrix} \end{array} \tag{9}$$

and the fuzzy performance matrix represents the overall performance of each alternative relative to all the subcriteria:

$$Z = \begin{bmatrix} \mathbf{x}\_{11} \otimes \mathbf{W}\_1 & \mathbf{x}\_{12} \otimes \mathbf{W}\_2 & \cdots & \mathbf{x}\_{1K} \otimes \mathbf{W}\_K \\ \mathbf{x}\_{21} \otimes \mathbf{W}\_1 & \mathbf{x}\_{22} \otimes \mathbf{W}\_2 & \cdots & \mathbf{x}\_{2K} \otimes \mathbf{W}\_K \\ \cdots & \cdots & \cdots & \cdots \\ \mathbf{x}\_{N1} \otimes \mathbf{W}\_1 & \mathbf{x}\_{N2} \otimes \mathbf{W}\_2 & \cdots & \mathbf{x}\_{NK} \otimes \mathbf{W}\_K \end{bmatrix} \tag{10}$$

*(f)* The ultimate values of the alternatives are calculated in the form of triangular fuzzy numbers:

$$F\_i = \sum\_{j=1}^{K} \mathfrak{x}\_{ij} \otimes \mathcal{W}\_j \tag{11}$$

*(g)* The final step includes defuzzification [32], ranking of alternatives and, in parallel, sensitivity analysis [33,34]. The optimal alternative is the one with the greatest weight. The sum of the weights of all the alternatives is equal to zero:

$$defuxzy(A) = \frac{(d-l) + (s-l)}{3} + l \tag{12}$$

The sensitivity analysis is performed by introducing the optimization index λ. The "total integral"—*I* is calculated, to express the expert's risk assessment (0—pessimistic, 1—optimistic, and 0.5—moderate):

$$I = \frac{(d\lambda + s + (1 - \lambda)l)}{2}, \quad \lambda \in [0, 1] \tag{13}$$

where: *l*, *s* and *d* are elements of the triangular fuzzy number.

A special-purpose application, Fuzzy-GWCS2 based on Microsoft Excel, was developed for the above mathematical optimization calculations. The objective was to provide clearer insight into the results and facilitate monitoring of changes in the final calculations during the sensitivity analysis.

### **4. Results and Discussion**

FAHP-aided selection of underground mining methods enables efficient decision-making and facilitates solving of complex problems that involve vagueness and multiple uncertainties, like in the case of the Borska Reka copper mine. The first uncertainty was associated with the identification of all lithostratigraphic units. This copper deposit is highly specific, with occurrences of numerous minerals. On the other hand, there was a lack of information on its geometry and certain physical indicators and parameters pertaining to the ore and surrounding rocks.

Calculations were made using the MCDM model (or the algorithm shown in Figure 2), to determine the best mining alternative for the Borska Reka copper mine. The procedure was gradual, following the above steps *a–g* and using the specially developed Fuzzy-GWCS2 application.

The selection of a team of experts and "exploitation" of their knowledge and experience play a key role in decision- making and underground mining method selection. Teamwork ensures technically sustainable, economically viable and, above all, safe mining of copper ore.

Successful underground mining method selection requires substantial knowledge about the geology of the mineral ore deposit. In addition, the depth of the mine necessitates exploratory drilling experience. Because of the specific features of copper deposits, petrologists, mineralogists and geochemists contribute key knowledge and analysis of the minerals and their physical parameters. Hydrogeologists examine groundwater flow to the mine and ways of protecting the mine. Experts in economic geology and mine management assess the technoeconomic viability of mining. Geologists define the characteristics of the ore deposit. As a result, mining experts gain insight into applicable underground mining methods and develop alternative solutions. Engineers then synthesize the information and define and asses the criteria than affect the selection of the preferred underground mining method. The quality of the identification of mining conditions and the experts' knowledge and experience directly influence the selection of the optimal method.

The given problem—selection of the optimal underground mining method—was examined in step *a*. Then the criteria and subcriteria that influence the selection were defined. Based on literature sources that address the selection of underground mining methods and the governing factors [18,35], the following three criteria were identified: technical, production, and economic.

The criteria were subdivided into subcriteria, in this case 18, as shown in Table 2. Given the different types of essentially opposed criteria, the MCDM approach was a reasonable and justifiable choice.

In addition, five different alternatives (underground mining methods) were defined, including:

Alternative 1—sublevel caving; Alternative 2—cut and fill; Alternative 3—shrinkage stoping; Alternative 4—block caving; Alternative 5—vertical crater retreat (VCR)

Mining methods depend on the shape, size and depth of the ore body, physical and mechanical properties of the ore and accompanying rocks, hydrologic conditions, sensitivity of ground surface to mining, mineral and chemical composition of the ore, mineral distribution, and value of ore. Consequently, all these characteristics are important and need to be taken into account when a decision is made about the optimal mining method.

In the case of ore bodies of irregular shape, such as that at Borska Reka, priority is usually given to a caving method. The ore body size is often the decisive factor, because it reflects ore reserves. In addition, the thickness of the ore body is important, as are its depth, angle, type of contact and tectonic circumstances. At large depths, the cut-and-fill method should be given priority over block caving.


### **Table 2.** Defining of criteria and subcriteria.

In the present case study, the mineral and chemical composition was important because of the presence of pyrite and pyrrhotite. Copper pyrite ore, with more than 40% of sulfur, as well as other sulfide ores with elevated concentrations of pyrite and pyrrhotite, are susceptible to oxidation, self-ignition and sticking. If the ore is left in a crushed state for a long time, it becomes oxidized and warm in contact with air and humidity. This reduces the ore utilization rate. As such, if the ore contains large amounts of pyrite and pyrrhotite, cut-and-fill methods are given priority over shrinkage stoping or block caving methods. If the ore is highly valuable, often the method of choice is less effective but with a much higher ore utilization rate than vice-versa.

Ore impoverishment is the reduction in metal content of the produced ore, relative to that of the excavated block. Shrinkage stoping and caving methods typically lead to greater impoverishment. Also, even cut-and-fill methods, where the ore is loaded by means of mechanical devices (scrapers or shovels) directly from the fill, tend to result in a higher level of impoverishment. In general, however, cut-and-fill and block caving, compared to other methods, cause less impoverishment. In the case of the room-and-pillar methods, secondary cutting invariably leads to impoverishment because the ore is mixed with side or roof gangue.

Physical and mechanical characteristics are also taken into account. Rock hardness and stability tend to be the most important parameters because they affect the span and surface of the tunnels. With regard to hydrologic circumstances, the amount and properties of groundwater need to be known, particularly its effect on "plastic" rocks such as clays. The presence of water-bearing rocks and stagnant groundwater hinder the shrinkage stoping method. The preferred technologies are cut-and- fill (with hydraulic or paste backfill) or room-and-pillar mining methods. If the ore mineral distribution is not uniform, cut-and-fill methods are given priority.

Safety at work is an extremely important factor. The economic advantages of a given mining method should not threaten people's lives, operations or the safety of mine installations. The selected mining method should not be capable of causing fire, inrush of groundwater or surface water, caving of underground or above-ground mine walls or other structural components, or endanger miners and mining. A healthy work environment requires good ventilation (fresh air supply and venting of harmful gases and dust), proper illumination and safe access to work stations, as well as machinery to relieve miners of heavy manual work and measures to ensure health protection. Certain underground mining methods degrade the environment by damaging the soil and potentially causing land subsidence. On the other hand, the productivity of the mining method is important in the technoeconomic

assessment of a mine, given that a higher productivity leads to greater output. The productivity of a method is based on the rate of mining of blocks or parts of the ore body, and the capacity of the ore body depends on the ability to mine all active levels. The costs of mining are also examined, because spending is required before there can be a return on investment (such as for shafts or declines, mining equipment, crushers, transporters, and venting and dewatering systems). Also, there are costs associated with excavation (e.g., materials for tunneling and consumables such as ANFO explosives, detonators, drilling tools, diesel fuel, oil, lubricants, loader and transporter tires, steel, steel cables and cement, as well as the preparation and distribution of backfill paste, along with associated labor) and maintenance costs (of machinery and installations, as well as depreciation, overhead, etc.).

In view of the above facts and given that Borska Reka belongs to the group of ore bodies with relatively high copper concentrations, that the ore body is deep and that there are structures and facilities on the land surface, the research warranted the consideration of five high-productivity mining methods, which would ensure economically viable mining.

The criteria, subcriteria and alternatives were evaluated and the scores were input parameters for the MCDA model. Their weights determined in the form of fuzzy numbers per steps *b* and *c*. Equation (1) was used to evaluate the criteria, Equation (2) the subcriteria, and Equation (3) the alternatives. Evaluation was based on pairwise comparison (of criteria, subcriteria and alternatives), using linguistic variables and their numerical values from FAHP scales (Table 1). Table 3 shows the criteria scores in the form of triangular fuzzy numbers and their relative importance. Equations (4) through (7) were used to calculate weights by fuzzy extent analysis.


**Table 3.** Evaluation of criteria.

With regard to the selection of the most suitable underground mining method in the present case study, the technical and production criteria were given a slight advantage over the economic criterion. Then the subcriteria were evaluated. Given that each criterion was subdivided into a number of subcriteria (Table 2), this step involved the determination of the importance of all the subcriteria in a group, relative to each of the criteria. Table 4 shows the technical subcriteria scores. The relative weights of the technical subcriteria are presented in Table A1 (Appendix A).

**Table 4.** Evaluation of technical subcriteria.



**Table 4.** *Cont.*

Among the technical subcriteria, the most important were *T*1—ore body depth and *T*2—ore body thickness, which were given a slight advantage over the other subcriteria, per the FAHP scale.

Table 5 shows the relative scores of the production subcriteria. The relative weight of each subcriterion in the form of a fuzzy number is presented in Table A2 (Appendix A). Among the production subcriteria, a slight advantage, per the FAHP scale, was given to *P*2—safety at work, *P*3—environmental impact, and *P*6—ventilation.


**Table 5.** Evaluation of production criteria.

Table 6 shows the relative scores of the economic subcriteria. The relative weight of each subcriterion in the form of a fuzzy number is presented in Table A3 (Appendix A).


**Table 6.** Evaluation of economic subcriteria.

The alternatives were evaluated in the next step, by pairwise comparison relative to each subcriterion. This involved 18 comparisons. The results are shown in Tables 7–9. The respective calculated weight vectors are presented in Tables A4–A6 (Appendix A), based on the fuzzy extent analysis applying Equation (4) through (7).

Table A7 (Appendix A) shows (per step *(d)*) the ultimate weights of the subcriteria, calculated applying the aggregation principle according to Equation (8). The triangular fuzzy numbers of the criterion weights were multiplied by the weights of their subcriteria calculated in the previous step **(***c***)**. Hence, one tier was eliminated from the criteria-subcriteria-alternatives hierarchy.

Then, using the equations described in step *e*, the fuzzy decision matrix was calculated for the five alternatives (Equation (5)), as was the fuzzy performance matrix (Equation (10)), which represented the overall performance of each alternative relative to all the subcriteria. It was a result of multiplying all the subcriteria weights by the elements of the decision matrix (Table A8, Appendix A).


**Table 7.** Evaluation of alternatives relative to technical subcriteria.


**Table 8.** Evaluation of alternatives relative to production subcriteria.


**Table 9.** Evaluation of alternatives relative to economic subcriteria.

Per steps *f* and *g*, Tables 10 and 11 show the ultimate scores of the five alternatives in the form of fuzzy numbers, obtained by adding the fuzzy numbers—elements of the fuzzy performance matrix, according to Equation (11). Then the ultimate weights of the alternatives are shown in the form of non-fuzzy numbers, after defuzzification employing Equation (12). The final ranking of the alternatives is based on the sensitivity analysis per Equation (13). Figure 4 shows the total integral value of moderate, pessimistic and optimistic experts' risk assessments, or the weights of the alternatives relative to the optimization index parameters. If the decision-maker's inclination is optimistic (α = 1), the weights of the alternatives vary over a very narrow range, compared to pessimistic (α = 0) and moderate (α = 0.5). Based on the sensitivity analysis, the average differences between the weights of the alternatives were in the 0.1–0.73% range for an optimization index of 0.5, and 0.75–7.8% for an optimization index of 0. Defuzzification yielded the weights of the alternatives in the form of "normal" or real numbers. The highest weight is the best score. Based on the results, alternative 5 is the optimal underground mining method, followed in descending order by Alternative 2, Alternative 1, Alternative 3 and Alternative 4.

According to the MCDA model, Alternative 5 (VCR) was proposed as the optimal underground mining method for the Borska Reka copper mine. This method does not require extensive preparations, its productivity is high and the costs of mining are relatively low. There are also other advantages, such as a high ore utilization rate, low ore impoverishment, and a high level of safety at work, which was one of the most important evaluation factors in the case study. For all these reasons, the method proposed for the given copper mine provides optimal mining conditions.


**Table 10.** Ranking and optimal alternative.


**Table 11.** Sensitivity analysis. **Table 11.** Sensitivity analysis.

Optimal alternative A5

*Symmetry* **2020**, *12*, 192 16 of 23

**Table 10.** Ranking and optimal alternative.

A1 0.020 0.196 1.955 0.1978 3 A2 0.023 0.222 2.019 0.2048 2 A3 0.016 0.173 1.814 0.1833 4 A4 0.014 0.133 1.407 0.1421 5

**TFN Real Number Ranking** 

**Figure 4.** Total integral values of moderate, pessimistic and optimistic expert's risk assessments. **Figure 4.** Total integral values of moderate, pessimistic and optimistic expert's risk assessments.

According to the MCDA model, Alternative 5 (VCR) was proposed as the optimal underground mining method for the Borska Reka copper mine. This method does not require extensive preparations, its productivity is high and the costs of mining are relatively low. There are also other advantages, such as a high ore utilization rate, low ore impoverishment, and a high level of safety at work, which was one of the most important evaluation factors in the case study. For all these reasons, the method proposed for the given copper mine provides optimal mining conditions. This procedure does not complete the analysis of the mining problem. Management support strategies are developed for the upcoming period of mining. Such strategies enable the management This procedure does not complete the analysis of the mining problem. Management support strategies are developed for the upcoming period of mining. Such strategies enable the management team to assume full professional responsibility for improving development plans. Additional future activities are defined to ensure mining efficiency and high productivity. This also includes the use of the latest technological achievements that help upgrade mining safety. On the one hand, the implementation of solutions and management team's commitment contribute to sustainable development of the entire process of mining operations, while on the other hand, they contribute to long-term stable technical, economic, and production conditions.

### activities are defined to ensure mining efficiency and high productivity. This also includes the use of **5. Conclusions**

the latest technological achievements that help upgrade mining safety. On the one hand, the implementation of solutions and management team's commitment contribute to sustainable development of the entire process of mining operations, while on the other hand, they contribute to The paper demonstrated that FAHP is an extremely useful technique in the mining industry, given that the criteria used in the case study were subjective and based on expert judgment (of mining engineers and geologists), which is an important consideration in underground mining.

team to assume full professional responsibility for improving development plans. Additional future

long-term stable technical, economic, and production conditions. **5. Conclusions**  The research indicated that an interdisciplinary approach connects underground mining with other areas of science. For example, it links mining with fuzzy logic (based on mathematics and psychology) and multiple-criteria decision-making.

The paper demonstrated that FAHP is an extremely useful technique in the mining industry, given that the criteria used in the case study were subjective and based on expert judgment (of mining engineers and geologists), which is an important consideration in underground mining. The paper described and analyzed in detail the factors that influence the selection of the optimal underground mining method, including (i) technical (ore body depth, ore body thickness, ore body shape, value of ore, ore body slope, rock hardness and stability, type of ore body and contact with neighboring rocks, and the mineral and chemical composition of the ore), (ii) production (productivity, capacity, safety at work, environmental impact, ore dilution, ore impoverishment, ventilation and hydrologic conditions), and (iii) economic (capital expenditure, costs of mining and costs of maintenance). These criteria, along with their subcriteria, are deemed to be universal and applicable to other underground mines.

The practical importance of the proposed methodology was demonstrated in a case study that included the evaluation of criteria, subcriteria and alternatives applying FAHP, and decision-making/selection of the optimal underground mining method. This reflects the primary academic contribution and implications for further research.

In addition, the approach implemented fuzzy logic in multiple-criteria optimization related to underground mining. On the one hand, the objective of applying the fuzzy approach to decision-making and problem-solving in cases where there are several alternatives and analyzing the relevant factors is to arrive at the optimal solution. On the other hand, expert intuition and experience play an important role in the assessment of the ore system and underground mining methods, while fuzzy logic in mathematical calculations enables such a heuristic approach to problem solving. Such an interdisciplinary approach contributes to the quality and sustainable management of underground mining.

**Author Contributions:** All the authors jointly contributed to the finalization of the paper: S.B. designed the MCDM model and conducted optimization calculations; B.G. and V.R.V., along with S.B., analyzed the criteria and alternatives, and D.B. critically reviewed the concept and design of the paper. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Acknowledgments:** We express our gratitude to the Serbian Ministry of Education, Science and Technological Development, for funding projects TR33029, TR33003 and OI176022.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

### **Appendix A**


**Table A1.** Weights of technical subcriteria.






**Table A4.** Weights of alternatives relative to technical subcriteria.

**Table A5.** Weights of alternatives relative to production subcriteria.



**Table A5.** *Cont.*

**Table A6.** Weights of alternatives relative to economic subcriteria.


**Table A7.** Ultimate weights of subcriteria.



**Table A8.** Elements of the performance matrix.

### **References**


© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

MDPI St. Alban-Anlage 66 4052 Basel Switzerland Tel. +41 61 683 77 34 Fax +41 61 302 89 18 www.mdpi.com

*Symmetry* Editorial Office E-mail: symmetry@mdpi.com www.mdpi.com/journal/symmetry

MDPI St. Alban-Anlage 66 4052 Basel Switzerland

Tel: +41 61 683 77 34 Fax: +41 61 302 89 18

www.mdpi.com ISBN 978-3-0365-1573-1