1. Introduction
Railway infrastructure safety is an important element of urban public safety [
1]. The safety of railway infrastructure is vital to secure goods and passengers from departures to destinations [
2]. Accidental events in railway transportation cause damages to human health, public property, environment, and the economy [
3]. In 2018, there were 1721 significant railway accidents in the European Union, with a total of 853 fatalities and 760 serious injuries. The total cost of railway accidents was approximately 5 billion EUR [
4]. Collisions and derailments account for about 200 accidents each year [
5].
To increase safety levels of assets, passengers, goods, and employees as well as prevent significant accidents it is essential to prioritize risk factors and assess railway infrastructure. However, the risk assessment of railway infrastructure is an emerging problem for railway planners. It is considered as a multi-criteria decision-making (MCDM) problem, since there are a finite set of available railway infrastructures, which safety levels need to be assessed, and numerous risk factors influencing railway infrastructure (i.e., evaluation criteria). Additionally, to mitigate information loss and prevent erroneous decisions, it is critical to efficiently represent and handle risk-related uncertain information when solving this complex MCDM problem.
Solving the railway infrastructure risk assessment problem is critical for railway transport worldwide. The key question is how to identify unsafe railway infrastructure so relevant authorities can undertake safety improvement projects on time. However, (i) only deterministic numbers and type-1 fuzzy sets have been used in the available studies to evaluate risk factors in railway transport; (ii) there is no methodological framework that can represent and handle uncertain, incomplete, and inconsistent risk-related information; and (iii) the available prioritization methods are unable to consider the subjective and objective importance of risk factors under a highly uncertain environment. As a result, the main objective of this paper is to introduce a picture fuzzy group MCDM approach for risk assessment of railway infrastructure.
Recently, Cuong and Kreinovich [
6] and Cuong [
7] proposed picture fuzzy sets (PFSs). PFSs represent an advanced generalization of intuitionistic fuzzy sets for efficient uncertainty modeling and solving real-life decision-making problems [
8]. They can efficiently describe fuzzy, uncertain, incomplete, and inconsistent information [
9]. PFSs are characterized by four functions, namely, the degree of positive membership, the degree of neutral membership, the degree of negative membership, and the degree of refusal membership [
10,
11]. Therefore, they are especially suitable for decision-making situations that require answers of the type yes, no, abstain, and refusal [
12]. This particular type of information is dominantly used in surveys and voting systems where experts are divided into previously mentioned four categories [
13]. As a result, PFSs can efficiently represent railway planners’ preferences and handle risk-related information. Unfortunately, railway planners are unable to naturally express their preferences by voting, since no previous research has applied a PFS-based MCDM approach for railway transport.
The Measurement of Alternatives and Ranking according to Compromise Solution (MARCOS) is one of the most recent MCDM methods developed by Stević et al. [
14]. It synthesizes the ratio and reference point sorting approaches to obtain ranking results [
15]. Compared with other MCDM methods, the MARCOS method is simple, effective, and easy to sort and optimize a decision-making process [
14,
16]. Its features are [
14,
15] consideration of ideal and anti-ideal solutions at the very beginning, the definition of utility degree to both solutions, ability to process large data sets, flexibility to analyze expert preferences, and a simple algorithm. However, the MARCOS method has not been extended before using PFSs, so it cannot reflect neutral/refusal information of decision-makers in the railway industry.
Based on the highlighted research gaps, the aims of this paper are (1) to utilize PFSs, which are superior in handling fuzzy, uncertain, incomplete, and inconsistent risk-related information, and help railway planners to naturally express their risk preferences by voting; (2) to determine the importance of each risk factor influencing railway infrastructure by using the novel picture fuzzy hybrid method for risk factor prioritization, which is developed by hybridizing the picture fuzzy direct rating and Tsallis–Havrda–Charvát entropy methods; (3) to rank railway infrastructures by employing the developed picture fuzzy MARCOS method for railway infrastructure ranking, which is for the first time extended under the picture fuzzy environment; and (4) to apply the formulated picture fuzzy group MCDM approach for risk assessment of railway infrastructure in the Czech Republic.
The rest of the paper is organized as follows:
Section 2 surveys the state-of-the-art research.
Section 3 reviews some definitions of PFSs.
Section 4 presents the introduced picture fuzzy group MCDM approach for risk assessment of railway infrastructure. A real-life case study is described in
Section 5.
Section 6 presents the case study results and discussions.
Section 7 gives the conclusions of the work and indicates possible extension areas.
2. Literature Review
The literature review is organized into three sub-sections. The first sub-section surveys existing MCDM approaches for railway transport. The second sub-section overviews applications and extensions of the MARCOS method. The last sub-section presents identified research gaps.
2.1. MCDM Approach for Railway Transport
The application of the MCDM approach for solving diverse problems in railway transport attracted a large interest of researchers, since it is important for the economy and society. The comprehensive summary of the state-of-the-art contributions is presented in
Table 1.
Mohajeri and Amin [
17] coupled the Analytical Hierarchy Process (AHP) with Data Envelopment Analysis (DEA) to evaluate sites for a railway station. Sivilevičius and Maskeliūnaite [
18] applied the AHP method to assess quality criteria of passenger railway trips in international routes. Brunner et al. [
19] utilized the AHP method to compare fixed-rail transit alignment locations plotted in the Geographic Information System (GIS) environment. Liu et al. [
20] integrated the analytic network process (ANP) and the technique of multiple regression to elicit the environmental impact of high-speed railway systems. Poorzahedy and Rezaei [
21] applied several MCDM methods to rank alternative configurations of a light railway transit network.
Celik et al. [
22] developed an interval type-2 fuzzy (VIšeKriterijumska Optimizacija i kompromisno Rešenje) VIKOR method to rank railway transit lines. They hybridized the statistical analysis, service quality (SERVQUAL) model, and direct rating method to classify and acquire criteria weights of passenger satisfaction criteria. Lu et al. [
23] provided a simple fuzzy framework to determine critical risk factors in railway reconstruction projects. Nathanail [
24] applied the Delphi-AHP approach to evaluate quality indicators of a railway network.
Aydin et al. [
25] coupled the statistical analysis, fuzzy AHP method, and Choquet integral to compare railway transit lines from the passenger satisfaction perspective. Dhir et al. [
26] used the cost–benefit analysis (CBA) and AHP method to identify the most suitable manufacturer for high-speed rail vehicles. An et al. [
27] employed the fuzzy AHP method to obtain contributions of hazard groups in a shunting depot. Zhao et al. [
28] used the fuzzy AHP method to prioritize fire emergency alternatives of an unattended train operation metro system.
Hamurcu and Eren [
29] utilized the AHP and Technique for the Order Preference by Similarity to Ideal Solution (TOPSIS) methods to select the best monorail technology. Montesinos-Valera et al. [
30] exploited the ANP method with benefits, opportunities, costs, and risks (BOCR) analysis to order projects for maintenance, renewal, and improvement of railway lines. Polat et al. [
31] applied the fuzzy AHP-TOPSIS approach to compare rail suppliers in intercity construction projects. Song et al. [
32] utilized the AHP method and interpretative structural modeling (ISM) approach to ascertain the importance and interrelations of vulnerability factors of urban rail transit operations, respectively. Stević et al. [
33] combined the best–worst method (BWM), simple additive weighting (SAW) method, and rough numbers to select used rail wagons for carrying out internal transport.
Hamurcu and Eren [
34] used the ANP method to prioritize railway construction projects. Mlinarić et al. [
35] employed the Delphi-AHP approach to prioritize a common set of key performance indicators of railway intelligent transportation systems. Sharma et al. [
36] combined the rough AHP and Multi-Attributive Border Approximation Area Comparison (MABAC) methods to rank railway stations from the service quality perspective. Wang et al. [
37] utilized the AHP method and Fuzzy Reasoning Approach (FRA) to assess safety performances of railway track systems. Yadegari et al. [
38] used the Delphi-ANP approach to prioritize factors affecting the development of an industry cluster of the rail industry.
Đorđević et al. [
39] formulated a rough FUll Consistency Method (FUCOM) to compute subjective weights of service quality criteria in railway transport from the position of passengers with disabilities. Jasti and Ram [
40] used the AHP method to appraise performance indicators and benchmark metro systems. Krmac and Djordjević [
41] utilized the AHP method to obtain weights of strengths, weaknesses, opportunities, and threats of implementing rail traffic management systems. Norouzi and Namin [
42] applied the fuzzy BWM-TOPSIS approach to sort risks of constructing railway megaprojects. Phanyakit and Satiennam [
43] employed the AHP-SAW approach to select a rehabilitation alternative for railway tracks structure. Wu et al. [
44] used the fuzzy ANP method to evaluate risks of constructing metro stations.
Recently, Blagojević et al. [
45] coupled the DEA and SAW methods to rank railway line sections in terms of safety. Brumercikova and Sperka [
46] used the AHP method to obtain an entrance order of freight carriers in railway stations. Huang et al. [
47] utilized two objective criteria weighting methods and historical data to elicit risks of transporting dangerous goods by railway. Kumar and Anbanandam [
48] applied the grey decision making trial and evaluation laboratory (DEMATEL) method to quantify relationships of prominent inhibitors to the implementation of railroad freight transport. Li et al. [
49] presented a hybrid Pythagorean MULTIMOORA method to order rail transit lines according to the passenger satisfaction level. Majumdar et al. [
50] exploited the TOPSIS method to prioritize key factors affecting metro rail infrastructure based on commuter perception. Sangiorgio et al. [
51] integrated the AHP method and linear programming approach to estimate the number of accidents and assess the safety of railway networks. Stoilova [
52] coupled the AHP and sequential interactive model of urban systems (SIMUS) methods to select railway passenger plans. Stoilova et al. [
53] used the SIMUS method to compare railway transport systems of different countries. Vesković et al. [
54] combined the pivot pairwise relative criteria importance assessment (PIPRECIA) and evaluation based on distance from average solution (EDAS) methods to rank options for achieving business balance in organizing railway traffic. Wang and Shi [
55] applied an interval-valued intuitionistic fuzzy weighted average operator and maximum entropy method to prioritize the traffic service quality of urban railway transit lines. Zhang and Sun [
56] formulated a Delphi-DEMATEL-ANP-TOPSIS approach to assess shunting derailment risk response strategies in railway stations. Tavassolirizi et al. [
57] used the DEMATEL-ANP approach to prioritize key delay factors in implementing railway development projects and identify their interactions.
2.2. Applications and Extensions of the MARCOS Method
The MARCOS method is one of the latest MCDM approaches. Its applications and extensions are summarized in
Table 2.
Badi and Pamucar [
58] coupled the grey direct rating and MARCOS methods to evaluate suppliers for the iron and steel industry. Ilieva et al. [
59] utilized the fuzzy MARCOS method to order cloud storage systems. Chakraborty et al. [
60] ranked suppliers for the iron and steel industry based on D numbers and the MARCOS method. Puška et al. [
61] used the MARCOS method to rank project management software solutions for facilitating business operations. Stanković et al. [
62] developed a fuzzy MARCOS method to identify a road section with the highest traffic risk for all participants. They applied the fuzzy PIPRECIA method to determine the importance of criteria affecting the traffic risk. Stević and Brković [
63] combined the FUCOM and MARCOS approach to compare performances of drivers engaged in international road transport. Stević et al. [
14] developed the MARCOS method and used it to solve the sustainable supplier selection problem in the healthcare industry. Ulutaş et al. [
16] integrated the correlation coefficient and standard deviation (CCSD), indifference threshold-based attribute ratio analysis (ITARA), and MARCOS methods to select manual stacker for performing logistics activities in warehouses.
2.3. Research Gaps
According to the performed comprehensive literature review, the research gaps are as follows: (i) No previous research has applied a PFS-based MCDM approach for railway transport. Hence, railway planners are unable to naturally express their preferences by voting; (ii) there is no risk assessment framework that can handle ambiguous, uncertain, and vague information, since only deterministic numbers and type-1 fuzzy sets have been used in the available studies; (iii) the direct rating and Tsallis–Havrda–Charvát entropy methods have not been hybridized before to determine the importance of risk factors; and (iv) the MARCOS method has not been extended previously using picture fuzzy sets.
3. Preliminaries
This section provides some definitions of picture fuzzy sets.
Definition 1 [
6,
7]
. Let PFS A on a universe X is an object in the form of:
whereis called the degree of positive membership of x inis the degree of neutral membership of x inis the degree of negative membership of x inandand satisfy the following condition:
The degree of refusal membership of x in the PFS A can be calculated as follows:
If X has only one element, thenis called a picture fuzzy number (PFN) in whichandFor convenience, a PFN is denoted by[64,65]. Definition 2 [
66]
. As a generalization of the algebraic product and sum, and Einstein T-conorm and T-norm, the Hamacher T-conorm and T-norm is more general and more flexible. Let
,
and be three PFNs, operational parameter ζ > 0, and λ > 0. The Hamacher T-norm and T-conorm operations of PFNs are defined as follows:
(c) Scalar multiplication Definition 3 [
67]
. Let and be two PFSs in X. A function En: PFSs(X) → [0, 1] is an entropy on PFS, if En satisfies the following axiomatic requirements: (a) Sharpness:if and only ifis a crisp set.
(b) Maximality:iffor all t = 1, …, f.
(c) Resolution:ifsatisfy eitherwhenorwhenfor all t = 1, …, f.
(d) Symmetry:wheredenotes the complement of
Definition 4 [
68]
. Tsallis–Havrda–Charvát entropy [69,70] is a generalized form of Shannon entropy [71]. This well-known one-parametric entropy is suitable for real-world applications, since its strength lies in properties and applications. Let be a PFS in X and χ > 0 (≠1) be the information measure parameter. The picture fuzzy Tsallis–Havrda–Charvát entropy measure of the PFS A is defined as follows: Definition 5 [
66]
. Let (l = 1, …, s) be a collection of PFNs, the operational parameter ζ > 0, and φ = (φ1, …, φs)T be the weight vector of the collection of PFNs with φl > 0 and The picture fuzzy Hamacher weighted average (PFHWA) operator is defined as follows: Definition 6 [
72,
73]
. Let be a PFN. A two-step defuzzification method to obtain a crisp value of the PFN is:
Step 1. Distribute the neutral degree to the positive and negative degrees as follows: Step 2. Calculate the crisp value y by: 4. Picture Fuzzy Group MCDM Approach for Risk Assessment
This section presents the developed picture fuzzy group MCDM approach for risk assessment of railway infrastructure (
Figure 1). The approach has three phases. In the first phase, linguistic importance evaluations towards railway infrastructures and risk factors are collected from invited experts and expressed as PFNs. In the second phase, subjective, objective, and hybrid importance of risk factors are computed. This phase represents the formulated picture fuzzy hybrid method for risk factor prioritization. Firstly, the picture fuzzy direct rating method is utilized to compute the subjective importance of risk factors. Secondly, the picture fuzzy Tsallis–Havrda–Charvát entropy method is applied to calculate the objective importance of risk factors. Thirdly, the subjective and objective importance are hybridized. In the last phase, the formulated picture fuzzy MARCOS method is used to rank railway infrastructures.
Let A = {A1, …, Am} (m ≥ 2) be a finite set of railway infrastructures that experts have to choose from, D = {D1, …, Dk} (k ≥ 2) be a set of invited experts, and C = {C1, …, Cn} (n ≥ 2) be a finite set of risk factors. The phases and encompassed steps of the developed group MCDM approach are given in the following:
Phase 1. Information collection and representation.
Step 1.1. Construct the linguistic evaluation matriceswhere m, k, and n are the number of railway infrastructures, invited experts, and risk factors, respectively;(i = 1, …, m; e = 1, …, k; j = 1, …, n) is the linguistic evaluation (i.e., yes, abstain, no, or refusal) given by the invited expert De towards the railway infrastructure Ai with respect to the risk factor Cj. Step 1.2. Determine the picture fuzzy evaluation matrixwhereis a PFN, which represents an evaluation of the railway infrastructure Ai with respect to the risk factor Cj given by the invited experts. Importance evaluations are represented as PFNs by computing the share of each item (i.e., yes, abstain, no, or refusal) in the related voting results of the invited experts. Step 1.3. Construct the linguistic risk factor importance matriceswhereis the linguistic importance evaluation (i.e., yes, abstain, no, or refusal) given by the invited expert De towards the risk factor Cj. Step 1.4. Determine the picture fuzzy risk factor importance matrixwhereis a PFN, which represents the importance evaluation of the risk factor Cj given by the invited experts. It is computed as the share of each item in the corresponding voting results of the invited experts. Phase 2. Picture fuzzy hybrid method for risk factor prioritization.
Step 2.1. Compute subjective importance of risk factors:whereis a PFN, which represents the importance evaluation of the risk factor Cj given by the invited experts; andrepresents the subjective importance vector of the risk factors, withandThe subjective importance of risk factors are computed by using the picture fuzzy direct rating method. Step 2.2. Compute objective importance of risk factors:whereis a PFN which represents an evaluation of the railway infrastructure Ai with respect to the risk factor Cj given by the experts; χ is the information measure parameter; andrepresents the objective importance vector of the risk factors, withandThe objective importance of risk factors are computed by using the picture fuzzy Tsallis–Havrda–Charvát entropy method. Step 2.3. Compute hybrid importance of risk factors:where γ ∈ [0, 1] represents the trade-off parameter;is the hybrid importance vector of the risk factors, withandThe risk factors have hybrid importance when γ ∈ (0, 1). They have exclusively subjective and objective importance when γ = 1 and γ = 0, respectively. Phase 3. Picture fuzzy MARCOS method for railway infrastructure ranking.
Step 3.1. Construct the picture fuzzy extended evaluation matrixwhereis a finite set of railway infrastructures;is a PFN which represents an evaluation of the railway infrastructure Ai with respect to the risk factor Cj given by the experts: In this step, the picture fuzzy evaluation matrix is extended by adding the anti-ideal and ideal railway infrastructures.
(1) Anti-ideal railway infrastructurewhere(j = 1, …, n) is a collection of PFNs, which represent anti-ideal values for each risk factor. (2) Ideal railway infrastructurewhere(j = 1, …, n) is a collection of PFNs which represent ideal values for each risk factor. Step 3.2. Determine the picture fuzzy additive relative importance of each railway infrastructure as follows:where ζ is the operational parameter. Step 3.3. Calculate utility degrees to the anti-ideal and ideal solutions of each railway infrastructure.
(1) Utility degree to the anti-ideal solution: (2) Utility degree to the ideal solution: Step 3.4. Calculate utility functions to the anti-ideal and ideal solutions.
(1) Utility function to the anti-ideal solution: (2) Utility function to the ideal solution: Step 3.5. Calculate the utility function of each railway infrastructure as follows: Step 3.6. Rank railway infrastructures.
Rank the railway infrastructures according to the increasing values of utility functions. The lowest value is the safest railway infrastructure.
5. Case Study
According to the 2020 census, the Czech Republic has 10.7 million inhabitants. More than 50% of the population use railway transportation. A real-life case study of risk assessment of railway infrastructure in the Czech Republic context is presented to evaluate railway infrastructure safety and demonstrate the utility of the formulated approach. The approach is used to prioritize twenty-three relevant risk factors (
Table 3) and assess six railway infrastructures (
Figure 2).
The investigated railway infrastructures in the Czech Republic context are:
(A1) Pardubice (
Figure 3a). The railway infrastructure of Pardubice is an important network hub. It is served by international and local trains. This railway infrastructure has connections to Prague, Brno, Bratislava, Budapest, and other large cities. The main railway station, operated by Czech Railways, is vital for passenger traffic and freight transport. Its passenger hall became a place where homeless people concentrate. This railway infrastructure was bombed several times during World War II and many reconstruction projects had been done. Additionally, few significant accidents occurred in the past.
(A2) Hradec Kralove (
Figure 3b). It is well-connected with Pardubice. These two railway infrastructures are 21 km away. Hradec Kralove is also an important hub for both passengers and freight. Several significant accidents occurred. Since lots of people pass through its station, the risk of a terroristic attack is possible as well.
(
A3) Kolin (
Figure 3c). This railway infrastructure is located between the capital of the Czech Republic Prague and Pardubice. It is of huge importance for railway transportation in the Czech Republic as a whole, due to a high frequency of trains on the way to Prague. Several years ago, on the railway infrastructure in Kolin, illegal driving was noticed, but luckily there were no injuries and damage declared.
(A4) Usti nad Labem (
Figure 3d). It is located near the border with Germany. This railway infrastructure counts many departures and arrivals of the local and international trains. Its railway station has been operable for more than sixty years. Due to its location, this railway infrastructure is exposed to various risks such as epidemics, terroristic attacks, and similar accidental events. In the past, some traffic accidents happened. For instance, a freight train passed a signal prohibiting driving and entered the built train path for a departing passenger train running along the adjacent track in the same direction. The drivers of both trains stopped in front of the switch, with the sides of the trains approaching each other less than a meter.
(A5) Chomutov (
Figure 3e). The railway infrastructure in Chomutov is located in the north-west part of the Czech Republic. It is 86 km away from Prague and 50 km away from Usti nad Labem. Recently, this railway infrastructure had several significant accidents. In 2010, a freight train knocked down a woman and seriously injured her. In 2016, two people died after being hit by a freight train. The accident occurred at a railway crossing secured by a light warning signal without barriers. According to the information of the Railway Inspectorate, it was in operation at the time of the accident.
(A6) Jaromer (
Figure 3f). This railway infrastructure connects Hradec Kralove with Pardubice. Some causes of accidents that happened in the past are walking on the track, passing a signal prohibiting driving, derailment of a train on the ground that was washed up on the track after torrential rain, and track defects.
A systematic approach is carried out to identify relevant risk factors influencing railway infrastructure in the Czech Republic from the literature [
74,
75]. Risk factors that do not apply to the Czech Republic context are not listed in
Table 3. Each risk is briefly described in
Table 3.
6. Results and Discussion
6.1. Experimental Results
Six assessed railway infrastructures in the Czech Republic context are “Pardubice” (
A1), “Hradec Kralove” (
A2), “Kolin” (
A3), “Usti nad Labem” (
A4), “Chomutov” (
A5), and “Jaromer” (
A6). Twenty-three risk factors influencing railway infrastructure are considered (
Table 3). Ten railway planners participated in the case study (
Table 4).
The online questionnaire approach via Google Forms was utilized to collect evaluations of risk factors and railway infrastructures due to the COVID-19 outbreak. In the first part of the online questionnaire, the invited railway planners are asked to enter their e-mail address, name, occupation, qualifications, total years of experience, and gender. In the second part, the experts are asked to evaluate twenty-three risk factors influencing railway infrastructure in the Czech Republic. In the last part, the participants are requested to provide linguistic evaluations of six investigated railway infrastructures with respect to the risk factors.
Phase 1. Information collection and representation.
Step 1.1. The invited railway planners provided evaluations of six railway infrastructures in the Czech Republic with respect to the risk factors via the online questionnaire. Importance evaluations can be yes, abstain, no, and refusal. Group refusal of voting either is invalid voting papers or does not take the vote. As a result, six linguistic evaluation matrices are constructed with the help of Equation (13). Linguistic evaluations given by 10 relevant railway planners towards the railway infrastructures with respect to twenty-three risk factors are presented in Table 5. Step 1.2. A picture fuzzy evaluation of railway infrastructure with respect to a risk factor is calculated as the proportion of each out of four items (i.e., yes, abstain, no, or refusal) in the corresponding online voting results of 10 questioned railway planners. The picture fuzzy evaluation matrix is given in Table 6. It is determined based on six linguistic evaluation matrices (Table 5) by using Equation (14). Step 1.3. Twenty-three risk factors influencing railway infrastructure are evaluated by the relevant experts (Table 7). A sample of the online questionnaire to determine the importance of the risk factors is given in the Appendix A. Ten linguistic risk factor importance matrices are constructed by using Equation 15. Step 1.4. A picture fuzzy evaluation of a risk factor is computed as the proportion of each out of four items in the corresponding results of the online questionnaire. The picture fuzzy risk factor importance matrix is provided in Table 8. It is determined based on ten linguistic risk factor importance matrices (Table 7) with the help of Equation (16). Phase 2. Picture fuzzy hybrid method for risk factor prioritization.
Step 2.1. The subjective importance and ranks of twenty-three risk factors are given in Table 9. The values are determined based on the railway planners’ voting on the risk factors influencing railway infrastructure (Table 8) by using the picture fuzzy direct rating method defined in Equation (17). Step 2.2. The objective importance and ranks of the risk factors are presented in Table 9. The values are calculated based on the railway planners’ voting on six railway infrastructures (Table 6) by using the picture fuzzy Tsallis–Havrda–Charvát entropy method defined in Equation (18). In the base case scenario, it is adopted that the information measure parameter χ of the developed picture fuzzy group MCDM approach for risk assessment of railway infrastructure is 1.4. Step 2.3. The hybrid importance and ranks of the risk factors influencing railway infrastructure are provided in Table 9. The values are computed by using Equation (19). In the base case scenario, it is assumed that the value of the trade-off parameter γ of the formulated picture fuzzy group MCDM approach is 0.5. This value equally appraises the picture fuzzy direct rating method and the picture fuzzy Tsallis–Havrda–Charvát entropy method. Phase 3. Picture fuzzy MARCOS method for railway infrastructure ranking.
Step 3.1. Anti-ideal and ideal railway infrastructures are determined by employing Equations (22) and (23), respectively. The picture fuzzy evaluation matrix (Table 6) is extended by adding the anti-ideal and ideal railway infrastructures. The resulting matrix is given in Table 10. Step 3.2. The picture fuzzy additive relative importance of the anti-ideal and ideal railway infrastructures, as well as six real-life railway infrastructures, are calculated by applying the picture fuzzy Hamacher weighted average operator defined in Equation (24). The obtained values can be found in Table 11. In the base case scenario, it is assumed that the operational parameter ζ of the developed picture fuzzy group MCDM approach for risk assessment of railway infrastructure is 3. Step 3.3. The utility degrees to the anti-ideal and ideal solutions of railway infrastructures are given in Table 11. They are calculated based on the picture fuzzy additive relative importance of railway infrastructures with the help of Equations (25) and (26). Step 3.4. According to Equation (27), the utility function to the anti-ideal solution is 0.708. On the other hand, the utility function to the ideal solution is 0.292. This value is obtained by using Equation (28).
Step 3.5. Table 11 provides the utility function for the assessed railway infrastructures. The values are calculated with the help of Equation (29). Step 3.6. Six railway infrastructures in the Czech Republic context are ranked according to the increasing values of the utility functions (Table 11). The ranking order is A6 (Jaromer) ≻ A2 (Hradec Kralove) ≻ A5 (Chomutov) ≻ A3 (Kolin) ≻ A1 (Pardubice) ≻ A4 (Usti nad Labem). As a result, the presented picture fuzzy group MCDM approach for risk assessment of railway infrastructure ranked “Jaromer” (A6) as the safest. Additionally, it is found that “Usti nad Labem” (A4) has the highest risk in the Czech Republic context. 6.2. Sensitivity Analyses
The sensitivity analyses to changes in the trade-off and operational parameters of the presented picture fuzzy group MCDM approach are performed to check the robustness of the obtained railway infrastructure risk assessment results.
The first sensitivity analysis thoroughly explores how the trade-off parameter γ of the novel picture fuzzy hybrid method influences the ranking order of the railway infrastructures in the Czech Republic context. This method based on the direct rating and Tsallis–Havrda–Charvát entropy prioritizes the risk factors influencing railway infrastructure. More detailed, the subjective importance of a risk factor is based on importance evaluations of that risk factor by the invited railway planners. Its value is computed by using the picture fuzzy direct rating method. On the other hand, the objective importance of a risk factor is based on the railway planners’ voting on six railway infrastructures. Its value is calculated by using the picture fuzzy Tsallis–Havrda–Charvát entropy method. Since different input data are used to obtain the subjective and objective ranking orders as well as the number of considered risk factors is large, significant differences in the ranking order of the risk factors are expected.
The parameter
γ allows practitioners to make a trade-off between the picture fuzzy direct rating and Tsallis–Havrda–Charvát entropy methods; i.e., a compromise amongst subjective and objective importance of the risk factors. In this regard, the parameter
γ is changed in the range [0, 1] with an increment value of 0.1 (
Figure 4). When
γ = 1, only the picture fuzzy direct rating method is applied to subjectively prioritize the risk factors. When
γ = 0, the picture fuzzy Tsallis–Havrda–Charvát entropy method is solely used to objectively evaluate the risk factors. Therefore, in the base case scenario,
γ was set to 0.5 to equally appraise both methods and generate hybrid risk factor importance. According to
Figure 4, “Jaromer” (
A6) is the best railway infrastructure under all γ values, since its utility function has the lowest value. Moreover, there is no change in the ranks of any railway infrastructure in all 10 new test cases; i.e., the railway infrastructure ranking order is
A6 ≻
A2 ≻
A5 ≻
A3 ≻
A1 ≻
A4. As a result, it is identified that the railway infrastructure risk assessment results of the investigated real-life context are very stable to changes in the trade-off parameter
γ.
The second sensitivity analysis is performed to check the robustness of generated solutions to changes in the operational parameter
ζ of the developed picture fuzzy group MCDM approach for risk assessment of railway infrastructure (
Figure 5). When
ζ = 1, the Hamacher T-norm and T-conorm reduce to probabilistic T-norm and T-conorm, respectively. When
ζ = 2, the Hamacher T-norm and T-conorm reduce to Einstein T-norm and T-conorm, respectively. As a result, in the base case scenario, it was assumed that
ζ is 3. Ten additional test cases are created to systematically analyze the influence of the operational parameter; i.e.,
ζ ∈ {0.5, 1, 1.5, 2, 2.5, 3.5, 4, 4.5, 5, 5.5}. In all additional test cases, the ranking order of the railway infrastructures in the Czech Republic context is
A6 ≻
A2 ≻
A5 ≻
A3 ≻
A1 ≻
A4 (
Figure 5). This result was obtained in the base case scenario. It is found that the ranks of the railway infrastructures are very stable to changes in the operational parameter
ζ.
According to the results of both sensitivity analyses, it can be outlined that the formulated picture fuzzy group MCDM approach for risk assessment of railway infrastructure is highly robust.
6.3. Comparative Analysis
The comparative analysis is performed to check the reliability of the presented results of the developed picture fuzzy group MCDM approach for risk assessment of railway infrastructure. The railway infrastructure risk assessment problem in the Czech Republic context is solved with five available state-of-the-art picture fuzzy approaches: (1) picture fuzzy TOPSIS [
76], (2) picture fuzzy EDAS [
77], (3) picture fuzzy grey relational analysis [
78], (4) picture fuzzy grey relational projection [
79], and (5) picture fuzzy cross-entropy [
80].
The comparison results are presented in
Table 12. The developed picture fuzzy group MCDM approach and picture fuzzy grey relational projection generate the identical ranking order of the railway infrastructures; i.e.,
A6 (Jaromer) ≻
A2 (Hradec Kralove) ≻
A5 (Chomutov) ≻
A3 (Kolin) ≻
A1 (Pardubice) ≻
A4 (Usti nad Labem). On the other hand, the picture fuzzy TOPSIS, grey relational analysis, and cross-entropy produce the same ranking order; i.e.,
A6 (Jaromer) ≻
A5 (Chomutov) ≻
A2 (Hradec Kralove) ≻
A3 (Kolin) ≻
A1 (Pardubice) ≻
A4 (Usti nad Labem). Six compared approaches assess “Jaromer” (
A6) as the safest railway infrastructure of the analyzed real-life context. The proposed approach, picture fuzzy grey relational projection and picture fuzzy EDAS, rank “Hradec Kralove” (
A2) as the second-best railway infrastructure. The other three approaches put this railway infrastructure in the third position and evaluate “Chomutov” (
A5) as the second-best. “Kolin” (
A3) holds the fourth position in five out of six approaches. Only the picture fuzzy EDAS puts this railway infrastructure in the third place. “Pardubice” (
A1) is the second-worst railway infrastructure in five out of six approaches. Finally, “Usti nad Labem” (
A4) is the worst-ranked in all approaches.
According to the results of the comparative analysis, it can be outlined that the developed approach for risk assessment of railway infrastructure is highly reliable.
The ranking similarity between six compared approaches is examined by Spearman’s rank correlation coefficient. This quantitative metric reveals the strength of the relationship between compared approaches (
Table 13).
The picture fuzzy group MCDM approach for risk assessment of railway infrastructure has 93.2% of ranks matched. Additionally, there is a perfect relationship between the formulated approach and picture fuzzy grey relational projection. Finally, a very strong correlation exists between the introduced approach and picture fuzzy TOPSIS, EDAS, grey relational analysis, and cross-entropy.
According to the results of Spearman’s rank correlation analysis (
Table 13), it can be outlined that the introduced approach for risk assessment of railway infrastructure produces highly consistent results.
7. Conclusions
The picture fuzzy group MCDM approach for risk assessment of railway infrastructure is introduced in this paper. Its major contributions are: (i) For the first time, PFSs are employed for representing railway planners’ preferences and handling risk-related information; (ii) the novel picture fuzzy hybrid method based on the direct rating and Tsallis–Havrda–Charvát entropy is provided to prioritize risk factors influencing railway infrastructure; (iii) the new picture fuzzy MARCOS method is developed to rank railway infrastructures; and (iv) the formulated approach is implemented in the Czech Republic context.
The merits of the presented real-life case study are: (a) The utility of the introduced approach is demonstrated; (b) the high robustness of the formulated approach is verified by two sensitivity analyses, since the ranks of the railway infrastructures are very stable to changes in the trade-off and operational parameters; (c) the high reliability of the developed approach is approved by the comparative analysis with the picture fuzzy TOPSIS, EDAS, grey relational analysis, grey relational projection, and cross-entropy; and (d) the high consistency with five state-of-the-art picture fuzzy approaches is confirmed by Spearman’s rank correlation analysis, since 93.2% of railway infrastructure ranks are matched.
The picture fuzzy group MCDM approach for risk assessment of railway infrastructure generates the following ranking order: A6 (Jaromer) ≻ A2 (Hradec Kralove) ≻ A5 (Chomutov) ≻ A3 (Kolin) ≻ A1 (Pardubice) ≻ A4 (Usti nad Labem). This novel approach identified “Jaromer” as the safest. On the other hand, the worst-ranked railway infrastructure is “Usti nad Labem”. It has the highest risk in the Czech Republic context. As a result, it is strongly recommended to undertake a safety improvement project for this railway infrastructure.
The picture fuzzy grey relational analysis and grey relational projection are one-parametric approaches with build-in distinguishing coefficients. The picture fuzzy TOPSIS, EDAS, and cross-entropy have no parameters. On the other hand, the formulated approach involves three intrinsic parameters, which is highly desirable for solving the complex railway infrastructure risk assessment problem. They are the trade-off parameter γ, the operational parameter ζ, and the information measure parameter χ. As a result, compared to the state-of-the-art picture fuzzy approaches, the provided three-parametric approach has superior flexibility in assessing the risk of railway infrastructure.
Limitations of this paper can indicate its possible extension areas. The limitations are: (1) The risk factors influencing railway infrastructure are not filtered; (2) interrelationships between the risk factors are mainly ignored. A well-known technique for filtering influential factors is the Delphi method. It is traditionally used to obtain a consistent flow of answers through the results of questionnaires. Its major features are anonymous response, iteration and controlled feedback, and finally statistical group response. In a future study, an online multi-round questionnaire approach, in line with the Delphi method, could be performed to collect risk-related information and overcome the first limitation. On the other hand, the DEMATEL method is mostly applied for modeling interrelationships of influential factors. This method has not been extended before into the picture fuzzy environment. One of the future researches may integrate a picture fuzzy DEMATEL method into the introduced methodological framework for risk assessment to handle the second limitation. Additionally, a comparison of the risk factors influencing railway infrastructure in the Czech Republic with risk factors proposed by other researchers can also be seen as an interesting topic that deserves a future research effort. Finally, the introduced methodological approach can be used not only for risk assessment of railway infrastructure but also in macro-level issues. Indicatively, various risks exist in different frames, due to climate change, like fires, floods, droughts, hurricanes, and tornadoes, which affect social and economic life seriously.