周浩霖，朱常龍，蔣軍成

(1.南京工業(yè)大學(xué) 安全科學(xué)與工程學(xué)院，江蘇 南京 211816；2.江蘇省危險化學(xué)品本質(zhì)安全與控制技術(shù)重點實驗室，江蘇 南京 211816)

0 引言

保護層分析(Layer of Protection Analysis，LOPA)是1種簡化的半定量分析理論，其應(yīng)用方法簡潔有效，在化工過程行業(yè)中已被廣泛運用，在事故預(yù)防中發(fā)揮著重要作用。相關(guān)學(xué)者圍繞LOPA開展了各類改進策略的研究并在定量化方面獲得了較大成功[1-2]。然而，罐區(qū)等重大危險源發(fā)生事故時往往影響范圍廣，常常會引發(fā)次生事故甚至導(dǎo)致災(zāi)難性連鎖事故，而LOPA更適合評估既定事故場景下的獨立安全屏障性能[3]，而非對所有安全屏障阻止事故傳播擴展的性能作定量結(jié)論，在利用LOPA對量化安全屏障在防止事故傳播的作用研究上仍有不足。

化工事故的傳播過程，通常被形象地稱為“多米諾效應(yīng)”。相較于單一事故場景，多米諾事故發(fā)生頻率低，但影響范圍更大，后果更嚴(yán)重[4]。國內(nèi)外曾發(fā)生了多起導(dǎo)致人員大量傷亡，嚴(yán)重經(jīng)濟損失的多米諾事故，造成了極其惡劣的社會影響：2009年齋普爾印度石油公司油庫的1起轉(zhuǎn)輸作業(yè)中汽油泄漏引發(fā)的連鎖火災(zāi)爆炸事故造成11人死亡，45人受傷，當(dāng)?shù)鼐用翊笠?guī)模撤離[5]；2019年響水天嘉宜公司硝化廢料火災(zāi)引發(fā)了連鎖爆炸事故，造成了78人死亡及巨量財產(chǎn)損失[6]。這些事故引發(fā)了人們對重大危險源多米諾效應(yīng)的廣泛關(guān)注，大量學(xué)者開展了多米諾效應(yīng)中安全屏障定量評估及優(yōu)化管理方面的研究[7-9]并取得了較大進展。

本文借鑒LOPA分析思路及表現(xiàn)形式，將LOPA中獨立保護層引申為多米諾擴展過程中安全屏障，通過引入可用性及有效性概念量化安全屏障在阻止多米諾擴展時的性能，實現(xiàn)LOPA及多米諾效應(yīng)定量風(fēng)險評估的有機結(jié)合，擴展LOPA邏輯適用范圍并使多米諾效應(yīng)定量風(fēng)險評估更加科學(xué)合理。

1 基于LOPA邏輯的多米諾效應(yīng)定量評估流程及分析思路

多米諾效應(yīng)定量風(fēng)險評估在傳統(tǒng)框架上[10]將獨立保護層引申為阻止多米諾擴展的安全屏障，考慮安全屏障在罐區(qū)多米諾擴展過程中的作用并將其以LOPA分析過程呈現(xiàn)。如流程圖1所示，步驟如下：

圖1 多米諾效應(yīng)定量風(fēng)險評估流程圖Fig.1 Flow chart of quantitative risk assessment of domino effect

1)確定初始事故場景的擴展向量并計算其擴展向量強度。

2)計算結(jié)果與表1[11]比較識別多米諾事故易損容器。

表1 池火災(zāi)擴展閾值表Table 1 Escalation thresholds of pool fire

3)計算未考慮安全屏障時多米諾擴展概率Pd1，識別初始事故與多米諾事故組合并計算未考慮安全屏障時多米諾事故發(fā)生頻率fd1。多米諾事故發(fā)生頻率fd即初始事故場景發(fā)生頻率fp與多米諾擴展概率Pd之積。相關(guān)研究表明Probit模型對一級多米諾擴展過程計算較為準(zhǔn)確[12]，Probit模型及其正態(tài)分布式見表2[13]及式(1)。

(1)

式中：P為概率單位值對應(yīng)概率；u為積分變量。

表2中Y為概率單位值；ttf為無安全屏障下容器的失效時間，min；I為容器所受熱輻射強度，kW/m2；概率單位值Y服從正態(tài)分布，將Y代入式(1)求得多米諾擴展概率。

表2 Probit模型Table 2 Probit model

4)由于擴展向量類型及作用時間不同，主動保護系統(tǒng)僅對擴展向量長時間作用的多米諾事故有效[14]，需識別對應(yīng)場景下阻止多米諾擴展的安全屏障。

5)安全屏障量化評估，步驟見后文。

6)列出考慮安全屏障時多米諾擴展事件樹。

7)計算考慮安全屏障時多米諾擴展概率Pd2及事故發(fā)生頻率fd2。

8)人體脆弱性模型計算區(qū)域內(nèi)個人風(fēng)險值IR并繪制個人風(fēng)險等值線圖。表3為人體脆弱性Probit模型[10]，求得個人死亡概率Pf后使用式(2)求得個人風(fēng)險值。

表3 人體脆弱性概率模型Table 3 Probabilistic model of human vulnerability

表3中：Y為概率單位值；D為劑量；I為人體所受熱輻射強度，W/m2；te為暴露時間，min。個人風(fēng)險值計算式為：

IRi=fi×Pf,i

(2)

式中：IRi為第i個事故造成的個人風(fēng)險；fi為第i個儲罐發(fā)生事故的頻率，次/a；Pf,i為第i個儲罐事故個人死亡概率。

9)若區(qū)域內(nèi)出現(xiàn)高于容許值的等值線，需另提出安全屏障將個人風(fēng)險降低至容許值以下。

2 安全屏障對多米諾擴展概率的影響

2.1 可用性及有效性

Landucci等[7]引入安全屏障的有效性及可用性?？捎眯员欢x為安全屏障的需求時的失效概率(probability of failure on demand,PFD)；有效性被定義為啟用安全屏障時阻止擴展的概率，以η表示。

2.2 安全屏障定量評估

2.2.1 可用性及有效性的量化

表4為各安全屏障對應(yīng)的可用性與有效性參考值及定義[7]，表中：ttfp為考慮防火層時容器失效時間，min；ttfc為防火層延長容器失效時間，min；tfm為控制事故所需時間，min；IWDS為考慮噴淋冷卻水時容器表面所受熱輻射強度，kW/m2；θ為熱輻射強度減緩系數(shù)；IHL為未考慮噴淋冷卻水時容器表面所受熱輻射強度，kW/m2。

表4 安全屏障對應(yīng)的可用性，有效性參考值Table 4 Reference values of availability and effectiveness of safety barriers

2.2.2 可用性及有效性計算邏輯

根據(jù)安全屏障的類型不同引入3種邏輯門用于定量評估安全屏障性能，見表5。

表5 安全屏障的邏輯門Table 5 Logic gates of safety barriers

3 案例分析

3.1 罐區(qū)基本情況

苯乙烯屬乙A類可燃液體，閃點31 ℃，沸點146 ℃，相對密度0.91，分子量104.1，飽和蒸汽壓0.7 kPa(20 ℃)。設(shè)2×2 000 m3苯乙烯罐組如圖2，參數(shù)見表6。

表6 案例分析：苯乙烯儲罐參數(shù)設(shè)置Table 6 Case study:parameters settings of styrene storage tanks

圖2 案例分析：苯乙烯儲罐布局Fig.2 Case study:layout of styrene tank farm

3.2 初始事故場景及頻率

苯乙烯罐常因儲存不當(dāng)發(fā)生物料聚合導(dǎo)致容器爆裂，泄漏大量物料形成可燃液池[15]。據(jù)統(tǒng)計分析，80%的初始火災(zāi)事故類型為池火災(zāi)[16]，苯乙烯罐泄漏的發(fā)生頻率假設(shè)為5×10-5次/a，參考如圖3的可燃液體泄漏事件樹[17]，可得出初始事故場景(池火災(zāi))發(fā)生頻率fp為3.47×10-6次/a(池火+閃火階段晚期池火)。

圖3 案例分析：可燃液體泄漏事件樹Fig.3 Case study:event tree of flammable liquid leakage

3.3 未考慮安全屏障時的多米諾擴展概率

據(jù)第1章所述步驟1)～3)，確定初始事故場景：假設(shè)初始事故場景為泄漏全部物料導(dǎo)致池火災(zāi)，池火焰為圓柱形，使用Mudan模型計算TK1罐池火災(zāi)對TK2罐表面受到的熱輻射，得出TK2儲罐接受熱輻射值為51.06 kW/m2，大于多米諾效應(yīng)擴展閾值(見表1)，則確定TK2儲罐為唯一多米諾易損單元。儲罐所受擴展向量強度、失效時間及未考慮安全屏障時多米諾事故發(fā)生頻率，見表7。

表7 案例分析：未考慮安全屏障的多米諾效應(yīng)評估結(jié)果Table 7 Case study:results of domino effect assessment without considering safety barriers

3.4 考慮安全屏障時的多米諾擴展概率

據(jù)第1章步驟4)～7)步識別安全屏障并計算多米諾擴展概率及事故發(fā)生概率。TK1初始事故場景為池火災(zāi)，主動保護系統(tǒng)有足夠時間啟動，可考慮主動保護系統(tǒng)。常用的安全屏障有泡沫滅火系統(tǒng)、噴淋冷卻水系統(tǒng)、防火層。對于泄漏引發(fā)的池火災(zāi)，泡沫滅火裝置作為安全屏障無法有效降低初始事故擴展，因此識別噴淋冷卻系統(tǒng)及防火層為安全屏障(邏輯門a)，1.89×10-2與4.33×10-2，有效性定義、數(shù)值及邏輯計算模式見表4及表5。案例中ttfc取15 min，θ取70%。圖4為考慮安全屏障時的池火災(zāi)多米諾事件樹；表8為考慮安全屏障前后的池火災(zāi)多米諾擴展概率及多米諾事故發(fā)生頻率。

表8 案例分析：考慮安全屏障后多米諾效應(yīng)評估結(jié)果Table 8 Case study:results of domino effect assessment considering safety barriers

圖4 案例分析：考慮安全屏障下的池火災(zāi)事件樹Fig.4 Case study:event tree of pool fire considering safety barriers

對比表7與表9，考慮安全屏障后的多米諾事故發(fā)生頻率為減輕與未減輕多米諾事故發(fā)生頻率之和(3.236×10-7次/a)，較未考慮安全屏障的多米諾事故發(fā)生頻率(3.36×10-6次/a)少了1個數(shù)量級，即安全屏障可以有效減少多米諾發(fā)生概率。

表9 案例分析：考慮安全屏障后事故發(fā)生頻率匯總Table 9 Case study:summary of accident frequencies considering safety barriers

3.5 罐組個人風(fēng)險

根據(jù)第1章步驟8)同樣使用Mudan池火災(zāi)模型計算熱輻射強度，據(jù)表3及式(2)可計算個人風(fēng)險值IR。

為方便討論，設(shè)IRprimary為初始事故個人風(fēng)險值，IRdomino為多米諾事故個人風(fēng)險值。根據(jù)表10個人風(fēng)險疊加情況，得到對應(yīng)個人風(fēng)險等值線圖5。

表10 案例分析：3種假設(shè)情況下的多米諾效應(yīng)個人風(fēng)險計算表達式Table 10 Case study:expressions for calculating individual risk of domino effect under three presumptions

通過圖5(a)～(b)對比可知：1)同一事故中，多米諾效應(yīng)發(fā)生時個人風(fēng)險等值線范圍較廣。2)多米諾效應(yīng)下，有安全屏障時代表1×10-5的風(fēng)險等值線從區(qū)域中消失。3)有安全屏障時，多米諾效應(yīng)個人風(fēng)險等值線范圍變小，幾乎與無多米諾效應(yīng)時一致。以上表明安全屏障顯著地降低了多米諾事故個人風(fēng)險。

圖5 池火災(zāi)個人風(fēng)險等值線Fig.5 Isogram of pool fire individual risk

3.6 基于LOPA邏輯的多米諾效應(yīng)定量風(fēng)險評估表

根據(jù)第1章步驟9)檢驗區(qū)域內(nèi)個人風(fēng)險是否超過不可接受值；若大于不可接受值，則更換有效性及可用性更高的安全屏障或增加安全屏障數(shù)量?；谝陨纤悸?，遵循LOPA邏輯，以表格形式直觀展示考慮安全屏障的多米諾效應(yīng)定量風(fēng)險評估結(jié)果，如表11所示，其可作為企業(yè)開展多米諾效應(yīng)定量風(fēng)險評估、檢驗所識別的安全屏障阻止多米諾效應(yīng)擴展性能的參考。

表11 基于LOPA邏輯的多米諾效應(yīng)定量風(fēng)險評估Table 11 Quantitative risk assessment of domino effect based on LOPA logic

4 結(jié)論

1)將保護層等安全屏障按功能劃分為降低罐體失效可能性、降低初始事故的影響、外部事故減緩措施3個類型，對不同類型的安全屏障指定相應(yīng)的可用性、有效性屬性及邏輯關(guān)系，采用事件樹分析計算罐區(qū)多米諾擴展概率，改進對罐區(qū)多米諾效應(yīng)定量風(fēng)險評估。提出基于LOPA邏輯的多米諾效應(yīng)定量風(fēng)險評估表，拓展LOPA方法的應(yīng)用并形成較為直觀的分析結(jié)果，為化工企業(yè)開展多米諾效應(yīng)定量風(fēng)險評估、檢驗安全屏障性能提供參考。

2)基于改進的罐區(qū)多米諾效應(yīng)定量風(fēng)險評估方法對某苯乙烯罐區(qū)開展案例分析，得出：多米諾效應(yīng)增加了區(qū)域內(nèi)個人風(fēng)險水平。相較完全未考慮安全屏障的情境，識別防火層，噴淋冷卻水系統(tǒng)作為安全屏障下的多米諾事故頻率及個人風(fēng)險降低了1個數(shù)量級，添加安全屏障可顯著降低罐區(qū)多米諾事故發(fā)生頻率及個人風(fēng)險。

3)主要針對罐區(qū)多米諾場景，所識別的安全屏障種類有限，后續(xù)需要進一步研究生產(chǎn)、運輸環(huán)節(jié)等場景下適用的安全屏障量化屬性。此外，僅考慮了1級多米諾擴展，而多級擴展過程對安全屏障可用性和有效性的影響較為復(fù)雜，還有待進一步的研究。