收稿日期：2022-08-12

基金項(xiàng)目：國家重點(diǎn)研發(fā)計(jì)劃（2019YFB2006305）

通信作者：王海清（1974—），男，博士、教授，主要從事新能源風(fēng)險(xiǎn)防控、工藝災(zāi)害分析等方面的研究。wanghaiqing@upc.edu.cn

DOI：10.19912/j.0254-0096.tynxb.2022-1206 文章編號：0254-0096（2023）11-0459-06

摘 要：針對不同領(lǐng)域國家標(biāo)準(zhǔn)中安全距離計(jì)算的差異性問題對某典型加氫站配置模型進(jìn)行定量風(fēng)險(xiǎn)分析（QRA），計(jì)算出加氫站設(shè)施與站外防護(hù)目標(biāo)安全距離和個人風(fēng)險(xiǎn)，實(shí)現(xiàn)多領(lǐng)域（氫能領(lǐng)域與危險(xiǎn)化學(xué)品領(lǐng)域）標(biāo)準(zhǔn)交叉驗(yàn)證并通過揭示加氫站氫能標(biāo)準(zhǔn)和?；窐?biāo)準(zhǔn)的適配性差異，為加氫站的設(shè)計(jì)、建設(shè)和風(fēng)險(xiǎn)控制提供多維度參考。

關(guān)鍵詞：氫；安全；加氫；定量風(fēng)險(xiǎn)分析；多標(biāo)準(zhǔn)交叉驗(yàn)證

中圖分類號：X937""""""""""" """""""""""" """"""文獻(xiàn)標(biāo)志碼：A

0 引 言

中國對氫能源需求日益增長，促進(jìn)加氫站高速發(fā)展，保證加氫站建設(shè)安全是關(guān)鍵的問題。目前已有文獻(xiàn)對加氫站（hydrogen refueling station，HRS）進(jìn)行研究，Ham等［1］總結(jié)出加氫站多種風(fēng)險(xiǎn)評估方法（包括定量風(fēng)險(xiǎn)分析（quantitative risk analysis，QRA））之間的差異和相同之處；Suzuki等［2］進(jìn)一步詳細(xì)構(gòu)建了加氫站內(nèi)外部模型并進(jìn)行了定量風(fēng)險(xiǎn)分析，得出組件的個人風(fēng)險(xiǎn)值；Kwon等［3］用Hy-KoRAM和Phast/Safeti兩個風(fēng)險(xiǎn)分析軟件對加氫站進(jìn)行了定量風(fēng)險(xiǎn)分析，并對結(jié)果進(jìn)行了對比驗(yàn)證。

為了合理建設(shè)加氫站，并滿足《加氫站技術(shù)規(guī)范》（GB 50516—2010）［4］、《危險(xiǎn)化學(xué)品生產(chǎn)裝置和儲存設(shè)施風(fēng)險(xiǎn)基準(zhǔn)》（GB 36894—2018）［5］和《危險(xiǎn)化學(xué)品生產(chǎn)裝置和儲存設(shè)施外部安全防護(hù)距離確定方法》（GB/T 37243—2019）［6］等標(biāo)準(zhǔn)，需計(jì)算加氫站的安全距離與個人風(fēng)險(xiǎn)值并比較評估不同標(biāo)準(zhǔn)的差異性、適配性和標(biāo)準(zhǔn)性。

1 加氫站模型及風(fēng)險(xiǎn)量化

1.1 加氫站內(nèi)部模型

加氫站內(nèi)部布局模型如圖1所示，圖中展示出傳統(tǒng)加氫站內(nèi)部模型的布局并描繪了主要設(shè)施具體位置。本文只考慮加氫站中壓縮機(jī)、氣瓶和分配器這3組設(shè)施的氫氣泄漏情景，不考慮加氫站內(nèi)部設(shè)施的3D結(jié)構(gòu)對風(fēng)險(xiǎn)的影響。假設(shè)每個設(shè)施的代表性泄漏點(diǎn)為其中心點(diǎn)，例如壓縮機(jī)節(jié)點(diǎn)的代表點(diǎn)在其中心。以加氫站的左下角為仿真研究的坐標(biāo)原點(diǎn)（0，0），相應(yīng)的壓縮機(jī)、儲氫氣缸和分配器的坐標(biāo)點(diǎn)分別為（12.2，14.95）、（16.23，15.4）和（14.5，5）。

1.2 加氫站周邊模型

由于一般防護(hù)目標(biāo)環(huán)繞加氫站而建，具體分類見表1，防護(hù)目標(biāo)與加氫站之間的距離采用GB 50516—2010中第4部分的站址選擇中加氫站的氫氣工藝設(shè)施與站外建筑物、構(gòu)筑物的防火距離，見表2。

表3為規(guī)定的加氫站各設(shè)施與防護(hù)目標(biāo)安全距離與坐標(biāo)（以加氫站左下角為原點(diǎn)），圖2描述了加氫站周邊防護(hù)目標(biāo)到加氫站中各設(shè)施的間距。

1.3 加氫站內(nèi)氫氣流程及設(shè)施組件

從外部供應(yīng)的氫氣由氫氣拖車運(yùn)送到加氫站。用于提供20 MPa的氫氣，然后通過壓縮機(jī)將氫氣增加到82 MPa以后，將氫氣儲存在氣缸中。將氫氣從氣缸輸送到分配器并充入燃料汽車的儲氫罐（氫氣壓力為82 MPa），在約3 min內(nèi)將氫氣壓力從10 MPa提升到82 MPa。

在每個設(shè)施中有各種類型的組件，例如氣瓶、閥門和管道等。通過組件數(shù)量和管道長度估計(jì)分析事故場景頻率和風(fēng)險(xiǎn)。表4中提供了每個設(shè)施中的組件數(shù)量。

1.4 QRA方法

QRA可定量計(jì)算出事故的后果，加氫站的QRA操作步驟如圖3所示。

1）定義HRS模型并收集滿足QRA分析的數(shù)據(jù)。

2）危險(xiǎn)識別并選擇泄漏情景，例如導(dǎo)致噴射火或爆炸的設(shè)施氫氣泄漏情景。

3）風(fēng)險(xiǎn)分析，包括頻率分析、個人風(fēng)險(xiǎn)計(jì)算。后果分析，例如噴射火或爆炸超壓，超壓事故后果對站外更具殺傷力。

4）結(jié)合頻率和后果影響計(jì)算個人風(fēng)險(xiǎn)。假設(shè)人員長期處于防護(hù)目標(biāo)場所且無保護(hù)，發(fā)生危險(xiǎn)化學(xué)品事故而導(dǎo)致的死亡頻率，通常由設(shè)施周圍的風(fēng)險(xiǎn)等值線來體現(xiàn)。

5）將獲得的個人風(fēng)險(xiǎn)值與國家標(biāo)準(zhǔn)進(jìn)行比較。

2 氫能領(lǐng)域與危化品領(lǐng)域相關(guān)標(biāo)準(zhǔn)的風(fēng)險(xiǎn)管控要求

氫氣同時(shí)屬于氫能領(lǐng)域和?；奉I(lǐng)域標(biāo)準(zhǔn)涵蓋的范疇，因此有必要對加氫站基于定量風(fēng)險(xiǎn)分析對多領(lǐng)域標(biāo)準(zhǔn)交叉驗(yàn)證，以確定每個標(biāo)準(zhǔn)的特點(diǎn)和標(biāo)準(zhǔn)之間的適配性和準(zhǔn)確性。

危險(xiǎn)化學(xué)品領(lǐng)域標(biāo)準(zhǔn)：

GB 36894—2018規(guī)定了危險(xiǎn)化學(xué)品生產(chǎn)裝置和儲存設(shè)施個人風(fēng)險(xiǎn)和社會風(fēng)險(xiǎn)的可接受風(fēng)險(xiǎn)基準(zhǔn)值。

GB/T 37243—2019規(guī)定了危險(xiǎn)化學(xué)品生產(chǎn)裝置和儲存設(shè)施外部安全防護(hù)距離確定方法。

氫能領(lǐng)域標(biāo)準(zhǔn)：

GB 50516—2010適用于新建、改建、擴(kuò)建的加氫站工程的設(shè)計(jì)、施工和建造。

2.1 不同領(lǐng)域多標(biāo)準(zhǔn)交叉驗(yàn)證

把氫氣視為新能源：采用新能源領(lǐng)域標(biāo)準(zhǔn)數(shù)據(jù)，將QRA計(jì)算結(jié)果（安全距離和個人風(fēng)險(xiǎn)值等）或直接與?；窐?biāo)準(zhǔn)對比。例如：1）采用GB 50156—2010的安全距離并使用定量風(fēng)險(xiǎn)評價(jià)得出個人風(fēng)險(xiǎn)再與GB 36894—2018的一般防護(hù)目標(biāo)的個人風(fēng)險(xiǎn)基準(zhǔn)對比。2）采用GB 50516—2010的安全距離與GB/T 37243—2019的外部安全防護(hù)距離公式計(jì)算得出各類防護(hù)目標(biāo)的安全距離對比。

把氫氣視為?；罚翰捎梦；奉I(lǐng)域標(biāo)準(zhǔn)數(shù)據(jù)，將QRA計(jì)算結(jié)果（安全距離和個人風(fēng)險(xiǎn)值等）或直接與氫能標(biāo)準(zhǔn)對比。例如：采用GB 36894—2018的一般防護(hù)目標(biāo)的個人風(fēng)險(xiǎn)基準(zhǔn)通過定量風(fēng)險(xiǎn)評價(jià)得出安全距離再與GB 50516—2010的安全距離對比。

2.2 同領(lǐng)域多標(biāo)準(zhǔn)交叉驗(yàn)證

采用危化品領(lǐng)域標(biāo)準(zhǔn)數(shù)據(jù)，將QRA計(jì)算結(jié)果（安全距離和個人風(fēng)險(xiǎn)值等）或直接與其他?；窐?biāo)準(zhǔn)對比。例如：采用GB 36894—2018求得的安全距離，并使用GB/T 37243—2019的外部安全防護(hù)距離公式計(jì)算得出各類防護(hù)目標(biāo)的安全距離，將計(jì)算得出的安全距離進(jìn)行對比分析。

圖4為GB/T 37243—2019外部安全防護(hù)距離確定流程。因?yàn)槌鞘兄兄辉试S建造二級、三級加氫站，加氫站的站內(nèi)氫氣儲量不超過5 t，未達(dá)到危險(xiǎn)化學(xué)品重大危險(xiǎn)源的標(biāo)準(zhǔn)，故采用事故后果法判別外部安全防護(hù)距離。

事故后果法是根據(jù)最嚴(yán)重事故情景以及GB/T 37243—2019給出的空氣沖擊波超壓安全閾值，按式（1）計(jì)算出外部安全防護(hù)距離。

[Δp=14×QR+4.3×Q（2/3）R2+1.1×Q（1/3）R3]" （1）

式中：[Δp]——空氣沖擊波閾值超壓值，Pa；[Q]—— 一次爆炸的梯恩梯炸藥當(dāng)量，kg；[R]——爆炸點(diǎn)距離防護(hù)目標(biāo)的距離（視作加氫站最外軸線），m。

3 案例分析

案例中使用HyRAM軟件工具包來執(zhí)行QRA分析。HyRAM整合了公開可用的相關(guān)數(shù)據(jù)和模型來評估氫在使用、運(yùn)輸和儲存等基礎(chǔ)設(shè)施的安全性，并在風(fēng)險(xiǎn)評估計(jì)算中結(jié)合壓縮氣體和液化燃料的不同設(shè)施組件故障的概率［7］以及熱通量和超壓對人的影響的概率模型［8］。案例分析步驟為：

1）篩選出9種組件并統(tǒng)計(jì)各設(shè)備組件的數(shù)量，見表4。管道外徑、管道壁厚等其他數(shù)據(jù)見表5。

2）暴露人員位置和數(shù)量由用戶輸入定義。通過定義每個人員的每個維度（[x，y，z]）精確地分配人員位置或使用隨機(jī)采樣的概率分布（均勻或正態(tài)概率分布）在指定范圍定義用戶位置。為貼合實(shí)際情況，目標(biāo)防護(hù)物中暴露人員的分布類型應(yīng)為均勻分布。

3）輻射熱通量的危害和超壓危害通常用熱劑量單位（V）通過式（2）來表示熱通量強(qiáng)度和暴露時(shí)間的關(guān)系：

[V=I （4/3）×t]"" （2）

式中：[V]——熱劑量，（W·s）/m；[I]——輻射熱通量，W/m；[t]——暴露持續(xù)時(shí)間，s。

指定的熱概率模型為Tsao，計(jì)算熱劑量函數(shù)的死亡概率[Y1]：

[Y1=-36.38+2.56×lnV]"""""" （3）

式中：[Y1]——熱劑量函數(shù)的死亡概率。

超壓概率模型為lung（Eisenberg），計(jì)算暴露于超壓的死亡概率[Y2]：

[Y2=-77.1+6.91×lnPS]"""""" （4）

式中：[Y2]——暴露于超壓的死亡概率；[pS]——峰值超壓，Pa。

設(shè)計(jì)方案的所有參數(shù)見表5，其他未標(biāo)注的參數(shù)均采用HyRAM中自帶的默認(rèn)值。

3.1 基于不同領(lǐng)域相關(guān)標(biāo)準(zhǔn)的定量風(fēng)險(xiǎn)分析

如圖2所示，在加氫站四周布置GB 36894—2018危險(xiǎn)化學(xué)品生產(chǎn)裝置和儲存設(shè)施風(fēng)險(xiǎn)基準(zhǔn)中的一般防護(hù)物中各級別的防護(hù)目標(biāo)，防護(hù)目標(biāo)到加氫站中各設(shè)施的間距采用GB 50516—2010加氫站的氫氣工藝設(shè)施與站外建筑物、構(gòu)建物的防火距離。使用HyRAM軟件計(jì)算個人風(fēng)險(xiǎn)并與GB 36894—2018的個人風(fēng)險(xiǎn)基準(zhǔn)對比。

將表4和表5的數(shù)據(jù)輸入HyRAM計(jì)算得出個人風(fēng)險(xiǎn)值見表6。將總個人風(fēng)險(xiǎn)值與GB 36894—2018中防護(hù)目標(biāo)個人風(fēng)險(xiǎn)基準(zhǔn)中的個人風(fēng)險(xiǎn)進(jìn)行對比，結(jié)果見表7。各種防護(hù)目標(biāo)的個人風(fēng)險(xiǎn)總和均小于相應(yīng)的各類防護(hù)目標(biāo)的個人風(fēng)險(xiǎn)值，根據(jù)GB 36894—2018和GB 50516—2010的交叉驗(yàn)證可得出：采用GB 50516—2010加氫站技術(shù)規(guī)范規(guī)定的外部安全距離完全滿足GB 36894—2018危險(xiǎn)化學(xué)品生產(chǎn)裝置和儲存設(shè)施風(fēng)險(xiǎn)基準(zhǔn)。但直接依據(jù)GB 50516—2010規(guī)劃城市中加氫站的周邊建筑，無法合理規(guī)劃城市中加氫站周邊建筑的布局，根據(jù)GB 36894—2018中各類別防護(hù)目標(biāo)的個人風(fēng)險(xiǎn)值，使用HyRAM軟件計(jì)算出外部安全距離（采用保守的安全距離）。結(jié)合圖1繪制加氫站及其周邊個人風(fēng)險(xiǎn)圖，結(jié)果如圖5所示。最內(nèi)層區(qū)域代表個人風(fēng)險(xiǎn)大于3×10-5（次/a）是無防護(hù)目標(biāo)區(qū)，次內(nèi)層區(qū)域代表個人風(fēng)險(xiǎn)介于1×10-5（次/a）和3×10-5（次/a）之間，是三類防護(hù)目標(biāo)區(qū)，只可設(shè)三類防護(hù)目標(biāo)，次外層區(qū)域代表個人風(fēng)險(xiǎn)介于1×10-5（次/a）和3×10-6（次/a）之間，可設(shè)二類和三類防護(hù)目標(biāo)，最外層區(qū)域代表個人風(fēng)險(xiǎn)小于3×10-6（次/a），可設(shè)全部類別防護(hù)目標(biāo)。

如圖5所示，針對加氫站北方向的一類防護(hù)目標(biāo)，滿足其個人風(fēng)險(xiǎn)基準(zhǔn)為3×10-6次/a時(shí)，可將距離加氫站最外軸線的安全距離從24.95 m縮短到12 m。對于加氫站西方向的二類防護(hù)目標(biāo)，滿足其個人風(fēng)險(xiǎn)基準(zhǔn)為1×10-5次/a時(shí)，可將距離加氫站最外軸線的安全距離從17.8 m縮短到5 m。至于加氫站南方向的三類防護(hù)目標(biāo)，在滿足其個人風(fēng)險(xiǎn)基準(zhǔn)為3×10-5次/a時(shí)，可將距離加氫站最外軸線的安全距離從原來的15.05 m縮短到1 m。

通過圖5不同等級的風(fēng)險(xiǎn)基準(zhǔn)計(jì)算出加氫站周邊各類防護(hù)目標(biāo)的保守安全距離，結(jié)果如表8所示。

3.2 基于同領(lǐng)域相關(guān)標(biāo)準(zhǔn)的定量風(fēng)險(xiǎn)分析

將表8中的外部安全距離與表9中的外部安全距離進(jìn)行對比，由表10可知GB 36894—2018的安全距離在一類防護(hù)目標(biāo)的安全距離比GB/T 37243—2019的安全距離在一、二類防護(hù)目標(biāo)的安全距離長，GB 36894—2018更保守。而GB/T 37243—2019的安全距離在一類防護(hù)目標(biāo)的安全距離比GB36894—2018的安全距離在一、二類防護(hù)目標(biāo)的安全距離長，GB50516—20251的安全距離是最長的，是最保守的安全距離。

4 結(jié) 論

本文通過詳細(xì)比對和驗(yàn)證GB 36894—2018、GB 50516—2010和GB/T 37243—2019這3個國家標(biāo)準(zhǔn)對氫氣設(shè)施的安全距離和個人風(fēng)險(xiǎn)值要求，采用國際知名的HyRAM軟件定量對加氫站進(jìn)行QRA，對上述多領(lǐng)域標(biāo)準(zhǔn)的風(fēng)險(xiǎn)控制要求進(jìn)行了交叉分析驗(yàn)證，得到以下主要結(jié)論：

1）氫能領(lǐng)域的標(biāo)準(zhǔn)GB 50516—2010中規(guī)定的安全距離是3個標(biāo)準(zhǔn)中最保守的，采用該標(biāo)準(zhǔn)的安全距離將完全符合危險(xiǎn)化學(xué)品領(lǐng)域標(biāo)準(zhǔn)，但可能會由于過于保守而浪費(fèi)城市土地及建站與供氫成本激增。

2）針對傳統(tǒng)危險(xiǎn)化學(xué)品領(lǐng)域，通過GB 36894—2018和GB/T 37243—2019兩個?；奉I(lǐng)域的相關(guān)標(biāo)準(zhǔn)對比，GB/T 37243—2019規(guī)定的三類防護(hù)目標(biāo)的安全距離更保守，GB 36894—2018規(guī)定的一、二類防護(hù)目標(biāo)的安全距離更保守。加氫站可依據(jù)危險(xiǎn)化學(xué)品領(lǐng)域標(biāo)準(zhǔn)GB 36894—2018的各類防護(hù)目標(biāo)的個人風(fēng)險(xiǎn)值，通過HyRAM計(jì)算出各類防護(hù)目標(biāo)的安全距離。

3）危險(xiǎn)化學(xué)品標(biāo)準(zhǔn)適用范圍更廣，相關(guān)標(biāo)準(zhǔn)若用于加氫站的設(shè)計(jì)和建造，風(fēng)險(xiǎn)值會偏大。而氫能標(biāo)準(zhǔn)對氫安全更加細(xì)致保守，風(fēng)險(xiǎn)值更小，更加適用于氫能源化工企業(yè)。不同標(biāo)準(zhǔn)交互比對可讓優(yōu)勢最大化，使QRA結(jié)果更加可靠。為交叉分析方法在其他新能源加工領(lǐng)域標(biāo)準(zhǔn)的制定提供參考。

［參考文獻(xiàn)］

［1］"""" HAM K， MARANGON A， MIDDHA P， et al. Benchmark exercise on risk assessment methods applied to a virtual hydrogen" refuelling" station［J］." International" journal" of hydrogen energy， 2011， 36（3）： 2666-2677.

［2］"""" SUZUKI T， SHIOTA K， IZATO Y I， et al. Quantitative risk assessment using a Japanese hydrogen refueling station model［J］. International journal of hydrogen energy， 2021， 46（11）： 8329-8343.

［3］nbsp;""" KWON D， CHOI S K， YU C. Improved safety by crossanalyzing quantitative risk assessment of hydrogen refueling"" stations［J］." International" journal"" of" hydrogen energy， 2022， 47（19）： 10788-10798.

［4］"""" GB 50516—2010， 加氫站技術(shù)規(guī)范［S］.

GB 50516—2010， Technical code for bydrogen fuelling station［S］.

［5］"""" GB 36894—2018， 危險(xiǎn)化學(xué)品生產(chǎn)裝置和儲存設(shè)施風(fēng)險(xiǎn)基準(zhǔn)［S］.

GB 36894—2018， Risk criteria for hazardous chemicals production unit and storage installations［S］.

［6］"""" GB 37243—2019， 危險(xiǎn)化學(xué)品生產(chǎn)裝置和儲存設(shè)施外部安全防護(hù)距離確定方法［S］.

GB 37243—2019， Method for determining the external safety protection distance of hazardous chemical production and storage facilities［S］.

［7］"""" GROTH K M， HECHT E S. HyRAM： a methodology and toolkit for quantitative risk assessment of hydrogen systems［J］. International journal of hydrogen energy， 2017， 42（11）： 7485-7493.

［8］"""" EHRHART B， HECHT E. Hydrogen plus other alternative fuels risk assessment models（HyRAM+） version 4.1 technical reference manual［R］. Sandia national lab.（SNL-NM）， albuquerque， NM （United States）； sandia national laboratories， SNL california， 2022.

［9］"""" TORVI D A， HADJISOPHOCLEOUS G V， HUM J. A new method for estimating the effects of thermal radiation from fires on building occupants［C］//Proceedings of ASME 2000 International Mechanical Engineering Congress and Exposition， Orlando， Florida， USA， 2021： 65-72.

QUANTITATIVE RISK ASSESSMENT OF URBAN HYDROGEN FUELLING STATIONS BASED ON MULTI-STANDARD CROSSANALYSIS

Wang Haiqing1，Zuo Honge1，Zheng Wei2，Ma Jiawen1

（1. Department of Safety Science amp; Engineering， China University of Petroleum（East China）， Qingdao 266580， China；

2. Hangzhou Hande Quality Certification Service Co.， Ltd.， Shanghai Branch， Shanghai 200072， China）

Abstract：Addressing the variability of safety distance calculations in national standards for different fields quantitative risk analysis （QRA） was carried out on a typical hydrogen fuelling station configuration model， and the safety distance and personal risk between the hydrogen fuelling station facilities and the protection targets outside the station were calculated. The cross-validation of standards in many fields （hydrogen energy field and hazardous chemicals field） was realized， and the adaptability difference between hydrogen energy standards and hazardous chemicals standards in hydrogen fuelling stations was revealed， which provided a multi-dimensional reference for the design， construction and risk control of hydrogen fuelling stations.

Keywords：hydrogen； safety； hydrogenation； quantitative risk assessment； multi-standard crossanalysis