劉金玉，王殿生，胡玉琴

(中國石油大學(xué)(華東) 理學(xué)院，山東 青島 266580)

引 言

靜電是油庫火災(zāi)爆炸事故的主要成因[1-3]，合理評估油庫區(qū)的靜電風(fēng)險具有重要的現(xiàn)實(shí)意義。于格非[4]研究了輕質(zhì)油品在裝卸和運(yùn)輸過程中靜電產(chǎn)生的主要因素以及靜電放電和靜電點(diǎn)火的條件等，但沒有采用定量的安全評價方法對油庫靜電安全性進(jìn)行評估，無法為運(yùn)營中的油庫提供行之有效的安全評價模型；紹虎[5]利用事故樹評價方法分析了儲罐靜電安全的影響因素，但僅適用于原油儲罐安全評價，無法運(yùn)用于整個油庫；易俊[6]、李凌峰[7-8]等人構(gòu)建了二階綜合模糊評估模型，該模型雖將靜電作為油庫火災(zāi)風(fēng)險的一個影響因素進(jìn)行研究，但并未深入分析靜電事故誘因。因此，目前關(guān)于油庫區(qū)靜電安全性的系統(tǒng)化評價研究還很少。

1 油庫靜電事故評價指標(biāo)體系

油庫系統(tǒng)中，石油因為流動、沉降、晃動、噴射、沖擊、過濾等運(yùn)動，時刻在產(chǎn)生電荷[9-10]，由于油品是電的不良導(dǎo)體，油品中電荷逸散的速度遠(yuǎn)遠(yuǎn)小于靜電產(chǎn)生的速度，靜電荷不斷積累，積累到一定程度就易發(fā)生靜電放電，導(dǎo)致油庫系統(tǒng)發(fā)生靜電火災(zāi)爆炸事故[11]。除此之外，人體靜電也是油庫的一大安全隱患，易產(chǎn)生靜電的衣物、鞋子會因為人體運(yùn)動摩擦起電，甚至人體會與其它帶電體發(fā)生感應(yīng)帶電[12]。

在全面分析油庫區(qū)靜電事故誘發(fā)因素的基礎(chǔ)上，參考國內(nèi)外油庫區(qū)靜電安全評估標(biāo)準(zhǔn)，結(jié)合油庫區(qū)日常管理經(jīng)驗，綜合考慮影響油庫區(qū)靜電安全的因素，建立適用于油庫系統(tǒng)的靜電安全評價指標(biāo)體系。該體系確定靜電產(chǎn)生、可燃油氣和靜電放電3個一級指標(biāo)，選取了25個二級評價指標(biāo)，分別為油品易產(chǎn)生靜電、管壁粗糙、油品流速過快、油品過濾產(chǎn)生靜電、油水混合、鞋子與地面摩擦、化纖品與人體摩擦、油液內(nèi)雜質(zhì)、絕緣材質(zhì)容器裝油、油霧與空氣摩擦、油品撞擊罐壁、噴濺式裝卸油品、油品揮發(fā)、閥門泄漏、設(shè)備腐蝕漏油、油品混裝、通風(fēng)不良、裝卸油品后未靜置、空氣濕度低、作業(yè)中與導(dǎo)體接觸、取樣和測量操作錯誤、未靜電接地、靜電接地電阻不達(dá)標(biāo)、接地裝置損壞和設(shè)備漆層導(dǎo)電性差。

2 油庫區(qū)靜電安全的模糊綜合評價模型

針對油庫區(qū)靜電安全問題的多準(zhǔn)則、多目標(biāo)、模糊性的特點(diǎn)，聯(lián)用層次分析和模糊評價法建立多層次模糊綜合評價模型(圖1)。該模型的層次結(jié)構(gòu)分為3層：目標(biāo)層A，即油庫區(qū)靜電事故；準(zhǔn)則層B，即影響靜電事故發(fā)生的3大主要因素，靜電產(chǎn)生B1、可燃油氣B2和靜電放電B3；指標(biāo)層C，即靜電安全評價體系中選擇的25個二級指標(biāo)，依次分別記為C1、C2、C3、…、C25?？紤]到靜電產(chǎn)生、積累到放電過程難以定量區(qū)分，因此在所構(gòu)建的層次模型中，每個二級指標(biāo)不一定只作用于一個一級指標(biāo)，有的二級指標(biāo)會作用于多個一級指標(biāo)。例如，噴濺式裝卸油品，不僅會產(chǎn)生靜電，還會產(chǎn)生可燃的油氣混合物，甚至?xí)l(fā)生火花放電，因此該指標(biāo)同時作用于3個準(zhǔn)則層指標(biāo)。

2.1 建立評價指標(biāo)集合和評語集

因素集就是決策系統(tǒng)各評價指標(biāo)的集合，用U表示，U=[u1,u2,…,un]。在所建立的油庫區(qū)靜電安全多層次模糊評價模型中，因素集分為準(zhǔn)則層因素和指標(biāo)層因素兩層，準(zhǔn)則層因素即3個一級指標(biāo)；指標(biāo)層因素即25個二級指標(biāo)。

評價集是給予評價的專家學(xué)者或工作人員對所評價對象可能做出的各種總的評估結(jié)果所組成的集合，通常用V表示，V=[v1,v2,…,vm]，從v1到vm每個評估等級都對應(yīng)著1個模糊子集。實(shí)際安全評價中，評價等級數(shù)m一般取3到7中的整數(shù)。如果m取值較大，難以給出準(zhǔn)確分明的評估等級描述語言，且給出的判定等級歸屬有效性低；反之，又無法達(dá)到模糊評價的目的。m一般取奇數(shù)，具體的等級劃分可以根據(jù)評價對象自定義[13]。

2.2 確定指標(biāo)權(quán)重

復(fù)雜的多因素油庫區(qū)靜電安全評價問題往往梳理成多層次單目標(biāo)問題，通過構(gòu)建兩兩比較矩陣來定量判斷因素的相對重要性[14-15]。按照標(biāo)度法給予一定的分值，建立比較判斷矩陣；通過求判斷矩陣的最大特征值及其對應(yīng)的特征向量來確定單層次權(quán)重；對求取的最大特征值進(jìn)行一致性檢驗，確保所求權(quán)重值的正確性[16]。n個因素的判斷矩陣C一般可表示為

C=[cij]n×n。

(1)

式中：Ci表示處于同一層次上的隸屬于同一因素X的各指標(biāo)，cij矩陣中表示Ci與Cj相對于因素X的重要程度；i，j=1，2，…，n。求解矩陣的最大特征值和對應(yīng)特征向量的公式為

CW=λmaxW。

(2)

式中：λmax為判斷矩陣C的最大特征值；W為λmax對應(yīng)的特征向量，W=[W1,W2,…,Wn]，Wi代表其對應(yīng)元素的單權(quán)重。

通過歸一化處理判斷矩陣C中的列向量得到新矩陣，再計算新矩陣行向量的算術(shù)平均值，即為權(quán)重向量的近似值，最后求解最大特征值，亦即

(3)

(4)

在兩兩比較判斷的過程中，可能會出現(xiàn)所構(gòu)造的判斷矩陣并不一定具有一致性。為保證求得特征向量的正確性，判斷矩陣需要進(jìn)行一致性程度檢驗。一致性指標(biāo)為

(5)

式中：λmax為最大特征值，n為矩陣的階。引入平均隨機(jī)一致性指標(biāo)RI，則隨機(jī)一致性比率

CR=CI/RI。

(6)

CR作為衡量判斷矩陣一致性的標(biāo)準(zhǔn)，當(dāng)CR<0.1時，一致性檢驗即可通過。

計算綜合權(quán)重，綜合權(quán)重又稱為層次總排序，指的是同一層次中所有因素對目標(biāo)層的相對重要性標(biāo)度。從上而下逐層求解綜合權(quán)重，假設(shè)已經(jīng)計算出了第k-1層m個指標(biāo)相對于最高層的權(quán)重向量為

(7)

第k層上有n個因素，其對于k-1層的某個指標(biāo)j的單層權(quán)重向量為

(8)

(9)

2.3 確定評價隸屬度矩陣

通過單因素的模糊評價來確定隸屬度矩陣R，單因素模糊評價是指評判某個指標(biāo)因素對于評價集的獨(dú)立隸屬度。因素ui的單因素評判方法：給出ui對評價集合V中每一個評價等級的隸屬度，將ui對V中第j個元素的隸屬度記為rij，ui，則單因素評價模糊集合Ri可表示為

Ri=[ri1,ri2,…,rim]。

(10)

式中：i=1，2，…，n；j=1，2，…，m。Ri為評價集V上的模糊子集，rij為關(guān)于因素ui具有評價Vj的程度。

n個因子的單因素評價集組成矩陣，得到一個總的評價矩陣即為綜合模糊評價矩陣R，可表示為

(11)

2.4 確定模糊綜合評判結(jié)果

單因素的模糊評價反映的是某個因素對評價指標(biāo)的影響效果，想要得到全面的評價結(jié)果，需要綜合考慮所有因素對評價目標(biāo)的影響，即模糊綜合評價。確定了指標(biāo)因子的權(quán)重集合W和綜合模糊評判矩陣R之后，利用R作模糊線性變換，將W轉(zhuǎn)化成評價集V上的模糊子集B，則代表系統(tǒng)評判集諸因素的相對權(quán)重的集合為

B=W*R=[b1,b2,…,bm]。

(12)

式中，*代表廣義模糊矩陣的合成運(yùn)算。針對油庫區(qū)靜電安全評價體系的特點(diǎn)，選用取小上界和型的評價模型M(·,?)來處理綜合隸屬度運(yùn)算，即

(13)

根據(jù)所建立的油庫區(qū)靜電安全評價層次結(jié)構(gòu)模型，需要進(jìn)行兩級模糊綜合評判，最終求解出綜合隸屬度B=[b1,b2,…,bm]；再按照最大隸屬度原則，確定出模糊綜合評價集B中最大的元素bj所對應(yīng)的評語vj，從而得出模糊綜合評價結(jié)果。

3 實(shí)例分析

以某石油公司油庫為例，采用所建立的模糊綜合評價模型和方法分析評估靜電風(fēng)險。該油庫裝備20臺大型設(shè)備，其中輸油泵機(jī)組12套，輸油管線12條，輕油儲罐、重油儲罐兩組儲罐共12個，儲油能力3×104m3，儲存油品主要有柴油、汽油、燃料油等。油庫系統(tǒng)中還裝備了輸油泵監(jiān)測儀、儲油罐液位計、在線石油含水分析儀以及管線監(jiān)測系統(tǒng)等輔助設(shè)備。整個油庫系統(tǒng)是根據(jù)“石油庫設(shè)計規(guī)范”(GBJ 74-84)建設(shè)的，按照“石油化工靜電接地設(shè)計規(guī)范”(SH3097-2000)和“液體石油產(chǎn)品靜電安全規(guī)程”(GB 13348-92)采取了相應(yīng)的靜電預(yù)防措施。調(diào)研時，該油庫處于正常運(yùn)營狀態(tài)。

3.1 評價指標(biāo)權(quán)重的計算分析

根據(jù)圖1所示層次結(jié)構(gòu)模型，準(zhǔn)則層對目標(biāo)層的判斷矩陣記為矩陣A，其特征向量記為W；指標(biāo)層對“靜電產(chǎn)生”、“可燃油氣”、“靜電放電”這3個準(zhǔn)則的判斷矩陣及其特征向量分別記為B1和W1、B2和W2、B3和W3；指標(biāo)層對目標(biāo)層的綜合權(quán)重記為WZ。兩兩比較各層因素的相對重要性，采用標(biāo)度法賦值給出的判斷矩陣A及B1、B2、B3分別見表1—表4，特征向量的計算結(jié)果分別為：

W=[0.25，0.25，0.50]T；

(14)

W1=[0.078，0.129，0.203，0.129，0.078，0.021，0.030，0.030，0.048，0.048，0.078，0.129]T；

(15)

W2=[0.043，0.130，0.216，0.074，0.074，0.074，0.043，0.216，0.130]T；

(16)

W3=[0.078，0.078，0.047，0.133，0.029，0.133，0.133，0.213，0.029，0.078，0.047]T。

(17)

根據(jù)已計算出的單層排序和式(9)，可計算出指標(biāo)層對目標(biāo)層的綜合權(quán)重

WZ=(0.019，0.032，0.051，0.032，0.019，

0.005，0.008，0.008，0.051，0.022，0.091，0.110，0.018，0.018，0.018，0.011，0.054，0.099，0.015，0.067，0.067，0.107，0.015，0.039，0.023)T。

(18)

由此可得出評價指標(biāo)層對目標(biāo)層的綜合權(quán)重排序(表5)。 因此，噴濺式裝卸油品C12、未靜電接地C22、裝卸油品后未靜置C18、油品撞擊罐壁C11這4個因素的權(quán)重皆大于0.09，對該油庫區(qū)的靜電安全具有非常重要的影響。

表1 準(zhǔn)則層對目標(biāo)層A的判斷矩陣Tab.1 Judgment matrix of criteria layer to target layer A

表2 指標(biāo)層對準(zhǔn)則層因子B1的判斷矩陣Tab.2 Judgment matrix of index layer to criteria layer factor B1

表3 指標(biāo)層對準(zhǔn)則層因子B2的判斷矩陣Tab.3 Judgment matrix of index layer to criteria layer factor B2

表4 指標(biāo)層對準(zhǔn)則層因子B3的判斷矩陣Tab.4 Judgment matrix of index layer to criteria layer factor B3

表5 指標(biāo)層因素的綜合權(quán)重排序和隸屬度估計Tab.5 Comprehensive weight and fuzzy membership degree estimation of index layer factors

3.2 模糊評價的隸屬度矩陣

采取安全、較安全、輕微危險、較高危險、高危險5個等級的評價方式建立評語集V=[v1,v2,v3,v4,v5]，采用專家打分法來確定隸屬度矩陣。邀請10位專家對評價模型中各指標(biāo)對應(yīng)各個評價等級的隸屬度給出估計值，分值取0～1，且每個指標(biāo)對應(yīng)各個評價等級的隸屬度之和為1。對10位專家給出的估計值進(jìn)行統(tǒng)計，取平均值作為各指標(biāo)的最終隸屬度估計值(表5)。

3.3 模糊綜合評價結(jié)果

由表5中C1～C12的隸屬度行組成評價矩陣R1，即為B1準(zhǔn)則的綜合隸屬度矩陣；C10～C18的隸屬度行組成評價矩陣R2，即為B2準(zhǔn)則的綜合隸屬度矩陣；C9、C11、C12、C18—C25的隸屬度行組成評價矩陣R3，即為B3準(zhǔn)則的綜合隸屬度矩陣。式(15)—(17)的權(quán)重向量W1、W2、W3分別對應(yīng)地轉(zhuǎn)化成評價集V上的模糊子集B1、B2、B3。由式(13)則有：

B1=W1*R1=[0.1300.4900.2420.0940.044] ；

(19)

B2=W2*R2=[0.1910.4570.2270.0890.036] ；

(20)

B3=W3*R3=[0.1660.5680.1660.0770.023]。

(21)

(22)

根據(jù)式(13)、式(14)和式(22)，計算出的最終模糊評價隸屬度為

B=W*R=[0.1630.5210.2000.0850.031]。

(23)

由式(23)可知，參照最大隸屬度原則，可知該油庫為較安全等級，這與該油庫的實(shí)際情況一致。盡管該油庫現(xiàn)階段可以正常運(yùn)行，但是存在一定的靜電安全隱患。需要加強(qiáng)管理，根據(jù)表5中的影響因素分析，采取對應(yīng)措施提高預(yù)防靜電事故的安全防范力度。

4 結(jié) 論

(1)采用層次分析和模糊分析相結(jié)合的方法，建立了油庫區(qū)靜電安全的多層次模糊綜合評價模型。該模型的評價指標(biāo)體系由3個一級指標(biāo)和25個二級指標(biāo)組成。

(2)應(yīng)用所建立的多層次模糊綜合評價模型，對一個實(shí)例油庫進(jìn)行了靜電安全評價，靜電安全等級的評定結(jié)果為較安全，即危險性一般，與實(shí)際情況相符。

(3)實(shí)例分析表明，所建立的油庫區(qū)靜電安全評價模型能夠有效地評價油庫區(qū)靜電安全的現(xiàn)狀，可解決靜電安全評價的復(fù)雜性、多層次性和不確定性難題。