Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2021 / Article
Special Issue

Operations Research in Engineering Problems: Potential Applications and Perspectives

View this Special Issue

Research Article | Open Access

Volume 2021 |Article ID 5557320 |

Mehdi Jahangiri, Abolfazl Moghadasi, Mojtaba Kamalinia, Farid Sadeghianjahromi, Sean Banaee, "Influence of the Inherent Safety Principles on Quantitative Risk in Process Industry: Application of Genetic Algorithm Process Optimization (GAPO)", Mathematical Problems in Engineering, vol. 2021, Article ID 5557320, 11 pages, 2021.

Influence of the Inherent Safety Principles on Quantitative Risk in Process Industry: Application of Genetic Algorithm Process Optimization (GAPO)

Academic Editor: Mohammad Yazdi
Received20 Feb 2021
Revised21 Apr 2021
Accepted27 Apr 2021
Published08 May 2021


Inherent safety (IS) refers to a set of measures that enhance the safety level of processes and equipment, rendering additional equipment and/or add-ons. The early design phase of processes is suited best for implementation of IS strategies as some of such strategies either are impossible to be implemented at the operation phase or substantially increase costs. The purpose of this study is to present a new approach called genetic algorithm process optimization (GAPO), by which processes can be made inherently safer even at the operation phase. This study simulates the IS principle, assessing its impact on quantitative risk and the possible consequences of process incidents identified by Hazard and Operation Study (HAZOP). The principle of intensification was simulated through GAPO, and feasibility of implementation was approved by HYSYS. Moreover, the integrated inherent safety index (I2SI) was used to evaluate and quantify the level of IS following implementation of GAPO compared to the initial design. Our result shows that GAPO substantially reduced the risk of consequences and quantitative risks and concomitantly improved the I2SI. The proposed GAPO can be applied to process operation as an approach to enhance IS at no cost and without decrease in production.

1. Introduction

Engineers and safety activists have always tried to optimize and increase the level of safety in industry through utilizing knowledge and advanced technologies to reduce the probability and severity of human and financial consequences [13]. Meanwhile, given the inherent nature of processes, reactions, raw materials, and products, process industries contribute a special role in the occurrence of incidents [46]. In process industries including petrochemical facilities, raw materials are converted into intermediate or final products using a physical and/or chemical chain of circumstances. In such industries, production, storage, transportation, use, and disposal of chemicals are inherently dangerous, and the potential for catastrophic accidents is a significant concern. Since simple negligence in these industries may lead to loss of life, damage to equipment, economic losses, and environmental pollution, great efforts have been made to minimize the chance for occurrence of such incidents.

In a general classification, approaches towards achieving higher levels of safety can be divided into two groups: conventional and inherent [710]. Conventional approaches include control of hazards by safety add-ons such as active or passive engineering strategies, procedural corrective actions, and preventive laws and regulations [3]. In passive safety systems, hazards are controlled through process or equipment design features without additional active functioning of any device. Active safety systems include process controls, safety instrumented systems (SIS), and mitigation systems. Such preventative measures mitigate hazards through controls and systems designed to monitor and maintain specific conditions, which may be triggered by an event. Procedural safety systems use personnel education and management, including standard operating procedures (SOPs), safety rules and procedures, training, emergency response plans, and management systems to control the hazards [10].

Many safety activists applied the inherent safety principles to reduce or eliminate risks and make processes or plants safer without knowing or categorizing the techniques as inherent safety [3, 11, 12]. In 1970, Klets came up with the idea that “What you do not have, cannot leak, or burn” [11]. Ten years after the explosion of Flixborough, and just a few weeks before the Bhopal tragedy, Klets introduced the four basic principles for inherent safety including intensification, substitution, attenuation, and limitation of the effects of failures. During the following years, these principles were reviewed and modified by various researchers [11]. For example, the I2SI described the inherent safety principles (ISPs) in five elements as limitation of effects, minimization, substitution, attenuation, and simplification [3].

Inherent safety refers to a set of measures that enhance the processes and equipment safety level without the need for further equipment and/or add-ons and prevents or reduces the severity of possible incidents [13]. Inherent safety seeks to eliminate the risk with different approaches rather than ignoring or controlling them [14, 15]. Inherent safety design strategies eliminate or reduce hazards to avoid or mitigate the consequences of incidents through principles such as substitution, moderation, and simplification of conventional safety approach. Although implementation of ISPs seems logical, simple, and obvious at first glance, it is challenging because of the required modification in process design, tools, and layouts. Moreover, enforcement of ISPs in a simple process is more difficult due to limited options for achieving inherent safety. For instance, chemicals may not be toxic, and processes are not complicated. These challenges are greater when the process is at the operation stage of the lifecycle, because any modification in the process must be logically feasible, justifiable, and cost-effective [16]. Many studies have emphasized that implementation of ISPs in the early stages of the system lifecycle is more economical and practical due to the lack of need to change the equipment and process layout [1719].

In the 1970s, quantitative risk assessment (QRA) was initially introduced to the nuclear industry and then integrated into the process industry as legal criteria. During the last four decades, substantial emphasis has been placed on QRA [20]. QRA can be used to numerically quantify the existing risks, determine the unacceptable risks according to a numerical criterion, and prioritize the risk control measures [21]. QRA is based on identifying the incident scenarios and evaluating the risk by defining the probability of failure, the prospect of various consequences, and the potential impact of those consequences. The QRA risk is defined as a function of probability, frequency, and consequence of a particular incident scenario. The QRA is commonly presented as individual and societal risk. The individual risk expresses the likelihood of experiencing fatal effects at a given location and is not affected by the distribution of population in the area. In other words, the term “individual risk” is used to calculate the risk of fatality for someone at a specific location, assuming an employee is always present at the location [22, 23]. The societal risk is a measure of the risk that the incidents pose to the population and takes into account population distribution in the area. The societal risk is expressed in terms of the likelihood of incident outcomes that affects a given number of people in a single incident. Societal risk is expressed using an F-N curve, which indicates that the expected frequency (F) of incident scenarios resulted in the number of N or more fatalities. The x-axis of the F-N curve represents the number of fatalities (N), which is depicted on a logarithmic axis with a minimum value of 1. The y-axis of the F-N curve represents the cumulative frequency of events with the number of fatalities equal to N or more [6, 9, 20]. PHAST software package is a tool that helps in calculating QRA. In 2020, Wu et al. used PHAST to study quantitative analysis of ground flare [24]. Also, Wang and Ma (2021) applied PHAST to calculate quantitative risk for hydrogen refueling station [25]. In 2019, Shuang et al. used QRA to conduct fire and explosion risk assessment in urban natural gas pipeline [26]. In 2016, Tong et al. studied long-distance oil and gas pipelines [23]. However, PHAST software can be used to develop consequence modeling and draw out valid data on QRA.

Genetic Algorithm is an approach, in which the machine can simulate the natural selection mechanism. This is done by searching the problem space to find the better answer, not necessarily the optimal one [2730]. GA is widely used in studies, in which the stochastic search algorithms are applied to find the optimal or the most immediate answer. In such algorithms, optimization strategies are defined by mixed continuous discrete variables and discontinuous and/or nonconvex system spaces [30]. As with early random generations, initial responses produced randomized generations of children by modifying and combining initial parent answers [31, 32]. This cycle will continue until the full recognition of specified process conditions. Approaching the optimal solution is a necessary condition for the process completion [30]. In optimization problems with various processes, variables, and conditions, GA can be used to find the best/optimum operating conditions for equipment and process to enhance inherent safety without reducing the production of final product (s) at industries.

Numerous studies have been conducted into the assessment of ISPs in reducing both the frequency (more influenced by the vessels construction material) and severity (mostly affected by nature and volume of chemical and process condition) [25, 710]. Most of these studies reach the conclusion that the best chance to make process inherently safer is at the design phase of process life cycle. Implementation of these approaches in operation phase would be costly or, in the best situation, need to change the design layout that is not practically applicable in all processes. For instance, Chen (2011) cited the purpose of the Explosion-Proof Technology in Oil Storage and Transportation Devices [5]; Palaniappan et al. (2004) suggested that layout design is playing a key role on making process inherently safer [15]; Rathnayaka et al. (2014), Tugnoli et al. (2008), and Kossoy et al. (2012) suggested the same point of view as Palaniappan [7, 8, 18]. Also, Syaza et al. (2016) purposed a graphical approach to make process inherently safer through research and development phase of process design. Although the earlier phase of process design is the greater opportunity to make the process inherently safer, there is a chance to implement the inherent safety in other phases of process life cycle. The present study aims to simulate the genetic algorithm process optimization (GAPO) and evaluate its effects on quantitative risk and inherent safety status at the operation phase of a methane recovery unit.

2. Methodology

The present study was carried out at the methane recovery unit of a petrochemical plant in the south of Iran to simulate the GAPO and evaluate its effects on quantitative risk. In GAPO, the level of produced methane was kept constant compared to the existing state, providing a higher level of safety without reducing methane separation or any add-on. The study was carried out in the following steps (Figure 1).

Figure 2 represents the process flow diagram of the studied methane recovery unit. The process consists of two drums for the refrigeration cycle, a feeding drum, a demethanizer tower, and two heat exchangers. First, the solid particles and the natural gas moisture are removed by filters and dryers. The feed gas (with the composition listed in Table 1) is then cooled to −94°C by the propane refrigeration (blue in Figure 2) and pumped into different parts of the demethanizer unit. The demethanizer tower is heated with a reboiler, which is embedded at the bottom (first tray) of the tower. Eventually, methane is separated from natural gas based on its bubble point and is exhausted from the top of the tower. The liquid gas, which is known as C2+, is exhausted from the bottom of the tower.

CompositionMole fraction


In this study, GAPO was used to optimize propane refrigerant consumption. For this purpose, the initial process was optimized using the GA Module in MATLAB. The implementation feasibility was investigated in terms of compatibility of the stability and energy equilibrium equations with fluids and thermodynamics rules using Aspen HYSYS (Figure 3) process simulator. The statistical population in GA was 10, and calculations are completed for 10 generations.

The role of 8 variables in reducing propane refrigerant consumption was investigated in this study. The initial value and the range of their changes are presented in Table 2.

Input variableValue

X1: temperature of stream 22-3 (°C)−59.94
X2: molar flow of stream 22-8 (Kgmol/h)293306
X3: molar flow of stream 22-15 (Kgmol/h)129306
X4: temperature of stream 22–31 (°C)−17.22
X5: temperature of stream 22–34 (°C)−49.75
X6: boil up ratio0.8012
X7: outlet pressure of expender32.01
X8: pressure of top stage of column31.51

Expander is assumed to work in the range of 95%∼105% of its normal capacity. Reflux Ratio is considered to be in the range of 65%∼110% of its normal flow rate.

The intended range of variables for optimization of propane consumption was determined according to the sensitivity analysis performed in the PHAST-Risk software for quantitative risk reduction, as well as variables changes versus the objective function. The objective function of the study was to reduce the propane refrigerant consumption and modify operational conditions to reduce the quantitative risk. Therefore, the optimal answer with minimum quantitative risk has been achieved. In addition, during the initial demethanizer process, the temperature of lines 22–26 (Figure 1), which plays the role of a reboiler in the system, was changed so that the LP steam consumption is approaching zero. The constraints of the problem are shown in Table 3; the objective function has been defined as the ratio of final to initial propane mass flow (kg/hr).

ConstrainsBase case

C1 recovery % (from demethanizer (T-2211)) > 8082
Minimum approach temperature of 10-E-2211 A-F > 4.1°C4.3
Minimum approach temperature of 10-E-2212 A-F > 0.54°C0.54

Next, the feasibility of GAPO was confirmed by Aspen HYSYS simulation to check whether the initial process was simulated for Aspen-HYSYS validation. The equation of the state of mixture fluid was PRSV.

The HAZOP method was used to identify the process hazards in the studied process unit. To this end, the process and equipment were divided into separate sections and nodes according to the nature of operation, including process lines, process vessels, process equipment, offsite systems, emergency shutdown systems, isolation at battery limits, and interface with other facilities. Members of the HAZOP team, consisting of 17 process, safety, and mechanical engineers, identified the hazards and possible incident scenarios for each node (Table 4).

ScenarioCompositionDescriptionLeak size (mm)Frequency (occurrence/year)Phase to releaseOutcome consequence

Sc01-T01-SComposition from Table 1Small leak of demethanizer tower251.38E − 0.03GasJet fire
flash fire
vapor cloud

Sc02-T01-MComposition from Table 1Medium leak of demethanizer tower1005.4E − 5GasJet fire
flash fire
vapor cloud

Sc03-T01-CRComposition from Table 1Large leak of demethanizer tower6094.8E − 5GasJet fire
flash fire
vapor cloud

Sc04-D32-SC3H8 (propane)Small leak of MP REFR circulation drum251.38E − 0.03LiquidJet fire
flash fire
vapor cloud

Sc05-D32-MC3H8 (propane)Medium leak of MP REFR circulation drum1005.4E − 5LiquidJet fire
flash fire
pool fire
vapor cloud

Sc06-D32-CRC3H8 (propane)Large leak of MP REFR circulation drum6094.8E − 5LiquidFireball
flash fire
vapor cloud
jet fire
pool fire

Sc07-D31-SComposition from Table 1Small leak of feed HP separator251.38E − 0.032 phaseJet fire
flash fire
vapor cloud

Sc08-D31-MComposition from Table 1Medium leak of feed HP separator1005.4E − 52 phaseJet fire
flash fire
vapor cloud

Sc09-D31-CRComposition from Table 1Large leak of feed HP separator6094.8E − 52 phaseJet fire
flash fire
vapor cloud

Sc010-D33-SC3H8 (propane)Small leak of LP REFR. Circulation drum251.38E − 0.03LiquidJet fire
flash fire
pool fire
vapor cloud

Sc011-D33-MC3H8 (propane)Medium leak of LP REFR. Circulation drum1005.4E − 5LiquidJet fire
flash fire
pool fire
vapor cloud

Sc012-D33-CRC3H8 (propane)Large leak of LP REFR. Circulation drum6094.8E − 5LiquidJet fire
flash fire
pool fire
vapor cloud

The consequences of possible scenarios in the studied process unit were modeled using the PHAST software package. The information needed for consequence modeling was derived from the HAZOP study and the operational conditions contained in the P and IDs and PFDs (Table 5). Leakage size and their frequencies were selected based on the guidelines of OGP (Risk Assessment Data Directory, Release Failure Frequency, 2010) and DNV. The required climatic information was obtained from the hourly records during the last five consecutive years by the Meteorological Organization (Table 6). Due to weather difference, consequence modeling was carried out in two climatic conditions of spring/summer (Weather 1) and fall/winter (Weather 2).

EquipmentCapacity (m3)Tem. (°C)Pressure (bar)Inventory (kg)

MP REFR circulation drum49.1−13.53.826840
LP REFR circulation drum54.4−43.81.231800
Feed HP separator82.6−6055.68119

PeriodAverage temperature (°C)Relative humidity (%)Wind velocity (m/s)Prevailing wind directionStability class

Weather 131.7546.832.71WNWA
Weather 221.1352.252.69WNWB

In this stage, quantitative risks (societal and individual) were calculated using PHAST-Risk software for the initial process design (IPD) and proposed GAPO. The quantitative risk was assessed based on process information derived from simulated processes (Aspen-HYSYS), consequence modeling, and population dispersion. The aim was to assess the effectiveness of GAPO in reducing the risks.

This study implemented I2SI, presented by Faisal Khan et al. (2004), to assess the process inherent safety status in GAPO compared to IPD. I2SI is calculated by dividing the inherent safety potential index (ISPI) to the hazard index (HI) for each subset or equipment (equation 1)), and finally, for the entire system (equation (2)). An increase in the score of I2SI reflects the improvement of the inherent safety status [33].

In this study, after determining the scenario and identifying the high-risk equipment in the process, I2SI was calculated for optimized process with GA compared to the IPD. Inherent safety status of process and equipment was calculated through I2SI approach according to the field expert’s opinion and information from consequences modeling (for each equipment and for the whole system).

3. Results and Discussion

3.1. GAPO and Aspen-HYSYS

Figure 4 shows mean GA fitness generation for 10 populations in each module compared to the best fitness. Best fitness refers to the fitness of the best individual in the population compare to the goal and condition of study. Each generation provides a new average population fitness that is called mean fitness. In this study, as can be seen, after 10 generations, the mean fitness approached the best fitness. Meanwhile, not only the difference between the mean fitness and best fitness is as low as acceptable (about 0.017934), but also the mean fitness after seven generations has been stabilized, showing that the obtained final objective function is optimum.

The optimized values of the parameters selected as genetic algorithm generations are shown in Table 7. As indicated, by increasing the temperature of stream 22–31 (variable X4) from −17.22°C to −5.25°C, the required propane refrigerant mass flow rate increased due to refrigeration reduction in 10-E-2211 cold box. Finally, the refrigerant mass flow rate was decreased by increasing propane temperature in stream 22–26. As the temperature of stream 22–3 (inlet stream to 10-D-2231) (variable X1) increased from −59.94°C to −56°C, propane refrigeration decreased due to increasing conversion of liquid to gas. The increase in gas was caused partly by the stream entering into the 10-TE-2271 expander. Since the pressure drop in the expander caused a cool down, the streams also reduced propane refrigeration flow rate. The pressure drops at the tower from 31.51 barg to 29 barg caused a reduction in the propane refrigeration flow rate. By reducing both the tower pressure (variable X8) and boil-up ratio (variable X6), the temperature of the outlet methane stream decreased, and consequently, the required cooling reduction was achieved. As the boil-up ratio dropped from 0.8012 to 0.7900, the gas flow in the tower also declined, which resulted in a flow rate reduction at tower side streams. The decrease in the rate at side streams eventually reduced the amount of propane refrigeration.


Population size10
Selection methodStochastic uniform
Probability of crossover0.8
Number of generations10

Input variableOptimized values
X1: temperature of stream 22-3 (°C)−56
X2: molar flow of stream 22-8 (Kgmol/h) 260000
X3: molar flow of stream 22-15 (Kgmol/h) 85000
X4: temperature of stream 22–31 (°C)−15.25
X5: temperature of stream 22–34 (°C)−48
X6: boil-up ratio0.79
X7: outlet pressure of expender29.3
X8: pressure of top stage of column29.00

Objective function0.8

Increase decrease (changing from main amount).

Table 8 shows the process information of methane recovery unit after simulation with genetic algorithm. As can be seen, after implementing the GAPO model, vessels have been changed to the lowest volume and safest status. The drop-in values of all parameters shown in Table 7 indicate a reduction in energy consumption. In other words, while the percentage of methane separation (production of methane and C2+) remained constant, the amount of propane consumption was reduced by 20٪, which reduces the volume of vessel, thereby reducing quantitative risk and increasing inherent safety of the process. Therefore, a feasible approach to improve inherent safety at no cost and without decrease in production would be the application of the proposed genetic algorithm process optimization.

Involved equipmentCapacity (m3)Tem. (°C)Pressure (bar)Inventory (kg)

Demethanizer 11043−96.372929220
MP REFR circulation drum42.1−13.53.823020
LP REFR. Circulation drum27.85−43.81.216280
Feed HP separator70.81−5655.66352

3.2. Hazard Identification and Consequence Modeling

Table 4 presents the identified scenarios for consequence modeling extracted from HAZOP study. This table reveals a total of 12 scenarios and details about the 4 most hazardous equipment devices including demethanizer, MP REFR circulation drum, Feed HP separator, and LP REFR. The circulation drum was extracted from HAZAOP study.

Table 9 represents the various consequences of scenarios in the process optimized with GA and the IPD. As can be seen, the vapor cloud radius decreased in all scenarios related to the demethanizer tower (scenarios 1 to 3). The highest reduction was observed in scenario 1, where the vapor cloud radius dropped to 44% at 23780 ppm in GAPO. This reduction is related to the 23% and 20% related to scenarios 2 and 3, respectively. But in scenarios 4 and 5, due to constant temperature and pressure in the drum, there is no decrease in the vapor cloud. However, in the worst-case scenario for this drum (Scenario 6), there is a 22% reduction in the vapor cloud radius. The radius of the vapor cloud in the feed drum is not significantly different due to a similar reason. In LP drum, due to reduction in the temperature and inventory, the vapor cloud radius decreased from 21% to 31% (Scenarios 10 to 12) in GAPO approach compared to the initial process.


Sc01-T01-SVapor cloud-23780 ppm (m2)5531
Jet fire-lethality of 100% (m)2321
Flash fire radius (m)12–3210–26

Sc02-T01-MVapor cloud-23780 ppm (m2)40933141
Jet fire-lethality of 100% (m)8275
Flash fire radius (m)125–186110–168

Sc03-T01-CRVapor cloud-23780 ppm (m2)160756129135
Jet fire-lethality of 100% (m)410373
Flash fire radius (m)668–850625–784

Sc04-D32-SVapor cloud-10000 ppm (m2)179174
Jet fire-lethality of 100% (m)3131
Flash fire radius (m)21–3720–36

Sc05-D32-MVapor cloud-10000 ppm (m2)50285049
Jet fire-lethality of 100% (m)106106
Flash fire radius (m)99–13899–138
Pool fire zone-lethality of 100% (m2)50285028

Sc06-D32-CRVapor cloud-10000 ppm (m2)8075962686
Jet fire-lethality of 100% (m)478426
Flash fire radius (m)353–480313–426
Pool fire zone-lethality of 100% (m2)2060115386
Fireball zone-lethality of 100% (m2)1367711304

Sc07-D31-SVapor cloud-22370 ppm (m2)118106
Jet fire-lethality of 100%(m)2827
Flash fire radius (m)17–4416–42

Sc08-D31-MVapor cloud-22370 ppm (m2)61615863
Jet fire-lethality of 100% (m)9694
Flash fire radius (m)152–220148–218

Sc09-D31-CRVapor cloud-22370 ppm (m2)4348137566
Jet fire-lethality of 100% (m)480468
Flash fire radius (m)406–496390–467
Fireball zone-lethality of 100% (m2)4004934618

Sc010-D33-SVapor cloud-10000 ppm (m2)7552
Jet fire-lethality of 100% (m)2016
Flash fire radius (m)11–209–17
Pool fire zone-lethality of 100% (m2)452254

Sc011-D33-MVapor cloud-10000 ppm (m2)1195939
Jet fire-lethality of 100% (m)4033
Flash fire radius (m)28–4625–42
Pool fire zone-lethality of 100% (m2)32151962

Sc012-D33-CRVapor cloud-10000 ppm (m2)68004961
Jet fire-lethality of 100% (m)141133
Flash fire radius (m)76–11768–104
Pool fire zone-lethality of 100% (m2)58053419

After optimization of the studied methane process, the operating process of equipment was not significantly changed, and therefore, intensity of jet fire was not highly different. For example, in all three scenarios related to the demethanizer tower in GAPO, the jet fire was reduced by 9%. In scenarios for MR drum and feed drum, the amount of jet fire was almost identical. But in the LP drum scenarios, the jet fire dropped from 6% to 20% due to the reduction of temperature and inventory in GAPO. For the other fire scenarios, the intensity of incidents has decreased; for instance, intensity of pool fire decreased by 44% for LP drum (Scenario 11). Finally, implementing GAPO, it can reduce both the hazardous materials in the process and the capacity of the vessels while keeping production constant.

3.3. Cumulative Quantitative Risk Assessment

The individual risk contour of 1 pmpy, for GAPO compared to the IPD, is shown in Figure 5. As illustrated, implementation of GAPO was associated with reduction in 1 pmpy individual risk from 726473.54 m2 in IPD to 596901.44 m2 in GAPO. This decrease was more tangible in the downwind, due to the vapor cloud tending to spread in the wind direction; therefore, the individual risk was higher in the downwind incidents.

Figure 6 shows the individual risk transect for GAPO compared to the IPD in 200 m from the incident location, where the workers are resting. As can be seen, the individual risk at the employee’s resting place reduced from 69.5E − 05/Avgeyear in IPD to 03.5E − 05/Avgeyear in GAPO.

In Figure 7, the results of societal risks are presented in GAPO and supported by Aspen HYSYS, compared to the IPD, which exists in the operation phase. As shown, the frequency of deaths is not significantly different in GAPO from the existing process. This uniformity is because the present study focused on the severity of incident consequences. As indicated, in the IPD, the number of deaths is 200 in almost frequency of 105/Avgeyear, which is reduced to 107/Avgeyear in GAPO. The number of deaths in integrated consequences of 12 scenarios in the IPD reduced from 600 to 400 in GAPO, which shows a 33٪ decrease.

3.4. Inherent Safety Assessment

Results of assessing process for inherent safety by I2SI for IPD and GAPO are present in Table 10. As the tables show, reducing hazard index for demethanizer tower has great impact on total inherent safety of initial process. Since the demethanizer tower is the most hazardous equipment on the process, implementing inherent safety into the process will improve the safety of the whole operation. This table also shows that the total I2SI of process increased significantly in the optimized process without a need to add safety equipment or use conventional safety methods and procedures.

EquipmentHazard index (HI)Inherent safety potential index (ISPI)Integrated inherent safety index (I2SI)Total I2SI of process

Demethanizer tower3.
MP REFR circulation drum1.471.750.531.180.320.68
Feed HP separator2.002.400.491.450.250.60
LP REFR. circulation drum1.321.852.392.501.811.35

4. Conclusion

The aim of this study was to investigate the effect of genetic algorithm process optimization on inherent safety enhancement and reduction of quantitative risks in operation phase of process. The results showed that, after optimizing the process with the genetic algorithm, the number of deaths decreased by one-third without reducing the amount of methane production. In addition, process optimization led to reduced energy consumption and improved efficiency. In conclusion, genetic algorithms can be implemented at no extra cost in all phases of the process life cycle to optimize processes and equipment, especially during the operation phase.


GAPO:Genetic algorithm process optimization
HYSYS:Hyprotech system
PFD:Process flow diagram
NG:Natural gas
PHAST:Process hazard analysis software tools
I2SI:Integrated inherent safety index
IS:Inherent safety
HI:Hazard index
ISPI:Inherent safety potential index
IPD:Initial process design
ISP:Inherent safety principle
P and ID:Piping and instrumentation diagram
GA:Genetic algorithm
RCY:Recycle operator
TEE:Flow splitter
FB:Fire ball
MC:Maximum concentration
QRA:Quantitative risk assessment
HAZOP:Hazard and operability study
ESDV:Emergency shut down valve
NRV:None return valve
VLE:Vapor-liquid equilibrium
ISDS:Inherently safer design strategy
D-2232:MP refrigerant circulation drum
D-2231:FEED HP-separator
D-2233:LP refrigerant circulation
Pmpy:Part million per year
VLE:Vapor-liquid equilibrium
PRSV:Peng–Robinson-Stryjek-Vera equations of state
SIS:Safety instrument system
SOP:Standard operation procedure
DNV:Det Norske veritas
OGP:Oil and gas producers
P and ID:Piping and instrumentation diagram.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This article was extracted from the thesis written by Abolfaz Moghadasi, M.Sc. student of Occupational Health Engineering, and was financially supported by Shiraz University of Medical Sciences, Shiraz, Iran (Grant No. 95-01-04–13890).


  1. H. M. Hassim and H. N. Hanafi, “Current status and future direction of inherently healthier design in Malaysia,” Energy Procedia, vol. 14, pp. 1939–1944, 2012. View at: Publisher Site | Google Scholar
  2. P. R. Amyotte, F. I. Khan, and T. A. Kletz, “Inherently safer design activities over the past decade,” in Proceedings of the IChemE Symposium Series, Manchester, UK, November 2009. View at: Google Scholar
  3. F. I. Khan and P. R. Amyotte, “I2SI: a comprehensive quantitative tool for inherent safety and cost evaluation,” Journal of Loss Prevention in the Process Industries, vol. 18, no. 4–6, pp. 310–326, 2005. View at: Publisher Site | Google Scholar
  4. F. I. Khan and P. R. Amyotte, “How to make inherent safety practice a reality,” The Canadian Journal of Chemical Engineering, vol. 81, no. 1, pp. 2–16, 2003. View at: Google Scholar
  5. C. Sining, “Application of inherent safety explosion-proof technology in oil storage & transportation device,” Procedia Engineering, vol. 15, pp. 4814–4818, 2011. View at: Publisher Site | Google Scholar
  6. H. Pasman and G. Reniers, “Past, present and future of quantitative risk assessment (QRA) and the incentive it obtained from land-use planning (LUP),” Journal of Loss Prevention in the Process Industries, vol. 28, pp. 2–9, 2014. View at: Publisher Site | Google Scholar
  7. A. Tugnoli, F. Khan, P. Amyotte, and V. Cozzani, “Safety assessment in plant layout design using indexing approach: implementing inherent safety perspective,” Journal of Hazardous Materials, vol. 160, no. 1, pp. 100–109, 2008. View at: Publisher Site | Google Scholar
  8. S. Rathnayaka, F. Khan, and P. Amyotte, “Risk-based process plant design considering inherent safety,” Safety Science, vol. 70, pp. 438–464, 2014. View at: Publisher Site | Google Scholar
  9. A. R. Carpenter, R. A. Ogle, S. J. Dee, and B. Cox, “Inherently safer design: lessons learned about the principle of simplification,” in Proceedings of the 2014 AIChE Spring National Meeting, New Orleans, LA, USA, March 2014. View at: Google Scholar
  10. E. J. Bernechea and J. Arnaldos, “Optimizing the design of storage facilities through the application of ISD and QRA,” Process Safety and Environmental Protection, vol. 92, no. 6, pp. 598–615, 2014. View at: Publisher Site | Google Scholar
  11. T. A. Kletz, “Inherently safer design: the growth of an idea,” Process Safety Progress, vol. 15, no. 1, pp. 5–8, 1996. View at: Publisher Site | Google Scholar
  12. A. M. Shariff and C. T. Leong, “Inherent risk assessment-a new concept to evaluate risk in preliminary design stage,” Process Safety and Environmental Protection, vol. 87, no. 6, pp. 371–376, 2009. View at: Publisher Site | Google Scholar
  13. M. Gentile, W. J. Rogers, and M. S. Mannan, “Development of a fuzzy logic-based inherent safety index,” Process Safety and Environmental Protection, vol. 81, no. 6, pp. 444–456, 2003. View at: Publisher Site | Google Scholar
  14. X. Lin, B. Chen, B. Yingli, and B. Zhu, “Study on effect of product liability to inherent safety,” Procedia Engineering, vol. 45, pp. 271–275, 2012. View at: Publisher Site | Google Scholar
  15. C. Palaniappan, R. Srinivasan, and R. Tan, “Selection of inherently safer process routes: a case study,” Chemical Engineering and Processing: Process Intensification, vol. 43, no. 5, pp. 641–647, 2004. View at: Publisher Site | Google Scholar
  16. P. R. Amyotte, A. U. Goraya, D. C. Hendershot, and F. I. Khan, “Incorporation of inherent safety principles in process safety management,” Process Safety Progress, vol. 26, no. 4, pp. 333–346, 2007. View at: Publisher Site | Google Scholar
  17. R. Rusli and A. Mohd Shariff, “Qualitative assessment for inherently safer design (QAISD) at preliminary design stage,” Journal of Loss Prevention in the Process Industries, vol. 23, no. 1, pp. 157–165, 2010. View at: Publisher Site | Google Scholar
  18. A. A. Kossoy, Y. G. Akhmetshin, and A. I. Benin, “Applying mathematical simulation for design of inherently safe chemical processes: I. providing inherent safety,” Russian Journal of Applied Chemistry, vol. 85, no. 9, pp. 1466–1474, 2012. View at: Publisher Site | Google Scholar
  19. S. Shah, U. Fischer, and K. Hungerbühler, “Assessment of chemical process hazards in early design stages,” Journal of Loss Prevention in the Process Industries, vol. 18, no. 4–6, pp. 335–352, 2005. View at: Publisher Site | Google Scholar
  20. J. Di Domenico, C. A. Vaz, and M. B. De Souza, “Quantitative risk assessment integrated with process simulator for a new technology of methanol production plant using recycled CO2,” Journal of Hazardous Materials, vol. 274, pp. 164–172, 2014. View at: Publisher Site | Google Scholar
  21. M. F. Milazzo, C. Vianello, and G. Maschio, “Uncertainties in QRA: analysis of losses of containment from piping and implications on risk prevention and mitigation,” Journal of Loss Prevention in the Process Industries, vol. 36, pp. 98–107, 2015. View at: Publisher Site | Google Scholar
  22. L. Ma, L. Cheng, and M. Li, “Quantitative risk analysis of urban natural gas pipeline networks using geographical information systems,” Journal of Loss Prevention in the Process Industries, vol. 26, no. 6, pp. 1183–1192, 2013. View at: Publisher Site | Google Scholar
  23. S.-J. Tong, Z.-Z. Wu, R.-J. Wang, and H. Wu, “Fire risk study of long-distance oil and gas pipeline based on QRA,” Procedia Engineering, vol. 135, pp. 369–375, 2016. View at: Publisher Site | Google Scholar
  24. T. Wu, Q. Kong, X. Wei, X. Gao, and X. Xiong, “Quantitative analysis of enclosed ground flare based on PHAST analysis,” Heat Transfer, vol. 49, no. 5, pp. 2746–2769, 2020. View at: Publisher Site | Google Scholar
  25. X. Wang and X. Ma, “Risk control analysis of safety accident in hydrogen refueling station based on PHAST software,” IOP Conference Series: Earth and Environmental Science, vol. 680, no. 1, p. 012119, 2021. View at: Google Scholar
  26. S. Li, C. Cheng, G. Pu, and B. Chen, “QRA-grid: quantitative risk analysis and grid-based pre-warning model for urban natural gas pipeline,” ISPRS International Journal of Geo-Information, vol. 8, no. 3, p. 122, 2019. View at: Publisher Site | Google Scholar
  27. P. Baybutt, “A critique of the hazard and operability (HAZOP) study,” Journal of Loss Prevention in the Process Industries, vol. 33, pp. 52–58, 2015. View at: Publisher Site | Google Scholar
  28. M. G. Jacobsen and S. Skogestad, “Active constraint regions for optimal operation of chemical processes,” Industrial & Engineering Chemistry Research, vol. 50, no. 19, pp. 11226–11236, 2011. View at: Publisher Site | Google Scholar
  29. M. V. Pathan, S. Patsias, and V. L. Tagarielli, “A real-coded genetic algorithm for optimizing the damping response of composite laminates,” Computers & Structures, vol. 198, pp. 51–60, 2018. View at: Publisher Site | Google Scholar
  30. R. Wang, Y. Wu, Y. Wang, and X. Feng, “An industrial area layout design methodology considering piping and safety using genetic algorithm,” Journal of Cleaner Production, vol. 167, pp. 23–31, 2017. View at: Publisher Site | Google Scholar
  31. A. Velásco-Mejía, V. Vallejo-Becerra, A. U. Chávez-Ramírez, J. Torres-González, Y. Reyes-Vidal, and F. Castañeda-Zaldivar, “Modeling and optimization of a pharmaceutical crystallization process by using neural networks and genetic algorithms,” Powder Technology, vol. 292, pp. 122–128, 2016. View at: Publisher Site | Google Scholar
  32. A. C. Caputo, P. M. Pelagagge, M. Palumbo, and P. Salini, “Safety-based process plant layout using genetic algorithm,” Journal of Loss Prevention in the Process Industries, vol. 34, pp. 139–150, 2015. View at: Publisher Site | Google Scholar
  33. F. I. Khan and P. R. Amyotte, “Integrated inherent safety index (I2SI): a tool for inherent safety evaluation,” Process Safety Progress, vol. 23, no. 2, pp. 136–148, 2004. View at: Publisher Site | Google Scholar

Copyright © 2021 Mehdi Jahangiri et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.