Next Article in Journal
Numerical Simulation of Vertical Well Depressurization with Different Deployments of Radial Laterals in Class 1-Type Hydrate Reservoir
Previous Article in Journal
A Novel Approach to Dynamic Line Rating Calculation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determining the Availability of Continuous Systems in Open Pits Using ANFIS and a Simulation Model

1
Mining and Metallurgy Institute Bor, Zeleni Bulevar 35, 19210 Bor, Serbia
2
Faculty of Mechanical Engineering Belgrade, University of Belgrade, Kraljice Marije 16, 11000 Belgrade, Serbia
3
Faculty of Mining and Geology Belgrade, University of Belgrade, Djusina 7, 11000 Belgrade, Serbia
*
Author to whom correspondence should be addressed.
Energies 2024, 17(5), 1138; https://doi.org/10.3390/en17051138
Submission received: 22 January 2024 / Revised: 20 February 2024 / Accepted: 23 February 2024 / Published: 28 February 2024
(This article belongs to the Section H: Geo-Energy)

Abstract

:
This paper presents a model for determining the availability of continuous systems at open pits using the neuro-fuzzy system. The concept of availability is divided into partial indicators (synthetic indicators and sub-indicators). The presented model in relation to already existing models for determining availability uses a combination of the advantages of artificial neural networks and fuzzy logic. The case study addressed the I ECC (bucket wheel excavator–conveyors–crushing plant) system of the open pit Drmno, Kostolac. In this paper, in addition to the ANFIS model for determining the availability of continuous systems, a simulation model was developed. The obtained results of the ANFIS model were verified with the help of a simulation model that uses certain assumptions about the distribution of failures. This paper was created as a result of several years of field and theoretical research into the availability of continuous systems in open pits, and completes a cycle that consists of several published articles on the subject of modeling the behavior of these systems in real time using a time picture of the state, expert assessment, simulation and AI models, while respecting the multidisciplinarity of the problem (mining technological, mechanical, and information technological aspects). The developed ANFIS model is a key instrument for improving operational efficiency and resource management in the mining sector. Its ability to accurately predict the availability of the ECC system brings not only operational benefits through reduced downtime and optimized maintenance, but also a potential reduction in overall costs at coal open pits. Such an innovative model marks a significant step forward in the mining industry, especially when it comes to continuous systems in coal open pits.

1. Introduction

The surface exploitation of coal deposits is carried out in difficult and complex conditions. The operation of continuous systems and their availability is of great importance for the stability of the energy system of the Republic of Serbia, because the largest part of the coal needed for the operation of the power plant is obtained precisely by exploitation with continuous systems. About 70 percent of the electricity in the Republic of Serbia comes from coal.
Continuous surface mining systems are systems where the flow of material is continuous. The mechanization that is applied is very complex and made according to special requirements, because these systems must be adapted to specific working conditions. The basic function of these systems when it comes to surface coal mine is to excavate, transport and deposit coal, which can be uniquely described as coal production.
The application of expert systems based on fuzzy logic and neural networks has played an increasing role in mining. Given that mining as an industrial sector has a major role in global economic development, many authors have tried to improve the operations applied in mining, increase efficiency, and reduce risks to people and the environment by applying expert systems. In the following chapter, an overview of articles concerning the application of expert systems in mining is given.

2. Literature Review

Expert systems based on fuzzy logic and neural networks use the knowledge and rules of experts from a specific domain to make decisions or solve problems. In mining, expert systems can have different applications, intended to improve efficiency, safety and productivity, to optimize the planning and execution of mining operations, to increase the safety of the operations themselves, and to enable the more efficient management of resources while reducing costs.
By implementing expert systems, it is also possible to see equipment’s conditions, predict equipment availability, predict maintenance needs, and provide recommendations for the repair or replacement of parts, which would directly affect downtime reduction and increase equipment efficiency.
In the paper “Adaptive neuro-fuzzy prediction of operation of the bucket wheel drive based on wear of cutting elements” [1], Miletic et al. define an ANFIS model that aims to determine how the wear of cutting elements affects the operation of a bucket wheel excavator.
In the paper “A Fuzzy Expert Model for Availability Evaluation” [2], Ivezic and others define the concept of the availability of auxiliary machinery, such as bulldozers. The formed expert fuzzy model analyzes and integrates the partial indicators of the availability of bulldozers working at surface mines (open pits) of the Electric Power Company of Serbia (Beograd, Serbia).
Petrović and others, in the paper “Fuzzy Model for Risk Assessment of Machinery Failures” [3], present a model for the implementation of negative risk parameters in the synthetic risk assessment model of the Lokotrack LT 1213S mobile crusher operating at the open pit Ladna Voda. The model shows that there is a high level of risk and that it is necessary to introduce the concept of risk-based maintenance.
In the paper entitled “Applying the Fuzzy Inference Model in Maintenance Centered to Safety: Case Study—Bucket Wheel Excavator”, Jovančić et al. [4] use the Fuzzy Inference Model to promote safety-centered maintenance, which has undergone online adaptation to work conditions. The model was tested on a case study of two SRs 1200 and SchRs 630 bucket wheel excavators.
Monjezi M. and others, in the paper entitled “Evaluation of effect of blast design parameters on flyrock using artificial neural networks” [5], applied the method of artificial neural networks to predict the flying of fragments during blasting at the Sungun copper mine (Sungun), Iran. Several ANN (artificial neural networks) models were run, and it was observed that a model trained with a back-propagation algorithm with a 9-5-2-1 architecture gave the best results. The flight of pieces was calculated side by side on the basis of available empirical models. Statistical modeling was also performed to compare the predictive ability of the artificial neural network against other methods. The comparison of the results showed the absolute superiority of the artificial neural network.
In the article by Qin J. and colleagues entitled “SVNN-ANFIS approach for stability evaluation of open-pit mine slopes” [6], an analysis of the stability of open pit mine slopes was conducted. The study utilized the SVNN-ANFIS model with a singular value for slope stability assessment. The findings from the applied methodology indicate a training accuracy of 99.20% and a testing accuracy of 97.62%.
The article “Predictive Model of Rock Fragmentation Using the Neuro-Fuzzy Inference System (ANFIS) and Particle Swarm Optimization (PSO) to Estimate Fragmentation Size in Open Pit Mining” [7] by Betty Vergara and others describes a predictive model used to estimate rock fragmentation size using the ANFIS in combination with Particle Swarm Optimization (PSO). Utilizing statistical measures such as the correlation coefficient (R2) and mean square error (RMSE), the study determined that the ANFIS-PSO model, with an R2 value of 0.85 and an RMSE value of 0.78, can be considered a dependable and satisfactory model for predicting rock fragmentation in the field.
In the article “Lost production costs of the overburden excavation system caused by rubber belt failure” [8], Bugarić and others examines the average costs incurred due to malfunctions (resulting in lost production) in the overburden excavation system at the Tamnava-east field open pit mine. The focus is on failures of rubber belts in the belt conveyor system, which operates on a bucket wheel excavator, belt wagon, and spreader. The study determines the unit cost of system malfunctions per hour of belt conveyor operation over the belt’s lifetime. The analysis is based on a proposed methodology that considers the working time to failure of rubber belts. This methodology accounts for sudden failures (tear, breakthrough), described by an exponential distribution, and gradual failures, described by a normal distribution. This approach aids in planning for malfunctions, determining spare rubber belt requirements, reducing operational costs, and suggesting optimal maintenance strategies.
In the article entitled “Reliability of rubber conveyor belts as a part of the overburden removal system-case study: Tamnava-east field open cast mine” [9], Bugaric U. and co-authors define the reliability function for belt conveyors in an ECS system (comprising bucket wheel excavator, belt conveyors, and spreader) at the Tamnava-East field open pit mine in Kolubara. The reliability function is established based on the length of conveyor belts and working hours. The methodology for determining the reliability function involves analyzing the operational time until the failure of belt conveyors. This operational time is characterized by a combination of exponential distribution (representing time until sudden failures) and normal distribution (representing time until gradual failures). Additionally, the study notes a linear relationship between the length of belt conveyors and the mean operational time until gradual failures.
In the paper “Safety Analysis and Synthesis Using Fuzzy Sets and Evidential Reasoning” [10], Wang J., Yang JB., and Sen P. introduce a distinctive approach to analyzing the safety of intricate technical systems. The authors advocate breaking down the system structure into hierarchical levels for a comprehensive assessment. Fuzzy logic is employed to characterize individual failures, and the information is synthesized using fuzzy inference, providing a nuanced understanding of system safety.
In the paper entitled “Reliability of Hydraulic Installation of Auxiliary Mechanization Machines—Applying the Theory of Fuzzy Sets and Factual Reasoning” [11], Tanasijević M. and colleagues conduct a reliability analysis of a technical system. They utilize fuzzy sets theory and integrate the obtained information through factual reasoning. Beyond the conventional quantitative assessment represented by a reliability function, the authors also consider the experiences and insights of maintenance and operational personnel. This approach provides a more realistic depiction of the reliability of the examined technical system.
In the paper entitled “Performance Evaluation of Hybrid FFA-ANFIS and GA-ANFIS Models to Predict Particle Size Distribution of a Muck-pile After Blasting” [12], Zhou J. and collaborators employed an ANFIS model to predict the particle size distribution resulting from the blasting process. The key contribution of the paper lies in optimizing the parameters of the ANFIS model, both initial and consequential, using the Firefly Algorithm (FFA) and Genetic Algorithm (GA).
The authors of the article “Predicting the Risk of Fault-Induced Water Inrush Using the Adaptive Neuro-Fuzzy Inference System” [13] developed a method for forecasting the likelihood of fault-induced water inrush in underground engineering. They utilized the Adaptive Neuro-Fuzzy Inference System (ANFIS) and identified six key parameters related to the aquifer, the water-resisting properties of the aquifuge, and mining-induced stresses. The ANFIS model was trained using twenty documented cases of fault-induced water inrush, and its predictive performance was tested with five additional cases. The conclusive outcomes demonstrated perfect alignment between the predicted results and the actual occurrences for the five test cases.
In the article “Study on ANFIS Application in Coal Mining Stray Current Security Prediction” [14], the authors introduced an Adaptive Neuro-Fuzzy Inference System model specifically tailored for predicting stray current security in coal mining workfaces. This model is capable of predicting workface stray current security using easily measurable parameters from non-production fields. If the stray current surpasses the defined standard, the system issues timely alarms. Additionally, the study compared the accuracy rates of security predictions using different membership functions. The findings reveal that ANFIS, based on subtractive clustering, achieves the highest prediction accuracy and faster computational speed.
In the article titled “Prediction of backbreak in open pit blasting by adaptive neuro-fuzzy inference system (ANFIS) model” [15] authored by Bazzazi A.A. and colleagues, the application of the ANFIS model was explored for predicting cracks, an undesirable outcome of the blasting process. The assessment of the ANFIS model’s performance relied on metrics such as the root mean squared error (RMSE), variance accounted for (VAF), and correlation coefficient (R2). The results presented in this study demonstrate the model’s outstanding predictive capabilities, indicating its excellent performance in predicting the occurrence of cracks during the blasting process.
The authors of the article [16] introduce a novel approach to enhance the efficiency of ANFIS by utilizing the Mine Blast Algorithm (MBA) for optimization. This marks the first instance of applying MBA to ANFIS learning. The optimized ANFIS, achieved through MBA, is then utilized for predicting the strength of Malaysian small and medium enterprises (SMEs). The findings demonstrate that the ANFIS rule-base optimized by MBA outperforms Genetic Algorithm (GA) and Particle Swarm Optimization (PSO) in terms of efficiency in parameter training and overall accuracy.
The idea for the ANFIS model shown in this article, in addition to the mentioned works, came from articles published in the previous period related to this model [17,18,19].
In the article “Analytical Determination of the Availability of a rotary excavator as a part of coal mining system—case study: Rotary Excavator SchRs 800.15/1.5 of the Drmno open pit” [17], authored by Bugaric U. and others, a model is presented for analytically determining the availability of the rotary excavator SchRs 800.15/1.5 within the ECC system at the Drmno surface mine. Utilizing this methodology enables the efficient identification of key operational factors over time. Through the application of statistical methods and modeling the work process as a function of time, the study defines functional relationships for parameters like availability, duration of failure, and work duration. The statistical analysis of parameter values provides insights into the current operational stage of the rotary excavator. In this specific instance, the excavator SchRs 800.15/1.5, based on the Bathtub curve, is identified to be in the exploitation phase, mirroring the actual operational context. The determined parameters play a vital role in assessing the availability of the rotary excavator.
In the paper entitled “Determining the Availability of Continuous Systems at Open Pits Applying Fuzzy Logic” [18], Gomilanovic M. and others present a model for determining the availability of continuous systems at open pits using fuzzy logic. The used model was constructed by synthesizing independent partial indicators of availability. This model relies on an expert system to assess the availability of continuous mining systems. The availability of the system, as a complex state parameter, is broken down into partial indicators, i.e., reliability and maintainability, and the fuzzy compositions used for the integration of partial indicators are max–min and min–max compositions. The advantage of this model compared to conventional models is that it takes into account the influence of partial indicators of availability, and it does not require long-term monitoring and records to yield a snapshot of the system’s state.
In the paper entitled “A Model for Determining Fuzzy Evaluations of Partial Indicators of Availability for High-Capacity Continuous Systems at Coal Open Pits Using a Neuro-Fuzzy Inference System” [19], Gomilanovic M. and others present a model for determining the fuzzy ratings of partial indicators of the availability of continuous systems in surface coal mines using a neuro-fuzzy inference system. The key advantage of this model lies in the fact that it relies on historical data of the specific system (I ECC system of the Drmno surface mine), instead of the general experience of experts and the usual assumed values for fuzzy ratings of partial indicators. As a result, the model can more accurately predict the availability of continuous systems based on expert ratings in a certain period of time. Another advantage of this model is that availability is estimated on a quarterly basis, thus providing a more accurate picture because it focuses on a smaller time period with similar characteristics, thus including certain external influences related to quarterly meteorological conditions. Based on the displayed values of MAE and RMSE statistics, we conclude that the model using the Gaussian function has better prediction capabilities compared to the other displayed models using the Sigmoid and Bell functions.
This article was created as a synthesis of the three mentioned works. In the first paper [17], the availability of one part of the system (the bucket wheel excavator) was calculated using assumptions about the distribution of time in failure and time in operation, which served to create a simulation model that uses theoretical assumptions about the distribution of time of the system, itself composed of four parts (excavator, beltwagon, conveyors and crushing plant). In the second paper, a further composition of the concept of availability was created based on several partial indicators, with the creation of a certain hierarchy between them. This is the availability decomposition used in the development of this ANFIS model. The advantage of the new model compared to the model in paper [18] is that it does not use pre-defined stages of assessment for the mentioned partial sub-indicators, but rather, their form is developed based on historical data of system availability. In the third paper, it is shown how, on the basis of historical data on the availability calculated on a quarterly level, the parameters of the fuzzy assessment of partial indicators can be estimated. The advantage of the new model compared to [19] is that the set of partial indicators has been expanded, and a hierarchical structure has been introduced among them with the application of the IF-THEN rule.
For more detailed information on new modern techniques, and learning expert systems based on fuzzy logic and neural networks, see [20,21,22,23,24,25,26,27,28,29,30,31,32,33].

3. Case Study: I ECC System

The evaluation of the model in this article was performed by focusing on the I ECC system of the Drmno open pit. Continuous surface mining systems usually consist of a bucket wheel excavator, multiple conveyor systems and a stacker or crushing plant.
The I ECC system of the Drmno open pit consists of a bucket wheel excavator (SRs 400.14/1.5), a beltwagon (BRs 2400), conveyors and a crushing plant. Figure 1 shows the Drmno open pit.
Figure 2 shows part I of the ECC system at the Drmno open pit. A more detailed description of the I ECC system is given below in Figure 2.
The bucket wheel excavator SRs 400.14/1.5 is a compact excavator with a relatively short boom in relation to the diameter of the rotor. It can perform selective excavation with a maximum digging height of 14 m. The theoretical capacity of the excavator is 2800 m3/h. The beltwagon BRs 2400 comprises the connection between the excavator and the conveyor. It allows coal blocks that are further away from the conveyor position to be mined. Its capacity is slightly higher than the capacity of the excavator, 3000 m3/h.
The transport system consists of seven belt conveyors. Five belt conveyors have a belt width of 1800 mm, with a conveying speed of 5.2 m/s and capacity of 7200 m3/h. The last two conveyors have a belt width of 2000 mm. The crushing plant is located at the end of the last conveyor. In it, coal is crushed using hammer crushers to a size below 30 mm, suitable for the needs of the thermal power plant.

4. Availability

Availability is a frequently employed term in maintenance engineering, indicating the service quality of an engineering system, such as machines. It involves the examination of vulnerable areas and asset management, and plays a crucial role in decision-making throughout the lifecycle management process [18].
On the basis of the time state picture [34], wherein the times in the correct state alternate with the times in failure, the availability can be displayed and calculated [34]. The time during which the technical system is operational (in working order) is divided into:
  • Time while the technical system waits to be put into operation ( t 11 );
  • Time when the technical system is in operation ( t 12 ).
When the technical system is down, the time is divided into:
  • Organizational time ( t 21 );
  • Logistic time ( t 22 );
  • Active repair time ( t 23 ) (time for corrective repairs ( t 231 ) and time for preventive repairs ( t 232 )).
Outage times t 21 and t 22 refer to breakdown, interventions, the procurement of spare parts, tools, trained workforce and administrative work, etc. Active repair time includes the processes of repair, assembly, disassembly, replacement, etc. [19,34].
Availability can be calculated using Equation (1). The equation is presented as the quotient of the total time during which the technical system is in a correct state (operational) and the total time consisting of the time in the correct state and the time in failure [19,34]:
A ( t ) = t 11 , t 12 t 11 , t 12 , t 21 , t 22 , t 231 , t 232

5. Methods and Material

5.1. Development of ANFIS Model

The ANFIS model shown and described below was developed in the Python programming language in the PyCharm 2023.2.1 editor (Community Edition, Jet Brains), which is open to user access.
ANFIS systems exhibit a synergy of artificial neural networks and fuzzy logic (fuzzy inference system). The advantage of these systems is reflected in the combinations of their positive features, namely, the ability to learn with artificial neural networks and the use of expert knowledge with fuzzy logic.
The architecture of the ANFIS system bears resemblance to that of artificial neural networks, where, based on the set of input–output data, a corresponding fuzzy inference system is formed, and the parameters of the membership functions that transform the input data are calculated. The general structure of the ANFIS model consists of five layers (Figure 3). Below is a brief description of the layers.
In the first layer, the input data are transformed into a system of appropriate fuzzy sets:
O 1 d , i = μ d , i x , i = 1 , 2 ,
where x is the input argument of the first layer, and μ d , i is the membership function of the corresponding linguistic variable d .
In the second layer of the ANFIS model, the output parameters from various variables in the preceding layer are integrated. The determination of output data involves:
O 2 , d 1 , d 2 , i , j = ω 2 , d 1 , d 2 , i , j = μ d 1 , i x · μ d 2 , j y , i , j = 1 , 2 ,
where d 1 and d 2 are two different variables.
In the third layer, the values derived from the second layer undergo a normalization process. The normalization procedure is conducted in the following manner:
O 3 , d 1 , d 2 , i , j = ω d 1 , d 2 , i , j ¯ = ω d 1 , d 2 , i , j d 1 , d 2 ω d 1 , d 2 , i , j , i , j = 1 , 2 .
The following layer involves the amalgamation of normalized values from the preceding layer with first-order polynomials:
O 4 , d 1 , d 2 , i , j = ω d 1 , d 2 , i , j ¯ f d 1 , d 2 , i , j = ω ¯ d 1 , d 2 , i , j p d 1 , d 2 , i , j x + q d 1 , d 2 , i , j y + r d 1 , d 2 , i , j , i , j = 1 , 2 .
where p d 1 , d 2 , i , j , q d 1 , d 2 , i , j and r d 1 , d 2 , i , j are the parameters of the fourth layer model.
In the fifth and final layer, the normalized values from the preceding layer are summed using the following formula:
O 5 d 1 , d 2 , i , j = i , j , d 1 , d 2 , ω d 1 , d 2 , i , j ¯ f d 1 , d 2 , i , j = d 1 , d 2 , i , j ω d 1 , d 2 , i , j f d 1 , d 2 , i , j d 1 , d 2 , i , j ω d 1 , d 2 , i , j
In Figure 4, the general architecture of the ANFIS model described in the previous part is shown.
The training of a neuro-fuzzy system is best done by applying a back-propagation process that uses the R M S E as the error function, defined by:
R M S E = 1 n i = 1 n ( y i y ^ i ) 2 ,
where y 1 , y 2 ,…, y n are actual values, and y ^ 1 ,   y ^ 2 ,…, y ^ n are values predicted by the ANFIS model.
When the input membership function parameters are set, the output from the ANFIS model is calculated as follows:
f = w 1 w 1 + w 2 f 1 + w 2 w 1 + w 2 f 2 = w 1 ¯ f 1 + w 2 ¯ f 2
Using f 1 = x · p 1 + y · q 1 + r 1 and f 2 = x · p 2 + y · q 2 + r 2 , the following equality is obtained:
f = w 1 ¯ x p 1 + w 1 ¯ y q 1 + w 1 ¯ r 1 + w 2 ¯ x p 2 + w 2 ¯ y q 2 + w 2 ¯ r 2
The training process, referred to as model training, centers around adjusting parameter values based on the provided training data. Essential to this process is the utilization of the back-propagation method, an algorithm designed to minimize the error between the network’s output and the desired output.
The determination of the availability of continuous systems and its partial indicators was processed using the results obtained through questionnaires related to the expert assessment of partial indicators of availability, and to historical data on downtime and work, which include the time period from 2016 to 2019.
The ECC system’s availability is contingent on specific factors, commonly classified into two groups: partial indicators, such as reliability and maintainability. These synthetic indicators further rely on a multitude of independent parameters (sub-indicators) (Figure 5), all treated as variables within the context of this ANFIS model.
Within this model, availability decomposes into partial sub-indicators that are assessed by experts via a questionnaire. Each component of the I ECC system, including the bucket wheel excavator, beltwagon, belt conveyors, and crushing plant, undergoes evaluation.
In the expert evaluation, 10 experts specializing in continuous systems within surface mining were interviewed. They offered assessments for the sub-indicators of availability during specific quarters, encompassing the timeframe from 2016 to 2019, for each component of the ECC system.
Data from 2016–2018 were used to train the ANFIS model (480 data points—training data set), while data from 2019 (160 data points—test data set) were used to test the obtained model. The experts gave grades in the questionnaire ranging from F (the worst grade) to A (the best grade). The layout of the questionnaire is shown in Figure 6; in this questionnaire, the expert was required to make assessments at the quarterly level in a predetermined period of time for each part I of the ECC system. The scores obtained in this way have been used as input data for this model.
Prior to model development, a database was established concerning the durations of mechanical, electrical, and other failures within the ECC system spanning four years (2016, 2017, 2018, 2019). Information from this database is employed for calculating historical availability on a quarterly basis, serving as the output data for the ANFIS model. The availability for each quarter was computed using the Formula (1).
In Table 1, part of the database is shown. The data were taken from the Electric Power Company of Serbia and contain information about downtimes on the specific system in the specified time period.
The system’s availability was assessed on a quarterly basis (based on the available data), and the resulting values are presented in Table 2.
The resulting ANFIS model was given the survey results for all nine partial sub-indicators for each part of the I ECC system as input parameters, while the output represents the corresponding availability in the quarter to which the survey results refer, which was obtained based on historical data taken from the Electric Power Company of Serbia.
In the first step of the model, fuzzification was performed, which represents the transformation of partial indicator scores, using membership functions, into the corresponding j -scale for j = 10. Predefined fuzzy sets are not used for probability functions, but membership functions are used instead, the parameters of which are estimated within the model training process. The utilized membership functions include the Bell-shaped membership function, the Gaussian membership function, and the Sigmoid membership function.
Using IF-THEN rules that are pre-defined, the synthetic indicator R is determined based on the partial sub-indicators o , c and b , and the synthetic indicator M is determined based on the partial sub-indicators t , e , u , d , m and s .
In the following, we will illustrate the determination of the synthesis indicator R using IF-THEN rules based on sub-indicators o , c and b . Let the IF-THEN rule be defined by IF o i AND c j AND b k THEN R l , where i, j and k are in the set { A , B , C , D , E , F }, and l is in the set { A ,   B , C , D , E }. Then, the fuzzy sets come together:
μ o i x · μ c j y · μ b k z ,
where x , y and z are the input values of grades i , j and k , respectively. For partial indicators o , c and b , we assign the value l . The fuzzy set corresponding to the rating l of the indicator R is the sum of all fuzzy sets assigned the value l . In a similar way, on the basis of sub-indicators t , e , u , d , m and s , the synthesis indicator M is calculated.
In the next step, using the IF-THEN rules, as described in the previous paragraph, the availability indicator A is determined by synthetic indicators R and M . Then, the Euclidean distance of the obtained fuzzy sets from the fuzzy sets assigned to the availability indicator A is determined based on the corresponding membership functions whose parameters we estimate within this ANFIS model. The distances d 1 ,   d 2 , d 3 , d 4 and d 5 determined in this way can be joined by the normalized reciprocal values of the relative distances, determined by:
μ i = d m i n d i d m i n d 1 + d m i n d 2 + d m i n d 3 + d m i n d 4 + d m i n d 5 , i 1 , 2 , 3 , 4,5 .
These values belong to the appropriate set of grades that determine the indicator of availability, i.e.,
A = μ 1 , E , μ 2 , D , μ 3 , C , μ 4 , B , μ 5 , A .
Finally, the linguistic description is transformed into a numerical designation:
1 μ 1 + 2 μ 2 + 3 μ 3 + 4 μ 4 + 5 μ 5 μ 1 + μ 2 + μ 3 + μ 4 + μ 5 .
Dividing by 5 gives the predicted value of availability, which is compared with the realized value of availability, calculated on a quarterly basis.
The IF-THEN rules used in this ANFIS model are shown in Table 3, Table 4 and Table 5. So, for example, the values shown in the first type of this table are interpreted as follows:
If the partial sub-indicator o is F (the working environment conditions typically do not align with the requirements for the equipment in use), and if the partial sub-indicator c is F (write-off machine, very high level of failure) and the partial sub-indicator b is F (underdeveloped basic engineering), then indicator R is unreliable, E .
A summary of the models considered for predicting availability is given in Table 6.

5.2. Development of Simulation Model

During the creation of the simulation model, all failures were classified into one of three types of failure (mechanical, electrical, and others). As in the case of the ANFIS model, the simulation model used data from three years (2016, 2017 and 2018) to obtain results.
In Table 7, the experimental and theoretical frequencies of machine failures by interval are given.
The distributions of mechanical failure times, considered in the 96th percentile of the data, conform to the Weibull distribution, with parameters γ = 5 , β = 0.9511 and η = 18.4311 . More precisely, the empirical distribution function is determined by:
F x = 1 e x p ( x 5 18.4311 ) 0.9511
The model was developed based on a total of 1238 instances of mechanical failures.
The testing of the hypothesis regarding the distribution of data was performed with the help of the Kolmogornov–Smirnov test, whose statistic value n   D n is equal to 1.7944, so with a significance level of 0.001 we cannot reject the null hypothesis that claims that the data are in accordance with the Weibull distribution. The Figure 7 shows the experimental and theoretical functions of the distribution of mechanical failures.
In Table 8, experimental and theoretical frequencies of electrical failures by interval are given.
The distribution of the duration of electrical failures, considered in the 98.5th percentile of the data, is in accordance with the Weibull distribution, with parameters γ = 5 , β = 0.9066 and η = 28.6022 . More precisely, the empirical distribution function is determined by:
F x = 1 e x p ( x 5 28.6022 ) 0.9066
The number of electrical failures on which this model was developed amounted to 908 failures. The testing of the hypothesis regarding data distribution was performed with the help of the Kolmogornov–Smirnov test, whose statistic value n   D n is equal to 1.2804, so with a significance level of 0.05, we cannot reject the null hypothesis that claims that the data are in accordance with the Weibull distribution. Figure 8 shows the experimental and theoretical distribution functions of electrical failure.
In Table 9, experimental and theoretical frequencies of other failures by interval are given.
The distribution of the durations of other failures, considered in the 100th percentile of the data, is in accordance with the exponential distribution, with parameters γ = 5 and λ = 0.0103 . More precisely, the empirical distribution function is determined by:
F x = 1 e x p 0.0103 · ( x 5 )
The number of other failures on which this model was developed amounted to 2030 failures. The testing of the hypothesis regarding data distribution was performed with the help of the Kolmogornov–Smirnov test whose value of the statistic n   D n is equal to 1.5761, so with a significance level of 0.01, we cannot reject the null hypothesis, which claims that the data are in accordance with the exponential distribution. Figure 9 shows the experimental and theoretical distribution functions of other failures.
In Table 10, the experimental and theoretical frequencies of duration between failures by interval are given.
The distribution of the durations between failures, considered in the 95th percentile of the data, is in accordance with the Erlang distribution, with parameters γ = 20 , k = 2 and λ = 0.0057 . More precisely, the empirical distribution function is determined by:
F x = 1 1 + 0.0057 · x 20 exp ( 0.0057 · x 20 )
The number of times between failures on which this model was developed was 5212. The testing of the hypothesis regarding the distribution of data was carried out with the help of the Kolmogornov–Smirnov test, whose value for the statistic n   D n is equal to 1.0192, so with a significance level of 0.2 we cannot reject the null hypothesis that claims that the data are in accordance with the Erlang distribution. Figure 10 shows the experimental and theoretical time distribution functions between failures.
During the duration of mechanical and electrical failures, the parameter (β) of the Weibull distribution is close to unity, but less than 1, which indicates that the equipment, parts, etc., both mechanical and electrical, have an approximately constant intensity of maintenance, which is already the case with other failures (exponential distribution), which means that the total intensity of maintenance of the ECC system is approximately constant when t→∞, that is, it can be considered a function of the convenience of keeping (maintainability) the entire ECC system roughly exponential, and it represents Poisson’s recovery process. Below (Figure 11), we show the maintainability function for different types of failures.
The obtained distribution of times between failures, of Erlang order k = 2, indicates that with time, the intensity of failure increases (which is shown in the Figure 12) i.e., the ECC system is at the end of the “exploitation” period and the beginning of the “obsolescence” period (periods II and III of the “bathtub” curve). For k = 1 (exponential distribution), the system is in period II of “exploitation”. The following figure (Figure 12) shows the intensity of the failure.
Figure 13 shows the frequency distributions of the considered failure types.
The algorithm of the developed simulation model is shown in Figure 14. The theoretical distributions obtained using the K-S test are used for generating failure duration times and times between failures in the following way. For mechanical and electrical failures, Weibull distribution is used (Figure 7 and Figure 8), while exponential distribution is used for other failures (Figure 9). For generating the types of failures, empirical distribution, shown in Figure 13, is used. Times between failures are generated using Erlang distribution, as shown in Figure 10.
The simulation experiment is performed for a time period of one year (tsim), while the simulation time (t) is calculated in seconds. The number of simulation (Nosim) experiments is one hundred.
At the beginning, for t = 0, the initial state of the system is defined as the state of the ECC system “running” (State = “1”). At this point, the time (tbf) and type (VRflr) of first failure are also generated.
While the simulation experiment is being performed, the model alternately compares simulation time (t) with the times of failure beginning (tbf) and failure end (tef).
If the simulation time (t) is equal to the time of failure beginning (tbf), the state of the ECC system is changed to “downtime” (State = “0”), and the time needed for repair is generated according to the current type of failure and the distribution of failure time duration. In other words, the time when failure ends (tef) is generated.
If the simulation time (t) is equal to the time when the failure ends (tef), the state of the ECC system is changed to “running” (State = “1”). Also, the type of failure (VRflr) and the time of next failure are generated, meaning that the time of the beginning of the next failure (tbf) is generated.
After that, the availability of the ECC system is checked. If the state of the ECC system is “1” (“running”), then the variable AECC is increased by one. Also, the current simulation time and appropriate state of the ECC system are written onto a file.
When all Nosim simulation experiments are executed, the average availability of the ECC system AECC and the stationary value of the ECC system availability ka are calculated, as are the changes in the ECC system’s availability in time.
Glossary:
tsim—duration of the simulation (s);
Nosim—number of simulations;
State—state ECC system (1—“running”; 0—“downtime”);
rnumber—random number generated by uniform distribution in the interval [0…1];
cs—current simulation;
tsim—simulation time;
TBF—time between failures (current);
DT—downtime failures (current);
tef—failure completion time (in simulation);
tbf—failure start time (in simulation);
VRflr—type of failure (1—mechanical; 2—electrical; 3—other);
AECC—system availability—ECC;
A (t)—availability of the system at a given time t;
ka—stationary availability value.
Figure 15 shows the dependence of availability on time, obtained as a result of the simulation model.
Based on the simulation model, the mean availability value is 0.8513, i.e., 85%, and the stationary value is 0.8489, i.e., 85%.

6. Conclusions

Conventional methods are sensitive to the nature of the data, as well as to the presence of noise and outliers. In addition, another problem faced by conventional methods is their behavior in response to unseen data, indicating robustness and adaptability beyond the training set. Machine learning techniques have a tendency to solve these problems, which is shown in this work, where the ANFIS model provides a more accurate prediction of availability compared to the conventional method (simulation model).
Bearing in mind that both in this model and in the model [19] related to the calculation of availability for the forms of fuzzy numbers of partial indicators, those that use the Gaussian function were selected. There is a belief that with the provision of a larger number of data related to different systems, a model could be developed for using transfer learning techniques for specific types of systems.
An ANFIS model that uses a Gaussian function has better predictive power than models that use a Sigmoid or Bell-shaped function.
Using the obtained ANFIS model and the scores of partial sub-indicators provided by the experts for the next time period, the prediction of availability is obtained, which amounts to 0.8090, that is, 81%. Based on the simulation model, the availability value is 0.8513, i.e., 85%. Bearing in mind that the average availability in the time period of three years that was considered within the training of ANFIS and the simulation model is equal to 0.8132, i.e., 81%, while the value of availability for the year 2019 is 0.7797, i.e., 78%, we conclude that the ANFIS model gives a closer picture of the state of availability of the I ECC system. An additional advantage of this model is its simplicity, as it does not include certain assumptions about data distribution. Figure 16 shows the obtained availability values.
The assessment of availability using the ANFIS model gives necessary information related to production planning in surface mines with continuous systems. The obtained availability value represents a limitation regarding the realization of coal mining, as well as transport and storage capacity, and defines the need for possible interventions regarding its increase if there is a demand for it from the aspect of the realization of the projected/required capacity.

Author Contributions

Conceptualization, M.G. and U.B.; methodology, (ANFIS model) M.G. and (simulation model) U.B.; writing—review and editing M.G., M.B., S.S. and N.S.; supervision, M.B., N.S. and S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Ministry of Science, Technological Development and Innovation of the Republic of Serbia, as well as the Contract on realization and financing of the scientific research work of the Mining and Metallurgy Institute Bor in 2024, contract number: 451-03-66/2024-03/200052.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to company privacy.

Acknowledgments

Gratitude is given to the Ministry of Science, Technological Development and Innovation of the Republic of Serbia; Mining and Metallurgy Institute Bor, Zeleni bulevar 35, Bor; PE Electric Power Industry of Serbia, Balkanska 13, 11000 Belgrade; PE Electric Power Industry of Serbia—Branch “TE-KO Kostolac”, Nikola Tesla 5-7, 12208 Kostolac; Ministry of mining and energy Republic of Serbia, Nemanjina 22-26, Belgrade.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Miletic, F.; Jovancic, P.; Milovancevic, M.; Ignjatovic, D. Adaptive neuro-fuzzy prediction of operation of the bucket wheel drive based on wear of cutting elements. Adv. Eng. Softw. 2020, 146, 102824. [Google Scholar] [CrossRef]
  2. Ivezic, D.; Tanasijevic, M.; Jovancic, P.; Đuric, R. A Fuzzy Expert Model for Availability Evaluation. In Proceedings of the 2019 20th International Carpathian Control Conference (ICCC), Krakow-Wieliczka, Poland, 26–29 May 2019; pp. 1–6. [Google Scholar] [CrossRef]
  3. Petrovic, D.; Tanasijevic, M.; Stojadinovic, S.; Ivaz, J.; Stojkovic, P. Fuzzy Model for Risk Assessment of Machinery Failures. Symmetry 2020, 12, 525. [Google Scholar] [CrossRef]
  4. Jovancic, P.; Tanasijevic, M.; Milisavljevic, V.; Cvjetic, A.; Ivezic, D.; Bugaric, U. Applying the Fuzzy Inference Model in Maintenance Centered to Safety: Case Study—Bucket Wheel Excavator, Applications and Challenges of Maintenance and Safety Engineering in Industry 4.0; IGI Global: Hershey, PA, USA, 2020; pp. 142–165. [Google Scholar] [CrossRef]
  5. Monjezi, M.; Mehrdanesh, A.; Malek, A.; Khandelwal, M. Evaluation of effect of blast design parameters on flyrock using artificial neural networks. Neural Comput. Appl. 2013, 23, 349–356. [Google Scholar] [CrossRef]
  6. Qin, J.; Du, S.; Ye, J.; Yong, R. SVNN-ANFIS approach for stability evaluation of open-pit mine slopes. Expert Syst. Appl. 2022, 198, 116816. [Google Scholar] [CrossRef]
  7. Vergara, B.; Torres, M.; Aramburu, V.; Raymundo, C. Predictive Model of Rock Fragmentation Using the Neuro-Fuzzy Inference System (ANFIS) and Particle Swarm Optimization (PSO) to Estimate Fragmentation Size in Open Pit Mining. In Advances in Manufacturing, Production Management and Process Control; Trzcielinski, S., Mrugalska, B., Karwowski, W., Rossi, E., Di Nicolantonio, M., Eds.; AHFE 2021. Lecture Notes in Networks and Systems; Springer: Cham, Switzerland, 2021; Volume 274. [Google Scholar] [CrossRef]
  8. Bugarić, U.; Tanasijević, M.; Polovina, D.; Ignjatović, D.; Jovančić, P. Lost production costs of the overburden excavation system caused by rubber belt failure. Eksploat. I Niezawodn. Maint. Reliab. 2012, 14, 333–341. [Google Scholar]
  9. Bugaric, U.; Tanasijevic, M.; Polovina, D.; Ignjatovic, D.; Jovancic, P. Reliability Of Rubber Conveyor Belts As a Part of the Overburedn Removal System—Case Study: Tamnava-East Field Open Cast Mine. Teh. Vjesn. 2014, 21, 925–932. [Google Scholar]
  10. Wang, J.; Yang, J.B.; Sen, P. Safety analyses and synthesis using fuzzy sets and evidential reasoning. Reliab. Eng. Syst. Saf. 1995, 47, 103–118. [Google Scholar] [CrossRef]
  11. Tanasijevic, M.; Ivezic, D.; Ivkovic, S. Reliability of hydraulic installation of auxiliary mechanization machines—Application of fuzzy set theory and factual reasoning. J. Appl. Eng. Sci. 2006, 4, 37–45. [Google Scholar]
  12. Zhou, J.; Li, C.; Arslan, C.A.; Hasanipah, M.; Amnieh, H.B. Performance evaluation of hybrid FFA-ANFIS and GA-ANFIS models to predict particle size distribution of a muck-pile after blasting. Eng. Comput. 2019, 37, 265–274. [Google Scholar] [CrossRef]
  13. Zhou, Q.; Herrera-Herbert, J.; Hidalgo, A. Predicting the Risk of Fault-Induced Water Inrush Using the Adaptive Neuro-Fuzzy Inference System. Minerals 2017, 7, 55. [Google Scholar] [CrossRef]
  14. Ma, C.Y.; Zhang, D.L.; Wang, Z.; Li, G.X.; Tang, J.J. Study on ANFIS Application in Coal Mining Stray Current Security Prediction. In Key Engineering Materials; Trans Tech Publications, Ltd.: Wollerau, Switzerland, 2010; Volumes 426–427, pp. 216–219. [Google Scholar] [CrossRef]
  15. Bazzazi, A.A.; Esmaeili, M. Prediction of backbreak in open pit blasting by adaptive neuro-fuzzy inference system (anfis) model. Arch. Min. Sci. 2012, 57, 865–876. [Google Scholar]
  16. Hussain, K.; Salleh, M.N.M.; Leman, A.M. Optimization of ANFIS using Mine Blast Algorithm for predicting strength of Malaysian small medium enterprises. In Proceedings of the 2015 12th International Conference on Fuzzy Systems and Knowledge Discovery (FSKD), Zhangjiajie, China, 15–17 August 2015; pp. 118–123. [Google Scholar] [CrossRef]
  17. Bugaric, U.; Tanasijevic, M.; Gomilanović, M.; Petrović, A.; Ilić, M. Analytic determination of the availability of a rotary excavator as a part of coal mining system-Case study: Rotary excavator SchRs 800.15/1,5 of the Drmno open pit. Min. Metall. Eng. Bor 2020, 3–4, 25–36. [Google Scholar] [CrossRef]
  18. Gomilanovic, M.; Tanasijevic, M.; Stepanovic, S. Determining the Availability of Continuous Systems at Open Pits Applying Fuzzy Logic. Energies 2022, 15, 6786. [Google Scholar] [CrossRef]
  19. Gomilanovic, M.; Tanasijevic, M.; Stepanovic, S.; Miletic, F. A Model for Determining Fuzzy Evaluations of Partial Indicators of Availability for High-Capacity Continuous Systems at Coal Open Pits Using a Neuro-Fuzzy Inference System. Energies 2023, 16, 2958. [Google Scholar] [CrossRef]
  20. Zhang, T.; Zhang, Y.; Katterbauer, K.; Al Shehri, A.; Sun, S.; Hoteit, I. Deep learning–assisted phase equilibrium analysis for producing natural hydrogen. Int. J. Hydrogen Energy 2024, 50 Pt B, 473–486. [Google Scholar] [CrossRef]
  21. Lara-Cerecedo, L.O.; Hinojosa, J.F.; Pitalúa-Díaz, N.; Matsumoto, Y.; González-Angeles, A. Prediction of the Electricity Generation of a 60-kW Photovoltaic System with Intelligent Models ANFIS and Optimized ANFIS-PSO. Energies 2023, 16, 6050. [Google Scholar] [CrossRef]
  22. AL-Qaysi, A.M.M.; Bozkurt, A.; Ates, Y. Load Forecasting Based on Genetic Algorithm–Artificial Neural Network-Adaptive Neuro-Fuzzy Inference Systems: A Case Study in Iraq. Energies 2023, 16, 2919. [Google Scholar] [CrossRef]
  23. Olayode, I.O.; Tartibu, L.K.; Alex, F.J. Comparative Study Analysis of ANFIS and ANFIS-GA Models on Flow of Vehicles at Road Intersections. Appl. Sci. 2023, 13, 744. [Google Scholar] [CrossRef]
  24. Ahmadianfar, I.; Shirvani-Hosseini, S.; He, J.; Samadi-Koucheksaraee, A.; Yaseen, Z.M. An improved adaptive neuro fuzzy inference system model using conjoined metaheuristic algorithms for electrical conductivity prediction. Sci. Rep. 2022, 12, 4934. [Google Scholar] [CrossRef]
  25. Ye, J.; He, X.L. A novel hybrid of ANFIS-based models using optimisation approaches to predict mine blast-induced flyrock. Int. J. Environ. Sci. Technol. 2022, 20, 3673–3686. [Google Scholar] [CrossRef]
  26. Vesović, M.; Jovanović, R. Heat Flow Process Identification Using ANFIS-GA Model. In Proceedings of the Sinteza 2023—International Scientific Conference on Information Technology and Data Related Research, Singidunum University, Belgrade, Serbia, 19–20 May 2023; pp. 44–51, Computer Science and Artificial Intelligence Session. [Google Scholar] [CrossRef]
  27. Kim, M.H.; Song, C.M. Prediction of the Soil Permeability Coefficient of Reservoirs Using a Deep Neural Network Based on a Dendrite Concept. Processes 2023, 11, 661. [Google Scholar] [CrossRef]
  28. Ganjidoost, H.; Mousavi, S.J.; Soroush, A. Adaptive Network-Based Fuzzy Inference Systems Coupled with Genetic Algorithms for Predicting Soil Permeability Coefficient. Neural Process. Lett. 2015, 44, 53–79. [Google Scholar] [CrossRef]
  29. Oladipo, S.; Sun, Y. Enhanced adaptive neuro-fuzzy inference system using genetic algorithm: A case study in predicting electricity consumption. SN Appl. Sci. 2023, 5, 186. [Google Scholar] [CrossRef]
  30. Lin, S.-S.; Song, J.-H.; Zhu, K.-Y.; Liu, Y.-C.; Chang, H.-C. Applying Adaptive Neuro-Fuzzy Inference System to Improve Typhoon Intensity Forecast in the Northwest Pacific. Water 2023, 15, 2855. [Google Scholar] [CrossRef]
  31. Dadios, E.; Baun, J.J.; Enriquez, M.L.; Janairo, A.G.; Concepcion, R., II; De Leon, J.A.; Francisco, K.; Mayol, A.P.; Bandala, A.; Vicerra, R.R. Adaptive Neuro-Fuzzy Inference System-Based GPS-IMU Data Correction for Capacitive Resistivity Underground Imaging with Towed Vehicle System; IntechOpen: Rijeka, Croatia, 2023. [Google Scholar] [CrossRef]
  32. Valencia, N.; Fotouhi, A.; Shateri, N.; Auger, D. Development of a Hybrid Adaptive Neuro-fuzzy Inference System with Coulomb-Counting State-of-Charge Estimator for Lithium–Sulphur Battery. Int. J. Fuzzy Syst. 2023, 25, 407–422. [Google Scholar] [CrossRef]
  33. Ivančan, J.; Lisjak, D.; Pavletić, D.; Kolar, D. Improvement of Failure Mode and Effects Analysis Using Fuzzy and Adaptive Neuro-Fuzzy Inference System. Machines 2023, 11, 739. [Google Scholar] [CrossRef]
  34. Djenadic, S.; Ignjatovic, D.; Tanasijevic, M.; Bugaric, U.; Jankovic, I.; Subaranovic, T. Development of the Availability Concept by Using Fuzzy Theory with AHP Correction, a Case Study: Bulldozers in the Open-Pit Lignite Mine. Energies 2019, 12, 4044. [Google Scholar] [CrossRef]
  35. Milovancevic, M.; Cirkovic, B.; Denic, N.; Paunovic, M. Prediction of Shear Debonding Strength of Concrete Structure with High- 312 Performance Fiber Reinforced Concrete. Structures 2021, 33, 4475–4480. [Google Scholar] [CrossRef]
Figure 1. Drmno open pit with the plotted position of the I ECC system (first photo photographed by M.G., second photo: source—Google Earth).
Figure 1. Drmno open pit with the plotted position of the I ECC system (first photo photographed by M.G., second photo: source—Google Earth).
Energies 17 01138 g001
Figure 2. I ECC system at the open pit Drmno with parts of system—(a) Bucket wheel excavator, (b) Beltwagon BRs 2400, (c) Belt conveyors, (d) Crushing plant (photographed by M.G.).
Figure 2. I ECC system at the open pit Drmno with parts of system—(a) Bucket wheel excavator, (b) Beltwagon BRs 2400, (c) Belt conveyors, (d) Crushing plant (photographed by M.G.).
Energies 17 01138 g002
Figure 3. ANFIS layers [19,35].
Figure 3. ANFIS layers [19,35].
Energies 17 01138 g003
Figure 4. ANFIS model architecture [13].
Figure 4. ANFIS model architecture [13].
Energies 17 01138 g004
Figure 5. Presentation of partial availability indicators [18] (synthetic indicators, sub-indicators).
Figure 5. Presentation of partial availability indicators [18] (synthetic indicators, sub-indicators).
Energies 17 01138 g005
Figure 6. Layout of the questionnaire.
Figure 6. Layout of the questionnaire.
Energies 17 01138 g006
Figure 7. Experimental and theoretical distribution functions of mechanical failures.
Figure 7. Experimental and theoretical distribution functions of mechanical failures.
Energies 17 01138 g007
Figure 8. Experimental and theoretical distribution functions of electrical failure.
Figure 8. Experimental and theoretical distribution functions of electrical failure.
Energies 17 01138 g008
Figure 9. Experimental and theoretical distribution functions of other failures.
Figure 9. Experimental and theoretical distribution functions of other failures.
Energies 17 01138 g009
Figure 10. Experimental and theoretical time distribution functions between failures.
Figure 10. Experimental and theoretical time distribution functions between failures.
Energies 17 01138 g010
Figure 11. Maintainability functions for different types of failures.
Figure 11. Maintainability functions for different types of failures.
Energies 17 01138 g011
Figure 12. Failure intensity.
Figure 12. Failure intensity.
Energies 17 01138 g012
Figure 13. Frequency distribution of the considered failure types.
Figure 13. Frequency distribution of the considered failure types.
Energies 17 01138 g013
Figure 14. Algorithm of the simulation model.
Figure 14. Algorithm of the simulation model.
Energies 17 01138 g014
Figure 15. Dependence of availability on time.
Figure 15. Dependence of availability on time.
Energies 17 01138 g015
Figure 16. Display of obtained availability values.
Figure 16. Display of obtained availability values.
Energies 17 01138 g016
Table 1. Presentation of part of the database on downtimes of the I ECC system.
Table 1. Presentation of part of the database on downtimes of the I ECC system.
DateMonthYearSystemObjectFailureStart of FailureEnd of FailureDowntimeTotal Downtime (min.)NoteShift
1 January 2016January2016I ECCBWE SRs-400Electrical10:00:0010:50:0000:5050/1
1 January 2016January2016I ECCCrush. plantOther13:00:0014:30:0001:3090/1
1 January 2016January2016I ECCBWE SRs-400Electrical19:00:0019:10:0000:1010/2
Table 2. Obtained results for the availability of the I ECC system.
Table 2. Obtained results for the availability of the I ECC system.
YearQuarterAvailabilityYearQuarterAvailability
201610.8300201810.8039
20.820320.8425
30.801830.8365
40.800840.7790
YearQuarterAvailabilityYearQuarterAvailability
201710.8079201910.8100
20.782520.7500
30.837030.7831
40.817740.7758
Table 3. IF-THEN rules for determining the indicator R —reliability.
Table 3. IF-THEN rules for determining the indicator R —reliability.
o c b R
FFFE
EEED
DDDC
CCCB
BBBA
BCBA
CBCB
CBBA
BCDB
CCBB
DCDC
EBEC
CDBA
EEDD
CCAA
DCBB
BBAA
BCDB
ABAA
DBBA
DEAB
ACCA
ACDB
ABBA
AEDB
ABCA
BCEB
BDDB
BEEC
BBAA
AAAA
FEED
FDDC
FCCB
FAAA
FEDD
EDCC
DCBB
CBAA
Table 4. IF-THEN rules for determining the indicator M —maintainability.
Table 4. IF-THEN rules for determining the indicator M —maintainability.
t e u d m s M
FFFFFFE
EEEEEED
DCCCBBB
DCCCCCB
DCACBBA
CBBBCBA
BCBBBBA
BDDCCCB
CCDBBBB
DCDCDBB
EDCBBCB
CBBBBBA
BDCCBBB
CCBBDCB
BBBBBAA
BACCCBA
CCBCBBA
CCCCBBB
CDCCDCB
CCBBCDB
CABCDCB
BCDDDBB
DDDDDDC
CCCCCCB
BBBBBBA
AAAAAAA
EEEDDDC
DDDCCCB
CCCBBBA
BBBAAAA
FEDCBAB
EEDCBAB
DEDCBAB
CEDCBAB
BEDCBAA
AEDCBAA
Table 5. IF-THEN rules for determining A —availability.
Table 5. IF-THEN rules for determining A —availability.
R M A
DDD
DCC
DBC
DAB
CDC
BDC
CCC
BBB
AAA
EDD
CBB
BAA
CAB
Table 6. A summarized overview of the models being examined for predictive availability.
Table 6. A summarized overview of the models being examined for predictive availability.
ANFIS Parameters123
Inputs (number)999
Function typeGaussian functionBell-shaped functionSigmoid function
Number of membership functions6 × 6 × 6 × 6 × 6 × 6 × 6 × 6 × 66 × 6 × 6 × 6 × 6 × 6 × 6 × 6 × 66 × 6 × 6 × 6 × 6 × 6 × 6 × 6 × 6
Training data set480480480
Test data set160160160
Number of iterations101010
Number of fuzzy rules39 ( R ) + 36 ( M ) + 13 ( A )39 ( R ) + 36 ( M ) + 13 ( A )39 ( R ) + 36 ( M ) + 13 ( A )
R M S E (Training data set)0.00000.00000.0000
R M S E (Test data set)0.00140.00150.0015
480 data—training data set, given as an attachment to this article (4 parts of the system × 3 years × 4 quarters × 10 experts). 160 data—test data set, given as an attachment to this article (4 parts of the system × 1 year × 4 quarters × 10 experts).
Table 7. Experimental and theoretical frequency of mechanical failures by interval.
Table 7. Experimental and theoretical frequency of mechanical failures by interval.
No.The Lower Bound of the IntervalThe Upper Bound of the IntervalExperimental pdfExperimental cdfTheoretical pdfTheoretical cdfKS Test
15140.26820.26820.22300.22300.0452
214220.29480.56300.28900.51200.0510
322310.16070.72370.17650.68850.0353
431400.03630.76010.11090.79940.0393
540480.07270.83280.07060.87000.0372
648570.03720.87000.04540.91530.0454
757660.04120.91110.02930.94470.0335
866740.01290.92410.01910.96370.0396
974830.01450.93860.01240.97610.0375
1083920.02100.95960.00810.98430.0247
11921000.01620.97580.00530.98960.0138
121001090.00400.97980.00350.99310.0133
131091180.01370.99350.00230.99540.0019
141181260.00320.99680.00150.99700.0002
Table 8. Experimental and theoretical frequency of electrical failures by interval.
Table 8. Experimental and theoretical frequency of electrical failures by interval.
No.The Lower Bound of the IntervalThe Upper Bound of the IntervalExperimental pdfExperimental cdfTheoretical pdfTheoretical cdfKS Test
15200.29960.29960.25710.25710.0425
220350.25000.54960.29560.55270.0031
335500.16080.71040.16880.72150.0111
450650.11450.82490.10200.82350.0014
565800.05180.87670.06330.88670.0101
680950.03850.91520.03990.92670.0115
7951100.01650.93170.02550.95220.0204
81101250.01870.95040.01640.96860.0182
91251400.00660.95700.01070.97930.0222
101401550.01100.96810.00700.98630.0182
111551700.00990.97800.00460.99090.0129
121701850.00880.98680.00300.99390.0071
131852000.00550.99230.00200.99590.0036
142002150.00771.00000.00130.99730.0027
Table 9. Experimental and theoretical frequency of other failures by interval.
Table 9. Experimental and theoretical frequency of other failures by interval.
No.The Lower Bound of the IntervalThe Upper Bound of the IntervalExperimental pdfExperimental cdfTheoretical pdfTheoretical cdfKS Test
15530.38030.38030.39100.39100.0107
2531010.28370.66400.23810.62910.0350
31011490.13840.80250.14500.77410.0284
41491980.08230.88470.08830.86240.0223
51982460.04330.92810.05380.91620.0119
62462940.02170.94980.03280.94900.0008
72943420.01720.96700.02000.96890.0019
83423900.00690.97390.01220.98110.0072
93904380.01130.98520.00740.98850.0032
104384860.00540.99060.00450.99300.0023
114865340.00150.99210.00270.99570.0036
125345830.00150.99360.00170.99740.0038
135836310.00150.99510.00100.99840.0033
146316790.00200.99700.00060.99900.0020
156797270.00200.99900.00040.99940.0004
167277750.00101.00000.00020.99960.0004
Table 10. Experimental and theoretical frequencies of duration between failures by interval.
Table 10. Experimental and theoretical frequencies of duration between failures by interval.
No.The Lower Bound of the IntervalThe Upper Bound of the IntervalExperimental pdfExperimental cdfTheoretical pdfTheoretical cdfKS Test
1201020.08190.08190.07800.07800.0040
21021840.15560.23750.15600.23390.0036
31842660.17650.41400.16600.39990.0141
42663480.14500.55910.14760.54750.0116
53484300.10610.66520.12020.66770.0025
64305120.09670.76190.09300.76070.0012
75125940.06700.82890.06960.83020.0014
85946760.05240.88120.05080.88100.0002
96767580.03510.91630.03640.91740.0011
107588390.01980.93610.02570.94310.0070
118399210.02110.95720.01800.96110.0039
1292110030.01400.97120.01240.97350.0023
13100310850.00920.98040.00860.98210.0017
14108511670.00540.98580.00580.98790.0021
15116712490.00480.99060.00400.99190.0013
16124913310.00520.99580.00270.99460.0012
17133114130.00230.99810.00180.99640.0017
18141314950.00191.00000.00120.99760.0024
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gomilanovic, M.; Bugaric, U.; Bankovic, M.; Stanic, N.; Stepanovic, S. Determining the Availability of Continuous Systems in Open Pits Using ANFIS and a Simulation Model. Energies 2024, 17, 1138. https://doi.org/10.3390/en17051138

AMA Style

Gomilanovic M, Bugaric U, Bankovic M, Stanic N, Stepanovic S. Determining the Availability of Continuous Systems in Open Pits Using ANFIS and a Simulation Model. Energies. 2024; 17(5):1138. https://doi.org/10.3390/en17051138

Chicago/Turabian Style

Gomilanovic, Miljan, Ugljesa Bugaric, Mirjana Bankovic, Nikola Stanic, and Sasa Stepanovic. 2024. "Determining the Availability of Continuous Systems in Open Pits Using ANFIS and a Simulation Model" Energies 17, no. 5: 1138. https://doi.org/10.3390/en17051138

APA Style

Gomilanovic, M., Bugaric, U., Bankovic, M., Stanic, N., & Stepanovic, S. (2024). Determining the Availability of Continuous Systems in Open Pits Using ANFIS and a Simulation Model. Energies, 17(5), 1138. https://doi.org/10.3390/en17051138

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop