陳國華，王永興，高子文

(1. 華南理工大學(xué) 安全科學(xué)與工程研究所，廣東 廣州 510640；2. 惠州大亞灣應(yīng)急管理有限公司，廣東 惠州 516081)

0　引言

化工園區(qū)主要是涉及化工行業(yè)的企業(yè)聚集區(qū)，其危險(xiǎn)物質(zhì)數(shù)量和潛在能量巨大[1]，企業(yè)內(nèi)部涉及的化學(xué)工藝復(fù)雜、設(shè)備裝置眾多、操作條件苛刻，稍有疏忽，極易造成嚴(yán)重的事故[2]。能否準(zhǔn)確判斷事故風(fēng)險(xiǎn)的發(fā)展演化過程，直接影響到應(yīng)急決策方案制定的科學(xué)性、有效性及應(yīng)急救援工作的順利開展[3]?；@區(qū)事故風(fēng)險(xiǎn)具有動(dòng)態(tài)性和不確定性，如何判斷事故風(fēng)險(xiǎn)演化狀態(tài)并及時(shí)采取有效措施，是園區(qū)發(fā)展中面臨的嚴(yán)峻現(xiàn)實(shí)。

熵(Entropy)是個(gè)源自熱力學(xué)領(lǐng)域的物理概念，是表示系統(tǒng)無序程度和混亂程度的狀態(tài)量[4]。風(fēng)險(xiǎn)熵是對系統(tǒng)風(fēng)險(xiǎn)程度不確定性的1種度量[5]，常用于系統(tǒng)風(fēng)險(xiǎn)評估、突發(fā)事件預(yù)測[6]，風(fēng)險(xiǎn)熵理論能夠從系統(tǒng)風(fēng)險(xiǎn)動(dòng)態(tài)變化的角度分析化工園區(qū)事故風(fēng)險(xiǎn)狀態(tài)和演化方向。國外學(xué)者Johnson等[7]研究了熵增與突發(fā)事件的關(guān)系，指出突發(fā)事件發(fā)生是系統(tǒng)熵增和紊亂的后果；Spartalis等[8]建立模糊熵的風(fēng)險(xiǎn)評估模型用于自然災(zāi)害的風(fēng)險(xiǎn)等級評估。國內(nèi)學(xué)者對熵理論的研究多與脆性理論、耗散結(jié)構(gòu)等結(jié)合，用于風(fēng)險(xiǎn)狀態(tài)的表征，陳偉珂等[9]通過分析熵增因素和負(fù)熵因素來說明突發(fā)事件的產(chǎn)生機(jī)理；汪送等[10]分析風(fēng)險(xiǎn)熵在事故演化網(wǎng)絡(luò)中的動(dòng)力學(xué)傳播過程來為事故分析和預(yù)防提供思路。

近年來國內(nèi)外化工事故不斷發(fā)生[11]，造成眾多嚴(yán)重后果。國內(nèi)外學(xué)者從事故多米諾效應(yīng)[12]、后果模擬[13]、事故發(fā)生機(jī)理等角度分析了化工事故的發(fā)展演化特點(diǎn)。化工園區(qū)災(zāi)害事故發(fā)生過程中系統(tǒng)可能出現(xiàn)崩潰突變的情況，風(fēng)險(xiǎn)熵理論能夠在初始事故發(fā)生后對系統(tǒng)風(fēng)險(xiǎn)進(jìn)行研判，以判斷是否可能會(huì)導(dǎo)致二次事故發(fā)生，當(dāng)前依據(jù)風(fēng)險(xiǎn)熵理論分析化工園區(qū)事故風(fēng)險(xiǎn)發(fā)展演化的研究則較少涉及。

基于以上分析，依據(jù)尖點(diǎn)突變模型和風(fēng)險(xiǎn)熵理論，結(jié)合化工園區(qū)實(shí)際情景，構(gòu)建化工園區(qū)系統(tǒng)風(fēng)險(xiǎn)熵突變模型，并對某?；穬?chǔ)存企業(yè)丁二烯泄漏情況下園區(qū)系統(tǒng)風(fēng)險(xiǎn)狀態(tài)進(jìn)行判斷，以期為此類事故引發(fā)的園區(qū)二次事故風(fēng)險(xiǎn)突變提供判別依據(jù)和參考。

1　化工園區(qū)事故風(fēng)險(xiǎn)躍遷機(jī)理分析

化工園區(qū)是由物質(zhì)、能量、信息等匯聚而成的復(fù)雜非線性系統(tǒng)，園區(qū)內(nèi)具有眾多危險(xiǎn)源且與外界環(huán)境有相互作用，如圖1所示。

圖1　化工園區(qū)物質(zhì)、能量、信息流動(dòng)示意Fig.1　Material, energy and information flow charts in chemical industry parks

系統(tǒng)風(fēng)險(xiǎn)熵與化工園區(qū)的風(fēng)險(xiǎn)狀態(tài)相對應(yīng)，它有2部分組成，一部分是園區(qū)系統(tǒng)內(nèi)部產(chǎn)生的風(fēng)險(xiǎn)熵增u，熵增會(huì)增加系統(tǒng)的紊亂程度；另一部分則是外界輸入改善風(fēng)險(xiǎn)的風(fēng)險(xiǎn)熵減v，熵減有利于降低系統(tǒng)風(fēng)險(xiǎn)。

用風(fēng)險(xiǎn)狀態(tài)參量λ表示化工園區(qū)系統(tǒng)風(fēng)險(xiǎn)偏離平衡態(tài)的程度，λ=0對應(yīng)于平衡態(tài)，風(fēng)險(xiǎn)狀態(tài)P是λ的函數(shù)，λc為臨界狀態(tài)，當(dāng)λ<λc時(shí)，園區(qū)處于系統(tǒng)風(fēng)險(xiǎn)演化初期，風(fēng)險(xiǎn)熵增較少，系統(tǒng)保持穩(wěn)定狀態(tài)。當(dāng)λ>λc，系統(tǒng)越過臨界狀態(tài)進(jìn)入非線性區(qū)，在系統(tǒng)熵增的作用下園區(qū)風(fēng)險(xiǎn)狀態(tài)發(fā)生隨機(jī)“微漲落”，在風(fēng)險(xiǎn)狀態(tài)臨界點(diǎn)(點(diǎn)c，d，e)，在突發(fā)因素或初始事故的誘發(fā)下園區(qū)風(fēng)險(xiǎn)狀態(tài)發(fā)生突變[14]，產(chǎn)生“巨漲落”，如圖2所示，會(huì)引起更為嚴(yán)重的后果。

圖2　化工園區(qū)風(fēng)險(xiǎn)狀態(tài)躍遷Fig.2　Chemical industry park risk state transition diagram

2　系統(tǒng)風(fēng)險(xiǎn)熵突變模型的原理及構(gòu)建

2.1　突變模型的原理

突變理論由法國數(shù)學(xué)家托姆(Renethom)于1972年創(chuàng)立，它通過勢函數(shù)來描述和預(yù)測事物發(fā)展過程中的突變。突變模型在社會(huì)學(xué)、經(jīng)濟(jì)學(xué)、水資源評價(jià)、能源等領(lǐng)域得到廣泛應(yīng)用，常見的有尖點(diǎn)突變模型、燕尾突變模型、蝴蝶突變模型等[15]。突變模型在分析系統(tǒng)的不連續(xù)變化過程中不需要涉及到系統(tǒng)的內(nèi)在機(jī)制，其特點(diǎn)是通過勢函數(shù)來表征系統(tǒng)可能出現(xiàn)的狀態(tài)，并通過研究參數(shù)變化時(shí)函數(shù)在臨界點(diǎn)附近的不連續(xù)特征來判斷系統(tǒng)是否發(fā)生突變。

表1　常用的突變模型

表1中，突變函數(shù)F(x1)表示1個(gè)系統(tǒng)中狀態(tài)變量x1的勢函數(shù)，狀態(tài)變量系數(shù)u1，u2，u3，u4表示該狀態(tài)變量的控制量。

2.2　化工園區(qū)系統(tǒng)風(fēng)險(xiǎn)熵突變模型構(gòu)建

化工園區(qū)是1個(gè)復(fù)雜的非線性生產(chǎn)系統(tǒng)，將導(dǎo)致系統(tǒng)紊亂程度增加的因素定義為熵增因子，降低系統(tǒng)風(fēng)險(xiǎn)的因素定義為熵減因子，熵增因子和熵減因子相互作用的過程中，風(fēng)險(xiǎn)狀態(tài)也會(huì)發(fā)生動(dòng)態(tài)變化，例如火災(zāi)爆炸中毒等事故是園區(qū)系統(tǒng)風(fēng)險(xiǎn)狀態(tài)突變的結(jié)果，這種狀態(tài)的突變不是瞬間完成的，而是系統(tǒng)結(jié)構(gòu)、功能朝某個(gè)方向演化的結(jié)果，當(dāng)園區(qū)有事故征兆或發(fā)生初始事故時(shí)，通過突變模型對系統(tǒng)風(fēng)險(xiǎn)狀態(tài)進(jìn)行判斷，以評估二次事故發(fā)生的可能性。

2.2.1突變數(shù)學(xué)模型的選擇

化工園區(qū)事故風(fēng)險(xiǎn)的動(dòng)態(tài)變化來源于熵增和熵減因子，風(fēng)險(xiǎn)熵增和風(fēng)險(xiǎn)熵減共同決定了事故風(fēng)險(xiǎn)的發(fā)展變化過程，由于事故狀態(tài)是唯一的狀態(tài)變量，風(fēng)險(xiǎn)熵增和熵減是2個(gè)控制變量，對于狀態(tài)變量個(gè)數(shù)為1個(gè)，控制變量個(gè)數(shù)為2個(gè)，應(yīng)對應(yīng)選擇表1中的尖點(diǎn)突變模型，用尖點(diǎn)突變理論進(jìn)行分析，如式1所示。

F(x)=x4+ux2+vx

(1)

式中：x為狀態(tài)變量，即用于反映園區(qū)的系統(tǒng)風(fēng)險(xiǎn)狀態(tài)；u=fu(u1,u2,…,un)為控制變量，是影響化工園區(qū)系統(tǒng)風(fēng)險(xiǎn)狀態(tài)，促使園區(qū)系統(tǒng)風(fēng)險(xiǎn)狀態(tài)發(fā)生突變的因素，即為熵增因子；v=fv(v1,v2,…,vm)為控制變量，是對化工園區(qū)緊急情況下采取的各種防控措施降低系統(tǒng)危險(xiǎn)程度、降低風(fēng)險(xiǎn)熵以減少事故損失的各種因素，即為熵減因子。

由F′(x)=0，得到風(fēng)險(xiǎn)熵突變模型平衡曲面M的方程為：

4x3+2ux+v=0

(2)

其對應(yīng)的平衡曲面如圖3所示。

圖3　系統(tǒng)風(fēng)險(xiǎn)熵尖點(diǎn)突變模型Fig.3　Cusp catastrophe model of system risk entropy

曲面由上、中、下3葉構(gòu)成，事故狀態(tài)隨狀態(tài)參量的變化而變化；隨著參量的增大，事故的平衡狀態(tài)由曲面下頁的穩(wěn)定狀態(tài)逐漸向上葉的事故激發(fā)狀態(tài)移動(dòng)。

2.2.2熵增和熵減估算模型

化工園區(qū)系統(tǒng)風(fēng)險(xiǎn)狀態(tài)的變化由熵增和熵減因子共同決定，設(shè)1個(gè)事故場景中共有n個(gè)風(fēng)險(xiǎn)熵指標(biāo)，整個(gè)事故的演化發(fā)展過程中需要計(jì)算的風(fēng)險(xiǎn)熵增(熵減)值的集合為T={A1，A2，…，Ak，…，Am}，Ak表示第k時(shí)刻的風(fēng)險(xiǎn)熵增(熵減)值(k=1，2，…，m)。分別構(gòu)建m×n階熵增矩陣Bg和m×n階熵減矩陣Bt,如式(3)所示。bg,ki表示第i個(gè)風(fēng)險(xiǎn)熵指標(biāo)第k時(shí)刻的風(fēng)險(xiǎn)熵增值(i=1,2,…,n)，bt,ki表示第i個(gè)風(fēng)險(xiǎn)熵指標(biāo)第k時(shí)刻的風(fēng)險(xiǎn)熵減值(i=1,2,…,n)。

(3)

設(shè)ri為事故的演化過程中第i個(gè)風(fēng)險(xiǎn)熵指標(biāo)的權(quán)重值(i=1,2,…,n)，其中0

(4)

潛在事故發(fā)展演化場景各個(gè)指標(biāo)的風(fēng)險(xiǎn)熵增和熵減值可以分別由式(5)，(6)求得。

(5)

(6)

式中：矩陣F為1×m階矩陣，且矩陣的每1個(gè)元素都為1；gi表示第i個(gè)風(fēng)險(xiǎn)熵指標(biāo)所有時(shí)刻熵增值之和，ti表示第i個(gè)風(fēng)險(xiǎn)熵指標(biāo)所有時(shí)刻熵減值之和(i=1,2,…,n)。

熵增和熵減值之和由各個(gè)指標(biāo)值根據(jù)權(quán)重計(jì)算得到，由式(7)，(8)分別計(jì)算系統(tǒng)整體的熵增u和熵減v值。

u=Gu×RT=F×Bg×RT

(7)

v=Tv×RT=F×Bt×RT

(8)

式中：RT為矩陣R的轉(zhuǎn)置矩陣。

2.2.3系統(tǒng)風(fēng)險(xiǎn)狀態(tài)的判斷

對式(2)中的x進(jìn)行求導(dǎo)，可得事故的臨界狀態(tài)。

F″(x)=12x2+2u=0

(9)

聯(lián)立式(2)，(9)可得系統(tǒng)的分叉集方程C，即平衡曲面上的折痕在控制平面上的投影。

Δ=8u3+27v2=0

(10)

分別根據(jù)式(7)和式(8)計(jì)算出的整體熵增和熵減，代入判別公式Δ=8u3+27v2判斷系統(tǒng)的狀態(tài)，若Δ>0，化工園區(qū)系統(tǒng)風(fēng)險(xiǎn)不被激發(fā)，保持穩(wěn)定狀態(tài)；若Δ<0，系統(tǒng)風(fēng)險(xiǎn)熵增和風(fēng)險(xiǎn)熵減構(gòu)成的坐標(biāo)(u,v)落在分叉集內(nèi)，系統(tǒng)風(fēng)險(xiǎn)發(fā)生突變，化工園區(qū)可能存在事故升級的風(fēng)險(xiǎn)；若Δ=0，化工園區(qū)系統(tǒng)風(fēng)險(xiǎn)進(jìn)入臨界狀態(tài)，此時(shí)若有不利因素影響可能促使風(fēng)險(xiǎn)狀態(tài)躍遷到突變狀態(tài)。

3　化工園區(qū)事故風(fēng)險(xiǎn)突變實(shí)例分析

3.1　某化工園區(qū)概況

某化工園區(qū)現(xiàn)規(guī)劃使用面積27.8 km2，園區(qū)內(nèi)32家企業(yè)構(gòu)成49個(gè)單元的危險(xiǎn)化學(xué)品重大危險(xiǎn)源，生產(chǎn)、

儲(chǔ)存、使用、經(jīng)營和運(yùn)輸過程中涉及到的危險(xiǎn)物質(zhì)共151種，各類危險(xiǎn)物質(zhì)的在線量為5.076 7 Mt，其中最多的危險(xiǎn)物質(zhì)為易燃液體，為4.928 Mt，占石化區(qū)危險(xiǎn)物質(zhì)總數(shù)的97.07%。這些危險(xiǎn)物質(zhì)廣泛存在于企業(yè)的生產(chǎn)經(jīng)營使用全過程中，其中大部分具有易燃易爆、有毒有害等危險(xiǎn)特性。

3.2　化工園區(qū)事故風(fēng)險(xiǎn)突變判斷思路

化工園區(qū)事故風(fēng)險(xiǎn)狀態(tài)判斷包含6個(gè)方面的內(nèi)容，主要包括突變模型研究、指標(biāo)參數(shù)確定，賦值等方面的內(nèi)容，具體判斷思路如圖4所示。

圖4　事故風(fēng)險(xiǎn)狀態(tài)判斷思路Fig.4　Accident risk state judgment

3.3　風(fēng)險(xiǎn)熵指標(biāo)體系構(gòu)建

根據(jù)調(diào)研統(tǒng)計(jì)，石化產(chǎn)業(yè)區(qū)儲(chǔ)罐數(shù)量共計(jì)682個(gè)，罐容達(dá)1 187 178 m3。該石化區(qū)危險(xiǎn)化學(xué)品庫房54個(gè)，共51 083.1 m2，儲(chǔ)存能力25 541.55 t，企業(yè)最大儲(chǔ)存量5 480.91 t。

3.3.1風(fēng)險(xiǎn)熵指標(biāo)建立

根據(jù)對園區(qū)事故風(fēng)險(xiǎn)來源分析可知，園區(qū)風(fēng)險(xiǎn)來源于企業(yè)內(nèi)部和企業(yè)外部的公共區(qū)域，考慮化工園區(qū)事故特點(diǎn)和危險(xiǎn)物質(zhì)情況，分為公共區(qū)域和企業(yè)端，確定風(fēng)險(xiǎn)熵監(jiān)測指標(biāo)，并針對一般工藝情況和環(huán)境情景給出風(fēng)險(xiǎn)熵指標(biāo)標(biāo)準(zhǔn)值、低限閾值、高限閾值和指標(biāo)符號，見表2。

表2　風(fēng)險(xiǎn)熵指標(biāo)

為了便于數(shù)據(jù)分析討論，選擇某?；穬?chǔ)存企業(yè)丁二烯泄漏情況下對數(shù)據(jù)進(jìn)行分析，在監(jiān)測平臺(tái)上采集風(fēng)險(xiǎn)監(jiān)測參數(shù)數(shù)據(jù)，數(shù)據(jù)采集過程中，分別監(jiān)測公共區(qū)域和企業(yè)的端的風(fēng)險(xiǎn)熵指標(biāo)變化情況，判斷泄漏情況下園區(qū)事故風(fēng)險(xiǎn)演化狀態(tài)。

3.3.2指標(biāo)權(quán)重確定方法

在風(fēng)險(xiǎn)熵指標(biāo)中，不同的指標(biāo)對于事故演化的重要程度不同，分配的權(quán)重系數(shù)也不同，并對最終的判斷結(jié)果有影響，因此要采用合理的方法賦予權(quán)重，采用改進(jìn)的層次分析法[16](AHP)處理專家調(diào)查問卷，并計(jì)算各風(fēng)險(xiǎn)熵指標(biāo)權(quán)重值。

3.3.3判斷矩陣的一致性檢驗(yàn)

在層次分析法的應(yīng)用中，判斷矩陣的一致性檢驗(yàn)是為了使得到的結(jié)論基本合理，根據(jù)9個(gè)風(fēng)險(xiǎn)熵指標(biāo)構(gòu)造了9階判斷矩陣。

1)計(jì)算一致性指標(biāo)CI。

CI=(λmax-n)/(n-1)

(11)

式中：λmax為判斷矩陣最大特征值；n為判斷矩陣階數(shù)。

2)根據(jù)表3查找平均隨機(jī)一致性指標(biāo)RI值。

3)計(jì)算一致性比率CR。

CR=CI/RI

(12)

當(dāng)CR<0.1時(shí)，認(rèn)為判斷矩陣具有滿意的一致性，可以用于評估，否則需要重新調(diào)整判斷矩陣。

表3　一致性指標(biāo)RI值

3.4　事故風(fēng)險(xiǎn)突變分析

石化區(qū)某?；反鎯?chǔ)企業(yè)球罐發(fā)生丁二烯泄漏事故，該企業(yè)采用DCS集散型控制系統(tǒng)，罐區(qū)儲(chǔ)罐主要工藝參數(shù)、全廠的可燃有毒氣體檢測報(bào)警器、火災(zāi)報(bào)警器引至全廠DCS中心控制室，并同時(shí)與園區(qū)應(yīng)急指揮平臺(tái)相連接。平臺(tái)監(jiān)測指標(biāo)出現(xiàn)明顯波動(dòng)，有可能由泄漏事故引發(fā)二次事故，根據(jù)企業(yè)內(nèi)和公共區(qū)域的監(jiān)測探頭從園區(qū)內(nèi)采集的數(shù)據(jù)，分為5個(gè)時(shí)刻分析園區(qū)內(nèi)系統(tǒng)風(fēng)險(xiǎn)熵情況，熵增值的判斷主要根據(jù)指標(biāo)檢測的波動(dòng)情況，熵減值的判斷主要根據(jù)應(yīng)急措施，系統(tǒng)自適應(yīng)特性，本質(zhì)安全設(shè)計(jì)等。

風(fēng)險(xiǎn)熵指標(biāo)評分值區(qū)間為(-1.00,1.00)，其中，熵增的取值區(qū)間為(-1.00,0.00)，劃分為5個(gè)等級：特別顯著(-1.00，-0.80]，比較顯著(-0.80,-0.60]，顯著(-0.60，-0.40]，一般(-0.40，-0.20]，較小(-0.20,0.00)。相應(yīng)熵減值的取值區(qū)間為(0.00,1.00)，也劃分為5個(gè)等級：較小(0.00, 0.20]，一般(0.20，0.40]，顯著(0.40，0.60]，比較顯著(0.60, 0.80]，特別顯著(0.80，1.00)，根據(jù)每位專家對每1指標(biāo)不同時(shí)刻的等級打分區(qū)間取中值再求取它們的平均值(即表4中指標(biāo)的風(fēng)險(xiǎn)熵增與熵減值)，經(jīng)過檢驗(yàn)，專家打分符合一致性要求，然后根據(jù)化工園區(qū)系統(tǒng)風(fēng)險(xiǎn)熵突變模型進(jìn)行判斷。

表4　園區(qū)風(fēng)險(xiǎn)熵值分析

1)權(quán)重計(jì)算。根據(jù)調(diào)查問卷，用改進(jìn)的AHP法計(jì)算得到各風(fēng)險(xiǎn)熵指標(biāo)的權(quán)重見表4，其一致性比率CR=0.000 07<0.10，滿足一致性要求。

2)事故風(fēng)險(xiǎn)狀態(tài)判斷。由表4數(shù)據(jù)，根據(jù)式(3)，(5)，(7)可得熵增值u=Gu×RT=P×Bg×RT=-1.503 84。同理根據(jù)式(3)，(6)，(8)可得熵減值v=Tv×RT=P×Bt×RT=0.816 84，帶入分叉集方程(10)，計(jì)算得：

Δ=8u3+27v2=-9.192 76<0

可知在丁二烯發(fā)生泄漏的情況下化工園區(qū)發(fā)生了風(fēng)險(xiǎn)狀態(tài)躍遷突變，系統(tǒng)處于無序崩潰的邊緣，此時(shí)極易由泄漏事故導(dǎo)致二次事故發(fā)生，因此要注意風(fēng)險(xiǎn)管控，及時(shí)降低系統(tǒng)風(fēng)險(xiǎn)，防止事故升級。

3.5　分析與建議

風(fēng)險(xiǎn)熵指標(biāo)數(shù)值如圖5所示。分析顯示，在L1(公共區(qū)域可燃?xì)怏w)、L2(公共區(qū)域有毒氣體)、L3(環(huán)境溫度)、L4(相對濕度)、L5(風(fēng)速)、L6(企業(yè)端可燃?xì)怏w)、L7(企業(yè)端有毒氣體)、L8(罐內(nèi)溫度)、L9(介質(zhì)液位)等9個(gè)指標(biāo)中， L1，L2，L6，L7指標(biāo)相對其他指標(biāo)熵增值較大，對系統(tǒng)風(fēng)險(xiǎn)突變影響較為顯著，這是由于丁二烯發(fā)生泄漏，企業(yè)端和公共區(qū)域的可燃?xì)怏w和有毒氣體含量明顯增加，監(jiān)測指標(biāo)數(shù)值發(fā)生波動(dòng)。L6，L7熵減值相對于其他指標(biāo)也較為明顯，可知在發(fā)生泄漏事故后企業(yè)采取了積極的應(yīng)急救援防護(hù)措施。

圖5　風(fēng)險(xiǎn)熵指標(biāo)數(shù)值Fig.5　Numerical value of risk entropy index

圖6　各因素占風(fēng)險(xiǎn)熵增比例Fig.6　The ratio of each factor to risk entropy increases

企業(yè)端和公共區(qū)域所占的風(fēng)險(xiǎn)熵增比例如圖6所示。分析表明，企業(yè)端是風(fēng)險(xiǎn)熵增的主要因素，可燃有毒氣體泄漏很容易導(dǎo)致化工園區(qū)系統(tǒng)崩潰，增加二次事故風(fēng)險(xiǎn)突變概率，因此要增加企業(yè)內(nèi)的監(jiān)測探頭，以做到事故的早期預(yù)警，配備應(yīng)急裝備物資，提升事故應(yīng)對能力，加強(qiáng)安全管理和事故隱患排查治理，針對企業(yè)可能發(fā)生的事故類型定期開展應(yīng)急演練，由于化工事故具有發(fā)展迅速的特點(diǎn)，因此事故的初始階段是最有利的處置時(shí)期。做好特殊情況下園區(qū)的應(yīng)急預(yù)案和風(fēng)險(xiǎn)防控措施，有利于將事故消滅在初始階段，防止二次事故的發(fā)生，降低事故損失。

4　結(jié)論

1)借鑒風(fēng)險(xiǎn)熵理論，將風(fēng)險(xiǎn)熵引入化工園區(qū)風(fēng)險(xiǎn)分析，提出化工園區(qū)事故風(fēng)險(xiǎn)躍遷概念，分析其形成機(jī)理。分析表明，化工園區(qū)事故風(fēng)險(xiǎn)的動(dòng)態(tài)性和不確定性來源于熵增和熵減因子，風(fēng)險(xiǎn)熵增和熵減共同決定了事故風(fēng)險(xiǎn)的突變演化過程。

2)根據(jù)化工園區(qū)事故風(fēng)險(xiǎn)特點(diǎn)，結(jié)合尖點(diǎn)突變模型，構(gòu)建化工園區(qū)系統(tǒng)風(fēng)險(xiǎn)熵突變模型，結(jié)果表明，該模型能對熵增(減)值進(jìn)行準(zhǔn)確的計(jì)算，可以為化工園區(qū)事故風(fēng)險(xiǎn)判斷提供依據(jù)。

3)針對某危化品儲(chǔ)存企業(yè)丁二烯泄漏事故進(jìn)行分析，并以監(jiān)測數(shù)據(jù)為基礎(chǔ)，結(jié)合問卷調(diào)查和專家打分，分析園區(qū)風(fēng)險(xiǎn)熵增(減)值，根據(jù)計(jì)算結(jié)果判斷化工園區(qū)事故風(fēng)險(xiǎn)狀態(tài)發(fā)生突變，有發(fā)生二次事故(火災(zāi)、爆炸、中毒等)的危險(xiǎn)。對各因素進(jìn)行分析，得出企業(yè)端是導(dǎo)致風(fēng)險(xiǎn)熵增的主要因素，據(jù)此可提出相應(yīng)的風(fēng)險(xiǎn)防控措施。

[1]邵輝, 楊麗丹, 李晶. 化工園區(qū)重大危險(xiǎn)源分級方法研究[J]. 工業(yè)安全與環(huán)保, 2011, 37(9):12-14.

SHAO Hui,YANG Lidan,LI Jing. The study on major hazards ranking for chemical industrial park [J]. Industrial Safety and Environmental Protection, 2011, 37(9):12-14.

[2]陳文濤, 葛世友. 基于三角模糊數(shù)的化工園區(qū)事故應(yīng)急救援能力評價(jià)[J]. 災(zāi)害學(xué), 2015,30(2):167-171.

CHEN Wentao, GEShiyou. Capability evaluation on emergency response ability for accidents in chemical industry parks based on triangular fuzzy number [J]. Journal of Catastrophology, 2015,30(2):167-171.

[3]陳國華, 安霆, 陳培珠. 危險(xiǎn)化學(xué)品多米諾事故應(yīng)急資源需求量估算模型[J]. 中國安全科學(xué)學(xué)報(bào), 2015, 25(4):87-93.

CHEN Guohua, AN Ting, CHEN Peizhu. Estimation model of emergency resource demands for chemical accidents involving Domino effect [J]. China Safety Science Journal, 2015, 25(4):87-93.

[4]Garrido P L, Goldstein S, Lebowitz J L. Boltzmann entropy for dense fluids not in local equilibrium.[J]. Physical Review Letters, 2003, 92(5):050602.

[5]胡瑞敏,呂海濤,陳軍.基于風(fēng)險(xiǎn)熵和Neyman-Pearson準(zhǔn)則的安防網(wǎng)絡(luò)風(fēng)險(xiǎn)評估研究[J]. 自動(dòng)化學(xué)報(bào), 2014, 40(12):2737-2746.

HU Ruimin, LYU Haitao, CHEN Jun.Risk evaluation model of security and protection network based on risk entropy and Neyman-Pearson criterion[J]. Acta Automatica Sinica, 2014, 40(12):2737-2746.

[6]于海峰, 王延章, 盧小麗,等. 基于知識(shí)元的突發(fā)事件風(fēng)險(xiǎn)熵預(yù)測模型研究[J]. 系統(tǒng)工程學(xué)報(bào), 2016, 31(1):117-126.

YU Haifeng, WANG Yanzhang, LU Xiaoli,et al. Emergency risk entropy forecasting model based on knowledge element[J]. Journal of Systems Engineering, 2016, 31(1):117-126.

[7]Johnson IV J J, Tolk A, Sousa-Poza A. A theory of emergence and entropy in systems of systems[J]. Procedia Computer Science, 2013, 20: 283-289.

[8]Spartalis S, Iliadis L, Maris F. An innovative risk evaluation system estimating its own fuzzy entropy[J]. Mathematical & Computer Modelling, 2007, 46(46):260-267.

[9]陳偉珂, 向蘭蘭. 基于熵及耗散結(jié)構(gòu)的公共安全突發(fā)事件的過程分析研究[J]. 中國軟科學(xué), 2007(10):149-154.

CHEN Weike, XIANG Lanlan. The process analysis on pubic safety emergencies based on entropy and dissipative structure theory [J]. China Soft Science, 2007(10):149-154.

[10]汪送, 王瑛, 杜純,等. 復(fù)雜系統(tǒng)風(fēng)險(xiǎn)熵的涌現(xiàn)與動(dòng)力學(xué)傳播分析[J]. 安全與環(huán)境工程, 2013, 20(2):118-120，134.

WANG Song,WANG Ying, DU Chun,et al. Emergence and dynamic transmission of risk entropy in complex system [J].Safety and Environmental Engineering, 2013, 20(2):118-120，134.

[11]李健, 白曉昀, 任正中,等. 2011—2013年我國危險(xiǎn)化學(xué)品事故統(tǒng)計(jì)分析及對策研究[J]. 中國安全生產(chǎn)科學(xué)技術(shù), 2014，10(6):142-147.

LI Jian，BAI Xiaoyun，REN Zhengzhong，et al. Statistic analysis and countermeasures of hazardous chemicals accidents occurring in China during 2011-2013 [J]. Journal of Safety Science and Technology, 2014，10(6):142-147.

[12]Darbra R M, Palacios A, Casal J. Domino effect in chemical accidents: main features and accident sequences.[J]. Journal of Hazardous Materials, 2010, 183(1-3):565-573.

[13]凌曉東, 黨文義, 于安峰,等. 環(huán)己烷泄漏事故模擬與風(fēng)險(xiǎn)控制措施研究[J]. 中國安全生產(chǎn)科學(xué)技術(shù), 2013, 9(8):33-39.

LING Xiaodong, DANG Wenyi,YU Anfeng,et al. Research on the simulation of cyclohexane leakage accidents and development of risk control measures[J]. Journal of Safety Science and Technology,2013, 9(8):33-39.

[14]凌復(fù)華. 突變理論及其應(yīng)用[M]. 上海：上海交通大學(xué)出版社, 1987.

[15]苗東升. 系統(tǒng)科學(xué)大學(xué)講稿[M]. 北京：中國人民大學(xué)出版社, 2007.

[16]陳國華, 梁韜, 張華文. 城域承災(zāi)能力評估研究及其應(yīng)用[J]. 安全與環(huán)境學(xué)報(bào), 2008, 8(2):156-162.

CHEN Guohua, LIANG Tao, ZHANG Huawen. Evaluation model for urban-regional disaster-bearing capacity and its application[J]. Journal of Safety and Environment, 2008, 8(2):156-162.