A Joint Stochastic / Deterministic Process with Multi-Objective Decision Making Risk-Assessment Framework for Sustainable Constructions Engineering Projects—A Case Study

: 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 ﬁxed-telecommunications technical projects of OTE SA (that is, the Greek Telecommunications Organization S.A.) by means of real accident data coming from two o ﬃ cial 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 di ﬀ erent conﬁgurations 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 ﬁnally with the Fault-Tree Analysis (FTA).


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 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. (2014aMarhavilas et al. ( , 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][50][51][52][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 (2009Marhavilas ( , 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. 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.

The Suggested Methodological Framework
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: Sustainability 2020, 12, 4280 and where B i (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).

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.

Working Environment & Processes 10
Deviation due to electrical problems, explosion, fire [e.g., (i) exposure to or contact with extreme temperature levels, (ii) exposure to or contact with electric current, etc.]

Working Environment & Processes 20
Deviation by overflow, overturn, leak, flow, vaporization, emission (e.g., exposure to or contact with hazardous substances or radiation, etc.)

Working Environment & Processes 30
Breakage, bursting, splitting, slipping, fall, collapse of material agent (e.g., slipping, collapse and being struck by falling objects)

Psychological/Human 40
Loss of control (total or partial) of machine, means of transport or handling equipment, hand-held tool, object, animal [e.g., (i) collision with an immobile object and falling against or being struck by moving objects, (ii) trapping, being crushed-inside or between objects, etc.]

Physical Activity 50
Slipping-stumbling and falling-fall of persons [e.g., (i) falling of person from a height, (ii) falling of person-on the same level]

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). 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.  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).
Sustainability 2020, 12, x FOR PEER REVIEW 9 of 24 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.
(a) (b)   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 (T 1 ) 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 Sustainability 2020, 12, 4280 8 of 21 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.

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  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 (T 2 ) 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.

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 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.

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].  [7,8,9] Extremely [8,9,9] 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:

Importance of Factor i over Factor j Fuzzy Number
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). Sustainability 2020, 12, x FOR PEER REVIEW 12 of 24 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. Other Deviations not listed above in this classification Figure 6. The hierarchy used for calculating the weights.

SAFETY
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. 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.

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.
Sustainability 2020, 12, x FOR PEER REVIEW 14 of 24 method of replacement, the following Boolean expression (which also constitutes the "minimal cut-set" expression) for the "top-event" was exported: 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.  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").

Base-Events
Description of "Base-Events"   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: 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.

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: 1.
The most significant hazard-source at the worksites of the C-S and OTE SA is the ES-C#50 one, according to the "SEPE" database (during 2009-2016).

2.
The most significant hazard-source at the worksites of the C-S and OTE SA is the ES-C#40 one, according to "IKA" database (during 2009-2016).

3.
The resulting hazard-sources ranking based on the calculated risk-value R is, on the one hand: There are other considerable hazard sources which present a risk value higher than 100.0, and according to the work of Marhavilas and Koulouriotis (2012b) [10], long-term actions are necessary for the extinction of their possible dangers. 6.
Besides, other hazard sources present a risk value smaller than 100.0, and according to the previous referenced paper [10], compulsory actions are not essential except for surveillance of the events.
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: The time-profiles of the ES-C#50 deviation show, in almost all cases, the appearance of a trend factor with a negative inclination in the curve of averaging i.e., the average risk-value decreases at the C-S workplaces (during 2009-2016).

2.
Instead, the profile of ES-C#50 presents a positive inclination (or a soft increase) in the curve of the risk-average, only at the workplaces of OTE SA.

3.
Besides, the analysis of the time-variations of the ES-C#40 hazard-source displays in its profiles the presence of a tendency factor with a negative slope in the curve of averaging, that means the average risk-value decreases at the workplaces of the C-S and OTE SA (during 2009-2016).

4.
Instead, the profile of ES-C#40 indicates a positive slope (or an enhancement) in the curve of the risk-average, only at the workplaces of the C-S.

5.
The analysis of the time-profiles of the magnitude R, concerning the ES-C#50, shows the existence of a periodic fluctuation with a periodicity of~2 years, which 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, according the "SEPE" and "IKA" database. 6.
Likewise, the time-profiles of R, concerning the deviation of ES-C#40, illustrate the appearance of a periodic fluctuation with a periodicity of~4 years (i.e., a harmonic of~2years).
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ček 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: • The AHP theory tries to measure the relative importance of alternatives with respect to each criterion by using pairwise comparisons. In order to do pairwise comparisons, experts use a conversion scale with crisp values, for expressing their ideas. So, an important limitation of AHP is the usage of crisp values to reflect human thinking.

•
Another limitation is coming from CR, which is a verification of the rational judgment performed by the experts.

•
Moreover, an additional limitation concerns the pairwise comparisons of criteria, which are carried out by several experts with a required significant working-experience in occupational safety.

•
As a classic multicriteria decision support tool, AHP is a subjective process, and the resulting results, either rankings or weights, are dependent on the way that the judgments are imposed and the criteria are compared to the rest. In other words, another decision maker (or a group of experts) could make "different" judgments regarding the relative importance of each pair of criteria, resulting in different resulting weights for the criteria and/or different rankings. This drawback leads to the need for applying group-decision-making processes for merging different judgments.

•
Besides, when there are slight differences between the criteria's weights, it is possible that the AHP could lead to the selection of a suboptimal alternative, instead of the optimal one. This phenomenon could be reduced by considering many criteria to be pairwise compared. A review of the AHP drawbacks can be found in the study of Whitaker (2007) [75].
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. Funding: This research received no external funding.

Conflicts of Interest:
We declare no conflict of interest.