昌偉偉，吳　兵

(中國礦業(yè)大學(xué)(北京) 資源與安全工程學(xué)院， 北京 100083)

0　引言

經(jīng)過幾十年的發(fā)展，我國已經(jīng)形成了規(guī)模龐大、門類齊全和配置完善的化學(xué)工業(yè)體系，化工園區(qū)作為化學(xué)工業(yè)體系的重要組成部分，對我國化學(xué)工業(yè)的發(fā)展起到了舉足輕重的作用。據(jù)中國石油和化學(xué)工業(yè)聯(lián)合會化工園區(qū)工作委員會所做的一次全國性調(diào)研統(tǒng)計(jì)，截至2015年底，全國重點(diǎn)化工園區(qū)或以石油和化工為主導(dǎo)產(chǎn)業(yè)的工業(yè)園區(qū)共有502家，其中，國家級化工園區(qū)(包括經(jīng)濟(jì)技術(shù)開發(fā)區(qū)、高新區(qū))47家，省級化工園區(qū)262家，地市級化工園區(qū)193家[1]。由于化工園區(qū)內(nèi)存在種類繁多的化學(xué)工藝過程和工藝設(shè)備，一旦發(fā)生事故，后果不堪設(shè)想。2010年，大連孤山化工園區(qū)“7·16”輸油管道爆炸事故，若不是及時關(guān)閉輸油管道臨近的儲油罐，勢必引發(fā)整個園區(qū)危險(xiǎn)化學(xué)品事故的多米諾效應(yīng)。2016年，天津港“8·12”特大火災(zāi)爆炸事故，亦是由于危險(xiǎn)品倉庫首先爆炸而引發(fā)園區(qū)危險(xiǎn)化學(xué)品事故多米諾效應(yīng)，現(xiàn)場形成6處大火點(diǎn)及數(shù)10處小火點(diǎn)，事故最后造成165人遇難，8人失蹤，798人受傷，304幢建筑物、12 428輛商品汽車、7 533個集裝箱受損[2-3]，是近年來我國最大的化學(xué)品爆炸事故。如此突出的安全風(fēng)險(xiǎn)，對于化工園區(qū)的安全發(fā)展和政府的安全監(jiān)管都帶來了極大的挑戰(zhàn)。

化工園區(qū)內(nèi)重大事故具有多米諾效應(yīng)，當(dāng)園區(qū)內(nèi)某單元(企業(yè))發(fā)生的事故能量足夠大時，其危害便會涉及園區(qū)內(nèi)其他單元，導(dǎo)致多米諾效應(yīng)。從事故發(fā)生的能量角度，應(yīng)用爆炸學(xué)、傳熱學(xué)、網(wǎng)絡(luò)理論等知識，可以找出重大事故發(fā)生時能量傳播的途徑，結(jié)合園區(qū)內(nèi)各單元現(xiàn)有的風(fēng)險(xiǎn)應(yīng)對措施，切斷事故的多米諾效應(yīng)[4]。從經(jīng)濟(jì)角度，園區(qū)內(nèi)各單位必然從自身利益出發(fā)，選擇最經(jīng)濟(jì)的方式避免事故多米諾效應(yīng)對自身企業(yè)產(chǎn)生不良影響。無論從何角度出發(fā)，預(yù)防多米諾效應(yīng)的有效手段便為園區(qū)內(nèi)各單元間相互合作[5]。

本文應(yīng)用博弈論理論，從經(jīng)濟(jì)角度分析化工園區(qū)內(nèi)各單元為了預(yù)防重大事故的多米諾效應(yīng)，彼此間的博弈過程，最終求解出均衡點(diǎn)，以期為預(yù)防化工園區(qū)事故多米諾效應(yīng)提供理論參考。

園區(qū)內(nèi)各單元的安全投入主要包括用于提升企業(yè)安全管理水平和預(yù)防事故發(fā)生的人員、設(shè)備、材料、改善作業(yè)環(huán)境和管理水平等方面的資金投入，其中，政府強(qiáng)制要求配置的各種人材機(jī)等費(fèi)用為安全投入的下限，各單元為進(jìn)一步提升安全水平而進(jìn)行的投入，可以更有效地預(yù)防事故，同時可以減少自身發(fā)生事故與相鄰單元發(fā)生事故造成本單元的損失。在本文中認(rèn)為企業(yè)高于政府要求部分的安全投入可以有效的預(yù)防事故多米諾效應(yīng)。

博弈論又稱對策論，可以定義為是對智能的理性決策者之間沖突與合作的數(shù)學(xué)模型的研究，博弈論為分析那些涉及兩個或更多參與者且其決策會影響互相間的利益的局勢提供了1種數(shù)學(xué)方法[6，7]。最初應(yīng)用于經(jīng)濟(jì)學(xué)范疇，后廣泛應(yīng)用于軍事、工程等學(xué)科。在安全領(lǐng)域，Nadeau[8]應(yīng)用博弈論分析不同社會成員間對于職業(yè)健康和安全目標(biāo)的策略選擇，通過建立博弈均衡點(diǎn)，找出不同社會成員的最佳安全策略；Reniers[9]應(yīng)用博弈論討論多企業(yè)間協(xié)同安全管理的策略選擇問題，對于同一個安全目標(biāo)，不同企業(yè)有不同的策略選擇，通過博弈分析，確定適合不同企業(yè)的最佳策略；宗翔等[10]采用博弈論原理分析隧道與地下工程的風(fēng)險(xiǎn)，將大自然與施工方作為博弈雙方，通過技術(shù)參數(shù)分析提出應(yīng)對地下工程風(fēng)險(xiǎn)的措施；李路等[11]將博弈論應(yīng)用于電網(wǎng)事故的應(yīng)急處理中，分析電網(wǎng)事故與調(diào)度之間的博弈關(guān)系，確定最佳的調(diào)度措施。采用博弈論對化工園區(qū)事故多米諾效應(yīng)進(jìn)行研究，提高各單位間的合作意愿，使園區(qū)內(nèi)各單元選擇最佳安全投入策略，避免產(chǎn)生重大事故的發(fā)生，對化工園區(qū)的安全管理將會產(chǎn)生積極的影響。

1　博弈模型建立

1.1　問題的提出

在化工園區(qū)內(nèi)存在生產(chǎn)不同類型化工產(chǎn)品的企業(yè)(本文稱之為單元)，各單元是構(gòu)成化工園區(qū)的基本元素，各單元不僅面臨本單元內(nèi)部發(fā)生事故對自身產(chǎn)生的損失，同時，也會面臨園區(qū)內(nèi)鄰近單元發(fā)生事故對其產(chǎn)生的影響，這種由相鄰單元發(fā)生事故造成本單元事故損失的現(xiàn)象，稱之為多米諾效應(yīng)。由于多米諾效應(yīng)的影響，各單元不僅要面臨本單元發(fā)生事故造成損失(本文稱之為直接損失)的風(fēng)險(xiǎn)，也面臨相鄰單元發(fā)生事故造成本單位損失(本文稱之為間接損失)的風(fēng)險(xiǎn)，因此各單元的安全投入既包括為了提升自身安全管理和預(yù)防事故水平的安全投入，又包括為了預(yù)防相鄰單元事故多米諾效應(yīng)的安全投入。由于各單元利益出發(fā)點(diǎn)不同，為預(yù)防多米諾效應(yīng)進(jìn)行安全投入的動力也不同，則會出現(xiàn)個別單元愿意為了預(yù)防多米諾效應(yīng)進(jìn)行安全投入，而個別單元不愿意為預(yù)防多米諾效應(yīng)進(jìn)行安全投入的情況，彼此間的博弈由此產(chǎn)生[12-14]。

1.2　條件假設(shè)

假設(shè)化工園區(qū)由n個基本單元(化工企業(yè))構(gòu)成，每個單元均為理性的個體，對博弈參數(shù)進(jìn)行定義如下：

1)將是否愿意為預(yù)防多米諾效應(yīng)進(jìn)行安全投入作為各單元的基本策略Ai，進(jìn)行安全投入為I，不進(jìn)行安全投入為NI，則Ai={I，NI}。

2)將單元i不進(jìn)行安全投入發(fā)生事故造成單元i自身的損失定義為Li≥0，其發(fā)生概率為0≤Pii≤1。

3)將單元j不進(jìn)行安全投入發(fā)生事故造成相鄰單元i的損失定義為Lj≥0，j≠i，其發(fā)生概率為0≤Pji≤1。

4)將單元i的為了預(yù)防多米諾效應(yīng)的安全投入定義為ci≥0。

5)定義Ψ為當(dāng)單元i選擇任意策略Ai時，園區(qū)內(nèi)采取策略I的單元數(shù)量，Φ為園區(qū)內(nèi)采取策略NI的單元數(shù)量。

6)定義li({Ψ},Ai)為當(dāng)單元i選擇選擇任一策略時，單元i的間接損失。

7)定義ui(i=1,…,n)為單元i無論選擇何種策略，其安全投入與事故損失之和。

基于以上的假設(shè)與參數(shù)定義，可以得出當(dāng)單元i選擇策略Ai=I時，其面臨的支付值為：

ui=ci+li({Ψ},I)

(1)

當(dāng)單元i選擇策略Ai=NI時，其面臨的支付值為：

ui=LiPii∏(1-Pji)+li({Ψ},NI)(1-Pii)(j≠i，j∈Φ)

(2)

公式(2)的前半部分表示單元i的直接損失，后半部分表示單元i的間接損失。通過各參與方不同策略選擇模型的建立，可以進(jìn)一步分析博弈模型的納什平衡條件和化工園區(qū)內(nèi)各單元間的合作是否處于博弈穩(wěn)定的情形。

1.3　支付矩陣建立

為了進(jìn)一步分析，假設(shè)1.2節(jié)建立的博弈模型參與方僅為2個單元，即單元1和單元2，則可以建立2個單元的支付矩陣，見表1。

表1　化工園區(qū)內(nèi)2單元的支付矩陣

假設(shè)企業(yè)為預(yù)防多米諾效應(yīng)進(jìn)行的安全投入對于預(yù)防化工園區(qū)重大事故多米諾效應(yīng)的是可靠并有效的，根據(jù)支付矩陣，會出現(xiàn)4種情況：

4)當(dāng)各單元均選擇策略NI時，各單元的支付值為自身發(fā)生重大事故的損失與對方發(fā)生重大事故對其的影響值之和，即單元1的支付值為P11L1(1-P21)+P21L1(1-P11)，單元2的支付值為P21L2(1-P12)+P12L2(1-P22)。

2　均衡分析

理性個體如何選擇策略與博弈模型內(nèi)的參數(shù)有關(guān)，不同的參數(shù)關(guān)系會產(chǎn)生不同的納什均衡點(diǎn)。針對表1建立的支付矩陣進(jìn)行分析，可以得出不同條件下的納什均衡點(diǎn)，見表2。

表2　不同參數(shù)條件下的納什均衡點(diǎn)

2.1　情形1

2.2　情形2

2.3　情形3

2.4　情形4

當(dāng)2個博弈單元中的1個單元選擇進(jìn)行合作策略(I，I)的支付值大于選擇不進(jìn)行合作的支付值時，必然會選擇代價最小的策略組合，(NI，NI)將會是穩(wěn)定的納什均衡點(diǎn)。

這種情形需對各單元進(jìn)行認(rèn)真仔細(xì)的分析，分析當(dāng)采取合作策略時安全投入的支付值與安全生產(chǎn)達(dá)到的水平之間的關(guān)系，根據(jù)具體情況選擇是否采取措施促使各單元選擇策略(I，I)。

3　實(shí)例分析

選擇某LPG園區(qū)的數(shù)據(jù)進(jìn)行實(shí)例分析，園區(qū)內(nèi)存在單元1和單元2，單元1包括12個LPG儲罐；單元2包括7個LPG儲罐。相鄰2個單元同時面臨著由于自身發(fā)生事故造成的損失和相鄰單元事故多米諾效應(yīng)而造成的損失。

3.1　參數(shù)選擇

Pii為由于單元內(nèi)部發(fā)生事故造成本單元損失的發(fā)生概率，對于本案例取P11為1.10×10-4，P22為0.64×10-4。

Pij為由于單元內(nèi)部發(fā)生事故造成相鄰單元損失的發(fā)生概率，對于本案例取P12為5.50×10-4，P21為4.34×10-4。

事故損失Li是通過發(fā)生事故后用于重建或修復(fù)生產(chǎn)單元發(fā)生的費(fèi)用來確定。對于單元1潛在的由于多米諾效應(yīng)造成的損失值L1=68 777 800元，單元2潛在的由于多米諾效應(yīng)造成的損失值L2=14 068 186元。

3.2　博弈分析

將數(shù)據(jù)代入表1，可以得到該LPG園區(qū)的支付矩陣，見表3。

表3　LPG園區(qū)內(nèi)兩單元支付矩陣

4)當(dāng)處于情形4時，由于各單元選擇合作策略(I，I)時的支付值要大于至少有1個單元選擇不合作的支付值，那么對于管理者，盡量選擇對1個單元或者全部單元采取激勵措施，使雙方均選擇策略(I，I)。

4　結(jié)論

1)建立了化工園區(qū)內(nèi)各單元為了預(yù)防事故多米諾效應(yīng)的博弈模型。博弈模型考慮了各單元自身發(fā)生事故的概率與損失，相鄰單元發(fā)生事故的概率與損失，各單元為了預(yù)防事故多米諾效應(yīng)而進(jìn)行的安全投入等因素。

2)應(yīng)用博弈模型建立相鄰2單元的博弈支付矩陣。討論2單元選擇不同策略時，各單元的支付值函數(shù)表明，各單元選擇的安全投入策略不同，為此付出的代價(支付值)不同。

3)討論不同的參數(shù)條件下博弈模型的納什均衡點(diǎn)。將博弈模型分為4種情形并確定每種情形下的納什均衡點(diǎn)，將每種情形對應(yīng)的納什均衡點(diǎn)與園區(qū)管理者希望的均衡點(diǎn)進(jìn)行比較，為園區(qū)管理者是否采取激勵措施提供依據(jù)。

[1]段雅楠. 保障園區(qū)安全 共享發(fā)展紅利——訪中國石油和化學(xué)工業(yè)聯(lián)合會化工園區(qū)工作委員會秘書長楊挺[J]. 現(xiàn)代職業(yè)安全, 2016(9):14-17.

DUAN Ya’nan. Safety of the park Sharing development dividends[J].Modern Safey, 2016(9):14-17.

[2]國家安全生產(chǎn)監(jiān)督管理總局.天津港“8·12” 瑞海公司危險(xiǎn)品倉庫特別重大火災(zāi)爆炸事故調(diào)查報(bào)告[EB/OL].2016.http://www.chinasafety.gov.cn/newpage/newfiles/201600812baogao.pdf.

[3]Fu G, Wang J, Yan M. Anatomy of Tianjin Port fire and explosion: Process and causes[J]. Process Safety Progress, 2016, 35(3):216-220.

[4]孫康,周武,柴瑞瑞. 化工園區(qū)危險(xiǎn)化學(xué)品事故多米諾效應(yīng)的網(wǎng)絡(luò)演化模型[J]. 電子科技大學(xué)學(xué)報(bào)(社科版),2015,17(5):37-41.

SUN Kang，ZHOU Wu，CHAI Ruirui. Network evolution model of hazardous chemical accident with domino effect in chemical industry park[J].Journal of UESTC(Social Sciences Edition),2015,17(5):37-41.

[5]李樹謙, 胡兆吉. 化工園區(qū)多米諾事故風(fēng)險(xiǎn)分析方法研究[J]. 安全生產(chǎn)與監(jiān)督, 2008(3):56-58.

LI Shuqian, HU Zhaoji. Study on the risk analysis method of domino accident in chemical park[J]. Work Safety & Supervision, 2008(3):56-58.

[6]羅杰·B·邁爾森.博弈論:矛盾沖突分析[M].北京:中國人民大學(xué)出版社,2015.

[7]謝識予.經(jīng)濟(jì)博弈論[M].上海:復(fù)旦大學(xué)出版社,2002.

[8]Nadeau S. Cooperation in health and safety:a game theory analysis[J].Law Rev (Risks),2003，1(3/4): 219-227.

[9]Reniers G. An external domino effects investment approach to improve cross-plant safety within chemical clusters[J]. Journal of Hazardous Materials, 2010, 177(1-3):167.

[10]宗翔,張冬梅,黃宏偉.博弈論在隧道及地下工程風(fēng)險(xiǎn)應(yīng)對中的應(yīng)用[J].地下空間與工程學(xué)報(bào)，2012，8(S2):1764-1765.

ZONG Xiang, ZHANG Dongmei, HUANG Hongwei. Application of game theory to risk response strategies in tunnel and underground engineering[J]. Chinese Journal of Underground Space and Engineering，2012,8(S2):1764-1765.

[11]李路, 周渝慧. 博弈論在處理大電網(wǎng)事故中的應(yīng)用[J]. 東北電力技術(shù), 2005, 26(4):25-27.

LI Lu，ZHOU Yuhui. The application of game theory to the handling of power grid collapse[J].Northeast Electric Power Technology， 2005, 26(4):25-27.

[12]Reniers G L L, Ale B J M, Dullaert W, et al. Designing continuous safety improvement within chemical industrial areas[J]. Safety Science, 2009, 47(5):578-590.

[13]Reniers G, Pavlova Y. Using Game Theory toimprove safety within chemical industrial parks[M]. London: Springer, 2013.

[14]Reniers G, Cuypers S, Pavlova Y. A game-theory based Multi-plant Collaboration Model (MCM) for cross-plant prevention in a chemical cluster.[J]. Journal of Hazardous Materials, 2012, 209-210(1):164-176.