Application of a Novel Hybrid f-SC Risk Analysis Method in the Paint Industry

: Occupational health and safety (OHS) activities aim to ensure the full mental and physical well-being of employees in the workplace. For this reason, it is essential to determine the precautions and classical Fine–Kinney methods were similar to the developed f-SC method. This research provides three contributions: (1) criteria must be weighted to determine risk scores, (2) using intuitive linguistic terms in scoring criteria made the risk analysis method more sensitive and appropriate, and (3) using MCDM methods instead of classical methods for the risk analyses in the OHS ﬁeld removes uncertainties.


Introduction
Despite the developments in technology, new difficulties have arisen in the field of occupational health and safety (OHS). For this reason, various accidents result in the death Sustainability 2021, 13, 13605 2 of 23 or injury of employees in the workplace [1]. OHS activities include hazards that may arise from processes carried out in the workplace as well as efforts to eliminate risks that may arise as a result of these hazards [2]. Each business can have its own sources of hazards, or hazards can be similar to those of other businesses. Identifying these hazards in risk assessment procedures is one of the most important steps to be carried out in the workplace. Various quantitative and qualitative risk analysis methods are used in the literature for risk assessment.
The most important stage of the risk analysis is the determination of the risk score. Necessary precautions are taken for the related OHS activities, and the time period required for the implementation of these measures is determined according to the obtained risk score. A mistake made at this stage may result in the late implementation of a precautionary measure. This can cause unwanted accidents to occur. For this reason, it is necessary to carefully identify potential hazards and sources of hazards while conducting risk analyses. In addition, the risk analysis used must be quite sensitive [3,4].
Nowadays, the use of multi-criteria decision-making (MCDM) methods has become quite common in decision-making processes. With these methods, the aim is that correct decisions are made from the simplest to the most complex state. In MCDM methods, the most important issue for a decision to be made precisely and accurately is the determination of the criteria. Making a mistake in the setting criteria can lead to irreversible results in the decisions made. There are many MCDM methods that have been used in different studies. In this context, the use of MCDM methods is becoming more widespread in the field of OHS in today's risk analysis [5].
In general, the risk score of the classical risk analysis methods used for risk assessment in the workplace is calculated by the mathematical multiplication of criteria given numerical values. For this reason, this classical approach creates disadvantages and uncertainties. Although the assigned numerical values have meaning, using linguistic terms in assigning the criteria and its alternatives is more intuitive for the decision maker, which will help to remove these disadvantages. Fuzzy logic appears to be a suitable method for this process. After expressing the criteria in terms of definitions, numerating these expressions with the help of fuzzy numbers will give more meaningful results and help to eliminate possible human errors [6].
In decision making analyses, every criterion is evaluated against each other, and these criteria are weighted because the importance level of each criterion is different, just as in classical risk analysis methods. For example, in the classical Fine-Kinney method, the risk score is obtained by multiplying the values of Probability (P), Severity (S), and Exposure (E), which are determined according to a relevant scale. However, in this method, a risk with P, E, and S values of 6, 6, and 3, respectively, is prioritized, compared to a risk with values of 1, 2, and 40. In other words, a situation that will cause a fatal accident or occupational disease in the case of an accident becomes less important than a situation that will result in a simple injury. For this reason, these disadvantages should be eliminated by weighting parameters with the help of decision-making methods.
In this study, a classical Fine-Kinney risk analysis conducted in the paint industry was modeled with a new hybrid MCDM method called f-SC, based on fuzzy logic. First of all, the criteria used in the risk score calculation was weighted with the fuzzy SWARA method by OHS experts. Afterwards, verbal expressions were used instead of the numerical values used in classical methods, and analysis was completed with the help of the fuzzy COPRAS method.
It is mandatory to make risk assessments and employ OHS experts in the workplace in Turkey because of the nation's 2012 enactment of the OHS law, numbered 6331. The using of risk analysis methods increased in workplaces in Turkey because of this. Although a risk analysis conducted for any industry has to be suitable for that industry, many companies use the classical Fine-Kinney method because of its simplicity. However, this method has limitations, e.g., the numerical data given the criteria cause uncertainties when calculating the risk score, and in this method, all the criteria have same weight. Thus, all of them affect Sustainability 2021, 13, 13605 3 of 23 the risk score evenly. Those can cause major faults when determining precautions according to the risk score. The aims of this paper are to (1) determine the effect of the criteria on the risk score with the f-SWARA method, (2) enable OHS experts to make independent scores, (3) minimize the uncertainties and errors that may arise from the numerical scoring in the classical Fine-Kinney method with the help of linguistic terms, and (4) create hazard clusters sensitively and rank hazards accurately using the fuzzy-based COPRAS method.
To the best of our knowledge, no study has been conducted with the f-SC method in the OHS field for a risk analysis. We want to develop the classical Fine-Kinney risk analysis with an integrated f-SC method because this method is used widely in Turkey in many industries. Different from other studies, we can explain why we developed this method, and our key points are as follows: There are a few studies in the literature using AHP for the weighting of the P, E, and S criteria for Fine-Kinney risk analysis in OHS, as opposed to other analysis methods such as FTA, FMEA, and HAZOP. We suggest the SWARA method for the criteria weighting process because of some of its advantages. SWARA requires considerably lower pairwise comparisons compared to AHP for the same number of criteria. Thus, SWARA is simpler than AHP [7,8]. As the number of criteria used in the AHP method increases, the number of comparisons significantly increases. This has a negative impact on comparisons. SWARA is not similar in this regard. SWARA, which contains significantly fewer pairwise comparisons than AHP, is an effective method for the use of less prepared participants [7][8][9][10]. There are studies using methods such as TOPSIS and VIKOR for risk analysis in OHS. We suggest the COPRAS method for risk analysis because it is quite simple in terms of use. It does not rely on a high number of long calculations compared to other methods [10,11]. In the COPRAS method, calculations can be made for criteria to be maximized and minimized, and these criteria are handled and evaluated separately [10][11][12]. The COPRAS method shows the degree of the benefit of the alternatives. By comparing the alternatives with each other, we can reveal how good or how bad it is compared to the other alternatives, and this helps to rank all of the alternatives [12]. The COPRAS method provides an opportunity to evaluate both the qualitative and quantitative criteria [12].
There are many studies that have used MCDM methods in the literature. Lately, some studies using MCDM methods have also applied them to risk analyses in the OHS field. We provide some examples in the following.
Wang and Xie studied the risk of fire occurring on an offshore platform. In this case, a fuzzy-logic-based fault tree method was used. In addition, they proposed a new method by integrating Bayes networks, which is a more comprehensive method than fuzzy logic. In their case study, they stated that fuzzy logic can be integrated directly into Bayes networks, and these networks can be used to predict marginal status. They stated that this integrated method is more flexible and useful than the traditional fault tree analysis method [13]. A study was carried out by Wang et al. in 2013 on oil tank fire and explosions, which occur frequently in refineries and cause life and economic losses and environmental pollution. In this study, the events that caused oil tank explosions were determined, and a fault tree was created. However, a new fuzzy logic-based analysis method was suggested by stating uncertainties existing in the traditional fault tree analysis method due to insufficient data. According to Wang et al., the new system helps to eliminate existing uncertainties [14]. Kou and others conducted a financial risk analysis in 2014. The TOPSIS, DEA, and VIKOR MCDM methods were chosen for that analysis. They proposed this method to rank a selection of popular clustering algorithms in the domain of financial risk analysis. They used a case study to validate their proposed approach. With the obtained results, they saw the effectiveness of MCDM methods in evaluating clustering algorithms [15]. Lui et al. studied decision making mechanisms in relation to emergency responses using a fault tree Sustainability 2021, 13, 13605 4 of 23 method in 2014. Emergency scenarios were predicted by this method, which was applied to H1N1 disease [16]. Lavasani et al. conducted a risk assessment study in 2015 on the petrochemical industry, where incidents such as fires, explosions, and toxic gas releases have been observed. In their study, the fuzzy logic fault tree method was used. For the quantitative fault tree analysis method, gathering enough data is significant. Therefore, a qualitative solution was quantitated with the help of fuzzy logic [17]. Rajakarunakaran et al. conducted a study on LPG tanks at gas stations in 2015 using the fault tree risk analysis method. The fault tree created in the study was quantitated with data received from experts and then scored with the help of fuzzy logic. They stated that this study contributed to the elimination of uncertainties in that method [18]. Guneri et al. conducted a study to choose the most suitable risk assessment method with a fuzzy-based AHP method for small and medium-sized enterprises (SMEs) in 2015. They chose the Turkish woven-printed-carton labelling company as a case study and determined four criteria: scope, practicality, cost, and sensitivity. They emphasized that decision-making methods can be useful for choosing suitable methods for companies [19]. In 2016, Othman et al. conducted a study to prioritize hazards with the AHP method in HAZOP analysis. The proposed method both prioritizes and quantitatively analyzes hazards. They applied two different case studies to show the efficacy of the proposed method. They emphasized that the proposed method helps to determine and prioritize hazards in the HAZOP method [20]. Dagsuyu et al. conducted a study on existing hazards and the assessment of their risks in sterilization chambers, one of the riskiest places in the health sector, in 2016. The Failure Mode and Effect Analysis method was used. Risk analysis was then modeled with fuzzy logic and presented as a new method. In addition, a 5 × 5 method risk analysis was conducted for the same location, and the two methods were compared. In conclusion, they emphasized that the proposed method was more effective than the other method [21]. Gül and Çelik conducted a study in 2017 on railroad transportation, the most important transportation system for crowded metropolises. Fuzzy logic based on the Fine-Kinney risk analysis method was used in their study. The numbers in the classical Fine-Kinney method were replaced with verbal definitions and then quantified with the fuzzy logic method [22]. Muhammet Gül et al. developed a hybrid risk analysis in ballast tanks used in the maritime industry in 2017. The Fine-Kinney risk analysis was modelled using the fuzzy logic method with Analytical Hierarchy Process (AHP) and Vlse Kriterijumska Optimizacija I Kompromisno Resenje (VIKOR) methods for this industry. They stated that the proposed method was useful and practical and stated that it could be adapted to different sectors [3]. Another study was conducted by Ali Kokangül et al. In 2017, the Fine-Kinney risk analysis was modeled using the AHP method to identify and score categorized hazards. The aim of that study was to minimize the interpretation differences in risk analysis [23]. Yazdi and Kabir, in 2017, proposed a method based on the classical fault tree analysis method, which is widely used in risk assessments in process industries. In their study, the fault tree analysis was modelled with fuzzy logic using Bayes networks. They emphasized that classical fault tree analysis has complexity in eliminating uncertainties in hazards and that this new proposed method works better in this regard. In addition, the proposed method was applied to an ethylene production facility as an example [24]. Gül and Ak, in 2018, compared two analysis methods: risk analysis modeled with a 5 × 5 matrix using the AHP method with the help of Pythagoras fuzzy clusters and a TOPSIS method via trapeze fuzzy clusters, a multicriteria decision-making method. They obtained similar results from both methods [25]. In 2018, Fattahi and Khalilzadeh developed a new method for a workplace in the steel sector using the Failure Mode and Effect Analysis method used in the risk assessment stage in various industries. They created a hybrid method. In addition to fuzzy logic based on AHP modeling for classical FMEA, they used fuzzy logic based on MULTIMOORA modeling for three parameters: time, cost, and profitability. They evaluated the data obtained by taking the weighted mean of the results of these two analyses [26]. Gul conducted a review study that evaluated 80 different studies based on MCDM risk analysis methods in 2018. He aimed to present these studies to researchers and determine the trend topics and methods in the OHS. It was found that many studies use the AHP method with the TOPSIS and VIKOR methods compared to others in the literature [27]. Our study will contribute to the literature because it is a novel hybrid method. Markowski and Siuta studied representative accident scenarios (RASs) with both traditional HAZOP and its fuzzy-based version. A liquefied natural gas (LNG) storing facility was used as the case study. They used the efficacy index (EI) and quality index (QI) to develop the risk correction index (RCI), which was used to determine the accident scenarios. They found that the fuzzy-based method was more sensitive than the traditional method. Moreover, they mentioned that the fuzzy-based method has some advantages compared to the traditional method [28]. Chen et al. conducted a case study to reduce risk levels at the Oilfield Gathering Station (OGS) integrating AHP and fuzzy comprehensive evaluation (FCE) in 2019. An evaluation index system of safety performance at the OGS was established, and it included a tank unit index, pipe unit index, digital monitoring unit index, and other systems. Each index was weighted with the AHP method. Combining the AHP and the FCE methods, the risk levels of representative enterprise S were validated. They found that the grade of S-OGS was low-risk. They emphasized that the study helped to improve the risk levels of the OGS [29]. In 2020, Simic et al. conducted a fuzzy-MCDM-based risk analysis study on railway infrastructure. They used fuzzy sets for express information about risks, and, with the help of Tsallis-Havrda-Charvát entropy, a hybrid method was conducted. Subsequently, fuzzy measurements of alternatives and a compromise solution method were developed to rank railway infrastructures. They demonstrated the substantiality of the proposed approach with sensitivity analyses [30]. In 2020, Marhavilas et al. developed a new risk analysis method with the joint collaboration of a deterministic (DET) and a stochastic (STO) process named MCDM-STO/DET for sustainable engineering projects. They used real work accident data from between 2009 and 2016 taken from "SEPE" (Labor Inspectorate, Hellenic Ministry of Employment) and "IKA" (Social Insurance Institution, Hellenic Ministry of Health). They chose the Greek construction sector and the fixed telecommunications technical projects of OTE SA (that is, the Greek Telecommunications Organization S.A.) as a case study. They aimed to help upgrade the OHS field and the sustainability of engineering projects with that study [31]. Marhavilas et al. proposed the HAZOP method integrated with a Decision Risk-Matrix (DRMA) and fuzzy-based AHP (f-AHP) to determine potential hazards in a crude-oil processing plant in Greece in 2020. They aimed to develop and generalize this method in that industry. The classical HAZOP was used to identify fault causes in abnormal circumstances, and DRMA was used for an assessment of the risks. Subsequently, f-AHP was used to extend the study and prioritize the hazards [32]. Marhavilas et al. conducted a HAZOP risk analysis integrated with classical and fuzzy-based AHP focused on the economic, health, and environmental impact of the sour-crude-oil process industry in 2020. They aimed to show the applicability of the classical HAZOP method and the usability of hazard analysis with the proposed integrated method arising from economic, health, and environmental circumstances. They also aimed to determine the potential destructive hazards with the classical HAZOP in abnormal circumstances in related industries. They emphasized that this proposed method helps to improve OHS levels and has certain advantages for decision makers [33]. In 2021, Shahri et al. studied the FMEA risk analysis method. They proposed the PF-VIKOR method, a Pythagorean fuzzy-based MCDM method, for a dehydration unit of an Iranian gas refinery. In that study, they first used Pythagorean fuzzy numbers to collect criteria: Severity (S), Occurrence (O), and Detection (D). Afterward, they used the k-means clustering algorithm to cluster failure modes, which was a novel method accepting PFNs as input. Finally, they used the PF-VIKOR approach to analyze the ordering of cluster risks. They compared the results with the classical FMEA method [34]. Boral and Chakraborty applied a MCDM-based risk analysis method for CNC machines in 2021. FMEA was chosen for the risk analysis in that study. They evaluated the failures of the CNC machines caused by human error. First, they identified major human errors with the FMEA method during the CNC machining process, which was performed by the experts. Afterwards, a fuzzy-based method, i.e., an interval type-2 fuzzy set (IT2FS) MCDM framework, was developed for prioritizing risks. In this process, they used the AHP, DEMATEL, and MARCOS as MCDM methods. Finally, they ranked the risks and compared these three fuzzy-based methods [35]. In 2021, Tabatabaeea and others developed a fuzzy-based MCDM method for identifying and evaluating risk factors for Building Information Modelling (BIM)-based industralised building system (IBS) projects. They used the Fuzzy Delphi Method to identify and determine risk factors. Afterwards, they performed the analyses with the help of DEMATEL and the Parsimonious-Fuzzy Analytic Network Process. The method allowed users to prioritize risk factors with an understandable procedure [36]. Wang et al. conducted a study to develop a new FMEA method combined with fuzzy-based MCDM methods. They aimed to eliminate the shortcomings of classical FMEA. For this reason, the risk factors were weighted as subjective and objective with social network analysis and the maximizing consensus method. The weights of these factors were integrated with the TOPSIS MCDM method to obtain an ideal solution. Finally, they made a case study to show the effectiveness of the proposed method [37].

Research Methodologies
Within the scope of this study, a risk analysis conducted with the classical Fine-Kinney method in a workplace in the paint industry was modeled by a new hybrid method called the f-SC risk analysis method. For this purpose, the criteria used in risk score evaluation were first weighted with fuzzy logic based on the SWARA method, and the analysis was then completed using fuzzy logic based on the COPRAS method. Moreover, the f-VIKOR method used for risk analysis in the field of OHS in the literature was applied to determine the differences between these methods.

Classical Fine-Kinney Risk Analysis Method
The Fine-Kinney risk analysis method is a quantitative method developed in 1976 in order to control hazards. This classical method includes 3 criteria: Probability (P), Exposure (E), and Severity (S) [38].

Probability (P)
Although probability values should be in the range of 0-1 mathematically, a different scale was developed for probability, similar to other factors in the Fine-Kinney method. In this scale, the probability of occurrence in the workplace is evaluated between "might well be expected" and "virtually impossible" risks [3,22,38].
Exposure (E) An exposure-oriented scale was determined for the exposure value. This scale ranges from continuous exposure to very rare, which is an annual exposure [3,22,38].
Severity (S) A classification was made for an accident that may occur; the severity value is determined according to the predicted injury or death status and accordingly the financial loss table that will arise [3,22,38].
Risk Score (RS) In the classical Fine-Kinney method, the risk score is determined by the multiplication of probability, severity, and exposure values. For these determined risk scores, there are five different risk situations and a priority of precautions [3,22,38].

Fuzzy Set Theory
The fuzzy set theory was developed by Zadeh to overcome the uncertainties of problems that can occur. In this method, problems are examined with the help of linguistic terms. These linguistic terms are then converted into fuzzy numbers with the help of decision makers. Triangular fuzzy numbers are shown as A = (l, m, u). Here, "l" indicates the lowest, "m" is the medium, and "u" is the highest value, respectively. Triangular fuzzy numbers are shown in Figure 1 [22,39].
The mathematical operations ( = ve = ) of any two triangular fuzzy numbers are shown in Table 1.

SWARA Method
The SWARA method was developed in 2010 by Kersuliene et al. Today, the use of this MCDM method is becoming increasingly widespread. In this method, the role of decision makers is vital compared to others. This method provides for the opportunity to choose one's own priorities for the determined criteria for decision makers [8,10].

COPRAS Method
The COPRAS (COmplex PRoportional ASsesment) method was developed in 1996 by Zavadskas and Kaklauskas. In this method, alternatives are ranked according to their importance and benefit values. While making an alternative assessment, the method aims to maximize beneficial values and minimize useless ones. This process is based on the Membership functions of triangular fuzzy numbers are shown in Equation (2).
The mathematical operations ( ) of any two triangular fuzzy numbers are shown in Table 1. Table 1. Arithmetic operations of triangular fuzzy numbers.

SWARA Method
The SWARA method was developed in 2010 by Kersuliene et al. Today, the use of this MCDM method is becoming increasingly widespread. In this method, the role of decision makers is vital compared to others. This method provides for the opportunity to choose one's own priorities for the determined criteria for decision makers [8,10].

COPRAS Method
The COPRAS (COmplex PRoportional ASsesment) method was developed in 1996 by Zavadskas and Kaklauskas. In this method, alternatives are ranked according to their importance and benefit values. While making an alternative assessment, the method aims to maximize beneficial values and minimize useless ones. This process is based on the importance of the decision options and the direct and proportional dependence of priority. In 2007, the fuzzy COPRAS method was developed by Zavadskas and Antucheviciene. This method was used to determine alternative decisions for the renewal of abandoned buildings in the rural areas of Lithuania [11].
In the classical COPRAS method, scoring is made by giving numerical values for alternatives. This situation may cause uncertainties and inadequacies in the analyses. For this reason, when scoring, it is more appropriate to use intuitive linguistic terms.

A Novel Hybrid f-SC Risk Analysis Method
Classical risk analysis methods used in workplaces assume that the weights of all criteria are equal. However, it is not possible for all criteria to have equal importance. For this reason, the criteria must be weighted with the help of expert decision makers. In this study, a Fine-Kinney risk analysis was conducted with the new hybrid f-SC method containing two MCDM methods: SWARA and COPRAS. SWARA was integrated with CO-PRAS since there is no criterion weighting in the COPRAS method. After determining the criterion weights with SWARA, a risk analysis was conducted with classical Fine-Kinney and modeled by COPRAS, and an alternative list for hazards was created. In this modeling, both SWARA and COPRAS methods were implemented based on fuzzy logic. In this way, the criteria of hazards were quantified with the help of fuzzy logic. In Figure 2, the stages of the f-SC risk analysis method are given.
In the classical COPRAS method, scoring is made by giving numerical values for alternatives. This situation may cause uncertainties and inadequacies in the analyses. For this reason, when scoring, it is more appropriate to use intuitive linguistic terms.

A Novel Hybrid f-SC Risk Analysis Method
Classical risk analysis methods used in workplaces assume that the weights of all criteria are equal. However, it is not possible for all criteria to have equal importance. For this reason, the criteria must be weighted with the help of expert decision makers. In this study, a Fine-Kinney risk analysis was conducted with the new hybrid f-SC method containing two MCDM methods: SWARA and COPRAS. SWARA was integrated with COPRAS since there is no criterion weighting in the COPRAS method. After determining the criterion weights with SWARA, a risk analysis was conducted with classical Fine-Kinney and modeled by COPRAS, and an alternative list for hazards was created. In this modeling, both SWARA and COPRAS methods were implemented based on fuzzy logic. In this way, the criteria of hazards were quantified with the help of fuzzy logic. In Figure  2, the stages of the f-SC risk analysis method are given. Our novel hybrid method consists of three main stages as seen in Figure 2. The steps of this method are as follows.
1st Stage: Determination of Criteria (Classical Fine-Kinney method) The classical Fine-Kinney method was used for determining the criteria for this study. In this method, the risk score (RS) is evaluated by multiplying Probability, Exposure, and Severity values [23,38]: In the Fine-Kinney method, separate tables have been defined by Kinney to determine Probability (P), Severity (S), and Exposure (E) values. The risk score is calculated by Our novel hybrid method consists of three main stages as seen in Figure 2. The steps of this method are as follows.
1st Stage: Determination of Criteria (Classical Fine-Kinney method) The classical Fine-Kinney method was used for determining the criteria for this study. In this method, the risk score (RS) is evaluated by multiplying Probability, Exposure, and Severity values [23,38]: In the Fine-Kinney method, separate tables have been defined by Kinney to determine Probability (P), Severity (S), and Exposure (E) values. The risk score is calculated by determining the numerical values corresponding to the statements in these tables separately for each parameter [38].
Probability, Exposure, and Severity values and the different risk situations arising from the Risk Score are shown in Table 2. The f-SWARA method was used for the weighting of the criteria (P, E, and S). The steps of the SWARA method are as follows [9,10]. 1st Step: Criteria are ranked from high to low according to the opinions of the decision makers. If there are multiple decision-makers, the ranking is made separately and a general ranking is determined by taking the geometric mean. 2nd Step: The relative significance of each criterion is determined. For comparing the j criterion with the (j + 1) criterion, how important the j criterion is in relation to the (j + 1) criterion is determined. This value (s j ) is called the "comparative significance of the mean value". 3rd Step: The k j coefficient is determined: 4th Step: The q j variable is determined: 5th Step: The relative weights of the evaluation criteria are determined as follows, where w j refers to the relative weight of the j criterion: In this study, the SWARA method was based on fuzzy logic. For this reason, decision makers expressed the importance of the criteria using the terms in Table 3. Afterwards, defuzzification was performed using fuzzy numbers corresponding to the linguistic terms.  (7,8,9) 3rd Stage: Analysis of the weighted criteria and hazards/alternatives (f-COPRAS method).
In the fuzzy COPRAS method, alternatives are evaluated via linguistic terms by decision makers and quantified with the help of fuzzy logic. In this way, uncertainties disappear [11]. The steps of this method are as follows. 1st Step: Definition of the linguistic terms. In this step, alternatives are determined based on fuzzy logic with linguistic terms according to the criteria. Linguistic terms based on fuzzy logic are shown in Table 4. Table 4. Linguistic terms for the hazards/alternatives.

Linguistic Terms Fuzzy Numbers
Very poor (VP) (0, 0, 1) Poor (P) (0, 1, 3) Medium poor (MP) (1, 3, 5) Fair (F) (3, 5, 7) Medium good (MG) (5,7,9) Good (G) (7, 9, 10) Very good (VG) (9,10,10) 2nd Step: Creation of the fuzzy decision matrix. Here, a fuzzy-based matrix with n criteria and m alternatives is created. This matrix, shown in Equation (6), is created by decision makers. 3rd Step: Defuzzification of the fuzzy decision matrix. The graded mean integration representation method was used for the defuzzification of the fuzzy decision matrix. This method is simple and useful [40]. The best non-fuzzy performance value for the fuzzy number A = (l i , m i , u i ) is calculated as follows: 4th Step: Normalization of the defuzzified decision matrix. The normalized value of this matrix is calculated with "w" data integrated into this equation with the help of f-SWARA (Equation (5)), as follows: After the calculation, the decision matrix is normalized as follows: 5th Step: Calculation of the weighted normalized decision matrix. Weighted normalized values are calculated as follows: After this process, the weighted normalized decision matrix is as follows: 6th Step: For each alternative, the sum of the attribute values (P i ) where large values are preferred is calculated: Here, k is the attribute value that needs to be maximized. 7th Step: For each alternative, the sum (R i ) of the attribute values in which small values are preferred is calculated: Here, (m-k) is the attribute value that needs to be minimized. 8th Step: The minimum R i value is determined: 9th Step: The relative weight of each alternative is calculated: 10th Step: K, the optimality criterion, is determined: 11th Step: Determination of the priority of the project. The alternative with the highest value is the best option for the project.
With this proposed method, we developed the classical Fine-Kinney method using the f-SWARA and f-COPRAS MCDM methods. We basically used the same equation (Equation (1)). First, the criteria of this equation (P, E, and S) were weighted using Equations (3)-(5) by 10 OHS experts. During the weighting process, all experts were asked to rank the criteria from effective to ineffective according to the f-SWARA method and to express the degree of the influence of each criterion in linguistic terms. All the experts who gave scores ranked the effect of the criteria identically. The criteria are listed from the most effective to the least effective in terms of severity, exposure, and probability. Table 5 contains the expressions used by the experts regarding the level of the influence for each criterion according to the linguistic terms in Table 3. Here, the effect of the severity (at the top) on the exposure and the effect of the exposure on the probability (at the bottom) are evaluated. Since there are no criteria on the effect on severity, there are no input values for this criteria's row. The criteria weighting calculations of the 1st decision maker are given in Table 6 as an example. Decision makers performed this process independently of one other. Afterwards, the geometric mean of the experts' data was used as the final value for the weighting value (w). After the weighting calculations for the criteria were determined by the decision makers with the help of Equations (3)-(5), the final weighting process was performed by taking the geometric mean of all values. The criterion weight values of all decision makers and the final values found with the geometric mean are shown in Table 7.
While evaluating the criteria, all decision makers chose the severity criterion as the most significant criterion with a common decision. By calculating the geometric mean of the values obtained as a result of the calculations made separately for each decision maker, the importance weights of the severity, exposure, and probability criteria were obtained as 0.518, 0.285, and 0.196, respectively. These obtained values were integrated into the f-COPRAS method, and the relevant weight ratio for each criterion was taken into account.
Subsequently, with the help of the f-COPRAS method, the hazards identified by the experts as suitable for the ISO 31000:2018 Risk Management Guidelines and the Risk Assessment for OHS Regulations depending on the OHS law in Turkey were scored using the data in Table 4, which contains the linguistic terms for valuing the decision matrix by the OHS experts. After this step, "w" data calculated by f-SWARA were integrated into f-COPRAS when normalizing the decision matrix with the equation  (17)) of the f-COPRAS method, N i values of each hazards were determined and ranked from low to high. According to the ISO 31000:2018 Risk Management Guidelines, "The purpose of risk identification is to find, recognize and describe risks that might help or prevent an organization [from] achieving its objectives. Relevant, appropriate and up-to-date information is important in identifying risks." Due to this, identifying hazards and risks is significant in the workplace. The risk identification process is made according to the risk score in the classical Fine-Kinney method. Using numbers in this method causes problems and uncertainties. The proposed model, which was developed according to IEC 31010:2019 Risk Management-Risk Assessment Techniques, helped to eliminate these shortcomings. In this method, decision makers used the intuitive linguistic terms instead of numbers. Figure 2 shows the stages of our proposed method corresponding to the relevant standards and regulations.

A Case Study for the f-SC Method in the Paint Industry
In this chapter, a Fine-Kinney risk analysis is implemented with the new hybrid f-SC method in one of the major businesses in the paint industry in Turkey. The paint industry is one of the sectors where chemicals are used extensively. For this reason, the measures taken in paint-related workplaces are classified as very dangerous and are essential in terms of the OHS. Therefore, necessary measures to prevent hazards or to minimize damage were proposed as a result of this risk analysis conducted with the new method. The facility where the application was conducted is very large, so the scope of this study was limited to the water-based paint production section and the area where the silos are located. In this context, 42 hazards were determined for this area.
In the conducted case study, first, the probability, exposure, and severity values determined by the classical Fine-Kinney method were expressed with linguistic terms by 10 decision makers who served as safety experts. In Table 8, the classical risk analysis conducted according to Equation (2) is given. The alternatives expressed with linguistic terms by the decision makers according to Table 4 are shown in Table 9 for the f-SC approach.   H15 H16 H17 H18 H19 H20 H21 H22 H23 H24 H25 H26 H27 After the criteria (P, E, S) for each hazard were expressed with linguistic terms according to Table 4, the fuzzy decision matrix was created using fuzzy numbers corresponding to the linguistic terms. The graded mean integration representation method given in Equation (7) was used in the created matrix, and the best non-fuzzy performance values for each criterion were calculated. They are shown in Table 10. Obtained defuzzified values were normalized using the f-SWARA method, and the criteria were weighted by 10 OHS experts. According to this process, possibility, exposure, and severity values were calculated as 0.196, 0.284, and 0.518, respectively, as shown in Table 7 in the previous section. After the calculations of the weights of the criteria with the f-SWARA method, by using Equation (10), the normalized values of alternatives were obtained. These values are shown in Table 11. Table 11. Normalized values of hazards/alternatives. After the normalization of alternatives, P i and R i values were determined with the help of Equations (12) and (13). Since there was no benefit criterion, the P i value was zero for all the alternatives. According to the calculated R i values, the lowest score was the H11 alternative with R min = 0.016. After the calculations, values were obtained. The ultimate values calculated using Equations (14)-(17) are given in Table 12. Furthermore, the rankings of the alternatives are given in Table 12. The alternatives are ranked from the most dangerous to the least dangerous in this table. Therefore, H11 with an N i value of 1.000 is in the 42nd row, and H7 with an N i value of 0.557 is in the 1st row.    The N i value for every alternative was calculated by COPRAS analysis. These values are given in Table 12 and Figure 3. In the calculated values, H11 with a N i value of 1 indicates the least dangerous situation. On the contrary, H7 with a N i value of 0.557 indicates the most dangerous situation. The priority for taking measures was determined for the alternatives/hazards between these values. The most dangerous hazard that should be dealt with immediately (H7) was in Cluster 1 only. The values in Table 13 were used for the determination of clusters, except for that containing H7. Clusters of the alternatives are shown in Figure 4.
Determining the precautions is an important step after the classification of hazards. Measures to be taken will be significantly helpful in eliminating or reducing risk. The suggested precautions for the hazards addressed in this study are given in Table 14   Determining the precautions is an important step after the classification of hazards. Measures to be taken will be significantly helpful in eliminating or reducing risk. The suggested precautions for the hazards addressed in this study are given in Table 14     As a result of the analysis performed with the f-SC method, the hazards are ranked in Table 12. A comparison of the hazard rankings obtained as a result of the classical Fine-Kinney method analysis, the f-VIKOR method analysis, and the f-SC method analysis is shown in Figure 5. In addition, the classes of hazards were determined for the f-SC, the classical Fine-Kinney, and f-VIKOR methods, and a comparison is shown in Figure 6.
supplied.  Materials on the roof must be fixed.


The weight of the hoses must not exceed 15 kg. During insertion and removal, the pipes must be transported by pulling them.  Maintenance and repair activities must be carried out periodically by the maintenance and repair team.  Employees must use helmets during work. As a result of the analysis performed with the f-SC method, the hazards are ranked in Table 12. A comparison of the hazard rankings obtained as a result of the classical Fine-Kinney method analysis, the f-VIKOR method analysis, and the f-SC method analysis is shown in Figure 5. In addition, the classes of hazards were determined for the f-SC, the classical Fine-Kinney, and f-VIKOR methods, and a comparison is shown in Figure 6. Figure 5 shows that, for all methods, the highest scored risk was Hazard 7 (H7). Moreover, the differences in the ranking of hazards can be seen in Figure 5.
When the hazard groups were examined, as shown in Figure 6, it can be seen that some hazard groups changed. This stage of the analysis was quite significant since it determined the timetable for the precautions.  Curved mirrors providing wide visibility must be used, and field of vision lighting must be used. Workers in the vicinity must wear reflective work clothes.

I
After the camlock clamps are installed, the necessary controls must be made before starting work.

I
A horizontal lifeline must be established outside. The warning sign and instructions must be hung, and whether the employees are following the instructions must be checked.

I
Employees must use personal protective equipment (PPE) suitable for the job, and training on the use of PPE must be provided to employees.
I MSDS forms of the chemicals used must be supplied to the employees, and the rules in these forms must be followed.

I
Employees must wear dust masks, and a ventilation system must be available.

I
Material packages must not exceed 25 kg. A forklift or pallet truck must be used during transportation.

I
Workers must use appropriate ear protectors, and regular ambient measurements must be performed.

I
Valves must be checked regularly, and malfunctioning valves must be reported to the maintainer.

I
The machine must have body grounding. In addition, the boilers must be grounded before the operation.

I
Drains must be checked, and appropriate signs and plates must be used.

I
Appropriate cleaning kits and PPE must be used when working in the area.

I
The working area must be tidy. The machine must not be used by anyone other than the operator.

I
Cleaning must not be performed while the mixer is running and washing tools must be used for cleaning. In addition, employees must be trained on boiler cleaning.

I
The boiler stretching operation must be performed slowly and carefully.

I
The materials must be placed according to the stacking rules. Hazardous materials must be transported closed and carefully. Pallets must be checked, and damaged and broken pallets must not be used.

I
Electric panels must be protected in accordance with the legislation.
I Pneumatic pumps must be fitted with a silencer. According to the noise measurement results, appropriate PPE must be used in the required areas.

I
The maximum speed limit within the facility must be 20 km/h, and all suppliers must be informed. In addition, vehicles entering the facility must be in a certain order.

I
The grounding of the vehicles must be performed properly, and periodic checks must be provided.

I
Employees must be given work clothes suitable for seasonal conditions. Cordless hand tools must be used.

I
Boilers must be used with wheels.

I
The filling line must have a suitable ventilation system, and the filling process must be performed with hoses. In addition, appropriate PPE must be used to prevent exposure to workers.

I
A scissor pallet truck must be used for lifting heavy materials.

I
A safety utility knife must be used in all areas.

I
The mixer must be stopped during the sampling process, and the bucket must be used during this process.
H10, H17, H31 IV I During the procedure, protective glasses must be used, work clothes must be worn, and hygiene rules must be followed.

I
All vehicles entering the facility must have an exhaust protector (flame/spark arrester). Unoperated vehicles must not be left running.

I
The boiler must be fixed with a seatbelt, and there must be an emergency stop button to cut energy in an emergency.

I
Pressure gauges and alarms in silos must be checked regularly. In addition, periodic inspections of silobas vehicles must be performed.

I
The driver must not be allowed to exit before the damper is closed. Periodic inspections must be performed.

Inspection reports
showing that the necessary periodic controls of the silos have been made must be supplied.

I
Materials on the roof must be fixed.
H3, H11, H28 VI I The weight of the hoses must not exceed 15 kg. During insertion and removal, the pipes must be transported by pulling them.

I
Maintenance and repair activities must be carried out periodically by the maintenance and repair team.

I
Employees must use helmets during work. Figure 5 shows that, for all methods, the highest scored risk was Hazard 7 (H7). Moreover, the differences in the ranking of hazards can be seen in Figure 5.
When the hazard groups were examined, as shown in Figure 6, it can be seen that some hazard groups changed. This stage of the analysis was quite significant since it determined the timetable for the precautions.

Discussion
It is essential to determine the activities to be implemented in order to ensure the health and safety of the employees in the workplace. While determining these activities, it is necessary to analyze the working environments in detail, and the hazards that may pose risks in these environments should be identified accurately and completely. For this reason, the chosen risk analysis methods should be appropriate for the work.
Traditional risk models are based on probability and classical set theory. In these methods, risk analyses are made using numerical data. However, it may not always be easy to obtain numerical data in the workplace. In such cases, existing risks will need to be evaluated based on the opinions and judgments of experts (decision makers) instead of numerical data. Thus, fuzzy logic models can be used where assessments can be made with insufficient data. While hazards are prioritized as in classical risk analysis methods, a controversial result that may lead to different interpretations can be obtained. For instance, high probability and low severity hazards can be classified at the same level as low-probability and high-severity hazards. These problems can be eliminated by using fuzzy logic-based MCDM risk analysis methods. In this study, the novel hybrid f-SC risk analysis method was applied in a workplace, i.e., a water-based paint manufacturing area in the paint industry. A classical Fine-Kinney risk analysis was carried out for this particular workplace.
In classical methods, a risk score obtained with numerical values is given to each criterion when performing risk analysis. If applied to such cases, the impact rates of all criteria are assumed to be equal. However, not all criteria have the same level of importance. Hence, in this study, before the risk analysis was performed with the f-COPRAS method, with the help of the f-SWARA method, 10 OHS experts were used as decision makers to weight the criteria. As a result of the study, the weights of probability, exposure, and severity criteria were determined as 0.196, 0.285, and 0.518, respectively.
Forty-two hazards were used as an alternative for the analysis and the criteria (P, E, and S) were weighted with the help of the f-SWARA method. N i values were then determined for each alternative using the relevant calculations of the f-COPRAS method. Since none of the criteria we used in our study were benefit criteria, this value was evaluated as "0". The f-COPRAS analysis makes a benefit-oriented ranking while listing the alternatives. In other words, the alternative with an N i value of 1, in the range of 0-1, is most useful. In this context, a smaller N i value corresponding to the risk score calculated for each alternative/hazard by the integrated method means a greater risk. As a result of the study, H7 corresponded to the most dangerous situation since it had the smallest N i value (0.557). On the other hand, H11, which had a value of 1, was the least dangerous situation.
After the risk analysis was completed, timetables for the necessary precautions were identified according to the classification, and the necessary precautions are given in Table 14 according to related standards and regulations. Health and safety precautions should be classified carefully and conducted on time; otherwise, great losses are likely to occur.
As a result of the analysis made within the scope of the study, different classifications were obtained for the classical Fine-Kinney method and the novel hybrid f-SC method combining the f-SWARA and f-COPRAS MCDM methods, since the criteria, without any weighting, do not have sensitivity and certainty in the classical Fine-Kinney method. With the new hybrid f-SC method, using MCDM methods based on fuzzy logic, the criteria were weighted, and the sensitivity was increased. Thus, the uncertainties of the classical Fine-Kinney method were minimized. There were differences between the risk classes arising from hazards in these two methods. Moreover, analyses were conducted with the f-VIKOR method, which has been used in previous studies for risk analysis in the field of OHS using existing data. In the analyses, the same criteria weighting values of the f-SC method were used. Using the results obtained from these analyses, a triple comparison was made. The results of the f-VIKOR method, which have been accepted in the literature, and the new hybrid f-SC method of this study were parallel to each other.

Conclusions
A decision-making process is needed for the determination of risks in the workplace. For this reason, our new hybrid f-SC multi-criteria decision-making method is suitable. Using MCDM methods based on fuzzy logic is more suitable since, compared to numerical values, linguistic terms are easier to understand. Thus, in this paper, considerable efforts were made to eliminate the uncertainties and lack of sensitivity of classical methods.
In this study, we used 42 hazards of a particular department in a chosen workplace for our proposed risk analysis method. It was chosen and limited in this way because increasing the number of alternatives causes additional mathematical operations for decision makers and causes difficulties in presenting the data. However, a risk analysis should be made in detail for any company to take significant precautions. In future work, the f-SC method can be used for risk analyses in other industries. Alternatives can be multiplied using different MCDM methods for risk analyses, and comparisons can thus be made with different methods for criterion weighting. In addition, an application can be developed for the novel f-SC method for use by OHS experts in various workplaces.
The results showed that all the criteria had different levels of importance in terms of the risk score, and changes in risk prioritizing and risk clusters can be made by weighting the criteria. It is easier for decision makers to score criteria when linguistic terms are used instead of numerical values. The case study made in this study indicated that using our proposed approach instead of the conventional method increases accuracy and sensitivity in risk analyses. Furthermore, similar results were obtained by applying the f-VIKOR method, a method that has been used in the literature, alongside our proposed approach, which demonstrates the validity of our study.
This paper provides four critical contributions to classical risk analysis methods: (1) unlike the traditional Fine-Kinney method, OHS experts weighted the criteria before the risk analysis in the proposed f-SC method. This changed the risk scores and their risk clusters; (2) OHS experts, independent from each other, considered linguistic terms, which made the valuing of the criteria more consistent and reliable; (3) By using the proposed method, hazards were analyzed and ranked accurately and with sensitivity; (4) The classical Fine-Kinney method was chosen in this study since it is widely used in the vast majority of industries in Turkey. The f-SC method eliminates the shortcomings and disadvantages of the classical Fine-Kinney method. Institutional Review Board Statement: Ethical review and approval were waived for this study because it was not required. The data is not considered to be sensitive or confidential in nature. Vulnerable or dependent groups were not included. The subject matter was limited to topics that are strictly within the professional competence of the participants.