This section follows the steps taken to collect the data for assessing the likelihood and consequences of complexity drivers and risk factors in relation to the cost overrun. In this study, experts’ subjective judgment is recorded through personal interviews from construction experts, who have been associated with metropolitan construction projects in under-developing countries. For the recording of subjective judgments, ten experts as DMs from a different area of expertise with more than fifteen years’ construction experience have been nominated to participate in the decision-making process by the project authorities. Subjective evaluation of complexity and risk related to cost overrun has been recorded in three rounds according to the requirement of the designed method. In the first round, data is collected to develop a hierarchical breakdown structure of potential complexity drivers and risk factors for cost overrun on a binary scale. In the second round, the linguistic scale is used to record the likelihood and occurrence of complexity and risk events. In the last round, prior and conditional probabilities are recorded for developing Bayesian inference of complexity-risk interdependencies along with the expected variation in project cost.
This research process contains the following steps in order to find complexity-risk breakdown structure and development of Bayesian inference for cost overrun.
Step 1: Identification of potential complexity and risk factors to form a hierarchical complexity-risk breakdown structure for cost overrun.
Step 2: Selection of suitable linguistic scales to record experts’ subjective judgments for both likelihood and consequences of complexity-risk events.
Step 3: Transformation of linguistic data into an appropriate trapezoidal fuzzy function in accordance with a recommended fuzzy scale set by the DMs. In the fuzzy decision matrix process, a collective opinion of DMs is obtained using fuzzy aggregation rule and fuzzy control.
Step 4: Development of BBN for complexity-risk interdependencies based on the fuzzy triangular distribution probability values of key prior complexity and conditional risk factors.
Step 5: For the analysis and re-evaluation process, project cost data against each complexity-driven risk criteria is collected across different construction projects, for sensitivity measurement. Three-point joint estimates of cost overrun (i.e., low, medium, and high) are determined against three estimates (i.e., pessimistic, most likely, and optimistic) of important risks that directly impact on cost overrun assuming complexity-risk interdependences.
The specific identity of experts involved in the decision-making process is not supposed to reveal due to the reason of anonymity. During the interview session, experts have been requested to record their judgments on a prescribed format of the questions. After getting an initial response on a binary scale, a structured framework of complexity drivers and risk factors is drawn for further work.
3.1. Identification of Complexity and Risk Breakdown Structure
After intensive literature review, different complexity elements and complexity-driven risk factors have been transformed into binary states of ‘True (T)’ or ‘False (F)’ for selection during interview process (Qazi et al. 2016). Afterwards, based on the experts’ aggregate opinion, a breakdown structure of complexity-risk elements is obtained for the prioritization process. Tables 1 and 2 present a total sixteen complexity elements and twenty risk factors from five different risk sources, i.e., engineering design, construction management, construction safety-related, natural hazards, and social and economic, as defined by Hastak and Shaked (2000), Dikmen et al. (2007), Fang et al. (2012), Qazi et al. (2016) and Samantra et al. (2017).
Table 1
List of Complexity Elements for Cost Chaos in Mega Construction Projects
Complexity dimensions | ID | Complexity drivers under specific dimensions | Reference |
Engineering design (D1) | CG1 | Lacking clarity and misalignment of goals | Zidane and Andersen, 2018 |
PS2 | Ambiguity in defining project scope | Eybpoosh et al., 2011; Ahmadi et al., 2017 |
QS3 | Quality standard | Kabir et al., 2016; Iqbal et al., 2015 |
CS4 | Conflicting standards | Shehu et al., 2014; Jato-Espino et al., 2014 |
IP5 | Inadequate control procedures | Cho and Eppinger, 2005 |
Construction management (D2) | TM6 | Uncertainty in technical methods | Arashpour et al., 2017; Samantra et al., 2017 |
IT7 | Innovative technology | Eybpoosh et al., 2011 |
LS8 | Lacking skill with the technology in use | Wang, 2011; Dikmen et al., 2010 |
LE9 | Lacking experience with project team involved | Floyd et al., 2017 |
CF10 | Cash flow during construction | Eybpoosh et al., 2011 |
BF11 | Recurring breakdowns of construction facilities and plan | Ou-Yang and Chen, 2017; Ou-Yang and Chen, 2017 |
Economic and social (D3) | SP12 | Number of stakeholders and their perspectives | Zayed et al., 2008; Afzal et al., 2018 |
PC13 | Political condition | Kim et al., 2009 |
MC14 | Market competition | Dikmen et al., 2010 |
CB15 | Corruption and bribery | Fazekas and Tóth, 2018; Tahir et al., 2011 |
GC16 | Geological conditions | Barakchi et al., 2017; Samantra et al., 2017 |
Table 2
List of Risk Factors for Cost Chaos in Mega Construction Projects
Risk dimensions | ID | Risk factors under specific dimensions | Reference |
Engineering design (D1) | PD1 | Inappropriate project designing and poor engineering process | Eybpoosh et al., 2011 |
DE2 | Design drawing errors | Dikmen et al., 2010; Eisenhardt and Graebner, 2007 |
IW3 | Inconsistency in work items | Jato-Espino et al., 2014 |
SI4 | Poor construction site inspections | Doloi et al., 2012 |
Construction management (D2) | CP5 | The poor construction planning process | Terstegen et al., 2016; Cho and Eppinger, 2005 |
ES6 | Insufficient experience and skill in construction works | Dikmen et al., 2010 |
DF7 | Delay in relocating existing facilities | Lazzerini and Mkrtchyan, 2011; Budayan et al., 2018 |
SM8 | Unstable supply of construction materials required | Boateng et al., 2015 |
Construction safety-related (D3) | PI9 | Inadequate protection of nearby infrastructure and facilities | Zhang et al., 2016 |
SP10 | Inadequate safety procedures of worker | Tahir et al., 2011; Shafiee, 2015 |
PE11 | Ineffective protection of the environment | Eybpoosh et al., 2011; Camós et al., 2016 |
TC12 | Mismanagement of traffic control | Satiennam et al., 2006 |
Natural hazards (D4) | HR13 | Heavy rainfall during construction | Wang et al., 2016 |
SC14 | Super cyclonic storm | Wang et al., 2016 |
EQ15 | Earthquake | Chang, 2014 |
WS16 | Ground water seepage problem | Samantra et al., 2017 |
Social and economic (D5) | PI17 | Political interference | Boateng et al., 2015; Valipour et al., 2015 |
PM18 | Increases in prices of critical construction materials | Senouci et al., 2016; Fazekas and Tóth, 2018 |
WS19 | Increases in workers’ salaries | Eybpoosh et al., 2011; Doloi et al., 2012 |
NR20 | Protest of nearby residents | Tsavdaroglou et al., 2018 |
[Insert Table 1]
[Insert Table 2]
3.2. Linguistic Scales and Fuzzy Membership Function
Under the condition of uncertainty in subjective risk information, getting exact data of likelihood of occurrence and consequences in relations with the complexity and risk events, the study necessitates the support of group decision to record experts’ subjective judgments on a linguistic scale (Kabir et al. 2016; Mehlawat and Gupta 2016). For the effectiveness of subjective judgments in the assessment process, a linguistic scale has been designed to measure the semantic strength of an event, as suggested by Rezakhani et al. (2014) and Samantra et al. (2017). In addition, by following Eq. 5, fuzzy numbers have been developed to represent a set of linguistic measurement values for each property of complexity and risk event (Fouladgar et al. 2012; Zhang et al. 2017).
Although many studies have followed different types of fuzzy membership functions for solving linguistic-based problems (Chang 2014; Russo and Camanho 2015), generalized fuzzy trapezoidal membership function has been selected for subjective assessments due to the possible optimized constraints of other fuzzy functions (Prascevic and Prascevic 2017). In this research, the likelihood and consequences for both complexity and risk events have been quantified by using linguistic scales of seven attributes (Samantra et al. 2017) and five attributes (Amiri and Golozari 2011), respectively, with corresponding fuzzy membership function (Prascevic and Prascevic 2017), as labelled in Tables 3 and 4, and in Fig. 1 as well.
Table 3
Fuzzy Linguistic Scale for Assessing Likelihood of Risk or Complexity
Linguistic scale | Description | Crisp values | Fuzzy membership function\({ ({\delta }}_{\tilde{{F}}}\left({Z}\right)\)) | Fuzzy numbers |
Very rare (VR) | It can be assumed that the possibility of occurrence is negligible | 1 | LVL = \(\tilde{F}\)(VRL, VRM, VRN, VRU) | (0, .1, .1, .2; 1) |
Rare (R) | Unlikely but possible to occur an event in operation | 2 | LR = \(\tilde{F}\)(RL, RM, RN, RU) | (.1, .2,.2, .3; 1) |
Occasional (O) | Likely to occur event in operation | 3 | LO = \(\tilde{F}\)(OL, OM, ON, OU) | (.2, .3, .3, .4; 1) |
Probable (P) | Likely to occur event several times in operation | 4 | LP = \(\tilde{F}\)(PL, PM, PN, PU) | (.3, .4, .5, .5; 1) |
Frequent (F) | Likely to occur frequently | 5 | LF = \(\tilde{F}\)(FL, FM, FN, FU) | (.5, .6, .6 .7; 1) |
Very frequent (VF) | Much frequent to occur | 6 | LVF = \(\tilde{F}\)(VFL, VFM, VFN, VFU) | (.6, .7, .7, .8; 1) |
Absolutely certain (AC) | Expected to occur with absolute certainty | 7 | LAC = \(\tilde{F}\)(ACL, ACM, ACN, ACU) | (.8,.9, .9, 1; 1 ) |
Table 4
Fuzzy Linguistic Scale for Assessing Occurrence of Risk or Complexity
Linguistic scale | Description | Crisp values | Fuzzy membership function\({ ({\delta }}_{\tilde{{F}}}\left({Z}\right)\)) | Fuzzy numbers |
Very low (VL) | It can be assumed that the magnitude if impact is negligible. | 1 | OVL = \(\tilde{F}\)(VLL, VLM, VLN, VLU) | (0, .1, .1, .2; 1) |
Low (L) | Possible to manage impact of event in operation | 2 | OL = \(\tilde{F}\)(LL, LM, LN, LU) | (.1, .2,.2, .3; 1) |
Moderate (M) | Likely to encounter an impact of risk in operation | 3 | OM = \(\tilde{F}\)(ML, MM, MN, MU) | (.2, .3, .3, .4; 1) |
High (H) | Significant impact of event in operation | 4 | OH = \(\tilde{F}\)(HL, HM, HN, HU) | (.3, .4, .5, .5; 1) |
Very high (VH) | Absolute certain impact of an event in operation | 5 | OVL = \(\tilde{F}\)(VHL, VHM, VLN, VHU) | (.5, .6, .6 .7; 1) |
[Insert Table 3]
[Insert Table 4]
[Insert Fig. 1]
3.3. Fuzzy Decision Matrix Process
-
During the decision-making process, using linguistic information of Tables 3 and 4, a fuzzy decision matrix, \({DM}_{C-RL/C}^{h}\) for an individual expert h (1, 2, 3,…., g) regarding the likelihood (L) and consequences (C) of complexity-risk elements is constructed separately as articulated in Eq. 6. Decision matrix addresses the states of n complexity influencing factors (C1, C2,…, Cn) and k risk factors (R1, R2,…, Rk) under m number of dimensions (D1, D2,…, Dm) by a group of g number of experts Hg. In complexity decision matrix\({DM}_{CL/C}^{h}\), \({\tilde{F}}_{ij}\) (i = 1, 2,…, m and j = 1, 2,…, n) is specified separately like comparing the states of L and C for each criterion Cm against a group of DMg. Similarly, in the risk matrix\({DM}_{RL/C}^{h}\), \({\tilde{F}}_{ij}\) (i = 1, 2,…, m and j = 1, 2,…, n) is specified separately for both states L and C for each criterion Cm against a group of DMg. Decision matrices \({DM}_{C-RL/C}^{r}\) for likelihood and consequence of complexity-risk criteria against individual DM have been expressed herein in Eq. 6, as recommended by Islam and Nepal (2016), Kabir et al. (2016) and Khodakarami and Abdi (2014).
\(\begin{array}{*{20}{c}} {D{M_1}}&{D{M_2}}& \cdots &{D{M_g}} \end{array}\)
$$DM_{{C - RL/C}}^{r}=\begin{array}{*{20}{c}} {{C_{1 \times 1}}} \\ {{C_{2 \times 1}}} \\ \vdots \\ {{C_{m \times 1}}} \end{array}\left[ {\begin{array}{*{20}{c}} {{{\widetilde {F}}_{1 \times 1}}}&{{{\widetilde {F}}_{1 \times 2}}}& \cdots &{{{\widetilde {F}}_{1 \times g}}} \\ {{{\widetilde {F}}_{2 \times 1}}}&{{{\widetilde {F}}_{2 \times 2}}}& \cdots &{{{\widetilde {F}}_{2 \times g}}} \\ \vdots & \cdots & \ddots & \vdots \\ {{{\widetilde {F}}_{m \times 1}}}&{{{\widetilde {F}}_{m \times 2}}}& \cdots &{{{\widetilde {F}}_{m \times g}}} \end{array}} \right]$$
6
-
Involving multiple experts create the decision-making process more complex and uncertain (Islam and Nepal 2016; Kabir et al. 2016). As the DMs belong to different designations, experience and qualification, therefore, their opinion holds different weight credibility in the decision process. The expert’s credibility factor is measured using \({W}_{h}\) as discussed by Kabir et al. (2016) and based on an expert’s experience \({E}_{h}\in \left[\text{0,1}\right]\), qualification \({Q}_{h}\in \left[\text{0,1}\right]\), and designation \({D}_{h}\in \left[\text{0,1}\right]\), particularly at risk management in the construction domain. In this paper, the formulation is modified using a normalization factor\({ \text{M}\text{A}\text{X}}_{h=1}^{H}{\left({D}_{h}{Q}_{h}{E}_{h}\right)}^{\beta }\). Table 5 describes the general profile of the experts and their weight criteria involved in the decision-making process. The normalization weight vector is computed such that the total credible weight of an expert is considered as 1 and remaining weights of DMs are\({ W}_{k}\le 1\). Therefore, for K number of experts, the credibility factor (\({W}_{k}\)) is derived as:
Table 5
General Profile of Experts and their Weight Criteria
Respondents | Designation | Experience | Qualification | DM’s weights | N-weights |
01 | Principal Engineer | 26 | Post-Graduate | 0.8 | 0.127 |
02 | Senior Engineer | 18 | Graduation | 0.5 | 0.079 |
03 | Principal Engineer | 18 | Graduation | 0.4 | 0.063 |
04 | Senior Engineer | 20 | Post-Graduate | 0.5 | 0.079 |
05 | Principal Engineer | 23 | Ph.D. | 0.6 | 0.095 |
06 | Principal Engineer | 20 | Post-Graduate | 0.5 | 0.079 |
07 | Senior Engineer | 24 | Graduation | 0.5 | 0.079 |
08 | Senior Engineer | 24 | Graduation | 0.7 | 0.111 |
09 | Principal Engineer | 28 | Post-Graduate | 0.8 | 0.127 |
10 | Chief Engineer | 30 | Post-Graduate | 1 | 0.159 |
$${W_h}=\frac{{{{\left( {{D_k}{Q_k}{E_k}} \right)}^\beta }}}{{MAX_{{k=1}}^{K}{{\left( {{D_k}{Q_k}{E_k}} \right)}^\beta }}}$$
7
Where β is a weight vector used to assign the weight of individual DM; therefore, the higher β value shows the dominance of DMg in assessment with higher \({D}_{h},{Q}_{h}\) and\({ E}_{h}\).
[Insert Table 5]
-
The fuzzy numbers in the decision matrix \({DM}_{C-RL/C}^{r}\)are multiplied by the weight score of the respective individual expert (wh). Afterwards, fuzzy multiplication rule as shown in Eq. 8 has been applied to get weighted elements of the matrix.
$${\left( {DM_{{C - RL/C}}^{h}} \right)_w}=~{w_i} \otimes {\left( {\bar {F}_{{ij}}^{h}} \right)_L}~,~~{w_i} \otimes {\left( {\bar {F}_{{ij}}^{h}} \right)_M}~~{w_i} \otimes {\left( {\bar {F}_{{ij}}^{h}} \right)_N}~~{w_i} \otimes {\left( {\bar {F}_{{ij}}^{h}} \right)_U}$$
8
Here, L, M, N and U mean the lowest, low moderate, moderate and highest possible number of fuzzy trapezoidal function, respectively, and the symbol ⊗ indicates fuzzy multiplication rule. All the matrices of individual DMs are transformed into one single matrix following the fuzzy arithmetic average (Islam and Nepal 2016) by using Eq. 9.
Elements of group matrix =
$$\bar {F}_{{ij}}^{G}=~\frac{1}{k}~\mathop \sum \limits_{{h=g}}^{{k=1}} \bar {F}_{{ij}}^{h} \otimes {w_h}=~{\left( {\frac{1}{k}~\mathop \sum \limits_{{h=1}}^{k} \bar {F}_{{ij}}^{h} \otimes {w_h}} \right)_L}~,~{\left( {\frac{1}{k}~\mathop \sum \limits_{{h=1}}^{k} \bar {F}_{{ij}}^{h} \otimes {w_h}} \right)_M}~{\left( {\frac{1}{k}~\mathop \sum \limits_{{h=1}}^{k} \bar {F}_{{ij}}^{h} \otimes {w_h}} \right)_N}~{\left( {\frac{1}{k}~\mathop \sum \limits_{{h=1}}^{k} \bar {F}_{{ij}}^{h} \otimes {w_h}} \right)_U}$$
9
-
The elements of the group matrix are defuzzified using the following Eq. 10. The defuzzification method–Center of Area (Hefei 2017; Yazdi and Kabir 2017)–is used to calculate the absolute values of likelihood and consequences.
$$Center{\text{ }}of{\text{ }}Area{\text{ }}\left( {COA} \right){\text{ }}={{{w}}_{{i}}}=\left( {\frac{{{{a}}+{{b}}+{{c}}+{{d}}}}{4}} \right)$$
10
-
Finally, fuzzy if-then rules between the likelihood and consequence presented in Table 6 are applied to finding the probability level of risk for Bayesian inference (Cárdenas et al. 2014; Kabir et al. 2016; Yazdi and Kabir 2017).
Table 6
Fuzzy Control Rules between Likelihood and Consequence
Level of severity | Consequences |
VL | L | M | H | VH |
Likelihood | VR | VR | VR | VR | R | R |
R | VR | R | R | R | O |
O | VR | R | O | O | P |
P | VR | O | O | P | F |
F | VR | O | P | P | VF |
VF | R | O | P | F | AC |
AC | R | P | F | VF | AC |
Furthermore, Bayesian inference is developed to find out the entire posterior probabilities of the complexity event and all intermediate-risk nodes of cost overrun.
3.4. Development of DAG Bayesian Belief Network
It is obvious that decision outcomes are dependent on the risk-bearing attitude of the domain experts. Though the linguistic probability scores of each identified complexity-risk factors may tend to alter when experts’ mental attitude is considered (Islam and Nepal 2016; Kabir et al. 2016). Therefore, considering experts’ different risk-bearing attitudes, i.e., pessimistic (P), most likely (ML) and optimistic (O), the integration of fuzzy scores into Bayesian inference has become an effective approach in order to analyze complexity-risk interdependencies (Islam et al. 2017; Yazdi and Kabir 2017). Tables 7 and 8 show the computation of linguistic probability scores of each identified complexity-risk factors with experts’ different risk-bearing attitudes.
Table 7
Prior Probabilities Scores of Complexity Elements along with Fuzzy Values
ID | Likelihood | Occurrence | Probability |
C.S | L.V | C.S | L.V | L.V | P | ML | O | C.S |
CG1 | 0.59 | F | 0.19 | VL | VR | 0 | 0 | 0.1 | 0.02 |
PS2 | 0.56 | F | 0.21 | L | O | 0.3 | 0.4 | 0.5 | 0.40 |
QS3 | 0.49 | P | 0.31 | H | P | 0.5 | 0.6 | 0.7 | 0.60 |
CS4 | 0.47 | P | 0.25 | L | O | 0.3 | 0.4 | 0.5 | 0.40 |
IP5 | 0.65 | F | 0.19 | VL | VR | 0 | 0 | 0.1 | 0.02 |
TM6 | 0.57 | F | 0.19 | VL | VR | 0 | 0 | 0.1 | 0.02 |
IT7 | 0.45 | P | 0.33 | H | P | 0.5 | 0.6 | 0.7 | 0.60 |
LS8 | 0.39 | O | 0.22 | L | R | 0.1 | 0.2 | 0.3 | 0.20 |
LE9 | 0.50 | P | 0.25 | L | O | 0.3 | 0.4 | 0.5 | 0.40 |
CF10 | 0.59 | F | 0.26 | M | P | 0.5 | 0.6 | 0.7 | 0.60 |
BF11 | 0.42 | P | 0.24 | L | O | 0.3 | 0.4 | 0.5 | 0.40 |
SP12 | 0.54 | F | 0.24 | L | O | 0.3 | 0.4 | 0.5 | 0.40 |
PC13 | 0.69 | VF | 0.28 | M | P | 0.5 | 0.6 | 0.7 | 0.60 |
MC14 | 0.56 | F | 0.28 | M | P | 0.5 | 0.6 | 0.7 | 0.60 |
CB15 | 0.73 | VF | 0.18 | L | O | 0.3 | 0.4 | 0.5 | 0.40 |
GC16 | 0.37 | R | 0.29 | M | R | 0.1 | 0.2 | 0.3 | 0.20 |
Table 8
Prior Probabilities Scores of Risk Factors along with Fuzzy Values
ID | Likelihood | Occurrence | Probability |
C.S | L.V | C.S | L.V | L.V | P | ML | O | C.S |
PD1 | 0.55 | F | 0.11 | VL | VR | 0 | 0 | 0.1 | 0.02 |
DE2 | 0.59 | F | 0.21 | L | O | 0.3 | 0.4 | 0.5 | 0.40 |
IW3 | 0.55 | F | 0.23 | L | O | 0.3 | 0.4 | 0.5 | 0.40 |
SI4 | 0.44 | P | 0.25 | M | O | 0.3 | 0.4 | 0.5 | 0.40 |
CP5 | 0.77 | VF | 0.16 | VL | R | 0.1 | 0.2 | 0.3 | 0.20 |
ES6 | 0.63 | F | 0.19 | VL | VR | 0 | 0 | 0.1 | 0.02 |
DF7 | 0.52 | F | 0.24 | M | P | 0.5 | 0.6 | 0.7 | 0.60 |
SM8 | 0.55 | O | 0.29 | M | O | 0.3 | 0.4 | 0.5 | 0.40 |
PI9 | 0.52 | O | 0.26 | M | O | 0.3 | 0.4 | 0.5 | 0.40 |
SP10 | 0.55 | P | 0.33 | H | P | 0.5 | 0.6 | 0.7 | 0.60 |
PE11 | 0.49 | O | 0.26 | M | O | 0.3 | 0.4 | 0.5 | 0.40 |
TC12 | 0.53 | P | 0.25 | M | O | 0.3 | 0.4 | 0.5 | 0.40 |
HR13 | 0.27 | R | 0.26 | M | R | 0.1 | 0.2 | 0.3 | 0.20 |
SC14 | 0.23 | R | 0.36 | H | R | 0.1 | 0.2 | 0.3 | 0.20 |
EQ15 | 0.30 | R | 0.27 | M | R | 0.1 | 0.2 | 0.3 | 0.20 |
WS16 | 0.39 | O | 0.41 | H | O | 0.3 | 0.4 | 0.5 | 0.40 |
PI17 | 0.68 | VF | 0.18 | L | O | 0.3 | 0.4 | 0.5 | 0.40 |
PM18 | 0.74 | VF | 0.15 | L | O | 0.3 | 0.4 | 0.5 | 0.40 |
WS19 | 0.58 | F | 0.21 | L | O | 0.3 | 0.4 | 0.5 | 0.40 |
NR20 | 0.47 | P | 0.38 | H | P | 0.5 | 0.6 | 0.7 | 0.60 |
In order to design a DAG-based Bayesian inference, the judgments of DMs are first transformed into fuzzy numbers, which provide a probability of risk occurrence. These probabilities are the input variables of BBN for representing the causal relationships among the complexity-risk elements of cost overrun.
[Insert Table 7]
[Insert Table 8]
The conditional probabilities of complexity-risk interdependencies have been recorded through the above-mentioned decision-making process along with the triangular distribution of cost data, i.e., low, medium and high. During the decision process, experts have been asked to first define complexity-risk interdependencies and then record conditional probability values directly into the network following the prior probability values of complexity elements.
[Insert Fig. 2]
The DAG in Fig. 2 presents interrelationships within complexity elements and risk factors, where complexity is considered as a parent node and risk as a child node, respectively. Consequently, risk factors, such as inappropriate project designing and poor engineering process (PD1), delay in relocating existing facilities (DF7) and increases in prices of critical construction materials (PM18), show high dependency in a network that directly impacts on the cost behaviour.
Bayesian inference in Fig. 3 shows normalized joint probability with complexity-risk interdependencies for cost overrun function calculated by using Equations 2 and 3 on three-point estimations. Utility node shows to address the cost overrun causes that can be controlled because of this complex relationship. Finally, three main risk factors, PD1, PM18 and DF7, are found to be important that reflect the direct impact on cost overrun depending on other posterior complexity elements and risk factors in a network.
[Insert Fig. 3]
[Insert Table 9]
Table 9
Cost Variation against the Probability States of Key Risk Sources (amount in million dollar)
PD1 | Pessimistic |
PM18 | Pessimistic | Most likely | Optimistic |
DF7 | Pessimistic | Most likely | Optimistic | Pessimistic | Most likely | Optimistic | Pessimistic | Most likely | Optimistic |
Low | 1.7 | 3 | 3 | 2.5 | 2.5 | 3 | 4 | 7.5 | 13.5 |
Medium | 3.5 | 4.5 | 5 | 4.5 | 5.5 | 6 | 6.5 | 7.5 | 22 |
High | 7.1 | 8 | 8 | 7.2 | 8 | 8 | 8 | 15 | 28 |
PD1 | Pessimistic |
PM18 | Pessimistic | Most likely | Optimistic |
DF7 | Pessimistic | Most likely | Optimistic | Pessimistic | Most likely | Optimistic | Pessimistic | Most likely | Optimistic |
Low | 7 | 8 | 8 | 8.5 | 11 | 14 | 14 | 20 | 20 |
Medium | 14.5 | 14.5 | 10 | 15 | 16 | 21 | 18 | 22 | 25 |
High | 21 | 21 | 22 | 23 | 23 | 25 | 25 | 28 | 30 |
PD1 | Pessimistic |
PM18 | Pessimistic | Most likely | Optimistic |
DF7 | Pessimistic | Most likely | Optimistic | Pessimistic | Most likely | Optimistic | Pessimistic | Most likely | Optimistic |
Low | 30 | 32 | 35 | 41 | 42 | 45 | 47 | 50 | 54.4 |
Medium | 35 | 35 | 35 | 45 | 60 | 60 | 60 | 65 | 71.5 |
High | 40 | 40 | 42 | 50 | 62 | 65 | 61 | 63 | 78.6 |
Three-point joint estimates of cost overrun (i.e., low, medium, and high) are determined against three estimates (i.e., pessimistic, most likely, and optimistic) of important risks that directly impact on cost overrun assuming complexity-risk interdependences. Table 9 illustrates the variation of cost in dollars against important risk factors that have been found through Bayesian inference. Additional cost required to manage risk within the complexity-risk network in ITPs has found between 1.7 million (in case of pessimistic approach of risk probability) to 78.6 million dollar (in case of optimistic approach of risk probability).
3.5. Simulation Modelling for Cost-Risk Re-evaluation
Risk analysis and re-evaluation
Right after taking the risk circulation conduct into view, the probability of risk is reconsidered and displayed as a numerical risk frequency (Afzal et al. 2020), for a thorough description of the methodology used for simulation modelling. The outcomes of Monte Carlo simulation are useful to estimate the cost-risk for the project, grounded on historic cost data and to calculate the entire costs. In order to re-evaluate the Fuzzy-BBN results, real cost data or each important risk factor is collected from different construction projects and there simulated values are used for further analysis. The overall cost of the project is calculated by the merge of base and risk costs of all several components. The supplementary cost vital for allay of the risk is estimated through a contingency model, for the timely completion of the project (Afzal et al. 2021).
The importance of the highlighted analysis is to know the risk factors of cost overrun that propagate a project into chaos and calculate the necessary amount of the additional cost needed to handle cot-chaos. In the next phase, specialists were requested to measure the total cost of risk for each identified risk in a network. It is clear that risk scores differ based on the experience of the specialist. In addition to this, the cost of a project is based on the risk level.
Furthermore, the present study aimed to evaluate the run over of cost while keeping the risk score into consideration that shows that how much risk is taken by a specialist like pessimistic, most likely, and optimistic. By applying a three-point calculation approach, cost variation has been designated for virtual decision-making (Afzal et al. 2020). Through revaluation process of model, it is validated that the maximum project cost flow is coming from the factors of inappropriate project designing and poor engineering process, delay in relocating existing facilities have caused improper cost management and increases in prices of critical construction materials.
Table 10 summarises cost-benefit contingency values created by each risk factor's simulated scoring and the corresponding cost to ameliorate each risk. The contingency index determines the budgetary allocation needed to reduce the project's risk impact despite of concluding prior to its accomplishment.
Table 10
Simulated Results of Cost-Risk Comparative Analysis
Di | Risk Factors IDs | Pessimistic (P) | Most likely (M) | Optimistic (O) | Re-evaluated simulation results |
Risk α= 0 | Cost α = 0 | Risk α = 0.5 | Cost α = 0.5 | Risk α = 1 | Cost α = 1 | Simulated risk | Simulated cost | Risk SD | Cost SD |
D1 | PD1 | 0.05 | 6.75 | 0.08 | 9.10 | 0.08 | 11.45 | 0.07 | 9.74 | 0.01 | 11.50 |
DE2 | 0.04 | 5.55 | 0.06 | 6.90 | 0.07 | 8.42 | 0.06 | 6.72 | 0.01 | 8.51 |
IW3 | 0.04 | 5.64 | 0.07 | 7.40 | 0.07 | 8.55 | 0.06 | 7.45 | 0.00 | 8.62 |
SI4 | 0.03 | 2.10 | 0.06 | 2.55 | 0.06 | 3.12 | 0.05 | 2.61 | 0.01 | 3.14 |
D2 | CP5 | 0.03 | 2.82 | 0.06 | 3.52 | 0.05 | 4.36 | 0.05 | 3.41 | 0.00 | 4.37 |
ES6 | 0.04 | 5.52 | 0.05 | 6.22 | 0.08 | 7.71 | 0.06 | 6.42 | 0.01 | 7.76 |
DF7 | 0.05 | 6.52 | 0.07 | 8.25 | 0.08 | 10.23 | 0.07 | 8.72 | 0.01 | 10.25 |
SM8 | 0.03 | 4.96 | 0.06 | 6.62 | 0.06 | 8.39 | 0.06 | 6.26 | 0.01 | 8.41 |
D3 | PI9 | 0.01 | 0.36 | 0.03 | 0.55 | 0.05 | 0.68 | 0.02 | 0.38 | 0.01 | 0.72 |
SP10 | 0.02 | 0.17 | 0.04 | 0.26 | 0.04 | 0.27 | 0.03 | 0.44 | 0.00 | 0.29 |
PE11 | 0.01 | 0.09 | 0.03 | 0.15 | 0.04 | 0.12 | 0.02 | 0.20 | 0.00 | 0.14 |
TC12 | 0.02 | 0.22 | 0.03 | 0.35 | 0.05 | 0.36 | 0.02 | 0.40 | 0.01 | 0.37 |
D4 | HR13 | 0.04 | 1.89 | 0.05 | 2.54 | 0.04 | 3.18 | 0.04 | 2.41 | 0.00 | 3.22 |
SC14 | 0.02 | 0.09 | 0.03 | 0.18 | 0.04 | 0.11 | 0.01 | 0.20 | 0.00 | 0.14 |
EQ15 | 0.04 | 0.17 | 0.04 | 0.24 | 0.06 | 0.33 | 0.02 | 0.31 | 0.01 | 0.36 |
WS16 | 0.03 | 0.58 | 0.05 | 0.66 | 0.02 | 0.74 | 0.03 | 0.71 | 0.00 | 0.76 |
D5 | PI17 | 0.04 | 3.83 | 0.07 | 4.53 | 0.06 | 5.65 | 0.04 | 4.66 | 0.01 | 5.64 |
PM18 | 0.05 | 7.58 | 0.08 | 10.12 | 0.09 | 12.05 | 0.08 | 9.95 | 0.01 | 12.11 |
WS19 | 0.02 | 1.82 | 0.03 | 2.55 | 0.08 | 3.08 | 0.05 | 2.48 | 0.00 | 3.14 |
NR20 | 0.03 | 0.12 | 0.04 | 0.19 | 0.02 | 0.14 | 0.01 | 0.15 | 0.00 | 0.16 |
[Insert Table 10]