Human Reliability Analysis (HRA) Using Standardized Plant Analysis Risk-Human (SPAR-H) and Bayesian Network (BN) for Pipeline Inspection Gauges (PIG) Operation: A Case Study in a Gas Transmission Plant


Mostafa Mirzaei Aliabadi ORCID 1 , Reza Esmaeili ORCID 1 , * , Iraj Mohammadfam ORCID 1 , Maryam Ashrafi ORCID 2


1 Center of Excellence for Occupational Health Engineering, Occupational Health and Safety Research Center, School of Public Health, Hamadan University of Medical Sciences, Hamadan, Iran

2 Department of Industrial Engineering and Management Systems, Amirkabir University of Technology, Tehran, Iran

How to Cite: Mirzaei Aliabadi M, Esmaeili R, Mohammadfam I, Ashrafi M. Human Reliability Analysis (HRA) Using Standardized Plant Analysis Risk-Human (SPAR-H) and Bayesian Network (BN) for Pipeline Inspection Gauges (PIG) Operation: A Case Study in a Gas Transmission Plant, Health Scope. Online ahead of Print ; 8(3):e87148. doi: 10.5812/jhealthscope.87148.


Health Scope: 8 (3); e87148
Published Online: July 28, 2019
Article Type: Research Article
Received: December 7, 2018
Revised: March 9, 2019
Accepted: April 13, 2019




Background: Pigging operation is one of the maintenance activities that is used to check pipeline functionality in operational conditions using a PIG device and high pressure of liquid/gas, which is potentially hazardous.

Objectives: The present study aimed at customizing SPAR-H methodology for the pigging operation using Bayesian networks (BNs). It also aimed at identifying and analyzing human errors in pigging operation in a gas transmission company.

Methods: The current article was composed of two main steps. In the first step, the SPAR-H BN model was developed using expert-elicited prior probabilities of pigging operation applied to Bayesian network. In this step, CPTs of PSF nodes are constructed using prior probabilities, which are achieved from expert opinion. The CPT of error node is developed using coding process of SPAR-H formula in a simulation node. In the second step, a descriptive study was carried out to estimate the probability of human errors in pigging operation in a gas transmission plant in Iran. First, hierarchical task analysis (HTA) was conducted by walking through the pigging operation and interviews with workers. Next, the SPAR-H BN model was utilized for estimation of human error probability.

Results: The developed model was tested on the pigging operation subtasks. In the considered case study, the mean probability of human error was estimated as 0.184. The highest probability of human error was related to “opening the kicker valve for enhancing pressure” subtask.

Conclusions: The BNs were helpful to adapt the SPAR-H methodology to the pigging operation using dedicated prior probabilities. Beside that, the probabilities of human error can be updated taking into account the more realistic operational and environmental conditions.


Human Error Standardized Plant Analysis Risk-Human (SPAR-H) Bayesian Network (BN) Pipeline Inspection Gauges (PIG)

Copyright © 2019, Author(s). This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 International License ( which permits copy and redistribute the material just in noncommercial usages, provided the original work is properly cited.

1. Background

The safety of all industries was examined by a group of continuous and systematic infrastructures, which are specified as “Critical Infrastructures” (1). These infrastructures can be systems, resources, a process, or a combination of these elements; the failure or unavailability of these infrastructures can interrupt the normal operation and safety of an industry (1). Maintenance procedures in complex industries are vulnerable to human performance (2, 3). Also, defects in human performance is one of the main reasons for accidents in maintenance operations of natural gas and oil industries. According to various studies in the context of human reliability, human error is the main cause of many accidents in various industries (4-7). For instance, more than 80% of failures in the chemical and petro-chemical industries (7) have been caused by human error. Therefore, it is required to analyze human reliability as well as system reliability. Over the years, various techniques have been developed for assessment of human reliability called human reliability analysis (HRA) methods. The HRA methods are a group of applicable tools for identification, modeling, qualification, and quantification of human error. Furthermore, HRA is a generic title for any method that estimates human error probability (HEP) for vital activities, to which human performance contributes (8). These methods are sometimes used to reduce human error in the context of probabilistic risk analysis (PRA) (9). At the moment, there are more than a dozen HRA methods for HEP in different systems, such as THERP (10), ASEP (11, 12), ATHEANA (13), SHERPA (14), CREAM (15), SPAR-H (16), and new methods are being developed. Today, various HRA methods have been utilized in multiple engineering applications (17-22). For example, Mandal et al. (18) identified probability of human error through overhead crane operation tasks using SHERPA methodology. Recently, Akyuz and Celik (19) applied a hybrid of HEART techniques and Interval type-2 fuzzy sets for assessment of human reliability for cargo operation.

Pipeline inspection gauge (PIG) operation, known as “pigging operation”, is one of the available maintenance activities for the inspection of oil and gas pipelines. Employing PIG devices are considered as a preventive maintenance procedure. Therefore, PIGs are applied for online monitoring of pipeline functionality in an operational condition. This online monitoring helps in taking corrective actions at the proper time and is useful for increasing the lifespan of pipelines (23). The most important identified risks in pigging process, such as high gas/liquid pressure in pipelines, excessive and uncontrolled speed of a PIG, the PIG striking equipment or operators and explosion, may lead to catastrophic accidents. Each year numerous incidents and accidents occur in the pigging process globally, some of which have been reported in the recent years. Table 1 describes some accidents in the current operation. These accidents are a small part of accidents that occur in pigging operations around the world. Although rigorous safety procedures are being applied in oil and gas industries, it seems that these procedures are not secure enough to prevent accidents and factors, such as human error, play a decisive role in these incidents.

Table 1. Examples of Accidents in the Pigging Operation
Place (Location)YearDescriptionsConsequence
Iran, Zahedan (24)February 6, 2017A pig device struck two operators.Two workers were killed
Iran, Fars (25)August 7, 2016In pressure decreasing procedure of receiver in pigging process of gas transmission pipeline an explosion happened.Two workers with burns covering 90 percent and 45 percent of their bodies and one operator were killed by shock wave from the blast.
India, Telangana (26)April 4, 2015The accident occurred when a PIG device, which is used to clean the pipeline, was ejected out of the pipeline.At least two operators were burnt alive in a gas pipeline explosion, four other workers were also injured in the blast.
Nigeria (27)March 13, 2015During an internal inspection project, two batch pigs were ejected under pressure from the barrel pig receiving trap.PIG device struck two operators and two workers were seriously injured.

In the present study, SPAR-H methodology and Bayesian Network (BN) were selected for analyzing human error. Although the SPARH technique was not originally developed for petrochemical industries, it has recently been applied for analysis of human reliability in risk assessment of a major accident in the Norwegian petroleum industry, offshore drilling operations in the oil and gas industry, and permit to work (PTW) system in a chemical plant (28-30).

A BN is a probabilistic graphic-based model that represents variables as nodes and also demonstrates the relationship between nodes using directed acyclic arcs. A BN employs a practical manner to utilize information and knowledge from various sources and handle missing data (31). The BN has been implemented in different contexts to estimate the probability of human error knowing that the model can conduct predictive analysis as well as diagnostic analysis. For instance, Islam et al. (3) applied BN for assessment of HEP during maintenance activities of marine systems. Musharraf et al. (32) applied BN during offshore emergency conditions.

2. Objectives

The main objective of the current study was to develop a SPAR-H- BN model to identify and analyze human errors in pigging operation in a gas transmission company. The application of the developed model assists the gas transmission company to analyze the probability of errors for pigging operation more accurately.

3. Methods

3.1. SPAR-H Methodology

In 1999, standardized plant analysis risk-human reliability (SPAR-H) was developed to succeed the accident sequence precursor (ASP) program (16, 33). Moreover, SPAR-H is the most practical technique that can be used by practitioners in any industry and this method calculates the likelihood of risks in post-initiator scenarios. The SPAR-H technique includes eight performance shaping factors (PSFs) that could influence human error probability (HEP), which is tabulated in Table 2 (34). There are two types of activity available in plant operation in the SPAR-H method: diagnosis and action. Moreover, for calculation of probability of human error, this methodology uses a dedicated multiplier for each PSF level and different formulas, according to negative levels of PSFs.

3.2. Fundamentals of Bayesian Networks (BN)

A BN is constituted from a group of unsystematic variables and their casual dependencies indicating probabilistic dependencies and interdependencies among variables through directed acyclic graph (DAG) (35). Besides that, BNs are developed from nodes and arcs (36). Nodes represent a random variable whereas probabilistic dependencies among variables are represented through arcs. Also, arcs directions show cause-effect dependencies among variables (37). A simple Bayesian network has been shown in Figure 1. Node C is a consequence or child node of A and B as parent nodes of C. Generally, in a BN, a node without any parent nodes and incoming arcs is called a “root node” and a node without child nodes and outgoing nodes is called a “leaf node” (38).

The BNs are used to specify a joint probability distribution through variables and DAGs (39, 40). Root nodes have a marginal probability distribution (MPD), and all other connected nodes have conditional probability distributions (CPD), which are dependent on the root nodes. Furthermore, a CPD indicates the effect of parent nodes on child node using quantifying process (40). For representing the quantitative effects of discrete parent nodes on their child node, a conditional probability table (CPT) has been assigned to each parent node and in the same way for each discrete root node or each continuous root node, a prior probability table or a function has been defined, respectively (41). According to the chain rule, a BN explains the joint probability distribution over all the variables, which are available in the DAG and in the following estimates the marginal and the conditional probabilities for each node of the network (42). Therefore, chain rule determines joint probability distribution through encoding the BNs process between variables (38) and calculates entries in joint probability distribution on the basis of BNs information. Joint probability distribution P(U) of variable U is:

P(X1, , Xn=i=1nPXiPUXi)

Where Пi is the parent set of variable Xi. Based on this, the marginal probability is calculated from the following equation:

PXi= except XiP(U)

In the present study, for assessment of HEP, a study group was organized consisting of a health and safety officer, resident supervisor (supervisor 1), gas transmission plant supervisor (supervisor 2), test-men operator, compressor operator, and researcher (as an interviewer). Note that all participated operators had more than four years of experience. In the first step, to capture the prior probabilities about PSFs levels, 45 questionnaires (SPAR-H worksheet) in pigging operation subtasks (three different pigging operations) were filled by experienced operators (experts). These prior probabilities were revised by gas transmission plant senior supervisor (with more than 15 years of experience). In the second step, researchers (as model builders) used prior probabilities to develop the final model. The process of developing the SPAR-H BN model is as follows:

3.3. Combine SPAR-H with Bayesian Network for Pigging Operation

When analyzing HFEs by the SPAR-H method, information and evidence of all PSFs are not available, thus, HRA practitioners have to deal with “insufficient information” for one or more PSFs in these cases. In SPAR-H method, value 1 is assigned to “insufficient information” level, which is the same as a “nominal” level value. This may lead to a HEP, similar to one that is calculated from cases with perfect evidence being generated. In fact, lack of evidence about some PSFs does not mean the evidence is at “nominal” level, which is one of the important issues in the SPAR-H method. Besides, BNs are beneficial for application of prior probabilities in conditions with “insufficient information”. In the present article, the SPAR-H BN model was developed using Agenarisk version 5.0.0 (43). The BN structure was constructed according to the SPAR-H methodology. The conditional probability table (CPT) was developed using expert judgment about performance shaping factors (PSFs) states in pigging operation. In this study, expert judgment has been incorporated with observational data and information contained in the literature (44) to capture the prior probability of each state of each PSF for action tasks. This information was used to extend marginal probability table (MPT) for the PSFs. The ultimate goal of SPAR-H BN model is estimation of MPT for the error node. Therefore, SPAR-H BN model used Equation 2 to capture marginal probability.

3.3.1. Arrangement of SPAR-H in a BN (Qualification Phase)

According to SPAR-H methodology, eight PSFs effect on the probability of error mathematically as in the following function:


Therefore, to construct a qualified model, two types of node were chosen: PSFs nodes and error node. In total, nine nodes exist: eight nodes of eight PSFs and one node of human error. Based on SPAR-H methodology instruction, PSFs in this method are independent and each of these PSFs directly effects on the probability of human error. Hence, arcs were used to demonstrate conditional independence among PSFs. Initial structure of SPAR-H BN is depicted in Figure 2. In the created model, the levels of each PSF were the same as each PSF level in the original SPAR-H methodology (except “insufficient information” level). The PSFs of SPAR-H method and their related levels and prior probabilities are listed in Table 2. The error node calculates human failure events probability based on PSFs nodes level and has two discrete levels: 1 (error occur = P (error)) and 0 (error not occur = 1 - P (error)). With regards to the original SPAR-H method, PSFs are independent so there is no casual arc between PSFs nodes and only the casual arcs exists between each PSF and error node.

Naïve SPAR-H BN model for pigging process
Figure 2. Naïve SPAR-H BN model for pigging process
3.3.2. Development of Conditional Probability Table (CPT) CPT of PSFs Nodes

For development of CPT of each PSF node, the expert-elicited probabilities of PSFs levels in pigging operation were used. The probability values for each PSF level are listed in Table 2. Bayesian networks have the ability to use prior possibilities (40). In SPAR-H methodology, in some cases, information about PSFs levels are incomplete. This methodology assigned “insufficient information” level with multiplier of one for these cases, which is the same as assigning “nominal” level. However, lack of access to information about PSFs does not mean the PSFs are in “nominal” level. Besides, “insufficient information” level in SPAR-H methodology was excluded in order to use available prior probabilities in Table 2.

Table 2. SPAR-H Method PSFs, Their Related State, Multiplier, and Prior Probability of Each State
PSF, StateMultiplier (Diagnosis Tasks)Multiplier (Action Tasks)Prior Probabilities (Action Tasks)
Available time
Time available ≥ 50 × the time required0.010.010.005
Time available ≥ 5 × the time required0.10.10.15
Nominal time110.80
Time available is ~ the time required10100.055
Inadequate timeP (failure) = 1P (failure) = 10.000001
Clear diagnosis0.1--
Moderately complex2236.6
Highly complex554.9
Available, but poor5530.6
Not available50502.5
Ergonomics/human machine interference (HMI)
Fitness for duty
Degrade fitness559.3
UnfitP (failure) = 1P (failure) = 10.000001
Work processes
Poor553.3 CPT of an Error Node

The SPAR-H methodology estimates probability of human error using predefined expression based on PSF composite. The value of HEP depends on the PSFs levels (P (error | PSF1-8). Considering this, to develop a CPT the number of probabilities will increase exponentially according to the number of parent nodes (45). In a case, such as SPAR-H method with the eight parent node, the number of probabilities is more than a thousand. Therefore, to develop the CPT of the error node, simulation node and coding process were applied to build error node CPT more conveniently and more accurately. It follows that:

Agenarisk software benefits from a node probability table (NPT) editing mode, which makes it possible to use a mathematical expression to build a CPT of nodes. Given that to construct the error node CPT, the SPAR-H formula (which is mentioned in the next section explicitly) was coded in the model directly. In the event that the “available time” and “fitness for duty” levels are “inadequate” and “unfit” respectively, the probability of human error is appointed equivalent to 1 notwithstanding the other PSFs levels (like original SPAR-H methodology). Based on the original SPAR-H method instruction, this method uses a correction factor for three or more negative levels of PSFs so a dummy node was used for counting the negative nodes in the Agenarisk model. At the end, a coding was added for testing if estimated value of error exceeded 1.0, rounding down the value of error to 1.0.

3.3.3. Estimating the Probability of Human Error

The BN model of SPAR-H uses the marginal probabilities of CPT of PSF composite to estimate the probability of human error. It is found that:


In cases with complete information about PSFs states, the SPAR-H BN model estimates HEP the same as the original SPAR-H methodology yet in cases with insufficient information, the SPAR-H BN model uses prior probabilities to estimate HEP. It demonstrates the strength of this model to quantify human error. Figure 2 exhibits the SPAR-H BN model with no new information from pigging operation. This means no evidence entered the model (prior model).

3.4. Application of Developed Model: Identifying and Analyzing Human Errors Using the SPAR-H BN Model

This descriptive cross-sectional study was done to analyze human reliability in pipeline inspection gauge (PIG) operation in a gas transmission company in Iran. First, hierarchical task analysis (HTA) was conducted via a walk through in pigging operation and interviewing workers along with assessment procedures of various tasks in the studied company (Table 3). Accordingly, investigating human error probability (HEP) in PIG operation tasks/subtasks set was a goal of the present study. Afterwards, the BN model of SPAR-H (for action tasks) and SPAR-H (for diagnosis tasks) method (16) were utilized to estimate the probability of human error.

According to SPAR-H, step-by-step guidance (46), the following steps were done for human reliability analysis in PIG operation:

Step 1: Categorizing the HFE as Diagnosis and/or Action

In this step, PIG operation recognized HFEs, were classified as either diagnosis tasks or action tasks or both. For quantification purpose of SPAR-H method, diagnosis tasks represent the complete spectrum of cognitive processing and action tasks represent tasks that have been executed. According to Gertman et al.’s study (16), SPAR-H method dedicated nominal human error probability (NHEP) values of 0.01 and 0.001 for diagnosis tasks and action tasks, respectively. These values refer to error rate in cognitive processing for diagnosis tasks and also refer to error rate on simple action implementation, such as pressing a button or turning a dial, and simple slips or lapses for action tasks.

Step 2: Assessment and Rating Eight Performance Shaping Factors (PSFs)

The next step was to analyze HFEs in pigging operation based on eight PSFs in Table 2. Each PSF was investigated with respect to HFEs contexts. Therefore, numerical values were allocated according to SPAR-H instruction to each HFE. To achieve this goal, a series of interviews with experts and operators, who participated in PIG operation along with observation of work procedure, was done. Note that in the action tasks, when there was some doubt about choosing the appropriate PSF levels, the multiplier was shown by a question mark (Table 4) and then, the BN model of SPAR-H was used to estimate HEP. It should be noted in some subtasks of action tasks, some doubt existed with regards to rating PSFs. Therefore, in these cases, SPAR-H BN model was used to estimate HEP.

Step 3: Estimating PSF-Modified HEP

After the PSFs levels were determined, the PSF composite was estimated through multiplying these PSFs levels. Then, the final HEP was calculated by multiplying NHEP by the PSF composite (Equation 5).


For condition with three or more negative PSFs levels (levels with multiplier greater than 1), HEP was calculated using the correction factor, according to Equation 6:


Where HEP is human error probability, NHEP is nominal human error probability for diagnosis task and action task activity are 0.01 and 0.001, respectively. PSF is performance shaping factors multiplier.

It should be noted that the prior probability of states of PSFs, which have been mentioned in previous sections, is only for PSFs of action task, moreover, the prior probability of diagnosis tasks is being studied. Therefore, in this step, Bayesian model of SPAR-H methodology was used to estimate probability of human error for action tasks. Original SPAR-H methodology was used for estimation of probability of human error in diagnosis tasks, which means that in a condition with partial information (“insufficient information”), the SPAR-H BN model and original SPAR-H methodology were applied for action tasks and diagnosis tasks, respectively.

Step 4: Calculation of Independent HEP

For a single HFE with composed diagnosis and action tasks, two HEPs were determined separately and then summed to produce the composite HEP. Human error probability for independent tasks was calculated according to Equation 7:


Where PW/OD is human error probability (HEP) for Independent task, HEPA is human error probability (HEP) for action tasks, HEPD is human error probability (HEP) for diagnosis tasks.

Step 5: Accounting for Dependency

In the SPAR-H method, dependency between two tasks refers to the extent, to which performance on one task effects performance of a subsequent task (47). Dependency may be recognized between two tasks carried out by the same crew or among multiple tasks conducted by different operators (16). As a result, the HEP on one task increases the probability of human error on a subsequent task (47). Given the dependence of some tasks, dependency among tasks in the studied process is calculated from Table 4. Accordingly, there were five types of dependency (complete, high, moderate, low, and zero) among tasks that were determined according to crew (same or different), time (close in time or not close in time), location (same or different), and cause (additional or not additional).

Table 3. Hierarchical Task Analysis (HTA) of Tasks in Pigging Operation
OperatorTaskSub Task (Symbol)
SupervisorChecking the gas flow rate in the pipelineEstimation of the average flow rate of the pipeline (X1)
Calculation of gas velocity based on gas flow rate (X2)
Calculation of the required average gas velocity for crossing the pig (X3)
Comparison of the average gas flow rate in the pipeline and the required average velocity for PIG crossing (X4)
Ensuring the PIG dimensionsControlling the changes in a pig that had been used previously (X5)
Comparing the diameter of the PIG with the barrel-shaped section of pipeline and trap (X6)
Ensuring the openness of the valves in linesPaying attention to the relief valve gauge (X7)
Checking valves manually (X8)
Issuing permit to workHandover the work to the operators (X9)
Signing the permit (X10)
Validation & revalidation after shift handover (X11)
Compressor operator and semi-skilled workerInjection of air to the back of PIGConnecting the hose to the pipeline (X12)
Turning on the compressor (X13)
Enhancing the pressure (X14)
Pressure decreasing (X15)
Turning off the compressor (X16)
Test manPIG launchingImporting the PIG to header and pipeline (X17)
Closing the header (X18)
Activation of PIG signaler (X19)
Checking vent and drain valves (X20)
Opening the kicker valve for enhancing pressure (X21)
Checking pressure gauges (X22)
PIG receivingOpening bypass valve for decreasing pressure (X23)
Opening drain and vent valve for evacuation of pipeline (X24)
Checking pressure gauges and ensuring full discharge of pressure in pipeline (X25)
Taking of the pig (X26)
Cleaning the garbage that remains from pigging operation for fire and explosion prevention (X27)
Table 4. An Example of the SPAR-H Worksheet for Dependency Determination in “Opening the Kicker Valve for Enhancing Pressure” Task
Condition NumberCrew (Same or Different)Time (Close in Time or not Close in Time)Location (Same or Different)Cause (Additional or not Additional)DependencyHEP Calculation Formula
1ScsnacompleteFor complete dependence the probability of failure is 1
3dnahighFor high dependence the probability of failure is (1 + Pw/od)/2
6amoderateFor moderate dependence the probability of failure is (1 + 6 × Pw/od)/7
8alowFor low dependence the probability of failure is (1 + 19 × Pw/od)/20
9DcsnamoderateFor moderate dependence the probability of failure is (1 + 6 × Pw/od)/7
13ncsnalowFor low dependence the probability of failure is (1+19 × Pw/od)/20
17zeroFor zero dependence the probability of failure is Pw/od
Table 5. Rating of PSFs and Estimation of HEPA and HEPD in Pigging Operation Tasks in the Studied Gas Transmission Plant
Task/sub taskOperatorDiagnosis or Action or BothPSFHEP
Available TimeStress/StressorsComplexityExperience/ TrainingProcedureErgonomicsFitness for DutyWork ProcessesHEPAHEPD
X1Supervisor 1D11215111-0.1
X2Supervisor 1D11215111-0.1
X3Supervisor 1D11215111-0.1
X4Supervisor 1D111111110.01
X5Supervisor 1D11111111-0.01
X6Supervisor 1D11111111-0.01
X7Supervisor 1D11111111-0.01
X8Supervisor 1A1111(?,?,0,0)1110.0023941-
X9Supervisor 2D12111111-0.02
X10Supervisor 2D11111111-0.01
X11Supervisor 2D11111111-0.01
X12Compressor operatorA111111110.001-
X13Compressor operatorA121111110.002-
X14Compressor operatorA1511(?,?,0,0)1110.01197-
X15Compressor operatorA1111(?,?,0,0)1110.0023941-
X16Compressor operatorA111111110.001-
X17Test man 1A1221110110.038499-
X18Test man 1A1111(?,?,0,0)10110.023941-
X19Test man 1A1111(?,?,0,0)1110.0023941-
X20Test man 1D11111111-0.01
X21Test man 1A(0,0,?,?,0)521(?,?,0,0)10110.20389-
X22Test man 1D11110.51110.005
X23Test man 2A(0,0,?,?,0)521(?,?,0,0)10110.20389-
X24Test man 2A1221(?,?,0,0)10110.083216-
X25Test man 2D12210.51110.02-
X26Test man 2A15(?,?,0)1(?,?,0,0)10110.1305-
X27Test man 2A1211(?,?,0,0)10110.044742-

4. Results

Application of the BN model of SPAR-H to the pigging operation is summed up in Figures 3 and 4. In total, five operators participated in pigging operation in the studied gas transmission company, and human reliability analysis was conducted for all tasks that were done in this operation. These operators included two supervisors (with seven and four years work experience, respectively). One compressor operator (with six years’ work experience) and two test-men (with eight and four years of work experience). Tasks and subtasks as a result of HTA are tabulated in Table 3. As the table illustrates, the pigging operation is comprised of seven tasks and also 27 subtasks. Table 5 presents the results of evaluation of each of the eight PSFs with regards to the context of subtasks, including diagnosis, action or both, along with the participated operators and HEPA or HEPD for all of them. Table 6 presents the values of HEPA, HEPD, PW/OD, PW/OD, the level of dependency, and final HEP for each subtask in pigging operation. According to the developed model outputs, the mean of final HEP in the studied pigging operation was 0.1724. The highest and the lowest final HEP was relevant to “opening the kicker valve for enhancing pressure” and “connecting the hose to the pipeline”, with values of 0.602 and 0.05, respectively. Figure 3 presents the comparison of final HEP among pigging operation tasks, and comparison of the final HEP in occupational groups that participated in pigging operation has been illustrated in Figure 4. These figures show the considerable differences between probabilities of human error in “test-man” operators with mean HEP.

Comparison of the mean human error probability among tasks of pigging operation
Figure 3. Comparison of the mean human error probability among tasks of pigging operation
The mean probability of human errors in the operators involved in pigging operation in the studied gas transmission company
Figure 4. The mean probability of human errors in the operators involved in pigging operation in the studied gas transmission company
Table 6. Final Human Error Probabilities in Pigging Operation Tasks
Task/Sub TaskOperatorHEPDHEPAPW/ODDependencyPW/DFinal HEP Using SPAR-H BN ModelFinal HEP Using SPAR-H Methodology
X1Supervisor 1-0. 10.1Low0.1450.1450.145
X2Supervisor 1-0. 10.1Moderate0.2280.2280.228
X3Supervisor 1-0. 10.1Moderate0.2280.2280.228
X4Supervisor 10.010.01Moderate0.1510.1510.151
X5Supervisor 1-0.010.01Low0.0590.0590.059
X6Supervisor 1-0.010.01Moderate0.1510.1510.151
X7Supervisor 1-0.010.01Low0.0590.0590.059
X8Supervisor 10.0024-0.0024Moderate0.1440.1440.1471
X9Supervisor 2-0.020.02Low0.0690.0690.069
X10Supervisor 2-0.010.01Moderate0.0590.1510.151
X11Supervisor 2-0.010.01Moderate0.0590.1510.151
X12Compressor operator0.001-0.001Low0.0500.0500.050
X13Compressor operator0.002-0.002Moderate0.1440.1440.144
X14Compressor operator0.012-0.012Moderate0.1530.1530.153
X15Compressor operator0.0024-0.0024Moderate0.1440.1440.147
X16Compressor operator0.001-0.001Moderate0.1430.1430.143
X17Test man 10.039-0.039Low0.0870.0870.0865
X18Test man 10.024-0.024High0.5120.5120.525
X19Test man 10.0024-0.0024Low0.0520.0520.0574
X20Test man 1-0.010.01Moderate0.0590.0590.059
X21Test man 10.204-0.204High0.6020.6020.642
X22Test man 10.0050.005Moderate0.1470.1470.147
X23Test man 20.204-0.204Low0.2430.2430.841
X24Test man 20.084-0.084Moderate0.2140.2140.285
X25Test man 20.02-0.02Moderate0.1600.1600.160
X26Test man 20.13-0.13High0.5650.5650.666
X27Test man 20.045-0.045Moderate0.1810.1810.220

5. Discussion

The current study was performed in order to model SPAR-H methodology in BN in the context of pigging operation and also to characterize HFEs along with determining HEP in PIG operation in a gas transmission company through implementation of SPAR-H BN model and original SPAR-H methodology for action tasks and diagnosis tasks, respectively. First, SPAR-H methodology was modeled using BN and expert-elicited probabilities. In the second stage, for analyzing human error, firstly, a hierarchical task analysis was performed. In a total of seven main tasks and also twenty-seven subtasks were recognized in the studied company. According to the final HEPs in Table 6, HEP intervals for pigging operation were in the range of 0.05 and 0.602 and also the mean probability of human error in pigging operation was 0.184. Furthermore, the highest probability of human error was related to “opening the kicker valve for enhancing pressure” subtask (0.602), which may increase the risk of launcher rupture or hose rupture of pressure compressor and as a result increases the risk of exposure to high pressure of gas/liquid (400 - 1200 PSI) and operator mortality. This error is due to high dependency on the previous task, no availability of time (doubt exists), extreme stress level, moderate complexity in subtask, poor ergonomic condition, and also lack of (doubt exists) applicable procedure. The “taking of the pig” subtask also had high HEP value (0.565), which is due to probability of PIG device being struck in pipeline and also risk of PIG device striking the operator that can lead to extreme stress level. Furthermore, task complexity (doubt existed between levels), poor ergonomic condition, lack of a strong procedure (doubt existed between levels), and high dependency on the previous subtask made it a critical subtask. The “closing the header” subtask was another critical subtask in pigging operation (0.512). This error can lead to the launcher header staying open and cause the operator to be exposed to high pressure (gas/liquid). “Calculation of gas velocity based on gas flow rate” and “calculation of the required average gas velocity for crossing the pig” subtasks also had HEP (0.228) higher than mean HEP. Error in these subtasks can have catastrophic consequences, such as pipeline rupture (also launcher/receiver), explosion or jet fire. Moreover, X23, X24, and X27 subtasks have the high HEP value, which are due to poor ergonomic condition, no availability of strong procedure (doubt exists) and high stress level. Given the above, the mean HEP value in present study was 0.184, which was close to HEP value in the “sweep the line by pigging in order to drain all cargo residues” sub task (0.186) in Akyuz and Celik’s (48) study. In the current study, 37% of subtasks had high or extreme stress level and 50% of these cases with stress had complexity. This shows high correlation between these PSFs, which is in agreement with Abreu et al.’s study (21).

Given Figures 3 and 4, the subtasks that have been done by test man operators had the highest probability of human error. Because of the high interaction of “test-man” operators by pigging operation, they are faced with more stressful conditions, more complex tasks, and awkward ergonomics conditions. Supervising operators had a vital role in the pigging operation yet because of appropriate level of PSFs (the nominal level of stress, nominal ergonomic condition, etc.), the HEP values were lower than mean HEP. However, in some supervisor subtasks, the complexity was high and also the procedure was poor.

According to Peterson’s (49) study, human errors caused by weakness in the procedures are responsible for 24% of industrial accidents. Furthermore, using improper procedures can lead to rule-based human error (50). As mentioned above, the available procedures in pigging operation are not applicable and strong, which has a considerable effect on probability of human error. Therefore, by revising the procedures the HEP values can be decreased.

5.1. Evaluation of the Final Model

Based on models purpose, different evaluation methods are available. If the model was developed to carry out inference, the evaluation process of the model can be done via stakeholders or expert panel opinion (51). If the model (the BN model) has a target node, a sensitivity analysis can be used to determine the effect of sensitivity nodes (parent nodes) on the target node (child node) (52).

5.2. Sensitivity Analysis

One of the most used features to validate BN models is sensitivity analysis. Therefore, in the present study, it was applied to test the developed model. Figure 5 illustrated the tornado graph for excepted values of error probability, which is resulted from the sensitivity analysis of the final model. According to tornado graph, the estimated final HEP are sensitive yet did not show abrupt changes to any minor change of the levels associated with the PSFs, except for levels of two PSFs (“time = not available” and “fitness for duty = unfit”, for these two level final HEP estimated equivalent to 1). Nevertheless, these changes are in accordance with the SPAR-H methodology completely. Therefore, sensitivity analysis approved the logically of the final model.

Tornado graph for excepted values of error probability
Figure 5. Tornado graph for excepted values of error probability
5.3. Expert Judgment

As discussed at the end paragraph of section 2.2, in the present article, a study group (an expert panel) was organized consisting of a health and safety officer, resident supervisor (supervisor 1), gas transmission plant supervisor (supervisor 2), test-men operator, compressor operator, and researcher (as an interviewer). After initial evaluation about PSFs levels, experts panel were requested to judge the final probability of human error in the action tasks. The results of expert judgment were in convergent with output results of the developed model (correlation value = 0.73). Hence, the logically of the final model was authorized by expert judgment.

5.4. Comparison of SPAR-H BN Model with SPAR-H Methodology

In Figure 3, results of final human errors were compared. According to the comparison of output values from both SPAR-H BN model and SPAR-H methodology, which are presented in Table 6, although the HEPs of diagnosis task (subtask X1 to X11) are equal (because SPAR-H BN model is developed for action tasks), the main difference between the two methods can be seen in action tasks (X12 to X27). It was found that due to more accurate calculation of HEP in SPAR-H BN model than SPAR-H methodology, outputs of SPAR-H BN model are more realistic and are lesser than SPAR-H methodology outputs. For instance, for “opening bypass valve for decreasing pressure” subtask the output of SPAR-H BN model and SPAR-H methodology was 0.243 and 0.841, respectively. This means that in these cases, SPAR-H methodology estimates HEP very pessimistically, also this problem can be seen in X18, X19, X21, X24, X26, and X27. Pessimistic calculation of HEP increases the costs of industry and also enhances the stresses in managerial level. In addition, SPAR-H BN model prevents pigging operation practitioners from the optimistic calculation of HEP. This can lead to a lack of attention to the system or human performance deficiencies.

5.5. Conclusions

The SPAR-H method demonstrates an acceptable spectrum of factors influencing human performance in assessment of industries operations through a framework of an acceptable number of PSFs. This method is easy to use, simple, and a useful tool for practitioners to utilize for calculating HEP. Even though some deficiencies exist in this methodology, through implementing SPAR-H in BNs, some of the issues were addressed. The BNs were helpful to adapt the SPAR-H methodology to the pigging operation using dedicated prior probabilities. One of the main advantages of the presented SPAR-H model in this study was that once doubt existed for choosing PSFs states, the prior probabilities were used to estimate HEP accurately. Therefore, the probabilities of human error can be updated taking into account the more realistic operational and environmental conditions, such that SAPR-H methodology did not benefit from these advantages. In total, SPAR-H BN is a practical model for managers to understand working condition more realistically than SPAR-H methodology to decrease risks contributed with human performance.




  • 1. De Felice F, Petrillo A, Zomparelli F. A hybrid model for human error probability analysis. IFAC-Papers OnLine. 2016;49(12):1673-8. doi: 10.1016/j.ifacol.2016.07.821.
  • 2. Aalipour M, Ayele YZ, Barabadi A. Human reliability assessment (HRA) in maintenance of production process: a case study. Int J Syst Assur Eng Manag. 2016;7(2):229-38. doi: 10.1007/s13198-016-0453-z.
  • 3. Islam R, Khan F, Abbassi R, Garaniya V. Human error probability assessment during maintenance activities of marine systems. Saf Health Work. 2018;9(1):42-52. doi: 10.1016/ [PubMed: 30363076]. [PubMed Central: PMC6111134].
  • 4. Wu S, Hrudey S, French S, Bedford T, Soane E, Pollard S. A role for human reliability analysis (HRA) in preventing drinking water incidents and securing safe drinking water. Water Res. 2009;43(13):3227-38. doi: 10.1016/j.watres.2009.04.040. [PubMed: 19493557].
  • 5. Helmreich RL. On error management: Lessons from aviation. BMJ. 2000;320(7237):781-5. doi: 10.1136/bmj.320.7237.781. [PubMed: 10720367]. [PubMed Central: PMC1117774].
  • 6. Ren J, Jenkinson I, Wang J, Xu DL, Yang JB. A methodology to model causal relationships on offshore safety assessment focusing on human and organizational factors. J Safety Res. 2008;39(1):87-100. doi: 10.1016/j.jsr.2007.09.009. [PubMed: 18325420].
  • 7. Kariuki SG, Lowe K. Integrating human factors into process hazard analysis. Reliab Eng Syst Saf. 2007;92(12):1764-73. doi: 10.1016/j.ress.2007.01.002.
  • 8. Swain AD. Human reliability analysis: Need, status, trends and limitations. Reliab Eng Syst Saf. 1990;29(3):301-13. doi: 10.1016/0951-8320(90)90013-d.
  • 9. Kirwan B. Human error identification in human reliability assessment. Part 2: Detailed comparison of techniques. Appl Ergon. 1992;23(6):371-81. [PubMed: 15676882].
  • 10. Swain AD, Guttmann HE. Handbook of human-reliability analysis with emphasis on nuclear power plant applications. Final report. United States: U.S. Department of Energy; 1983.
  • 11. Swain AD. Accident Sequence Evaluation Program: Human reliability analysis procedure. United States: Sandia National Labs; 1987.
  • 12. Kirwan B. Validation of human reliability assessment techniques: Part 1 - Validation issues. Safety Science. 1997;27(1):25-41. doi: 10.1016/s0925-7535(97)00049-0.
  • 13. Cooper SE, Ramey-Smith AM, Wreathall J, Parry GW. A technique for human error analysis (ATHEANA). United States: U.S. Department of Energy; 1996.
  • 14. Kirwan B. Human error identification in human reliability assessment. Part 1: Overview of approaches. Appl Ergon. 1992;23(5):299-318. doi: 10.1016/0003-6870(92)90292-4.
  • 15. He X, Wang Y, Shen Z, Huang X. A simplified CREAM prospective quantification process and its application. Reliab Eng Syst Saf. 2008;93(2):298-306. doi: 10.1016/j.ress.2006.10.026.
  • 16. Gertman D, Blackman H, Marble J, Byers J, Smith C. The SPAR-H human reliability analysis method. US Nuclear Regulatory Commission; 2005.
  • 17. Akyuz E, Celik M. Application of CREAM human reliability model to cargo loading process of LPG tankers. J Loss Prev Process Indust. 2015;34:39-48. doi: 10.1016/j.jlp.2015.01.019.
  • 18. Mandal S, Singh K, Behera RK, Sahu SK, Raj N, Maiti J. Human error identification and risk prioritization in overhead crane operations using HTA, SHERPA and fuzzy VIKOR method. Expert Syst Appl. 2015;42(20):7195-206. doi: 10.1016/j.eswa.2015.05.033.
  • 19. Akyuz E, Celik E. A modified human reliability analysis for cargo operation in single point mooring (SPM) off-shore units. Appl Ocean Res. 2016;58:11-20. doi: 10.1016/j.apor.2016.03.012.
  • 20. Ribeiro AC, Sousa AL, Duarte JP, Frutuoso e Melo PF. Human reliability analysis of the Tokai-Mura accident through a THERP–CREAM and expert opinion auditing approach. Safety Science. 2016;87:269-79. doi: 10.1016/j.ssci.2016.04.009.
  • 21. Abreu JM, Hardy N, Pereira FC, Zeifman M. Modeling human reliability in the power grid environment: An application of the SPAR-H methodology. Proc Hum Factors Ergon Soc Ann Meet. 2016;59(1):662-6. doi: 10.1177/1541931215591144.
  • 22. Ung ST. A weighted CREAM model for maritime human reliability analysis. Saf Sci. 2015;72:144-52. doi: 10.1016/j.ssci.2014.08.012.
  • 23. Canavese G, Scaltrito L, Ferrero S, Pirri CF, Cocuzza M, Pirola M, et al. A novel smart caliper foam pig for low-cost pipeline inspection - Part A: Design and laboratory characterization. J Petrol Sci Eng. 2015;127:311-7. doi: 10.1016/j.petrol.2015.01.008.
  • 24. Shana Petro Energy Information Network. 2017. Available from:
  • 25. Ta'adol Newapaper. 2018. Available from:
  • 26. Indian Express. 2018. Available from:
  • 27. Pigging products. 2018. Available from:
  • 28. Jahangiri M, Hoboubi N, Rostamabadi A, Keshavarzi S, Hosseini AA. Human error analysis in a permit to work system: A case study in a chemical plant. Saf Health Work. 2016;7(1):6-11. doi: 10.1016/ [PubMed: 27014485]. [PubMed Central: PMC4792918].
  • 29. Gould KS, Ringstad AJ, van de Merwe K. Human reliability analysis in major accident risk analyses in the norwegian petroleum industry. Proc Hum Factors Ergon Soc Ann Meet. 2012;56(1). doi: 10.1177/1071181312561421.
  • 30. van de Merwe K, Qie S, Gould KS. The application of the SPAR-H method in managed-pressure drilling operations. Proc Hum Factors Ergon Soc Ann Meet. 2016;56(1):2021-5. doi: 10.1177/1071181312561422.
  • 31. Chen SH, Pollino CA. Good practice in Bayesian network modelling. Environ Model Software. 2012;37:134-45. doi: 10.1016/j.envsoft.2012.03.012.
  • 32. Musharraf M, Hassan J, Khan F, Veitch B, MacKinnon S, Imtiaz S. Human reliability assessment during offshore emergency conditions. Saf Sci. 2013;59:19-27. doi: 10.1016/j.ssci.2013.04.001.
  • 33. Rasmussen M, Standal MI, Laumann K. Task complexity as a performance shaping factor: A review and recommendations in standardized plant analysis risk-human reliability analysis (SPAR-H) adaption. Saf Sci. 2015;76:228-38. doi: 10.1016/j.ssci.2015.03.005.
  • 34. Laumann K, Rasmussen M. Suggested improvements to the definitions of standardized plant analysis of risk-human reliability analysis (SPAR-H) performance shaping factors, their levels and multipliers and the nominal tasks. Reliab Eng Syst Saf. 2016;145:287-300. doi: 10.1016/j.ress.2015.07.022.
  • 35. Hu Z, Mahadevan S. Uncertainty quantification and management in additive manufacturing: current status, needs, and opportunities. Int J Adv Manuf Tech. 2017;93(5-8):2855-74. doi: 10.1007/s00170-017-0703-5.
  • 36. Ashrafi M, Davoudpour H, Khodakarami V. A bayesian network to ease knowledge acquisition of causal dependence in CREAM: Application of recursive noisy-OR gates. Qual Reliab Eng Int. 2017;33(3):479-91. doi: 10.1002/qre.2021.
  • 37. Fenton N, Neil M, Lagnado D, Marsh W, Yet B, Constantinou A. How to model mutually exclusive events based on independent causal pathways in Bayesian network models. Knowl Base Syst. 2016;113:39-50. doi: 10.1016/j.knosys.2016.09.012.
  • 38. Boudali H, Dugan JB. A continuous-time Bayesian network reliability modeling, and analysis framework. IEEE Trans Reliab. 2006;55(1):86-97. doi: 10.1109/tr.2005.859228.
  • 39. Zhang Q. Dynamic uncertain causality graph for knowledge representation and probabilistic reasoning: Directed cyclic graph and joint probability distribution. IEEE Trans Neural Netw Learn Syst. 2015;26(7):1503-17. doi: 10.1109/TNNLS.2015.2402162. [PubMed: 25781960].
  • 40. Pearl J. Fusion, propagation, and structuring in belief networks. Artif Intell. 1986;29(3):241-88. doi: 10.1016/0004-3702(86)90072-x.
  • 41. Li PC, Chen GH, Dai LC, Zhang L. A fuzzy Bayesian network approach to improve the quantification of organizational influences in HRA frameworks. Safety Science. 2012;50(7):1569-83. doi: 10.1016/j.ssci.2012.03.017.
  • 42. Trucco P, Cagno E, Ruggeri F, Grande O. A Bayesian Belief Network modelling of organisational factors in risk analysis: A case study in maritime transportation. Reliab Eng Syst Saf. 2008;93(6):845-56. doi: 10.1016/j.ress.2007.03.035.
  • 43. Bayesian network software for risk analysis and decision making. 2018. Available from:
  • 44. Hallbert B. The employment of empirical data and Bayesian methods in human reliability analysis: A feasibility study. US Nuclear Regulatory Commission, Office of Nuclear Regulatory Research; 2007.
  • 45. Tang Z, McCabe B. Developing complete conditional probability tables from fractional data for bayesian belief networks. J Comput Civ Eng. 2007;21(4):265-76. doi: 10.1061/(asce)0887-3801(2007)21:4(265).
  • 46. Whaley AM, Kelly DL, Boring RL, Galyean WJ. SPAR-H step-by-step guidance. Helsinki, Finland: Idaho National Laboratory (INL); 2012.
  • 47. Deng X, Jiang W. Dependence assessment in human reliability analysis using an evidential network approach extended by belief rules and uncertainty measures. Ann Nucl Energ. 2018;117:183-93. doi: 10.1016/j.anucene.2018.03.028.
  • 48. Akyuz E, Celik M. A hybrid human error probability determination approach: The case of cargo loading operation in oil/chemical tanker ship. J Loss Prev Process Indust. 2016;43:424-31. doi: 10.1016/j.jlp.2016.06.020.
  • 49. Petersen D. Human-error reduction and safety management. 1984.
  • 50. Barry K. A Guide to practical human reliability assessment. London: Taylor & Francis; 1994. doi: 10.1201/9781315136349.
  • 51. Aguilera PA, Fernandez A, Fernandez R, Rumi R, Salmeron A. Bayesian networks in environmental modelling. Environ Model Software. 2011;26(12):1376-88. doi: 10.1016/j.envsoft.2011.06.004.
  • 52. Martins MR, Maturana MC. Application of Bayesian Belief networks to the human reliability analysis of an oil tanker operation focusing on collision accidents. Reliab Eng Syst Saf. 2013;110:89-109. doi: 10.1016/j.ress.2012.09.008.