**Panagiotis K. Marhavilas 1,\*, Michael G. Tegas 2, Georgios K. Koulinas <sup>1</sup> and Dimitrios E. Koulouriotis <sup>1</sup>**


Received: 29 April 2020; Accepted: 21 May 2020; Published: 23 May 2020

**Abstract:** This study, on the one hand, develops a newfangled risk assessment and analysis (RAA) methodological approach (the MCDM-STO/DET one) for sustainable engineering projects by the amalgamation of a multicriteria decision-making (MCDM) process with the joint-collaboration of a deterministic (DET) and a stochastic (STO) process. On the other hand, proceeds to the application of MCDM-STO/DET at the workplaces of the Greek construction sector and also of the fixed-telecommunications technical projects of OTE SA (that is, the Greek Telecommunications Organization S.A.) by means of real accident data coming from two official State databases, namely of "SEPE" (Labor Inspectorate, Hellenic Ministry of Employment) and of "IKA" (Social Insurance Institution, Hellenic Ministry of Health), all the way through the period of the years2009–2016.Consequently, the article's objectives are the following: (i) The implementation and execution of the joint MCDM-STO/DET framework, and (ii) to make known that the proposed MCDM-STO/DET algorithm can be a precious method for safety managers (and/or decision-makers) to ameliorate occupational safety and health (OSH) and to endorse the sustainable operation of technical or engineering projects as well. Mainly, we mingle two different configurations of the MCDM method, initially the Analytical Hierarchy-Process (the typical-AHP), and afterwards the Fuzzy-Extended AHP (the FEAHP) one, along with the Proportional Risk Assessment Technique (PRAT) and the analysis of Time-Series Processes (TSP), and finally with the Fault-Tree Analysis (FTA).

**Keywords:** MCDM; AHP; fuzzy AHP; PRAT; FTA; time-series; risk assessment; OSH; telecom technical projects; constructions

## **1. Introduction**

Taking into consideration that one of the intentions of "sustainable engineering" is to create structures that validate evolution without jeopardizing the natural environment's quality, and also the coming generations' capability to meet their own necessities, modern engineering tactics concentrate (among other options) on the systems' design/operation in a manner that admits the sustainable usage of resources. Thus, decision-making processes have to be enhanced by techniques and methods which permit the decision-makers to use a broad variety of sustainable alternatives.

The protection of the employees' health and safety in the workplace is a universally vital aspect, in parallel with the effort of enhancing the productivity and upgrading the sustainable development of a firm. Moreover, the World Health Organization (WHO) defines the sustainable development as a

policy *"to meet the needs of today's global population without adversely a*ff*ecting its health and the environment and without destroying (or endangering) the global resource base, and thus without jeopardizing the ability of the coming generations to meet their needs"* [1–3].

Occupational accidents, on the one hand, have a significant impact upon human body-integrity, and affect the employees, business operation, and overall the sustainability performance of firms, and on the other hand, generate high expenses for the social health and insurance system of any country, and demote the society's sustainability. Each corporation is progressively more concerned with the improvement of sustainability and occupational safety and health (OSH) performance, and this is attained by controlling OSH risks, in accordance with their OSH policy, and in the context of applying risk assessment and analysis (RAA) techniques and the state's legislation. Besides, the Occupational Safety & Health Administration (OSHA) recognizes that actions to ensure occupational safety and health can be incorporated into sustainability efforts and take advantage of the dynamics of sustainability to make companies' workplaces safer and healthier [3–5].

In this sense, the conception of sustainability is used as a frame to categorize existing and/or develop new OSH standards (e.g., ISO 45001, OSHAS 18001, etc.), and also to generate multitudinous efficacious and modern RAA techniques [5].

In general, the variety of RAA assessment methods is such that there are numerous proper approaches for most occasions, and plentiful studies build up hybrid models merging separate techniques. For example, RAA techniques are classified, on the one hand, into quantitative/qualitative/ hybrid techniques, and on the other side, into deterministic and stochastic ones.

However, the examination of the scientific literature with reference to RAA methodologies unveils that every technique presents specific limitations or restrictions in its applicability to miscellaneous accidents, and consequently a separate (or solitary) technique: (i) Can't generate either a realistic risk prediction model or a sufficient RAA process in the worksites, and (ii) can't achieve the optimum risk-assessment result in the workplaces, but on the other hand, future prospects must focus: (i) On the joint estimation/forecasting process, and (ii) on the parallel utilization of a deterministic approach with a stochastic one [6–10].

The above topics constitute the research context of this study, and consequently, its foremost contribution is the under mentioned:(i) The implementation of an alternative RAA methodological framework for upgrading the OSH's situation and sustainability of engineering projects by the combination of a multicriteria decision-making (MCDM) process with the joint collaboration of a stochastic (STO) and a deterministic (DET) process, and (ii) the application of it on the worksites of the Greek construction sector and of fixed-telecommunications technical projects of OTE SA (the Greek Telecommunications Organization SA), by means of real accident data.

The paper is organized into the subsequent sections: (1) Introduction, (2) literature survey, (3) development of the suggested methodological framework, (4) application of the new method on the constructions and telecom sector and analysis of their results, and (5) discussion of the main points and conclusions.

#### **2. Literature Survey**

Plenty of research studies have been conducted in copious construction industries contributing to the scientific field of OSH by revealing the features (advantages and disadvantages) and also presenting the implementation of numerous RAA techniques. Below, we conduct a survey of the scientific literature concerning RAA techniques, which are used in the OSH field along with the constructions sector, and which are concentrated, on the other hand, on several deterministic (DET), stochastic (STO), and MCDM processes (like PRAT, FTA, TSP, AHP, FEAHP, etc.), that have been incorporated in the newly proposed integrated-method. More explicitly, our study develops a novel method which incorporates three different types of techniques (DET, STO, and MCDM), so the presented literature survey covers these three special classes of techniques.

On the whole, RAA techniques are categorized into three significant groups (e.g., quantitative/ qualitative/hybrid techniques) whereas another methodical categorization comprises the deterministic and the stochastic approaches [8–17].

More specifically, Marhavilas et al. (2011) [8], based on the literature's examination, determined the dominant RAA methods, and categorized them into three principal classes, i.e., the qualitative, the quantitative, and the hybrid (semi-quantitative or qualitative–quantitative) methods, wherein the first ones are grounded mainly on the safety-managers' capability, the second ones express the risk by a mathematical relation in conjunction with real accident data recorded in a workplace, and finally, the third ones are characterized by a high complexity due to their ad-hoc feature. A plethora of other works confirm the results of this study, and also the importance of its classification, like the articles of (i) da Costa et al. (2020) [18], Gul and Celik (2018) [19], AriamuthuVenkidasalapathy et al. (2018) [20], Ozkan and Uluta (2017) [21] (*for quantitative RAA*), (ii) Athar et al. (2019) [22], Sanmiquel-Pera et al. (2019) [23], Domínguez et al. (2019) [24] (*for qualitative RAA*), and (iii) Kharzi et al. (2020) [25], Cinar and Cebi (2020) [26], Mutlu et al. (2019) [27], Bora et al., (2019) [28], Kamsu-Foguem and Tiako (2017) [29], Zheng et al. (2017) [30], Mentes and Ozen (2015) [31] (for hybrid RAA).

One of the most practical and skillful numerical formula for evaluating the quantified risk (R) at workplaces is the proportional risk-assessment DET-technique (PRAT) which takes into consideration the severity (S) of a harm (i.e., the potential consequences of an accident), the occurrence likelihood (P) i.e., the probability factor of that harm, and the exposure factor [32–37]. This concept becomes operational, and hence many studies, like the papers of Gul and Celik (2018), Gul et al. (2018), Kokangül et al. (2017), Marhavilas (2015), Marhavilas et al. (2011b), Marhavilas and Koulouriotis (2008) [19,38–42] combined exposure-factors with likelihood-factors to arrive at the frequency (F), and to evaluate, on the one hand, the risk by

$$R = P \cdot S \cdot F \tag{1}$$

and consequently, unsafe OSH situations at worksites by means of exposure, which estimates the hazard's occurrence frequency.

Continuing on, fault-tree analysis (FTA) and event-tree analysis (ETA) are, according to Marhavilas et al. (2014a, 2014b) [43,44], well known DET-methods and valuable analytic tools worldwide, which have been used successfully by reliability experts in failure-analysis and for the reliability and safety of complex technical-systems, and also of OSH-systems. Nowadays, many studies continue incorporating FTAs/ETAs in OSH RAA, like: (i) The work of Mutlu et al. (2019) [27], which depicts a OSH RAA approach in the textile industry by integrating FTA with Failure-Mode and Effect-Analysis (FMEA) and BIFPET, (ii) the articles of Babaei et al. (2018) [45] and Yasli and Bolat (2018) [46], which include FTA in fuzzy OSH RAA methods, (iii) the paper of Gul and Ak (2018) [47], which illustrates the collaboration of MCDM with FTA, (iv) the work of Fuentes-Bargues et al. (2017) [48], which combines Hazard and Operability (HAZOP) and FTA in risk analysis of fuel storage-terminals, and (v) the paper of Marhavilas et al. (2014a) [43], which physically embody FTA in DET OSH-RAA techniques.

Over and above that, the stochastic (STO) processes, like the Time-Series Processes (TSP), Markov-chains, etc., play an important role in the recognition of potential hazardous sources and also in the methodologies of accident analysis, RAA, prediction and modeling, and additionally their contribution in the analysis of accident tendency, and to the forecasting of the safety's situation based on the real accident observations and/or data, is vital [7,17,49–53].

Moreover, there are many studies, in the relative scientific literature, that develop MCDM methods by hybridizing existing, very popular approaches, to improve results in risk assessments. Zheng et al. (2012) [54] employed a variation of fuzzy Analytical Hierarchy-Process (AHP) which uses trapezoidal fuzzy numbers, for risk evaluation under hot and humid environment conditions, and Fattahi and Khalilzadeh (2018) [55] introduced a framework for assessing risks in steel industry using Failure-Mode and Effect-Analysis (FMEA), fuzzy AHP, and Fuzzy Multi-Objective Optimization on the basis of Ratio Analysis (MULTIMOORA) methods. In the study of Mete (2019) [56], a hybrid fuzzy analytical

hierarchy process and data envelopment analysis-based framework were proposed for managing uncertainty, while assessing risks in pipeline construction projects. Recently, Oturakci (2019) [57] introduced a novel fuzzy sets and Analytical Hierarchy Process based approach for environmental risks evaluation, while Yucesan and Kahraman (2019) [58] applied the Pythagorean fuzzy analytical hierarchy process method for assisting a team of experts to evaluate risks in hydroelectric power plants. In the study of Boral et al. (2020) [59], a new approach was developed which combines the Fuzzy Analytical Hierarchy Process with the modified FMAIRCA (Fuzzy Multi-Attribute Ideal Real Comparative Analysis) in order to improve the disadvantages of the FMEA technique, used for proactive risk assessment.

It is worth mentioning that there are numerous other published papers, pertaining to OSH management systems and to RAA techniques, and focusing on construction sites. Therefore, Marhavilas (2009, 2015) [40,60], Ardeshir et al. (2014) [61], Anil Kumar et al. (2015) [62], Marhavilas and Vrountas (2018) [63], and Koulinas et al. (2019a; 2019b) [3,64] presented the implementation of various OSH RAA techniques at workplaces of the constructions section.

#### **3. The Suggested Methodological Framework**

Figure 1 displays the flow sheet of the suggested RAA generic algorithmic framework, which couples a multicriteria decision-making (MCDM) process with both a stochastic (STO) and a deterministic (DET) process. More especially, we join the MCDM-AHP/Fuzzy-Extended AHP (FEAHP) techniques with the DET-PRAT and DET-FTA processes, and also with the STO-TSP method.

Consistent with this diagram, the real sequential data A(t) (for instance accident data) recorded at any workplace can be utilized firstly in the phase of hazard-identification, and subsequently in the DET-PRAT process for quantitatively calculating the risk as a function of time R(t). Moreover, the DET-TSP process can be applied in parallel with DET-PRAT process in order to achieve a joint-evaluation of their results. Particularly, when the time-series data A(t) are sufficient (according to Nyquist-Shannon's sampling-theorem) with reference to a specific time-period, the generated by PRAT time-series R(t) could make-up inputs for TSP, which means they could be analyzed by miscellaneous substantial methods (like the moving-average one, the seasonal-adjustment, the auto-regression, the function-fitting one, etc.), unveiling in many cases significant features for the time-series variation of R(t) and also for the company's RAA situation (such as a long-term trend, a periodicity, etc.).

It is worth noting that the combination of PRAT and TSP techniques can give, by the calculation of risk (R), the essential information to safety-managers in order to (i) identify (taking into consideration the calculated risk level) the most important sources of danger at the companies' worksites, and also (ii) to achieve a preliminary ranking of them. Thus, the resulted ranking (1st ranking) can be precious for evaluating the urgency of investing a restricted budget in specific measures in order to attain the highest OSH-protection with the lowest cost.

To continue, the next segment in the above referred flow sheet is the MCDM one, which represents the application of a multicriteria decision-making technique, like the AHP and/or the FEAHP one, for accomplishing the final ranking (2md ranking) in the sources of danger defined at worksites, in order to achieve the finest (optimum) distribution of an investment budget (i.e., with the minimum cost) to appropriate measures reaching the best OSH-protection.

Ultimately, the proposed RAA generic algorithmic framework is completed by the application of the FTA process in order to: (i) Graphically illustrate possible occurrences, which can result in adverse events (the "base-events"), (ii) relate their sequences, which could lead to a "top-event", and (iii) to calculate the occurrence likelihood of the "top-event". Therefore, the occurrence-probability of the top-event is calculated by using the following equations:

$$P\left(\bigcup\_{i=1}^{n} B\_{i}\right) = \sum\_{i=1}^{n} P(B\_{i}) - \sum\_{i=1}^{n-1} \sum\_{j=i+1}^{n} P\left(B\_{i} \cap B\_{j}\right) + \sum\_{i=1}^{n-2} \sum\_{j=i+1}^{n-1} \sum\_{k=j+1}^{n} P\left(B\_{i} \cap B\_{j} \cap B\_{k}\right) - \dots + (-1)^{n-1} P(B\_{1} \cap B\_{2} \cap \dots \cap B\_{n}) \tag{2}$$

or

$$P\left(\bigcup\_{i=1}^{n} B\_i\right) = 1 - \prod\_{i=1}^{n} \left[1 - P(B\_i)\right] \tag{3}$$

and

$$P\left(\bigcap\_{i=1}^{n} B\_i\right) = \prod\_{i=1}^{n} \left[P(B\_i)\right] \tag{4}$$

where *Bi* (*i* = 1, ... , *n*) are the identified (by FTA) "base-events", and taking into account that the first (or equivalently the second) relation is used when the output-event is logically-valid (according to Boolean algebra) if one of the input-events applies (i.e., fulfilling an OR-Gate logical-operation), while the third relation is used when the output-event is logically-valid only if all input-events apply (i.e., through an AND-Gate Boolean-function).

**Figure 1.** The joint stochastic/deterministic with multi criteria decision-making (MCDM) risk-assessment methodological framework.

#### **4. Application of the Suggested Methodological Framework on the Constructions Sector**

We proceed to the application of the novel methodological framework on the worksites of the Greek construction-sector(C-S)and of fixed-telecommunications technical-projects of OTE SA (the Greek Telecommunications Organization SA), by using real accident data recorded in two official State data sources, namely of "SEPE" (Labor Inspectorate, Hellenic/Greek Ministry of Employment) and of "IKA" (Social Insurance Institution, Hellenic/Greek Ministry of Health) with reference to the years 2009–2016.

The hazard identification (for instance, by HAZOP) and the risk analysis of the sources of danger (or deviations) existing at the worksites of the Greek constructions section are unveiled in Table 1, which depicts the deviations according to the classification of ESAW (European Statistics on Accidents at Work) by EU (2013) [65]. More specifically, the column "c" illustrates the description of the various hazards (or deviations) determined by the usage of "SEPE" and "IKA" statistical databases, while columns "b" and "a" showing the ESAW codes of the deviations, and the general types of the deviations (1st level of categorization), respectively.

**Table 1.** The identified (throughout the period of2009–2016) deviations at the workplaces of the Greek constructions section (C-S) and of fixed-telecommunications technical-projects of OTE SA (the Greek Telecommunications Organization SA), in accordance with the classification of ESAW (European Statistics on Accidents at Work) by EU (2013) [65].


#### *4.1. Application of the PRAT and TSP Processes*

The application of PRAT process was achieved by having calculated the risk (with its implementation equation 1), regarding all the hazard sources of Table 1, and using accident data from the above referenced databases ("SEPE" and "IKA") throughout the period of 2009–2016.

In Figure 2, we display time-variations of the magnitude of risk (R) calculated by the PRAT and by means of real data of "SEPE", through the years 2009–2016, that concern the most important hazard-source (i.e., the one with the highest risk) at the workplaces of the construction-sector (C-S) and also of the fixed-telecommunications technical-projects of OTE SA, which has been proven to be the risk-source of *"slipping—stumbling and falling—fall of persons"* (i.e., with ES-C#50).

**Figure 2.** Time-variations of the risk (R) calculated by proportional risk-assessment technique (PRAT) utilizing accident data of "SEPE", during 2009–2016, and concerning the most considerable hazard-source of "slipping—stumbling and falling—fall of persons" (with ES-C#50) at the workplaces of the construction sector (C-S) (**a**) and also of OTE SA (**b**).

In Figure 3, we illustrate time-profiles of risk (R) calculated by the PRAT and via accident data of "IKA", during the same period of 2009–2016, that regard the most significant hazard-source at the worksites of the C-S and of the technical-projects of OTE SA as well, which has been proven to be the risk-source of *"Loss of control (total or partial) of machine, means of transport or handling equipment, hand-held tool, object, animal"*(i.e., with ES-C#40).

**Figure 3.** Time-variations of the risk (R) calculated by PRAT utilizing accident data of "IKA", during 2009–2016, and concerning the most considerable hazard-source of "loss of control (total or partial) of machine, means of transport, or handling equipment, hand-held tool, object, animal" (with ES-C#40) at the workplaces of (**a**) the C-S and (**b**) OTE SA.

It is worth noting that the analysis by the TSP process of the time-profiles of the magnitude of R, concerning the deviation with ESAW-code of ES-C#50, shows (by the blue triangular symbols in the graphs 'a' and 'b' of Figure 2) the existence of a periodic fluctuation with a periodicity (T1) that is equal to ~2 years. This seems to be a permanent feature in the "behavior" of the ES-C#50 deviation at the workplaces of the C-S and OTE SA, because apart from the "SEPE" database, this periodicity obviously exists in the graphs 'a' and 'b' of Figure 4, which illustrates the R(t) profiles of ES-C#50 that were recorded at the same workplaces, but grounded, on the other hand, on the data of the "IKA" database.

**Figure 4.** Time-profiles of the risk (R) calculated by PRAT utilizing accident data of "IKA", during 2009–2016, and concerning the hazard-source of "slipping—stumbling and falling—fall of persons" (with ES-C#50) at the workplaces of the C-S (**a**) and also of OTE SA (**b**).

Moreover, the time-profiles of the risk-value R, concerning the deviation of ES-C#40, illustrate (by the blue triangular symbols in the graphs 'a' and 'b' of Figure 3) the appearance of a periodic fluctuation with a periodicity (T2) equal to ~4 yrs, and accordingly the presence of a harmonic of T1-2years, i.e., equivalent to T1.

Hence, the presence of a dominant periodic-component (~2yrs) in the most significant hazard source (taking into account the risk-value R) constitutes a steady characteristic for the dynamic behavior of a working environment (or working system), like the workplaces of the C-S and/or of the OTE fixed-telecommunications technical-projects, relatively to OSH safety.

What is more, the time-domain analysis of the time-series of the ES-C#50 deviation shows in its time-profiles (according to Figure 2a, Figure 4a,b) the appearance of a trend factor (or slope) with a negative inclination in the curve of averaging i.e., the average risk-value decreases at the constructions-sector workplaces, during the years 2009–2016, and only the profile of ES-C#50 in drawing 2b, presents a positive inclination or a soft increase in the curve of the risk-average at the workplaces of OTE SA. Besides, the analysis of the time-variations of the ES-C#40 hazard-source displays in its profiles (according to Figure 3a,b and Figure 5b) the appearance of a tendency factor with a negative slope in the curve of averaging, which means the average risk-value decreases at the construction-sector workplaces, during the years 2009–2016, and only the profile of ES-C#40 in drawing Figure 5a indicates a positive slope or an enhancement in the curve of the risk-average.

In addition, the comparison between Figures 2a and 5a, and also between Figures 2b and 5b, reveals that the profile of the deviation ES-C#50 presents all over the period 2009–2016 higher risk-values R in comparison with deviation ES-C#40, confirming that the most significant hazard source at the worksites of the C-S and SA is the ES-C#50 one, according to the "SEPE" database, contributing in this way to the 1st ranking of deviations (as stated by Figure 1). On the other hand, the comparison between Figures 3a and 4a, and also between 3b and 4b, discloses that the profile of deviation ES-C#40 shows in the course of the period 2009–2016 greater risk-values in comparison with deviation ES-C#50, confirming that the most significant hazard-source at the worksites of the C-S and OTE SA is the ES-C#40 one, according to the "IKA" database, contributing to the 1st ranking of Figure 1, as well.

**Figure 5.** Time-variations of the risk (R) calculated by PRAT utilizing accident data of "SEPE", during 2009–2016, and concerning the hazard-source of "loss of control (total or partial) of machine, means of transport, or handling equipment, hand-held tool, object, animal" (with ES-C#40) at the workplaces of (**a**) the C-S and (**b**) OTE SA.

We note that the other hazard-sources (i.e., with ES-C #10–#30, and #60–#99) of Table 1 present noticeably lower risk-values R compared to the ones of ES-C#40 and ES-C#50, so the resulted deviation-ranking is: (i) ES-C#40 and ES-C#50, consistent with the "IKA" database, and (ii) ES-C#50 and ES-C#40, according to the "SEPE" database. Consequently, a second and more trustworthy ranking (denoted as 2nd ranking in Figure 1) of the hazard-sources (like the one applied by a MCDM method) is essential, which could be utilized by the safety-responsible (or managers) in the procedure of prioritizing investments in OSH-prevention measures. Inevitably, in the next section, we show the application of the MCDM process.

#### *4.2. Application of the Typical-AHP and the Fuzzy-Extended-AHP MCDM-Processes*

*The concept of the Typical Analytical Hierarchy Process (AHP):* The Typical-AHP, introduced by Saaty (1990) [66], is one of the most popular multicriteria decision making methods. While applying the AHP, the problem under study is modeled as a hierarchy, and the local and global weights are calculated with respect to each criterion and/or sub-criterion as an alternative of the problem. The method takes as input the judgments of the decision maker via linguistic variables, transforms them to numerical values using a scale, and applies pairwise comparisons for extracting the ranking of factors according to their importance. The major advantage of the AHP method is its integrated function for estimating the decision maker's consistency of judgments, which is a key characteristic for producing reliable results. In that manner, the consistency ratio (CR) is calculated to monitor if the decision maker has a consistent axiom system and is not making judgments "at random". Note that a pairwise comparison matrix is considered as consistent if its CR index is less than or equal to 10%.

*The concept of Fuzzy-Extended-AHP Process (FEAHP):* In the present study, the Fuzzy Extended AHP (FEAHP) proposed by Chang (1996) [67] is used for taking account of the uncertainty of the decision maker's judgments. The FEAHP is a popular extension of AHP because of its "ease of use" feature and of its proven efficiency. In our approach, the only prerequisite from the decision maker are the judgments, which are assigned to quantitative values using the scale proposed by Lamata (2004) [68] and illustrated in Table 2. The corresponding Triangular Fuzzy Number (TFN) value is illustrated in the last column. Note that the value of the standard Saaty's scale [66] is highlighted in bold letters.


**Table 2.** Linguistic variables and corresponding triangular fuzzy numbers (TFNs) by Lamata (2004) [68].

The process works as described by Chang (1996) [57], and accordingly, the output triangular fuzzy numbers are transformed into crisp numbers with the average value approach by Zimmermann (2001) [69] as the defuzzification process. Thus, for a triangular fuzzy number M = (l,m,u), the corresponding crisp number C as:

$$\mathbb{C}(M) = (l + m + u) / 3\tag{5}$$

*Application of the typical-AHP and FEAHP MCDM processes:* The problem of assessing nine hazard factors with OSH impact at the workplaces of the C-S and OTE SA is modeled as a hierarchy, in which the goal is preventing safety, and the 1st level consists of nine hazard factors (as in Figure 6).

**Figure 6.** The hierarchy used for calculating the weights.

The numbers filled in the comparison matrix (Table 3) by the decision maker correspond to linguistic variables (Table 2) and express his preferences regarding the importance of hazards. As for the consistency of the judgments, the typical-AHP is employed first, and the consistency ratio (CR) is computed. In the present case, CR was found to be less than 0.1, meaning that the resulting weights are reliable, since the judgments of the decision maker are consistent.


**Table 3.** The judgments of the decision maker.

The resulting rankings and weights computed by the typical-AHP and the fuzzy extended FEAHP for every hazard are illustrated in Table 4, with a descending order regarding their importance. Note that these results were extracted using the same decision maker's judgments as shown in the pairwise comparison matrix of the Table 3.



As for the results of the typical-AHP calculations, the factors "slipping—stumbling and falling—fall of persons" (ES-C#50) and "loss of control (total or partial) of machine, means of transport, or handling equipment, hand-held tool, object, animal" (ES-C#40) are more important, since their weights are 31.51% and 22.63%, respectively. The next more important group includes the factor "breakage, bursting, splitting, slipping, fall, collapse of material agent" (ES-C#30) with a weight of 16.06%, and the (ES-C#10) "deviation due to electrical problems, explosion, fire" with a 10.04% score. Regarding the less important hazards, the sum of their weights is about 19.77% of the total risk, which is considered quite important. More specifically, this group of hazards includes the "body movement under or with physical stress (generally leading to an internal injury)" (ES-C#70), "body movement without any physical stress (generally leading to an external injury)" (ES-C#60), "deviation by overflow, overturn, leak, flow, vaporization, emission" (ES-C#20), and "shock, fright, violence, aggression, threat, presence" (ES-C#80), and "other deviations not listed above in this classification" (ES-C#99).

Applying the Fuzzy Extended AHP (FEAHP) to the judgments of Table 3 results in a similar hazards' ranking, as illustrated in the last column of Table 4. But, notwithstanding this, the relative weights' distance between the first deviation (with ES-C#50) and the second one (with ES-C#40) is very much reduced, meaning that the FEAHP describes the uncertainty of the expert's choices, in a different manner than typical-AHP.

Additionally, it is observed that, except for the first ranked deviation of ES-C#50, which has clearly lower score than before, and the second deviation of ES-C#40, which has almost the same weight, all the other hazards' weights are higher, and they are responsible for 53.29% of the total risk, instead of the total 45.86% calculated with typical-AHP.

Given the fact that these conclusions are extracted from a consistent pairwise comparison matrix and the importance of hazards is close to one another, it is preferable to use FEAHP instead of the typical-AHP under these circumstances, as stated by Chan et al. (2019) [70]. The use of FEAHP ensures that the proposed approach is more efficient since it handles better the judgments' uncertainty, which can affect the health and safety of employees.

#### *4.3. Application of the FTA Process*

Taking into account the MCDM outcomes extracted (in Table 4) by the utilization of the typical-AHP and FEAHP, along with the result that the hazard-sources' ranking unveiled as the most important deviation at the workplaces of the C-S and of OTE SA, the ES-C#50 one, while as the second-ordered deviation the ES-C#40 one, we complete the proposed RAA generic algorithmic framework by the application of the FTA process, as far as the ES-C#50 at the worksites of the fixed-telecommunications technical-projects of OTE SA is concerned: (i) For graphically illustrating the "base-events" and their sequences to the "top-event", and (ii) for calculating the occurrence probability of the "top-event" (ES-C#50) of the most significant deviation.

Therefore, the resulting FTA drawing is illustrated in Figure 7, where the likelihood of the "top-event" ES-C#50 has been calculated to be P = 51.3% according to the estimated (by the safety manager) likelihoods of the "base-events" E1–E12 depicted in the Table 5 and the utilization of equations (2), (3), and (4) as well. The "intermediate events" G1–G4 (illustrated in Table 6) have also been inserted in this FTA graph. Particularly, the events G1, G2, and G4 constitute the outputs of three transitional OR-gates, while G3 comprises the output of an AND-gate.

**Figure 7.** The fault-tree analysis (FTA)Drawing, as far as the most significant hazards' source of ES-C#50 at the worksites of the fixed-telecommunications technical-projects of OTE SA is concerned.


**Table 5.** The estimated likelihoods of the "base-events" E1–E2 concerning the FTA graph of the most significant hazards' source of ES-C#50 ("top-event").

**Table 6.** The "intermediate-events" G1–G4 concerning the FTA graph of the most significant hazards' source of ES-C#50 ("top-event").


To boot, the "minimal cut-sets" were determined according to the theory (e.g., [71]), and the fault-tree was first "translated" to its equivalent Boolean equations, and then, by the "top-down" method of replacement, the following Boolean expression (which also constitutes the "minimal cut-set" expression) for the "top-event" was exported:

$$\text{"TOP\\_EVENT"} = \sum\_{i=1}^{4} E\_i + (E5 \cdot E6) + \sum\_{j=7}^{12} E\_j \tag{6}$$

Taking into consideration that (i) the probability of the "top-event" is greater than 50% (during the period of 2009–2016), and (ii) the work of Marhavilas and Koulouriotis (2012a) [9], we characterize the ES-C#50 deviation as a high-risky hazard source.

#### **5. Discussion and Conclusions**

The way that a technical project attains its sustainability goal is implemented by emphasizing, among other aspects, on sustainable living, for instance on the actualization of health, that includes among other things, the buildings' and infrastructure technology [72,73]. What's more, the OSHA accepts that actions to guarantee OSH can be involved in sustainability efforts and take advantage of the dynamics of sustainability to make each workplace safer and healthier [3,4].

Besides, the notion of sustainability is utilized as a frame to sort subsistent (or to grow new) OSH-standards, and also to develop novel and innovative OSH-RAA assessment techniques. Thus, a sizeable part of the scientific literature employs the topic of sustainability as a means for the evolution of OSH-RAA techniques [5]. In general, the diversity of OSH-RAA techniques is such that there are

plenty of efficient approaches for most cases, and many studies build up hybrid models merging distinct RAA approaches.

In this sense, this article concentrates on: (i) The implementation of a new RAA methodological framework (the MCDM-STO/DET one) for sustainable engineering projects by the amalgamation of a multicriteria decision-making (MCDM) process with the joint-collaboration of a deterministic (DET) and a stochastic (STO) process, and (ii) the application of MCDM-STO/DET at the workplaces of the Greek construction sector and also of the fixed-telecommunications technical projects of OTE SA (that is the Greek Telecommunications Organization S.A.), by means of real accident data coming from two official State databases, namely of "SEPE" (Labor Inspectorate, Hellenic Ministry of Employment) and of "IKA" (Social Insurance Institution, Hellenic Ministry of Health), all the way through the years of 2009–2016.

Hence, the article's objectives are next: (i) The presentation and execution of the joint MCDM-STO/DET framework at the workplaces of the Greek construction sector and also of the fixed-telecommunications technical projects of OTE SA, and (ii) to make known that the proposed MCDM-STO/DET scheme can be a precious method for safety managers (and/or decision-makers) to ameliorate OSH and also to support the sustainable operation of technical or engineering projects.

Mainly, we mingle two different configurations of the MCDM method, initially the typical AHP technique, and afterwards the fuzzy-extended AHP (FEAHP) one, along with the DET-PRAT and the STO-TSP methods, and finally with the DET-FTA one.

The results of the DET-PRAT utilization could be summarized on the subsequent topics:


The results of the DET-TSP application by the time-domain analysis of the time-series of accident data and risk values (R) could be summarized as follows:


In other words, the application of the DET-PRAT technique, along with the DET-TSP one, revealed (by the analysis of the time-profiles of the most significant hazard sources) two significant factors: (i) The presence of a dominant periodic-component (~2years), and (ii) the appearance of a trend factor in the magnitude of risk (R) with a decreasing tension (in almost all the cases), constituting (during 2009–2016) two permanent characteristics for the dynamic behavior of the working environment (or working system) at the workplaces of the C-S and also of the OTE fixed-telecommunications technical-projects, relative to OSH safety, and should be taken into consideration by any responsible (or manager, officer, etc.) for OSH issues, both at the business level and also at the State level.

Taking into consideration that the resulting rankings of the hazard-sources (at the workplaces of the C-S and of OTE SA) by using PRAT and the data of "SEPE" and "IKA" were dissimilar and/or contradictory, a second (final) and more trustworthy ranking of the hazard-sources was achieved by the application of a MCDM method (in its two different configurations, the typical-AHP and the FEAHP one) in order to attain the optimum distribution of an investment budget to appropriate measures reaching the best OSH-protection (i.e., in the procedure of prioritizing investments in OSH-prevention measures). The typical-AHP and the FEAHP were employed in the present study for assessing hazards with an efficient and "easy-to-use" multicriteria method and its fuzzy extension. The resulting rankings show: (i) The necessity of using the FEAHP in OSH applications, since it can better handle the uncertainty of the decision maker's judgments, and (ii) that the most significant hazard-source is ES-C#50 (while as second, the ES-C#40 one).

Finally, we completed the suggested RAA algorithmic framework by the application of the DET-FTA process, as far as the most significant deviation ES-C#50 at the worksites of the fixed-telecommunications technical-projects of OTE SA is concerned, for illustrating the "base-events" and their sequences to the "top-event", and for calculating the probability of the "top-event" (ES-C#50), which was P = 51.3% (i.e., a high-risk hazard source).

The partnership of PRAT with TSP technique can give, by the computation of risk (R), indispensable information to safety-managers in order to identify the most significant hazard-sources at the corporations' worksites, and also to realize a preliminary ranking (1st ranking) of them. This ranking can be worthy in evaluating the urgency of investing a restricted budget in particular measures for attaining the highest OSH-protection with the minimal cost. Moreover, the next module in the suggested algorithmic framework is the MCDM one, which represents the application of a multicriteria decision-making technique, for accomplishing the final (definitive) ranking (2nd ranking) of the hazardous sources at the worksites, in order to achieve the finest (optimum) distribution of an investment budget (i.e., with the minimal cost) for appropriate measures for reaching the best OSH-protection. In addition, the FTA technique could also efficiently facilitate the aforementioned 1st and 2nd rankings.

It is worth mentioning that to facilitate the calculations, we could use the following integrated software packages: (i) The GRAPHER by Golden Software Inc, and the SPSS one by SPSS Inc, for achieving the TSP analysis, and (ii) the Fault-Tree+ one by ISOGRAPH Inc for the FTA analysis.

Besides, the results of the MCDM calculations, for both Typical-AHP and FEAHP, have been extracted using MS-Excel spreadsheets developed by one of the authors. Our results, calculated by the Typical-AHP, could be confirmed using the software available at http://fuzzymcdm.upol.cz/. This software applies the Typical-AHP (as presented by Saaty; 1990) [66], but it supports a different extension of AHP to fuzzy logic (as presented by Holeˇcek and Talašová; 2016) [74] than the FEAHP performed in this study, i.e., calculations with FEAHP only include the application of the equations as

described in the work of Chang (1996) [67]. Because of this, the results with this software might be slightly different.

The main limitations of using this approach are listed below:


This multiparty approach can assist safety risk managers by the prioritization of hazards, to reduce the impact of the most important risks. In addition, it could be used by risk managers to allocate more effectively the constrained budget and to reduce mortgage costs, and thus, for reserving extra funds for spending on safety measures, and to further reduce other risk factors' impact. The main contribution of this study is the cooperative use of a proven efficient multicriteria decision-making method with a joint STO/DET one for identifying and prioritizing hazards.

Accordingly, as a future improvement, the FTA technique could be utilized for contributing more capably to the 1st and 2nd ranking of the hazard-sources, according to the configuration of Figure 8, which constitutes an alternative form of the proposed RAA framework (illustrated in Figure 1).

As a general conclusion, every RAA approach has specific limitations and/or restrictions in its applicability to miscellaneous accidents, and consequently, a separate technique can't accomplish the optimum risk-assessment result at the workplaces. On the other hand, collaborative methodologies, like the suggested MCDM-STO/DET one, would be efficient for sustainable engineering projects. Thus, the novelty of this article is fulfilled through the suggested newfangled risk assessment and analysis (RAA) methodological approach by the amalgamation of a multicriteria decision-making (MCDM) process with the joint-collaboration of a deterministic (DET) and a stochastic (STO) process in order to upgrade OSH and to validate the sustainable operation of technical or engineering projects, as well.

**Figure 8.** An alternative form of the proposed joint stochastic/deterministic/MCDM RAA methodological framework.

**Author Contributions:** Conceptualization, P.K.M.; data curation, P.K.M., M.G.T., G.K.K.; formal analysis, P.K.M.; investigation, P.K.M.; methodology, P.K.M., DEK; writing—original draft, P.K.M., G.K.K.; writing—review and editing, P.K.M., G.K.K., D.E.K. All authors have read and agreed to the published version of the manuscript

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

**Conflicts of Interest:** We declare no conflict of interest.
