李 陽，張明廣，錢城江，桂 陽，冷 源

基于理想優(yōu)基點法的原油儲罐區(qū)布局多目標(biāo)決策

李 陽1，張明廣1，錢城江2，桂 陽1，冷 源1

(1.南京工業(yè)大學(xué)江蘇省城市與工業(yè)安全重點實驗室 江蘇 南京 210009；

2.南京國泰協(xié)創(chuàng)安全科技有限公司，江蘇 南京 211500)

原油罐區(qū)是重大的火災(zāi)爆炸危險源，儲罐之間易產(chǎn)生多米諾效應(yīng)加之危險源與居民區(qū)的距離較近，使罐區(qū)的風(fēng)險較大?；诙嗄繕?biāo)決策理論對原油罐區(qū)的布局進(jìn)行研究，首先分別計算出每種布局方案下土地使用面積、罐區(qū)內(nèi)多米諾效應(yīng)發(fā)生的綜合失效概率、罐區(qū)對周邊居民區(qū)的個人風(fēng)險；然后運用多目標(biāo)決策理論中的理想優(yōu)基點法對各備選方案進(jìn)行綜合比較，選擇最優(yōu)的布局方案；最后通過對某居民區(qū)附近的原油罐區(qū)布局計算，從而驗證該方法的可行性。

原油罐區(qū)布局；理想優(yōu)基點法；多目標(biāo)決策；多米諾效應(yīng)；個人風(fēng)險

原油罐區(qū)是原油儲、運的重要場所，儲罐數(shù)量多、儲量大、生產(chǎn)設(shè)施多、工藝復(fù)雜，是火災(zāi)爆炸的危險源。國內(nèi)消防部門對煉油廠的調(diào)查結(jié)果表明，在全部油罐火災(zāi)中，原油罐占40%，居油罐火災(zāi)事故的第一位[1]。臺灣高雄大學(xué)[2]統(tǒng)計了1960—2003年間242起儲罐事故，其中原油儲罐事故66起，排名第一。原油儲罐發(fā)生火災(zāi)事故的特點主要表現(xiàn)在火焰溫度高、輻射熱強(qiáng)，相鄰儲罐在強(qiáng)輻射熱的影響下，液態(tài)油品產(chǎn)生相變，油罐內(nèi)部因壓力升高而導(dǎo)致油罐破損[3]，很容易引發(fā)二次事故，導(dǎo)致事故多米諾效應(yīng)，從而使得火勢擴(kuò)大，形成罐區(qū)火災(zāi)。

國內(nèi)對原油罐區(qū)布局的設(shè)計主要依據(jù)各種指令性的規(guī)范，如《石油化工企業(yè)設(shè)計防火規(guī)范》(GB 50160—2008)、《石油化工企業(yè)廠區(qū)總平面布置》(SH/T 3053—2002)、《石油庫設(shè)計規(guī)范》(GB 50074—2014)等。但隨著原油罐區(qū)大型化，居民區(qū)不斷向城市的郊區(qū)發(fā)展的趨勢，標(biāo)準(zhǔn)規(guī)范中的安全距離已無法在有限的空間下，能夠做到合理布局儲存裝置又不浪費一寸土地，且能確保裝置間具有足夠的安全距離、減少和消除多米諾效應(yīng)、減弱罐區(qū)對周邊人員密集區(qū)域的風(fēng)險[4]。目前，國內(nèi)雖然對原油罐區(qū)布局的研究較多[4-5]，但主要集中于單目標(biāo)優(yōu)化方法的研究，如安全距離法、基于后果的方法、基于風(fēng)險的方法和綜合性的方法等，而對罐區(qū)布局的多目標(biāo)優(yōu)化方法的研究較少[6-8]。國外對原油罐區(qū)布局的設(shè)計主要依據(jù)Seveso Directive III，且主要是從風(fēng)險、脆弱性的角度分析其合理性，忽略了儲罐間的多米諾效應(yīng)[9]。

罐區(qū)裝置布局實質(zhì)上是個多目標(biāo)決策問題，如何在有限的土地資源下，通過罐區(qū)的不同布局方案決策，選擇一個能使罐區(qū)內(nèi)儲罐間的多米諾效應(yīng)、對周邊設(shè)施的風(fēng)險達(dá)到最小的布局方案，則是本文主要的研究內(nèi)容。本文采用理想優(yōu)基點法對原油罐區(qū)布局方案進(jìn)行多目標(biāo)決策研究。首先，根據(jù)有限個評價對象與理想化目標(biāo)的接近程度進(jìn)行排序，在現(xiàn)有的對象中進(jìn)行相對優(yōu)劣的評價；然后將原始數(shù)據(jù)進(jìn)行同趨勢和歸一化處理，可從同一指標(biāo)與最優(yōu)指標(biāo)的距離之比計算中體現(xiàn)出同一指標(biāo)間的接近度，這在一定程度上反映了離散程度，排序結(jié)果充分利用了原始數(shù)據(jù)信息；最后運用各級標(biāo)準(zhǔn)和待評價對象與規(guī)范統(tǒng)一的正、負(fù)理想點的均方根距離來計算其對應(yīng)貼近度，進(jìn)而得出相應(yīng)的評價結(jié)果[10]。該方法在有限的數(shù)據(jù)情況下，能夠定量反映不同評價單元的優(yōu)劣程度，且直觀、可靠、操作性強(qiáng)。

1 原油罐區(qū)布局多目標(biāo)決策體系構(gòu)建

本文從土地使用面積、儲罐綜合失效概率以及對居民區(qū)的個人風(fēng)險三個方面，運用理想優(yōu)基點法選擇原油罐區(qū)最優(yōu)布局方案，圖1為建立的原油儲罐布局多目標(biāo)決策體系，整個體系分成計算層、屬性層、方案層和決策層。

圖1 原油罐區(qū)布局多目標(biāo)決策體系Fig.1 System of the layout of a crude oil tank farm based on multi-objective decision making

2 實例應(yīng)用與分析

本文以某居民區(qū)附近的原油罐區(qū)的布局設(shè)計為例，詳細(xì)說明建立的原油罐區(qū)布局多目標(biāo)決策體系的整個計算、決策過程。

2.1 原油罐區(qū)布局的基本情況

如圖2所示，在距離居民區(qū)150 m處計劃建造一個(座)原油庫，常壓儲罐儲存原油12 000 m3。本文的目標(biāo)是在考慮土地資源、罐區(qū)內(nèi)多米諾效應(yīng)以及儲罐對居民區(qū)的風(fēng)險的情況下，通過多目標(biāo)決策，得到最優(yōu)的原油罐區(qū)布局。

圖2 某原油罐區(qū)Fig.2 Hypothetical crude oil tank farm

在滿足相關(guān)規(guī)范[11-12]的前提下，本文設(shè)計了6種可能的原油儲罐區(qū)布局方案，詳見圖3，各布局方案下的儲罐參數(shù)見表1。

圖3 原油儲罐區(qū)的不同布局方案Fig.3 Alternatives for the layout of a crude oil tank farm

2.2 原油罐區(qū)多米諾綜合失效概率計算

Abdolhamidzadeh等[13]運用蒙特卡羅法(Monte Carlo Method)開發(fā)出一種在理論分析的基礎(chǔ)上計算裝置多米諾效應(yīng)概率的計算方法。本課題組在其研究成果的基礎(chǔ)上，對其算法做了相應(yīng)的改進(jìn)，如在傳播概率與隨機(jī)數(shù)比較時增加一個裝置間的內(nèi)循環(huán)，以降低運算時間，并根據(jù)改進(jìn)的算法開發(fā)出一款計算模擬軟件，稱之為“事故連鎖多米諾效應(yīng)概率計算方法”。該計算方法基于蒙特卡羅原理，運用多次試驗來模擬實際的多裝置系統(tǒng)發(fā)生事故連鎖多米諾效應(yīng)的情況，其最終的輸出結(jié)果就是系統(tǒng)中每個裝置的多米諾綜合失效概率(裝置自身失效概率以及在多米諾效應(yīng)作用下的失效概率)。

表1 不同原油罐區(qū)布局相關(guān)參數(shù)

原油儲罐初始失效概率為5×10-6[14]，對于儲罐間的傳播概率Pij的計算，Bagster等[15]提出儲罐間傳播概率與它們之間距離的平方呈反比的關(guān)系，即距離越遠(yuǎn)導(dǎo)致發(fā)生多米諾效應(yīng)的概率越小，因此給出如下公式：

(1)

(2)

式中:rth為初級事故能引起破壞的最大距離，即達(dá)到多米諾效應(yīng)破壞閾值的距離(m);rij為儲罐之間的距離(m);Qiy為池火災(zāi)下儲罐之間的熱輻射通量(kW/m2);Ith為常壓容器池火熱輻射強(qiáng)度的破壞閾值(kW/m2)，通常取15kW/m2[9]。

根據(jù)公式(2)，可計算得到各布局方案儲罐多米諾效應(yīng)破壞閥值的距離rth，詳見表2。

表2 各布局方案下儲罐多米諾效應(yīng)破壞閾值的距離

以布局方案X5為例，根據(jù)公式(1)分別計算出不同儲罐間的傳播概率Pij，并輸入軟件，最終計算得到布局X5的儲罐多米諾綜合失效概率，見圖4。

圖4 布局方案X5的模擬結(jié)果Fig.4 Simulation results of the layout scheme X5

同理，可依次計算出每種布局方案下儲罐的多米諾綜合失效概率，詳見表3。由表3可見，隨著罐區(qū)內(nèi)儲罐數(shù)目的增加，在考慮儲罐間的多米諾效應(yīng)的情況下，儲罐的多米諾綜合失效概率變大。

表3 不同布局方案下儲罐的多米諾綜合失效概率

2.3 原油罐區(qū)對居民區(qū)的個人風(fēng)險計算

2.3.1 后果分析

脆弱性模型是一種被廣泛認(rèn)可且用于評價人體對有毒物質(zhì)、熱輻射或超壓的劑量-反應(yīng)關(guān)系的方法，該模型由標(biāo)準(zhǔn)正態(tài)概率分布函數(shù)的累積表達(dá)推導(dǎo)而來，其計算公式如下：

(其中，D為劑量，D=6×10-3I1.33t[16];μ和σ為高斯分布的中位數(shù)和方差)；Y可通過下式計算:

Y=-14.9+2.56lnD

(4)

在考慮多個事故場景的情況下，總體的脆弱性Fd可以通過單個事故場景脆弱性組合進(jìn)行計算。Cozzani等[17]提出的應(yīng)用較為廣泛的方法是假設(shè)整體脆弱性為多米諾事故中所有場景的死亡概率之和(上限為1):

(5)

式中:Fp為初始事故造成的人體脆弱性；Fd，i為多米諾事故場景i造成的人體脆弱性。

2.3.2 多米諾場景的事故概率

多米諾場景的事故概率的計算考慮同一初始事故引發(fā)的1個以上的2級事故。若不考慮2級事故的進(jìn)一步升級，可認(rèn)為升級事故是獨立事件，則含k個發(fā)生2級事故的設(shè)備同時發(fā)生第m種多米諾場景的期望概率為[18]

(6)

(7)

2.3.3 考慮多米諾效應(yīng)的居民區(qū)個人風(fēng)險

基于以上分析，在考慮多米諾效應(yīng)的前提下居民區(qū)的個人風(fēng)險R[18]可由下式計算:

(8)

2.3.4 不同布局方案下考慮多米諾效應(yīng)的居民區(qū)個人風(fēng)險值

通過以上計算方法，最終可計算得出各布局方案下考慮多米諾效應(yīng)的居民區(qū)個人風(fēng)險值，見表4。

表4 不同布局方案下考慮多米諾效應(yīng)的居民區(qū)

2.4 多屬性決策的理想優(yōu)基點法

在對各布局方案進(jìn)行決策時，主要考慮了罐區(qū)土地面積、罐區(qū)儲罐綜合失效概率以及居民區(qū)的個人風(fēng)險值這三個屬性，各屬性的優(yōu)劣狀況最終都可以用定量的數(shù)值表示，因此在整個多屬性決策過程中可以通過數(shù)值的計算得出最優(yōu)解。

2.4.1 求方案的優(yōu)基點

采用式X={X1,X2,X3,X4,X5,X6}的約束集表示6個可行的布局方案，根據(jù)不同布局方案各屬性狀況信息(見表5)，將其轉(zhuǎn)換為對所有屬性狀況都用定量表達(dá)的如下數(shù)值矩陣：

表5 不同布局方案各屬性狀況

布局方案各屬性狀況土地面積A1/m2儲罐綜合失效概率A2/(10-6)個人風(fēng)險值A(chǔ)3/(×10-6)X116005.023.4X220445.428.6X3214210.321.8X424488.427.1X5270411.48.56X6285613.04.17

因為屬性A1、A2、A3這三對屬性值都是要求達(dá)到最優(yōu)極小化，而且都有不同的量綱單位，因此需要對矩陣中的值進(jìn)行無量綱化處理，即

(9)

2.4.2 求理想優(yōu)基點

由各優(yōu)基點對應(yīng)分量中的最大值組成的點，對于方案來說是最“理想”的點；相反地，由最小值組成的點，對于方案來說是最“不理想”的點，分別記[19]：

(10)

作m維屬性空間Rm中的點，有

(11)

(12)

式中：zΔ為理想優(yōu)基點；z為反理想優(yōu)基點。

2.4.3 計算綜合距離

優(yōu)基點在三維空間上的反映就是一個點，找最優(yōu)的方案，其實就是在某距離意義下尋找最接近理想優(yōu)基點而最遠(yuǎn)離反理想優(yōu)基點的點。如圖5所示，在三維空間畫出優(yōu)基點z1,z2,z3,z4,z5,z6以及理想優(yōu)基點zΔ和反理想優(yōu)基點z。

圖5 各優(yōu)基點在三維空間中的分布Fig.5 Distribution of each cardinal point in three- dimension space

各優(yōu)基點與理想優(yōu)基點和反理想優(yōu)基點的距離計算公式如下：

(13)

以di作為衡量zi到zΔ和z的綜合距離，其計算公式如下[19]:

(14)

根據(jù)公式(13)和(14)，可計算得到不同布局方案下各優(yōu)基點與理想優(yōu)基點和反理想優(yōu)基點的距離以及綜合距離，見表6。

2.4.4 綜合距離優(yōu)選排序

根據(jù)距離最小原則，設(shè)xi，xj∈X(i≠j)，若di

由表6的計算結(jié)果可得：d4>d3>d2>d1>d5>d6，根據(jù)距離最小原則，可知方案X6是最優(yōu)偏好解。

表6 不同布局方案下、di的計算結(jié)果

雖然方案X6的儲罐數(shù)量、儲罐間的多米諾效應(yīng)、土地使用面積是最多的，但罐區(qū)對周邊的居民區(qū)的個人風(fēng)險是最小的。由于理想優(yōu)基點法是用空間距離大小來表征方案的優(yōu)劣，出現(xiàn)最優(yōu)結(jié)果是方案X6的情況，則是因為在計算時個人風(fēng)險屬性的極差相對于其他屬性而言較大，所以導(dǎo)致風(fēng)險所占的比重較大。

3 結(jié)論與建議

(1) 本文建立了原油罐區(qū)布局多目標(biāo)決策體系，從罐區(qū)的土地面積、多米諾效應(yīng)、對周邊居民區(qū)的個人風(fēng)險三個方面綜合考慮原油罐區(qū)的布局，并以理想優(yōu)基點法作為罐區(qū)布局方案選擇的多目標(biāo)決策方法，最終得到最優(yōu)的布局方案。

(2) 理想優(yōu)基點法是以靠近正理想解和遠(yuǎn)離負(fù)理想解兩個基準(zhǔn)作為評價依據(jù)來確定方案的排序，具有直觀的幾何意義。該方法基于數(shù)據(jù)樣本本身，具有一定的客觀性，且對原始數(shù)據(jù)的利用比較充分，信息損失比較少，適用于化工裝置的多目標(biāo)決策。但本文的研究還存在數(shù)據(jù)量有限、數(shù)據(jù)波動較大、沒有典型的分布規(guī)律等缺陷，建議下一步可以通過擴(kuò)充數(shù)據(jù)容量、對樣本數(shù)據(jù)進(jìn)行內(nèi)在關(guān)聯(lián)分析等方法來解決此問題。

[1] 許學(xué)瑞，帥健，吳宗之.大型原油庫火災(zāi)定量風(fēng)險評價[J].油氣儲運，2015(5):482-487.

[2] Chang J I，Lin C C.A study of storage tank accidents[J].JournalofLossPreventionintheProcessIndustries，2006(19)：51-59.

[3] 苑靜.石油儲罐火災(zāi)爆炸危害控制的研究應(yīng)用[D].天津：天津理工大學(xué)，2008.

[4] 安君.多米諾效應(yīng)風(fēng)險評價在化工裝置優(yōu)化布局中的應(yīng)用[D].沈陽：沈陽航空航天大學(xué)，2013.

[5] 葉君樂，蔣軍成，張明廣，等.基于本質(zhì)安全目標(biāo)的化工廠裝置平面布局優(yōu)化[J].安全與環(huán)境學(xué)報，2011，11(1):167-171.

[6] 米紅梅，張悅，宋占兵.基于VEGA和DCGA的化工園區(qū)布局優(yōu)化研究[J].中國安全生產(chǎn)科學(xué)技術(shù)，2013，9(5):45-51.

[7] 張悅.基于風(fēng)險分析的化工園區(qū)布局優(yōu)化方法研究[D].北京：中國礦業(yè)大學(xué)(北京)，2013.

[8] 陳國華，施文松，趙遠(yuǎn)飛.基于風(fēng)險的化工園區(qū)布局優(yōu)化決策支持系統(tǒng)[J].中國安全科學(xué)學(xué)報，2012，22(7):141-146.

[9] Khakzada N，Reniers G.Risk-based design of process plants with regard to domino effects and land use planning[J].JournalofHazardousMaterials,2015，299：289-297.

[10]問國強(qiáng).基于指標(biāo)規(guī)范變換TOPSIS法的地表水環(huán)境質(zhì)量評價[J].安全與環(huán)境工程,2015(3):65-69.

[11]GB 50074—2014 石油庫設(shè)計規(guī)范[S].

[12]GB 50160—2008 石油化工企業(yè)設(shè)計防火規(guī)范[S].

[13]Abdolhamidzadeh B，Abbasi T，Rashtchian D，et al.A new method for assessing domino effect in chemical process industry[J].JournalofHazardousMaterials，2010，182(1/2/3):416-426.

[14]Cozzani V，Gubinelli G，Salzano E.Escalation thresholds in the assessment of domino accidental events[J].JournalofHazardousMaterials，2006，A129:1-21.

[15]Bagster D F，Pitblado R M．The estimation of domino incident frequencies—An approach[J]．ProcessSafetyEnvironmentalProtection,1991，69(4)：195-199．

[16]邢志祥，王云慧，楊扣華，等.化學(xué)流程工業(yè)多米諾效應(yīng)風(fēng)險評估的研究進(jìn)展[J].中國安全科學(xué)學(xué)報，2014(10):144-150.

[17]Cozzani V，Antonioni G，Spadom G.Quantitative assessment of domino scenarios by a GIS-based software tool[J].JournalofLossPreventionintheProcessIndustries，2006，19(5):463- 477.

[18]馬科偉.基于多米諾效應(yīng)的區(qū)域定量風(fēng)險評估方法研究[D].杭州：浙江工業(yè)大學(xué)，2010.

[19]胡毓達(dá).多目標(biāo)決策:實用模型和選優(yōu)方法[M].上海：上?？茖W(xué)技術(shù)出版社，2010.

Multi-objective Decision Making for Layout of Crude Oil Tank Farms Based on the Method of Ideal Outranking Cardinal Point

LI Yang1,ZHANG Mingguang1,QIAN Chengjiang2,GUI Yang1,LENG Yuan1

(1.JiangsuKeyLaboratoryofUrbanandIndustrialSafety,NanjingTechUniversity,Nanjing210009,China;2.NanjingGuoTaiXieChuangSafetyTechnologyCo.,Ltd,Nanjing211500,China)

A crude oil tank farm is a major hazard source of fires and explosions.As domino effect tends to occur between storage tanks and hazard sources and the residential areas are close,the risk becomes relatively high.This paper conducts a research on the reasonable layout of crude oil tank farms based on the method of multi-objective decision making.First,the paper calculates the land use area,the comprehensive failure probabilities of domino effect in tank farms,and the individual risk of tank farms to the surrounding residential areas under each layout scheme respectively.Then,the paper applies the method of outranking cardinal point from the theory of multi-objective decision making to compare various alternatives,and selects the optimal layout scheme.Last,the paper calculates the layout of a crude oil tank farm near a residential area to testify the feasibility of the method.Key words:crude oil tank farm layout;method of ideal outranking cardinal point;multi-objective decision making;domino effect;individual risk

韋朝陽(1965-)，男，研究員，主要從事重金屬污染的生態(tài)與健康風(fēng)險及環(huán)境修復(fù)等方面的研究。E-mail:weicy@igsnrr.ac.cn

1671-1556(2016)02-0066-05

2015-06-10

2016-01-06

江蘇省自然科學(xué)基金面上項目(BK2012824)

李 陽(1991—)，男，碩士研究生，主要研究方向為化工裝置安全評價。E-mail:LY_njtech@126.com

X928.03;TE88

A

10.13578/j.cnki.issn.1671-1556.2016.02.013