Multi-Attribute Technological Modeling of Coal Deposits Based on the Fuzzy TOPSIS and C-Mean Clustering Algorithms

Miloš Gligorić 1,*, Zoran Gligorić 1,*, Čedomir Beljić 1, Slavko Torbica 1, Svetlana Štrbac Savić 2 and Jasmina Nedeljković Ostojić 3 1 Faculty of Mining and Geology, University of Belgrade, Ðušina 7, 11000 Belgrade, Serbia; cedomir.beljic@rgf.bg.ac.rs (Č.B.); slavko.torbica@rgf.bg.ac.rs (S.T.) 2 The School of Electrical and Computer Engineering of Applied Studies, Vojvode Stepe 283, 11000 Belgrade, Serbia; svetlanas@viser.edu.rs 3 Department of Geodesy, Belgrade University College of Applied Studies in Civil Engineering and Geodesy, Hajduk Stanka 2, 11000 Belgrade, Serbia; jasminaon@vggs.rs * Correspondence: milos.gligoric@rgf.rs (M.G.); zoran.gligoric@rgf.bg.ac.rs (Z.G.); Tel.: +381-11-3219-178 (Z.G.)


Introduction
Thermal power plants around the world primarily rely on coal produced by surface and underground mines.Today's coal mining industry and coal-based thermal power plants are faced with economic uncertainties and increasingly stringent environmental requirements.In such an environment managing the upstream process composed of coal mine and plant must be improved constantly.We focus on the coal mine, i.e., on the coal deposit representing the source of feedstock for the thermal power plant.Technological requirements of the thermal power plant are primarily related to some predefined values of the heating value, sulfur and ash content of the mined coal.Such requirements protect the power plant from being penalized by the environmental authorities, but on the other hand push coal mines into economically unviable situations.
Energies 2016, 9, 1059 2 of 23 One way to increase the efficiency of energy production in the power plants based on fossil fuels is to supply raw materials with specific and relatively stable quality parameters [1].Naworyta et al. developed a model based on the availability of coal stock and a blending yard for bed blending and availability of a second excavator and the possibility to blend two lignite streams on the belt conveyor [1].
Deposit characteristics belong to the set of the internal (endogenous) conditions influencing the efficiency of the upstream process.Such parameters cannot be managed by planners.The principal contribution to the efficiency is that flexible approaches, whose task is to describe as really as possible the production planning environment, get included into mining operations by way of creation the artificial technological models of deposits.In order to create a real model of a deposit it is necessary to have alphanumeric and graphical databases containing relevant data obtained by geological exploration.These databases serve as a basis for planning and design of mining operations, aiming at the economically viable exploitation of the mineral-raw material [2].
Different clustering algorithms are used to solving many problems in geology and mining.Abedi et al. used clustering methods, including self-organizing map and fuzzy c-means to prepare mineral prospectivity maps [3].Fink et al. have introduced a clustering method using clocked objective functions and Softassign techniques to optimize an appropriately formulated objective function that allows clustering between mutually constraining heterogeneous features.The heterogeneous features are spatial and mineral features [4].Richards et al. introduced a weighted method of clustering the individual units of a segmented image to analyze geologic maps generated from experts' analysis of remote sensing images [5].Ren and Qian proposed a novel power-efficient and anti-fading clustering based on a cross-layer that is specific to the time-varying fading characteristics of channels in the monitoring of coal mine faces with wireless sensor networks [6].Weintraub et al. proposed cluster analysis approach based on K-mean algorithm which approximates reasonably well the original large scale, detailed models.The model has been developed for planning in the copper mines of the CODELCO Company [7].
There are several mine planning algorithms available that can deal with multiple attributes, uncertainty and a large number of blocks.Gilani and Sattarvand [8] developed a stochastic open pit mine production planning algorithm based on an ant colony optimization approach to integrate geological uncertainty described through a series of the simulated ore bodies.Lamghari and Dimitrakopoulos [9] applied a network flow algorithm to prepare production plans in multi-processor open pit mines.Bendorf and Dimitrakopoulos [10] developed a stochastic production scheduling model to integrate geological uncertainty described by sets of equally possible scenarios of the unknown body.
Interaction between production planning models and the developed technological model is supported by the following facts: every mine production planning model contains operational constraints concerning the processing capacity.Within this constraint there is a parameter indicating that the block is an ore block or not.Similarly, the technological model gives information if a coal block belonging to a defined technological mining cut (cluster) is suitable for customers or not.The attribute variability magnitude in existing planning models is defined by the lower and upper value constraints.This means that we have 2n a constraints where n a is the total number of block attributes.The technological model significantly reduces the number of constraints i.e., only one constraint is employed.This means that if the relative closeness of the mining cut is greater than a predefined value, all blocks belonging to the cut will be processed, otherwise they will be treated as waste.Most models are based on the maximization of the net present value and they are static models with respect to asset (metal) price and production costs.Including volatility of the asset price makes them impossible to solve.The technological model provides an opportunity to overcome this situation and makes the planning process more realistic.Forecasted price and costs are used to calculate the unit value of the each cluster for every year of the project time (V cluster ij , where i denotes the i-th cluster in the j-th year of the project time).When we calculate the cash flow, it is only need to explore which blocks of the clusters can be mined for a given year.In this way we add a dynamic dimension to the problem of production planning and make it more realistic.The model also gives information about how much the quality of the block deviates from the desired value and helps planners deal with the process of blending.
The main aim of the developed model is to define the quality zones within a deposit with respect to the given technological requirements.In that sense, aggregation of mining blocks does not decrease the freedom of the optimizer through the planning process with respect to the operational constraints such as precedence constraint, mining capacity and reserve constraint.These constraints are not a function of quality.
If the production planning concerns a bench by bench scenario then a technological model is made for every bench separately.In the case of multiple bench mining all blocks are subject to clustering and we obtain a 3D model.The algorithm does not take into account boundary conditions.Once economic values for all blocks have been defined, an open pit optimal boundary can be determined by inputting the values in optimization algorithms such as the Lerchs-Grossmann algorithm, Korobov algorithm and floating cone algorithm [11][12][13].
The block model for a real life mineable reserve can contain hundreds of thousands to millions of blocks.In the presence of such a number of blocks making an optimal production plan is a very hard combinatorial and time consuming task.For that reason we intend to reduce the initial large scale problem to the small scale one, i.e., we try to create a base which will make mine production planning much easier.The reduction process is focused on the combination of mineable blocks into block aggregates based on the measurement of the relative closeness of block attributes to the given technological requirements.Calculation of the relative closeness can be treated as an Alternatives, Criteria, Evaluations model; where alternatives are mineable blocks, criteria are distances between block attributes and required values and evaluations represent the rating of the alternative with respect to defined criterion.Obviously, this is a multi-criteria decision-making problem.The uncertainty related to the input data is treated by fuzzy set theory, more specifically by triangular fuzzy numbers.In order to calculate the relative closeness the fuzzy Technique for Order Preference by Similarity to the Ideal Solution (TOPSIS) method is used.
Afterwards, the obtained set of the relative closenesses is clustered into an adequate number of block aggregates called technological mining cuts.Fuzzy C-mean clustering algorithm is used as a way to partition the coal deposit into an adequate number of technological mining cuts.The optimal number of cuts is obtained by the Fukuyama-Sugeno validity functional.If we use a different number of criteria in the relative closeness calculation it indicates a different number of ways of clustering.Selection of the optimal way of clustering is based on the comparison of the obtained adjusted Rand indexes, entropies and Fukuyama-Sugeno validity functionals.

The Concept of the Model
Creating the technological model of a coal deposit represents the approach that documents variability within a coal deposit with respect to the multi-attribute technological requirements of the coal customer, primarily related to the coal quality, ash and sulfur content; in our case these are coal thermal power plant requirements.Building the model can be treated as technological mapping.It is used to assist with new mine project development or major operating mine expansion prior to significant capital investment.The approach allows the location and quantification of zones of quality in 3D space, efficient mine production planning, economic analysis forecasting, reduction of project risks and even better environmental planning.Technological mapping is undertaken during the feasibility planning stages for a new mine project development or mine expansion as a support function for production planning.At this point geostatistical block model of deposit is developed and attributes are estimated and we use geostatistical information to build up our model.Since the technological model represents the basis for production planning it can be used for both short and long-term planning.In the presence of multiple customers the model must be modified.In such case the number of clusters is strictly equal to the number of customers and the fuzzy TOPSIS method measuring the relative closeness must be replaced by the fuzzy multi-attribute distance function measuring the distance between each block and the requirements of each customer.Clustering is now performed over the set of obtained distances.This algorithm does not allow criteria to be divided.
The traditional block model of coal deposit means the deposit is divided into an adequate number of blocks having the same size.Such a model is created using the data obtained from exploration drilling and application of geostatistical methods.Supposing what portion of the coal deposit is defined to be mineable in economical way (mineable reserves), i.e., the ultimate number of mineable blocks is defined.
Let B = {b k } k≥1 be the set of mineable blocks in the coal deposit.Each block is characterized by a certain number of attributes, such as dimensions, location, density, tonnage, heating value, sulfur and ash content.Without loss of generality we suppose the density and tonnage are the same for every block.With respect to the thermal power plant's requirements, heating value, sulfur and ash content are the main modeling attributes of.Such attributes are used to determine the possible technological value of every block in B.
One of methods of reducing the size of the production planning problem is to first combine all mineable blocks into block aggregates based on relative closeness of their attributes to the given technological requirements.Definition 1.Let {b k } k≥1 be a sequence of mineable blocks.A subsecuence of {b k } k≥1 is a sequence b k i i≥1 with k i < k i+1 for all i ≥ 1 [14].Definition 2. A point c is called a cluster point or accumulation point of a sequence {b k } k≥1 if for any ε > 0 and any K ≥ 1 there is an k ≥ Considering Definitions 1 and 2 we can create cluster composed of mineable blocks.We refer to this cluster composed of blocks as a technological mining cut.Obviously, the mineable reserves can be represented by the set composed of all technological mining cuts.Formally, the clustering structure of the coal deposit is represented as a set of the following subsets Definition 5. Let S = { s i } i=1,2,...,k be a set of the triangular fuzzy numbers describing the relative closeness of every mineable block to the A tar .
Considering Definitions 3-5 we can define technological mining cut (cluster) as a three attributes spatial object, composed of mineable blocks, with relative closeness belonging to one and only one predefined cluster interval.Therefore, the problem statement is to divide the area of the coal deposit into N technological mining cuts.The procedure of the solving the problem is divided into two main phases: • calculation of the relative closeness of every mineable block to the target values based on the Technique for Order Preference by Similarity to the Ideal Solution (TOPSIS) method, • clustering the obtained values based on the fuzzy C-mean clustering method.

The Relative Closeness
In order to define the relative closeness of block to the target values we apply the concept of multi-criteria decision making.To decrease the uncertainty of the input attributes we apply fuzzy set theory i.e., triangular fuzzy numbers [15,16].Application of the fuzzy set theory has found wide use for solving problems in the mining industry [17][18][19][20][21].
Let X be a classical set of objects, called the universe, whose generic elements are denoted by x.The membership in a crisp subset X is often viewed as characteristic function µ A from X to {0, 1} such that: where {0, 1} is called a valuation set.If the valuation set is allowed to be real interval [0, 1], A is called a fuzzy set and denoted by A and µ A (x) is the degree of membership of x in A.
In order to calculate the relative closeness, s i = f b act i , A tar , between the i-th mineable block and requirements of the thermal power plant we apply the modified TOPSIS method.Modification refers primarily to the different sequence of steps and calculation of attribute weights than the ones used in the original method.The TOPSIS method is a technique for order preference by similarity to an ideal solution proposed by Hwang and Yoon (1981) [25].The basic concept of this method is that the chosen alternative should have the shortest distance from the positive ideal solution and the farthest distance from the negative ideal solution.
Application of fuzzy TOPSIS can be found in many scientific papers.Chen [26] extended the use of TOPSIS to the fuzzy environment and gave numerical examples of system analysis engineer selection for a software company.Chu [27] presented a fuzzy TOPSIS model under group decisions for solving the facility location selection problem.Yang and Hung [28] proposed the use of TOPSIS and fuzzy TOPSIS methods for plant layout design problem.
The modified fuzzy TOPSIS method is based on the following steps: Step 1: Normalization of the input data The space of the input data is defined by the union of the mineable block attributes and target values.It can be represented by the following input data matrix: where k is the total number of the mineable blocks in the coal deposit and m is the total number of the attributes.In our case the input data matrix is as follows: Energies 2016, 9, 1059 Note, in the input data matrix, requirement of the thermal power plant is treated as a fictitious mineable block.Each element of the matrix Y is transformed using the following equations: Equation ( 4) refers to the block attributes while Equation ( 5) refers to the targets.The normalized input data matrix is as follows: Step 2: Weights of the attributes and targets In the original TOPSIS method the global weight of each criterion is calculated, while in our model we calculate the local weight of each attribute within the mineable block and technological requirement (targets) separately.The weight of a block attribute is calculated as follows: The weight of target value is calculated as follows: The weights of target attributes are assumed to have equal importance.
Step 3: Construct the decision matrix The first phase in solving the problem of the coal deposit partitioning into finite number of technological mining cuts can be treated as an Alternatives, Criteria, Evaluations model.A finite set of alternatives is defined by the set B = {b k } k≥1 , i.e., each mineable block represents one alternative.A finite set of criteria is defined by the set F = f 1 , f 2 , . . ., f j , where each element represents the distance function between the j-th block attribute and the j-th target value.A set of evaluations of alternatives with respect to given criteria is defined by the set X = x ij , where each element represents the rating of the alternative with respect to defined criterion.
Heating value, sulfur and ash content are used as three main attributes in the calculation of distance functions, i.e., in the process of the block evaluation.Accordingly we have the set F composed of three basic distance functions, i.e., the set of criteria is defined as F = { f 1 , f 2 , f 3 }.Evaluation of each element of the set X is obtained as follows: where term w j r j refers to the weighted normalized target value of the j-th attribute while w ij r ij to the actual value of the mineable block.
Energies 2016, 9, 1059 7 of 23 The criterion based on the heating value (f 1 ) is divided into two criteria with respect to the sign of the distance function (f 1 > 0; f 1 < 0).Suppose there is a sequence of mineable blocks with ascending order of the heating value and required (target) value, see Figure 1.
distance functions, i.e., in the process of the block evaluation.Accordingly we have the set F composed of three basic distance functions, i.e., the set of criteria is defined as = , , .Evaluation of each element of the set is obtained as follows: where term ̃ refers to the weighted normalized target value of the j-th attribute while ̃ to the actual value of the mineable block.
The criterion based on the heating value (f1) is divided into two criteria with respect to the sign of the distance function (f1 > 0; f1 < 0).Suppose there is a sequence of mineable blocks with ascending order of the heating value and required (target) value, see Figure 1.According to Equation ( 9) we obtain descending order of the distance function values, see Figure 2. The sign of the distance function is defined by the following sign function of f 1 : Energies 2016, 9, 1059 7 of 21 According to Equation ( 9) we obtain descending order of the distance function values, see Figure 2. The sign of the distance function is defined by the following sign function of f1: By the sign of the heating value distance function we generally separate the mineable blocks into two subsets, where the first subset is composed of the mineable blocks having a heating value greater than the target ( = −), while the second are the blocks having a lower value ( = +).
Accordingly, we obtain the two following criteria and .It is very important to emphasize that one and only one of these two criteria exists in the i-th mineable block.It is defined as follows: Note, [0, 1] indicates only the existence of the criterion in the i-th mineable block not the evaluated values.Criterion should be maximized while is minimized.Accordingly, the set of criteria is transformed and is as follows; = , , → , → .Evaluation of each element of the set with respect to criteria based on the sulfur and ash content (f3, f4) is also evaluated by Equation (11).Both criteria should be maximized.Finally we obtain the following decision matrix:  By the sign of the heating value distance function we generally separate the mineable blocks into two subsets, where the first subset is composed of the mineable blocks having a heating value greater than the target (sign( f i1 ) = −), while the second are the blocks having a lower value (sign( f i1 ) = +).Accordingly, we obtain the two following criteria f 1 1 and f 2 1 .It is very important to emphasize that one and only one of these two criteria exists in the i-th mineable block.It is defined as follows: Note, [0, 1] indicates only the existence of the criterion in the i-th mineable block not the evaluated values.Criterion f 1  1 should be maximized while f 2 1 is minimized.Accordingly, the set of criteria is transformed and is as follows; Evaluation of each element of the set X with respect to criteria based on the sulfur and ash content (f 3 , f 4 ) is also evaluated by Equation (11).Both criteria should be maximized.Finally we obtain the following decision matrix: Step 4: Define the ideal and the negative-ideal solutions Let us suppose that A + identifies the ideal solution and A − the negative one.They are defined as follows: where: With benefit and cost attributes, we discriminate between criteria that the decision maker desires to maximize or minimize, respectively.
Step 5: Measure the distance between alternatives and ideal solutions To calculate the m-Euclidean distance from each alternative (mineable block) to A + and A − the following equations can be easily adopted: Step 6: Measure of the relative closeness to the ideal solution The elements of the set S = { s i } i=1,2,...,k are calculated by the following equation: A very important concept related to the applications of triangular fuzzy numbers is the process of defuzzification.It converts a triangular fuzzy number into a crisp value The most commonly Energies 2016, 9, 1059 9 of 23 used defuzzification method is the centroid defuzzification method, which is also known as center of gravity [29].The defuzzification formula of triangular fuzzy numbers (a,b,c) is: and it will be used to express the fuzzy relative closeness as crisp value.Finally, we obtain the set of the relative closeness S = {s i } i=1,2,...,k that should be partitioned into adequate number of clusters.
Applying the same concept to the sulfur and ash content attribute we can obtain the following decision matrix:

Coal Deposit Partitioning Model
In order to divide a coal deposit into an adequate number of technological mining cuts we apply the fuzzy C-mean clustering algorithm [30][31][32][33] over the set S = {s i } i=1,2,...,k .
Algorithm aims to determine cluster centers c n (n = 1, 2, . . ., N) and fuzzy partition matrix U by minimizing the following function: subject to: The exponent ω is used to adjust the weighting effect of membership values.A large ω will increase the fuzziness of the function J.The value of ω is often set to 2. Applying partial derivative to the function J with respect to variable u ω ni and c n the following update equations are obtained: Based on Equations ( 21)-( 26) we describe the model to coal deposit partitioning into technological mining cuts as follows: Step 1: select an integer number of technological mining cuts i.e., clusters (N) and threshold value ε; let ω = 2; Step 2: input a set of initial cluster centers [c 1 , c 2 , . . . ,c n ], composed of the increasing order values randomly chosen from the interval [min{s i }, max{s i }], i = 1, 2.., k; Step 3: compute all (c n − s i ) 2 and then all u ni according to Equation (26); Step 4: update the set of initial cluster centers according to Equation (25); Energies 2016, 9, 1059 10 of 23 Step 5: compute the value of the objective function J according to Equation ( 21) and compare J (t+1) with J (t) , where t is the iteration number.If |J (t+1) − J (t) | < ε then stop otherwise return to Step 2.
Determination of the optimal number of technological mining cuts (N) is based on the fuzzy validity criterion.The Fukuyama-Sugeno validity functional is used to define it [34][35][36]: Since we have two clustering approaches based on four and six criteria, respectively, it is necessary to compare them.There are numerous measures for comparing clustering results.In this paper the adjusted Rand index (ARI) is used for comparison [37][38][39][40].
The range of ARI is 0 ≤ ARI ≤ 1, with only extreme values below zero.A value 0 indicates that there is no similarity, whereas a value of 1 indicates a similarity.If the value of ARI is high (ARI > 0.85) then mine production planners can select any one of the obtained models, otherwise it is necessary to select the optimal partitioning.
The selection is based on the entropy as a measure of quality of the obtained clusters [41,42].This measure considers the overlaps between clusters P and E. The entropy of cluster P and E is H(P) and H(E).The optimal technological model of the coal deposit is selected according to: min[H(P), H(E)].If there is no significant difference between H(P) and H(E), (<20%), then some additional way of selection have to be employed.For that purpose we can use Fukuyama-Sugeno validity functional (FS).For compact and well separated clusters it is expected small value of FS.Optimal technological model of coal deposit is selected according to: min[FS(P), FS(E)].A model of multi-attribute technological model of a coal deposit is represented in Figure 3.

Numerical Example
The developed model is tested on a small hypothetical coal deposit.A block model of the coal deposit was created on the basis of exploration drilling and geostatistical methods.The input parameters required for the model testing are given in Figure 4, Tables 1 and 2.

Numerical Example
The developed model is tested on a small hypothetical coal deposit.A block model of the coal deposit was created on the basis of exploration drilling and geostatistical methods.The input parameters required for the model testing are given in Figure 4, Tables 1 and 2.
Membership degree The process of updating the cluster center vector is represented in Table 5. Calculation of the optimal number of technological mining cuts based on the Fukuyama-Sugeno validity functional is represented for the case of four criteria and five cluster centers.The set of initial cluster centers is [c1, c2, c3, c4, c5] = [0.38,0.42, 0.46, 0.50, 0.54].Calculation of the first block membership degree indicating with what degree the relative closeness s1 belongs to the initial cluster center vector [0.38, 0.42, 0.46, 0.50, 0.54] is represented in Table 4.The process of updating the cluster center vector is represented in Table 5.After nine iterations threshold requirement is satisfied and technological model of coal deposit for matrix (78 × 4) and five cluster centers is obtained, see Figure 5. Fuzzy characteristics of obtained technological mining cuts (clusters) are represented in Table 6.Fuzzy characteristics of obtained technological mining cuts (clusters) are represented in Table 6.The Fukuyama-Sugeno validity functional was analyzed for values of N between 2 and 5 and the obtained results are represented by Figure 6.
Energies 2016, 9, 1059 16 of 21 The Fukuyama-Sugeno validity functional was analyzed for values of N between 2 and 5 and the obtained results are represented by Figure 6.The optimal number of technological mining cuts is 4 and the technological model of the coal deposit composed of four mining cuts is represented by Figure 7.According to Equation (20) the input decision matrix composed of six criteria is as follows: is obtained by applying Equations ( 14)- (19) and represented in Table 7.The optimal number of technological mining cuts is 4 and the technological model of the coal deposit composed of four mining cuts is represented by Figure 7.The Fukuyama-Sugeno validity functional was analyzed for values of N between 2 and 5 and the obtained results are represented by Figure 6.The optimal number of technological mining cuts is 4 and the technological model of the coal deposit composed of four mining cuts is represented by Figure 7.According to Equation (20) the input decision matrix composed of six criteria is as follows: is obtained by applying Equations ( 14)- (19) and represented in Table 7.According to Equation (20) the input decision matrix composed of six criteria is as follows: Triangular fuzzy number is represented in the vertical form (*). The set of the relative closeness S = {s i } i=1,2,...,78 is obtained by applying Equations ( 14)-( 19) and represented in Table 7.The same procedure is performed for the input decision matrix (78 × 6).Fukuyama-Sugeno validity functional was analyzed for values of N between 2 and 5 and obtained results are represented by Figure 8.The same procedure is performed for the input decision matrix (78 × 6).Fukuyama-Sugeno validity functional was analyzed for values of N between 2 and 5 and obtained results are represented by Figure 8.The optimal number of technological mining cuts is four and the technological model of coal deposit composed of four mining cuts is represented by Figure 9.

Discussion
Table 8 summarizes the characteristics of the two models obtained with respect to technological requirements and two sets composed of four and six criteria, respectively.
Distribution of the technological mining cuts with respect to defined technological requirements is represented by Figure 10.The optimal number of technological mining cuts is four and the technological model of coal deposit composed of four mining cuts is represented by Figure 9.The same procedure is performed for the input decision matrix (78 × 6).Fukuyama-Sugeno validity functional was analyzed for values of N between 2 and 5 and obtained results are represented by Figure 8.The optimal number of technological mining cuts is four and the technological model of coal deposit composed of four mining cuts is represented by Figure 9.

Discussion
Table 8 summarizes the characteristics of the two models obtained with respect to technological requirements and two sets composed of four and six criteria, respectively.
Distribution of the technological mining cuts with respect to defined technological requirements is represented by Figure 10.

Discussion
Table 8 summarizes the characteristics of the two models obtained with respect to technological requirements and two sets composed of four and six criteria, respectively.Distribution of the technological mining cuts with respect to defined technological requirements is represented by Figure 10.

Discussion
Table 8 summarizes the characteristics of the two models obtained with respect to technological requirements and two sets composed of four and six criteria, respectively.
Distribution of the technological mining cuts with respect to defined technological requirements is represented by Figure 10.Technological models of the coal deposit obtained for the input data matrices (78 × 4) and (78 × 6) are denoted by P and E, respectively.The information on the overlap between P and E is summarized in form of the following P × E confusion matrix: The elements of the alternative matrix (AV) are calculated using the values of elements of the matrix V as follows:  The value of ARI indicates there is no high degree of agreement between two obtained technological models.The entropy of model P and E is calculated as follows: k log v i. k = − According to: min[FS(P), FS(E)] = min[−0.324128,−0.083654] = −0.083654, the optimal technological model of coal deposit is also model P (see Figure 7).

Conclusions
Having the ability to make flexible production plans is recognized as a crucial element for the long term success of any mining company, especially in the coal sector.The efficiency of a plan depends directly on the quality of the input data.For that purpose we have developed the technological model of coal deposits to help mining engineers create production plans in an easier way.In nature this model is heterogeneous because the characteristic technological mining cuts can be located at different parts of the coal deposit.This means the cut is not continuous.Future research will be focused on the possibility of creation of a homogeneous model where each technological mining cut is continuous.

Figure 1 .
Figure 1.Ascending order of heating value and target value.

Figure 1 .
Figure 1.Ascending order of heating value and target value.

Figure 2 .
Figure 2. Descending order of the distance function values.

Figure 2 .
Figure 2. Descending order of the distance function values.

Figure 3 .
Figure 3. Model of the technological coal deposit.

Figure 3 .
Figure 3. Model of the technological coal deposit.

Figure 3 .
Figure 3. Model of the technological coal deposit.

Figure 4 .
Figure 4. Block model of a coal deposit.

Figure 4 .
Figure 4. Block model of a coal deposit.

Figure 7 .
Figure 7. Optimal technological model of the coal deposit for a (78 × 4) matrix.
is represented in the vertical form (*). The set of the relative closeness = , ,…,

Figure 7 .
Figure 7. Optimal technological model of the coal deposit for a (78 × 4) matrix.
is represented in the vertical form (*). The set of the relative closeness = , ,…,

Figure 7 .
Figure 7. Optimal technological model of the coal deposit for a (78 × 4) matrix.

Consequently, any mineable block in B belongs to exactly one and only one technological mining cut. Definition 3. Let B act = b act i
Let A tar = h tar , u tar , p tar be a set of the triangular fuzzy numbers describing the target values of attributes required by the thermal power plant.

Table 1 .
Attributes of the mineable coal block.

Table 4 .
Membership degree of the first block to the initial cluster center vector.

Table 5 .
Updating the cluster center vector.

Table 4 .
Membership degree of the first block to the initial cluster center vector.

Table 5 .
Updating the cluster center vector.
= 14 v .2= 21 v .3= 27 v .4= 16 v .. = 78 According to: min[H(P), H(E)] = min[0.5326,0.5879] = 0.5326, the optimal technological model of coal deposit is model P. Since the values of H(P) and H(E) are very similar we used the FS validity functional as additional way of selection.