Hydrogen Leakage Simulation and Risk Analysis of Hydrogen Fueling Station in China

: Hydrogen is a renewable energy source with various features, clean, carbon-free, high energy density, which is being recognized internationally as a “future energy.” The US, the EU, Japan, South Korea, China, and other countries or regions are gradually clarifying the development position of hydrogen. The rapid development of the hydrogen energy industry requires more hydrogenation infrastructure to meet the hydrogenation need of hydrogen fuel cell vehicles. Nevertheless, due to the frequent occurrence of hydrogen infrastructure accidents, their safety has become an obstacle to large-scale construction. This paper analyzed ﬁve sizes (diameters of 0.068 mm, 0.215 mm, 0.68 mm, 2.15 mm, and 6.8 mm) of hydrogen leakage in the hydrogen fueling station using Quantitative Risk Assessment (QRA) and HyRAM software. The results show that unignited leaks occur most frequently; leaks caused by ﬂanges, valves, instruments, compressors, and ﬁlters occur more frequently; and the risk indicator of thermal radiation accident and structure collapse accident caused by over-pressure exceeds the Chinese individual acceptable risk standard and the risk indicator of a thermal radiation accident and head impact accident caused by overpressure is below the Chinese standard. On the other hand, we simulated the consequences of hydrogen leak from the 45 MPa hydrogen storage vessels by the physic module of HyRAM and obtained the ranges of plume dispersion, jet ﬁre, radiative heat ﬂux, and unconﬁned overpressure. We suggest targeted preventive measures and safety distance to provide references for hydrogen fueling stations’ safe construction and operation.


Introduction
In the context of carbon neutrality, the global energy structure is shifting from fossil fuels to renewable energy to reduce carbon dioxide emissions [1].Hydrogen is an important way to support renewable energy development due to its advantages of extensive source, clean and carbon-neutral, high calorific value, and abundant application scenarios [2].Countries around the world are accelerating the pace of hydrogen energy development.The US promulgated the "Hydrogen Roadmap U.S.", Japan proposed to create a "Hydrogen based society" in the "Basic Hydrogen Strategy", South Korea released the "Hydrogen Economy Development", in which it planned to establish a "Hydrogen Economy Demonstration City".China proposed the "Medium and Long-term Plan for the Development of Hydrogen Energy Industry (2021-2035)" in March 2022, which established hydrogen's strategic position and development goals.The planning goal is 50,000 fuel cell vehicles by 2025 [3].China is prioritizing the promotion of hydrogen fuel cell buses in several cities such as Beijing and Guangzhou.According to the "National Demonstration City Cluster of Hydrogen Fuel Cell Vehicle Statistic and Analysis Report (2022.4)",as of 30 April 2022, China's monitoring and management platform for new energy vehicles had accessed 8198 hydrogen fuel cell vehicles, of which 4241 were buses, accounting for 51% [4].The growth of hydrogen fuel cell buses has begun to slow down because of various reasons such as economic situation, technology, and supporting infrastructure.A hydrogen fueling station is a special place for filling hydrogen fuel for hydrogen storage bottles of hydrogen fuel cell vehicles, hydrogen internal combustion engine vehicles, or hydrogen and natural gas hybrid fuel vehicles.As of March 2022, with the support of Chinese policies, the number of hydrogens fueling stations reached 264 [5].Due to accidents such as the explosion of the Sandvika hydrogen fueling station in Norway [6] and the hydrogen tank explosion in Gangneung, South Korea [7], the safety of the hydrogen fueling station has become a major obstacle to hydrogen promotion.The main accidents of hydrogen fueling station are caused by hydrogen leaks, which can produce serious loss of life or property.For pushing hydrogen energy commercialization and hydrogen infrastructure construction, the risk assessment of hydrogen fueling stations has attracted extensive attention of many experts and scholars.
Risk analysis of hydrogen fueling station usually adopts the QRA method.The safety distance of hydrogen fueling stations means the fire separation distance for buildings, that is the interval distance to prevent a building on fire from igniting adjacent buildings within a certain period of time which is convenient for firefighting, and it is the minimum distance to avoid direct or indirect damage when a disaster occurs.Sakamoto et al. designed the physical model of the hydrogen fueling station, analyzed the scenario of multiple safety measures failure, and determined the combinations of failures that could cause accidents [8].Gye et al. analyzed main risks of the high-pressure hydrogen fueling station in a densely populated area [9].Tsunemi et al. conducted a consequence analysis and probability analysis on the hydrogen leakage scenarios in conjunction with failure trees, and concluded safety distances of 6 m, 11-12 m, and 11-12 m for the dispenser, accumulator, and compressor, respectively [10].Kikukawa et al. analyzed the risk of hydrogen leakage in a 70 MPa hydrogen fueling station by extrapolation from a 35 MPa hydrogen fueling station, identified 721 failure scenarios by HAZOP (Hazard and Operability Studies) and FMEA (Failure Mode and Effect Analysis), and determined the effectiveness of safety measures in risk reduction if applied accurately and with a safety distance of 6 m [11].Yuan et al. calculated the fatality probability, individual risk, and social risk of jet fire and explosion based on the actual situation of a hydrogen fueling station [12].Li et al. quantitively analyzed the consequences of physical explosion, flash, jet fire, and vapor cloud explosion of hydrogen fueling station under different hydrogen storage pressures, leak diameters, and wind velocities, and drew the effect rules [13].
With the continuous improvement of risk analysis software, quantitative risk analysis of hydrogen fueling station is increasingly integrated with HyRAM, Phast/Safeti, FLACS, etc. Kwon et al. concluded the hazards influence regions of hydrogen leakage from core equipment combining Hy-KoRAM with Phast/Safeti [14].Kim et al. concluded that the maximum distance of jet fire and heat radiation were 8.2 m and 10.6 m, respectively, which were caused by the 5 mm 2 hydrogen leakage of the compressor at 90 MPa using Hy-KoRAM, and the individual risk metrics were in the ALARP (as low as reasonably practicable) region and acceptable region [15].Kim et al. simulated several scenarios of hydrogen leakage through FLACS and concluded the maximum expansion distance of the explosion was 17 m [16].Zhou et al. used CASST-QRA to conclude the effect areas of the explosion consequences under the fire prevention distances given by the national standard and design institute [17].
Comparing software simulation results with experiments for risk analysis, Park et al. designed experiments by thermal imaging camera and Schlieren imaging system and compared them with the simulation results of HyRAM, which yielded relatively conservative software results, and used RISKCURVES for analyzing individual and social risk to determine that the safety distance is 29 m between hydrogen fueling station and surroundings [18].
The above studies lack the quantitative analysis of the hydrogen leakage risk from hydrogen fueling stations and fire separation distances designed under the current Chinese regulations and the hydrogen fueling demands of hydrogen fuel cell buses.Since the development of hydrogen fueling stations started late in China, the research on quantitative risk analysis of hydrogen fueling stations is insufficient.Although the number of hydrogens fueling stations has ranked among the top in the world, the safety problems are prominent.There is an urgent need for quantitative risk analysis of hydrogen fueling stations that is suitable for Chinese conditions.
This study designed the hydrogen fueling station scenario according to the Chinese standard specifications and hydrogen fuel cell buses demand and calculated and analyzed the accidents frequencies, risk indicators, components failure frequencies, and physical models of consequences, resulting in more targeted recommendations for component precautionary measures and safety precautionary distance combining Quantitative Risk Assessment method with HyRAM developed by Sandia National Laboratories in the US.

QRA Method
The QRA method is widely used in safety management, provides a technical means of quantitative risk assessment, and is suitable for hydrogen refueling stations [8][9][10][11][12][13][14][15][16][17][18]. Figure 1 shows a flow chart of our study.We combined the original QRA procedure (scenario definition, risk identification, probability analysis, and consequence analysis) with the physical models of hydrogen leakage consequences to assess risk and propose preventive measures for reducing risks.

Scenario Simulation Conditions
The simulation scenario is hydrogen leakage at a 35 MPa hydrogen fueling station.The station is 90 m in length and 50 m in width.We set a non-combustible solid wall with a height of 2.5 m, referring to "Hydrogen Fueling Station Technical Specifications" [19].Assuming that the hydrogen fueling station operates continuously throughout the year, the station is staffed by three employees at all times and operates for 8760 h all year round.
China has two types of hydrogen fueling stations, one is the on-site station and the other is the off-site station.The on-site station express hydrogen production equipment in the hydrogen fueling station can produce hydrogen by itself and realize compression and filling.The off-site hydrogen fueling station refers to the delivery of hydrogen from outside of the station to the hydrogen fueling station, and filling hydrogen after secondary pressurization.The transportation modes include high-pressure gaseous hydrogen longtube trailers, liquid hydrogen tankers, and pipeline.Due to immature technology and high cost of liquid hydrogen tankers and pipeline, most of the hydrogen fueling stations built in China are off-site stations.Therefore, we designed an off-site hydrogen fueling station in this manuscript to enhance the applicability of the results to the actual situation in China.
The hydrogen fueling station is grade three, with a capacity of less than 3000 kg and a single tank capacity of less than 800 kg.The equipment configuration is shown in Table 1.The total gaseous hydrogen mass of cylinder assemblies for storage of gaseous hydrogen is 464.41 kg at 20 °C.Assuming a hydrogen fuel cell bus's gaseous hydrogen storage capacity is 20 kg, this station can refill 18 hydrogen fuel cell buses every day.The annual demand for this hydrogen fueling station is 6570 vehicles.According to the "Hydrogen Fueling Station Technical Specification" [19], the connection of hydrogen pipelines adopts weld fittings or bite-type fittings that have passed the hydrogen compatibility assessment.The connection between hydrogen pipelines, equipment, and valves shall be flanged or threaded.According to "Compressed Hydrogen Dispenser for Automobiles" [20], the hydrogen dispenser needs to install a gas filter between the gas source interface and the intake valve.Referring to the specifications for outer diameter and thickness of "Dimensions, Shapes, Masses and Tolerances of Seamless Steel Tubes" [21] and "Seamless Stainless Steel Pipes for Fluid Transport" [22], we set the pipe's outer diameter to 10 mm and the pipe's wall thickness to 1.6 mm.The layout of the hydrogen fueling station strictly abides by the fire separation specified in the "Hydrogen Fueling Station Technical Specification" [19].The process flow diagram and layout are shown in Figure 2. The tube trailer for gaseous hydrogen transports the hydrogen to the hydrogen unloading area, then connects the unloading hoses.The hydrogen enters the 45 MPa hydrogen compressor through the hydrogen pipeline, and after being pressed, enters the 45 MPa cylinder assemblies for storage of gaseous hydrogen, which are divided into three types: high, medium, and low for graded filling.The ratio of the number of graded cylinders is 5:2:2, thereby improving the filling efficiency and reducing the cost of hydrogen storage [23].Where the hydrogen fueling cell vehicle is filling, the cylinder uses the pressure difference to send the gaseous hydrogen to the 35 MPa hydrogen dispensers along the pipeline.

Probability Models
The hydrogen leakage scenarios were divided into five sizes by HyRAM, which were 0.01%, 0.1%, 1%, 10%, and 100% of the pipeline's cross-sectional area (inner diameters of 0.068 mm, 0.215 mm, 0.68 mm, 2.15 mm, and 6.8 mm) [24].In the 0.01%, 0.1%, 1%, and 10% leakage scenarios, the frequency of each component is independent, while for the 100% leakage scenario, accidents and shutdown failure factors are added, and the vehicle hydrogenation demand is considered, the leakage frequency is calculated as the following equations [24].
Exposure to the effects of thermal radiation and overpressure from hydrogen leaks can cause harm or death.The harm models of thermal radiation and overpressure incorporate the probability formula to calculate the probability of harm or death in accidents.The probability formula is: where  is the normal cumulative distribution function [24].
Thermal radiation damages personnel, equipment, and structures listed in Table 2.Because the exposure time required for damage is more than 30 min, people can complete the evacuation before the structure is damaged.Therefore, we analyzed the influence of thermal radiation on personnel.Thermal radiation's damage degree is shown in Table 3.The formula of harm from radiant heat fluxes is: where V is the thermal dose unit, as an index of harm measuring, I is the radiant heat flux in W/m 2 , and t is the exposure time in seconds [24].Four thermal probability models are encoded in HyRAM, which are Eisenberg, Tsao and Perry, TNO, and Lees models.Hydrogen leakage will produce infrared radiation and ultraviolet radiation, but the infrared radiation in the hydrogen fires is far less than the ultraviolet radiation.The result of the Tsao and Perry model would be conservative for exposure to hydrogen fires.Although the Eisenberg model only considers the influence of ultraviolet radiation, the results have better accuracy and generality at present [28].Therefore, we adopted the Eisenberg model [29] in HyRAM.The formula is: The possible accident effect of overpressure has two types, including direct effect and indirect effect, as shown in Table 3.According to Table 3, we concluded that the pressure required to throw a person against obstacles or to generate missiles that can penetrate the skin is far less than the pressure required to cause fatal lung damage.In addition, persons in the building structure have a higher probability of being crushed to death.Therefore, we adopted Models (7) and (8) in HyRAM.
TNO-Head impact model [30]: where Ps is peak overpressure in Pa, i is the impulse of the shock wave in Pa×s.

Risk Indicator
According to the harm models selected above, we used the Potential Loss of Life (PLL) to calculate the risk level of thermal radiation and overpressure accidents and analyzed the results referring to the individual acceptable risk criteria in China.The PLL is calculated as follows [24]: where n is one of the possible accident scenarios, fn is the frequency of scenario n, and cn is the expected number of fatalities for scenario n.
The hydrogen fueling station is a low-density place since the number of staff and customers is less than 30.The individual acceptable risk criterion in China for new stations is 1 × 10 −5 fatalities/year in this place [31].

Gas Plume Dispersion
HyRAM uses the physical model of plume dispersion described by Houf and Winters [32].The model is a one-dimensional model that takes into account buoyancy effects but does not consider the aspect of wind [24].The jet plume models are Equations ( 10 ) where: v is the velocity of hydrogen, ρ is the density of hydrogen, the subscript amb denotes ambient, the subscript cl is the centerline, B is a characteristic half-width of the jet with half value of vcl,  is the ratio of density spreading relative to velocity, r is a radius that is perpendicular to the stream-wise direction, Y is the mass fraction of hydrogen.
We analyzed the diffusion distance and mole fraction of hydrogen leak in five leak size scenarios, including 0.068 mm, 0.215 mm, 0.68 mm, 2.15 mm, and 6.8 mm leak sizes in hydrogen storage vessels with a pressure of 45 MPa.We set the contour mole fraction to 0.04 in HyRAM, which is the lower explosion limit (LEL) of hydrogen.

Jet fire Temperature and Trajectory
The jet fire physical models [24,33,34] used in HyRAM are Equations ( 13)- (18). where: L * is a non-dimensional flame length, Lvis is the visible flame length, fs is the mass fraction of fuel in a stoichiometric mixture of fuel and air, dj is the orifice diameter, ρj and ρamb are the density of fuel at the orifice and air, qrad (x,r) is the radiant heat flux measured at a particular axial location, x, and radial position, r, Srad is the total emitted radiative power from a flame, C * is the normalized radiative heat flux, Xrad is the radiant fraction, m • fuel is the mass flow rate of fuel, ∆Hc is the heat of combustion (120 MJ/kg for hydrogen [35]), τf is the flame residence time, ap is the Plank-mean absorption coefficient for an optically thin flame (0.23 for hydrogen [36]), Tad is the adiabatic flame temperature, ρf is the flame density, uj is the velocity of the jet at the exit (orifice), Pamb is the ambient pressure, Wmix is the mean molecular weight of the stoichiometric products of hydrogen combustion in air, R is the universal gas constant.

Radiative Heat Flux
Heat radiation is the main heat transfer mechanism in hydrogen fires.The curved flame of radiative heat flux from the buoyancy corrected is calculated by a weighted multisource model [37], which is calculated as: where: VF is the view factor, proportional to the heat flux transmitted to the observer, τ is the transmissivity, Af is the surface area of the flame, the constraint is 1 nN  , , wi is the emitter strength weighting parameter, D and β are the distance and angle, respectively, between the observer and unit normal to the point emitter.

Unconfined Overpressure
The hydrogen storage area of the station is an open area, then the scenario of overpressure due to delayed ignition of leaking hydrogen is calculated by the TNT-equivalence method [27], which estimates the effects of an explosion relating to the trinitrotoluene explosion.The formulas are Equations ( 22) and (23). where: is the mass of TNT that contains the same energy as combusted gaseous hydrogen, Fequiv is the equivalence factor, which is set as 3% [38], ∆Hc,TNT = 4.68 MJ/kg [24], R is the distance of overpressure, TNT R  is the scaled distance of the equivalent mass of TNT.

Consequence Analysis
The main risks of hydrogen fueling station are fire and explosion triggered by gaseous hydrogen leak.The volume of the leak determines the degree of accident losses.Highpressure hydrogen storage vessels, compressors, dispensers, hoses, trailers, etc. all have the possibility of hydrogen leakage.Hydrogen leakage causes frostbite, suffocation, jet fire, flash fire, and explosion accidents.
In cases where a hydrogen leak is detected and isolated by safety devices, the leaks are shut down, and no accidents occur generally.When hydrogen leakage is not detected and isolated, if there is no ignition source, a non-ignition leakage accident will occur, which may cause frostbite or suffocation.If there is an ignition source to ignite immediately, a jet flame will occur causing a personal burn accident.If the ignition is delayed, the hydrogen will accumulate, in the case of obvious overpressure, a vapor cloud explosion accident will occur, causing shock injuries to personnel; in the case of hydrogen accumulation without overpressure, flash fire accidents will occur.When the hydrogen leakage rate is less than 0.125 kg/s, the probability of immediate ignition is 0.8% and the probability of delayed ignition is 0.4%.When the hydrogen leakage rate is between 0.125 kg/s to 6.25 kg/s, the probability of immediate ignition is 5.3% and the probability of delayed ignition is 2.7%.When the hydrogen leakage rate is greater than 6.25 kg/s, the probability of immediate ignition is 23%, and the probability of delayed ignition is 12% [39].
Table 4 shows the results of risk indicators calculated by HyRAM.The PLL value under the conditions of Eisenberg and TNO-Structure Collapse models is 3.175 × 10 −5 , which is higher than the standard (1 × 10 −5 ) and requires appropriate prevention measures to reduce the risk.The risk contribution refers to the ratio of each accident scenario to all accident scenarios.Among the hazards of collapse and thermal radiation, the highest risk contribution is 94.932% for the explosion scenario triggered by a 100% hydrogen leak size.The PLL value under the conditions of Eisenberg and TNO-Head Impact models is 4.414 × 10 −9 , which is far below the standard, and the risk contribution is concentrated in the jet fire scenario triggered by a 100% hydrogen leak size.

Probability Analysis
Based on the random leakage frequency data of components integrated into the Hy-RAM, we can get a probability rank of consequences, including timely shutdown, jet fire, explosion, and unignited leakage.The probability of leakage detected and isolated is 90%, which is default value of HyRAM [24].
As is shown in Figure 3, the accidents are mostly concentrated in leakages of small sizes, such as 0.068 mm and 0.215 mm.Among them, the "No Ignition" of 0.068 mm leakage scene occurred the most, reaching 0.257 times/year.Among the five leakage scenarios of 0.068 mm, 0.215 mm, 0.68 mm, 2.15 mm, and 6.8 mm, except for the consequence of leakage detected and isolated, the proportions of "No Ignition" are the largest.Due to the interaction between dispensers and human, the probabilities of jet fire and explosion are greatly increased in the 6.8 mm leakage scenario.In the small size leakage scenarios of 0.01% and 0.1%, the component with the highest failure frequency is the flange, which is far higher than other components, followed by the compressor, valve, and instrument.In the 1% leakage scenario, the highest failure fre-quency is still the flange.However, as the leakage size increases, the flange's failure frequency drops sharply.Meanwhile, the failure frequencies of compressors, valves, and instruments are prominent.In the 10% and 100% size leakage scenarios, the filter's failure frequency exceeds the flange's and is the highest frequency, followed by the instrument.In addition, in the leakage scenario of 100% size, not only the failure of each component needs to be considered, but also the following scenarios should be taken into account, namely the pipeline rupture caused by overpressure during filling, the customer driving away without completing the filling process, or the nozzle release and manual valve and solenoid valve shutdown failure.The frequencies of factors that may cause "accidents and shutdown failure" are introduced in Figure 4f, where nozzle release, manual valve, and solenoid valve shutdown failure all have a high frequency.According to Equations ( 1)-( 3), the frequency of accidents and shutdown failure is 1.8 × 10 −5 times/year, which is much lower than the other components.Therefore, the prevention of hydrogen leakage accidents in hydrogen fueling stations needs to consider not only the core equipment, but also the purchase, installation, and maintenance of flanges, valves, instruments, and filters.

Plume Dispersion
Figure 5 shows the hydrogen plume dispersion of five leak sizes in 45 MPa.At the same pressure, the hydrogen diffusion distance increased as the leak diameter increased.The maximum diffusion distance triggered by a 6.8 mm leak diameter is approximately 33 m, while in the leak diameter of 0.068 mm, the maximum diffusion distance is only approximately 0.35 m.

Jet Flame Temperature and Trajectory
Figure 6 shows the simulated jet flame shape of hydrogen for five leak sizes.We took into account the pressure of 45 MPa of the hydrogen storage vessels.The maximum distance of the jet flame for the 0.068 mm, 0.215 mm, 0.68 mm, 2.15 mm, and 6.8 mm leak diameters reached approximately 0.17 m, 0.6 m, 1.7 m, 5.5 m, and 17.5 m, respectively.Without considering the barrier of the partition wall, the jet flame has influenced the refilling area and unloading area, and there are risks of burn injury to staff and customers.

Radiative Heat Flux Effect
HyRAM is used to analyze the influence region of radiant heat flux.The heat flux rate is divided according to Table 2.The influence region of radiative heat flux for a 6.8 mm leak diameter at a pressure of 45 MPa is shown in Figure 7.In the simulation scenario of the most serious 100% pipe cross-sectional area leakage in five sizes, the radiative heat flux of 35 kW/m 2 was transmitted up to 15m, which has a 1% probability of death within 10 s of exposure in this area.Therefore, this region belongs to a high-risk region.The area of heat flux propagation (i.e., heat flow rate over 4 kW/m 2 ) where first-degree burns and more serious consequences may occur to the human is more than 30 m.The staff of the station should maintain a safe distance of at least 15 m and be able to escape more than 30 m within the 20 s to prevent the occurrence of first-degree burn hazards when a 6.8 mm diameter leak occurs at 45 MPa pressure.

Unconfined Overpressure
The result of the simulated overpressure is shown in Figure 8.We analyzed the leakage scenario of 6.8 mm diameter at 45 MPa, the most serious one among the five leak sizes.The radius of the overpressure area is more than 15 m, which is affected by the delayed ignition of the leaking gaseous hydrogen over 5 kPa.Among this area, the radius is approximately 7.5 m of the region affected by overpressure higher than 16 kPa, which has risks of structure collapse, possible fatality by being projected against obstacles, and skin lacerations by missiles.The diameter of the area higher than 70 kPa is approximately 5 m, which is a high-risk area with the risks of destruction of buildings, heavy machinery damage, and a 100% probability of fatality from missile wounds.We considered the above four hazards overall and integrated the specific effect regions of the above four models in Figure 9.The region of radius 17.5 m that the jet fire reached is the most harmful area for people with the risk of burns or even death.People should avoid staying here for a long time.The gaseous hydrogen plume diffused the farthest distance, reaching 33 m.This area must forbid ignition sources.The radius of the radiant heat flow area reached 30 m, and the damage caused in this area takes 20 s to accumulate.Therefore, people should escape from this area within 20 s when a leak occurs.Referring to the fire separation distance between the hydrogen storage vessel and other facilities in the station as specified in the "Hydrogen Fueling Station Technical Specification" [19], the fire separation distance of 8 m between the third-level hydrogen storage vessel and dispenser is less than the effect range of 17.5 m of jet fire.Due to the high frequency of interaction between the hydrogen dispenser and people, the risk caused by hydrogen leakage under the condition of 6.8 mm diameter and 45 MPa pressure will do terrible damage to the staff of the station and customers, without fire barriers.

Conclusions
This paper adopted the QRA method to analyze the risk of hydrogen fueling stations designed referring to Chinese specifications, combined with HyRAM.We calculated that the PLL value of thermal radiation and structure collapse accidents caused by hydrogen overpressure explosion exceeded the Chinese individual acceptable risk standard and the PLL value of thermal radiation and head impact accidents caused by hydrogen overpressure explosion was lower than the standard.We concluded the distribution of high-incidence consequences of five leak sizes and high-frequency failure components triggering hydrogen leak, including flanges, valves, compressors, instruments, and filters.Based on the above results, we proposed specific recommendations for safe hydrogen fueling station construction and operation.
The purchase of flanges, valves, compressors, instruments, and filters should strictly follow the specified parameters.The installation should strictly comply with the design torque, protect the joint surface, and ensure the welding quality, among which the installation of instruments should also abide by the "Technical code of construction of instrumentation engineering for petrochemical industry" [40].Pay attention to the anti-corrosion of pressure vessels by controlling the purity of gaseous hydrogen strictly.Set up a safe, effective, and reasonable control system to ensure smooth connection and control of monitoring, alarm, cut-off, and exhaust system.
We simulated the physical scenarios of plume dispersion, jet fire temperature and trajectory, radiative heat flux, and unconfined overpressure by the physic module of Hy-RAM, and concluded as follows: (1) As the leak size increases, the distance of hydrogen diffusion increases.The distance reached 33 m in the scenario of a 100% leak, 6.8 mm pipe cross-sectional area at 45 MPa.Ignition sources should be prohibited in this area.(2) In the simulated jet fire scenario, the jet distance of flame was positively correlated with the leak size.The flame jet distance exceeded 17.5 m and the central temperature exceeded 2000 K in the 6.8 mm leak size and 45 MPa.The staff and customers should avoid staying here for a long time without fire barriers.(3) The radiant heat flux of the high-risk area was more than 15 m horizontally along the leakage direction, reached 5 m perpendicular to the leakage direction, and rose more than 6 m under the condition of 6.8 mm and 45 MPa, which had a probability of death.The horizontal direction of the area that would cause first-degree burns to people reaching 30 m. (4) The radius of the affected area that would cause indirect injury to people was 7.5 m approximately in the accident simulation of the unconfined overpressure caused by delayed ignition and accumulation of leaking hydrogen, with 6.8 mm leak diameter and 45 MPa.The staff must wear safety helmets and stay away from buildings when they are in this area due to work requirements.
Considering the calculation of risk indicators (PLL value) and the hazard simulation, Figures 5-9, overall, the fire-proof and explosion-proof design of the buildings should be strictly controlled at the early stage of hydrogen fueling station construction, especially in the area where people stay for a long time, to reduce the occurrence of structure collapse fatal accidents.When a leakage accident occurs without being dealt with in a timely manner, staff and customers must keep a safe distance of more than 17.5 m from the leakage source, and should not stay more than 20 s within a 30 m radius of the leakage source.

Figure 1 .
Figure 1.QRA flow chart for the hydrogen fueling station.

Figure 2 .
Figure 2. Hydrogen fueling station process flow and layout.

Figure 3 .
Figure 3. Scenario outcome: (a) Frequencies of scenarios outcome, (b) accident distribution of 0.01%, 0.1%, 1%, or 10% release, (c) accident distribution of 100% release.The cut-set refers to the predicated failure frequency of components in each scenario.It is shown in Figure 4.

Figure 5 .
Figure 5. Simulated hydrogen plume dispersion in five leak sizes at 45 MPa.

Figure 6 .
Figure 6.Simulated jet flame temperature/trajectory of five leak sizes at 45 MPa.

Figure 7 .
Figure 7. Simulated radiative heat flux distribution for a 6.8 mm leak diameter at 45 MPa.

Figure 8 .
Figure 8. Simulated radiative heat flux distribution for 6.8 mm leak diameter at 45 MPa.

Figure 9 .
Figure 9.Effect area of accidents on the hydrogen fueling station.

Table 2 .
Radiant heat flux harm criteria for people.

Table 3 .
Damage to humans, structures, and equipment from overpressure events.