魏沁汝 ,李 杉 ,周永淳,韋孟佼
1.西南石油大學(xué)石油與天然氣工程學(xué)院 (四川 成都 610500)
2.中國石油西南油氣田川西北氣礦 (四川 邛崍 611500)
3.中國石化西南油氣分公司川西采氣廠 (四川 德陽 618000)
4.中國石油西南油氣田分公司輸氣管理處 (四川 成都 610213)
管道業(yè)飛速發(fā)展,相應(yīng)的事故也悄然而來。如今輸油管道趨向大管徑、高壓力,各種泄漏引發(fā)的爆炸、火災(zāi)、污染等連鎖事故逐漸引起人們高度關(guān)注。2013年11月22日10時25分,位于山東省青島經(jīng)濟技術(shù)開發(fā)區(qū)的中石化股份有限公司,下屬管道儲運分公司東黃輸油管道原油泄漏并發(fā)生爆炸,造成62人死亡、136人受傷,直接經(jīng)濟損失7.5億元[1]。2012年10月3日,印度尼西亞南蘇門答臘省一輸油管道發(fā)生爆炸并引發(fā)大火[2],事故造成5人死亡,18人受傷。 這一系列人員傷亡、財產(chǎn)損失、環(huán)境污染事故均是由多米諾效應(yīng)引起的,由此可見多米諾效應(yīng)研究是非常有必要的。
多米諾效應(yīng)是指事故的擴大和連鎖反應(yīng)。一旦發(fā)生多米諾效應(yīng),事故種類增加,危險等級提高,影響范圍也將迅速擴大。歐盟1982年和1996年兩度頒布塞維索法令[3-4],要求對工廠布局、土地使用方面考慮多米諾效應(yīng)。隨即歐美各國開始研究多米諾效應(yīng),但大多針對工廠與廠房。2000年以后我國開始研究多米諾效應(yīng),目前研究尚處于初級階段。胡小榮、龔時雨[5]等人分析了多米諾效應(yīng)初始事件和事故場景;劉茂[6]等人對多米諾效應(yīng)的風(fēng)險分析方法進行了探索;戴雪松,王艷華[7]等人對化工重大事故多米諾效應(yīng)做了相關(guān)研究;在城市燃氣管網(wǎng)事故多米諾效應(yīng)方面劉艷華[8]進行了相關(guān)分析;張青松[9]對油庫池火災(zāi)場景進行了研究;師立晨[10]等對重大危險源多米諾效應(yīng)進行了分析,并探究了熱輻射、超壓以及碎片引發(fā)的多米諾效應(yīng)閥值。在石油化工行業(yè),多米諾效應(yīng)多針對儲罐等化工園區(qū),并無輸油管道相關(guān)探索。且傳統(tǒng)的多米諾效應(yīng)分析立足于危險設(shè)備,運用火災(zāi)爆炸理論分析其發(fā)生途徑和后果,相對單一。所以統(tǒng)計典型輸油管道泄漏事故,判定事故多米諾效應(yīng),分析多米諾效應(yīng)產(chǎn)生的不同途徑與后果,有利于輸油管道設(shè)計選線,且為輸油管道泄漏事故處理提供理論指導(dǎo)。
近年來由于腐蝕破裂、疲勞損壞和第三方破壞(特別是打孔盜油)導(dǎo)致輸油管道泄漏事故時有發(fā)生,有些泄漏在管道工作人員及時發(fā)現(xiàn)情況下得到控制,有些泄漏事故則產(chǎn)生多米諾效應(yīng),造成不堪設(shè)想的后果。為總結(jié)輸油管道泄漏事故后果規(guī)律,搜集2000~2013年國內(nèi)外主要的輸油管道泄漏事故51起,列出事故主要后果,見表1。
由于上述51起事故均為較嚴重的輸油管道泄漏事故,有2個共同點:財產(chǎn)損失和管道停輸。除此之外,還主要造成了污染河流海域、污染土地(農(nóng)田、耕地等)、人員傷亡、影響居民用水、以及其他特殊情況如破壞動植物生存環(huán)境,旅游業(yè)、捕魚業(yè)受損等事故后果,事故后果所占比例如圖1所示。
表1 國內(nèi)外主要輸油管道泄漏事故
續(xù)表1
注:不包括財產(chǎn)損失、管道停輸。
從圖1中可以看出成品油管道泄漏事故主要產(chǎn)生了6種事故后果,其中人員傷亡、污染土地、污染河流海域情況最多,分別為26%、26%、28%。2種組合后果污染河流和人員傷亡、污染河流和影響居民用水則分別占4%,6%。單一影響到居民水源情況相對較少。一些特殊情況占8%。
根據(jù)Cozzani[46]等人給出的事故多米諾效應(yīng)3個存在條件:存在引發(fā)多米諾效應(yīng)的初始事故;事故蔓延開來必須導(dǎo)致至少1個二次事故發(fā)生;初始事故發(fā)生1個或者多個二次事故,產(chǎn)生了大于初始事故后果的嚴重事故。51個事故存在觸發(fā)多米諾效應(yīng)的初始事件即輸油管道泄漏,初始事件泄漏的蔓延觸發(fā)周圍環(huán)境發(fā)生二次以上事件,如污染、爆炸、甚至海邊旅游和漁業(yè)受損等,二次以上事件產(chǎn)生擴大了初始事件后果,所以判定51個事件產(chǎn)生了多米諾效應(yīng)。例如,2011年9月12日肯尼亞輸油管道事故,首都內(nèi)羅畢東部工業(yè)區(qū)的燃油庫一條地下輸油管線發(fā)生破裂泄漏(初始事件產(chǎn)生),隨即泄漏原油擴散流入居民區(qū)(初始事件蔓延),民眾搶油,有人吸煙(觸發(fā)臨界條件)引發(fā)火災(zāi)爆炸(二次事故),導(dǎo)致至少82人死亡,118人受傷 (泄漏事故后果擴大)。傷亡事故的產(chǎn)生加大了泄漏事故后果,滿足輸油管道多米諾效應(yīng)的3個條件,所以判定肯尼亞輸油管道泄漏事故產(chǎn)生了多米諾效應(yīng)。其余事故依次類推。
由于多米諾效應(yīng)是一個初始事件引發(fā)多級事件的連鎖反應(yīng),故運用事件樹來分析最為合理。事件樹分析法(Event Tree Analysis,ETA)[47]是一種邏輯的演繹法,在給定一個初因的事件后,分析初因事件可能導(dǎo)致各種不同的事件序列結(jié)果,可以定量和定性地評價系統(tǒng)特性。
根據(jù)統(tǒng)計分析得知,51個輸油管道事故產(chǎn)生了多米諾效應(yīng),其地點多發(fā)生在特殊自然環(huán)境如海洋河流、農(nóng)田、以及一些特殊人類活動場所;多米諾效應(yīng)事故種類有污染、火災(zāi)爆炸、停止人類生產(chǎn)活動等;多米諾效應(yīng)后果為人員傷亡、影響居民用水、農(nóng)田污染、財產(chǎn)損失等。由此建立輸油管道泄漏事故多米諾效應(yīng)事件樹模型,事件模型分為海底管道與陸地管道。
1)模型一:海底管道泄漏事故多米諾效應(yīng)事件樹模型。根據(jù)上文輸油管道泄漏事件統(tǒng)計分析結(jié)果,海底管道泄漏引發(fā)的事故后果有污染周圍海域,影響海洋生物的生存,污染海灘,影響當?shù)芈糜螛I(yè)、捕魚業(yè)等相關(guān)產(chǎn)業(yè)。由此結(jié)果建立海底管泄漏多米諾效應(yīng)事件樹模型,如圖2所示。
由圖2知,海底管道泄漏產(chǎn)生了3中不同結(jié)果,此處分析最復(fù)雜的第一種結(jié)果。由于海底管道泄漏產(chǎn)生1次事故,污染附近海域;接著泄漏擴散成功,觸發(fā)海水水質(zhì)變化產(chǎn)生2次事故,污染面積擴大,破壞海洋生物生存環(huán)境,部分海洋生物生存瀕危,如泰國海底管道泄漏事故中一些魚類和珊瑚的死亡;當擴散蔓延至海邊時,污染面積急速增加,魚類等生物死亡,有旅游勝地或者捕魚業(yè)受損,最終影響當?shù)亟?jīng)濟,使泄漏后果擴大,多米諾效應(yīng)產(chǎn)生。
圖2 海底輸油管道泄漏事故多米諾效應(yīng)事件樹模型
2)模型二:陸地管道泄漏事故多米諾效應(yīng)事件樹模型。根據(jù)輸油管道泄漏事件統(tǒng)計分析結(jié)果,陸地管道引發(fā)的事故后果有人員傷亡,鐵路停運,污染農(nóng)田、河流、水庫,影響居民飲用水,管道停輸,財產(chǎn)損失。陸地管道較海底管道面臨的周圍環(huán)境情況更多,如農(nóng)田、河流、公路、濕地、工廠等。綜合泄漏后果與管道周邊情況,建立陸地管道泄漏事故多米諾效應(yīng)事件樹如圖3所示。
此處研究的事故都是比較嚴重的泄漏事故,由圖3和表2可以看出當初始事件管道泄漏A1發(fā)生,在不產(chǎn)生多米諾效應(yīng)情況下則僅會產(chǎn)生財產(chǎn)損失、管道停輸S6。當發(fā)生多米諾效應(yīng)則可能產(chǎn)生5種情況造成更加嚴重的后果(S1至S5)污染土地、污染河流、影響居民用水、人員傷亡和影響人類其他活動。
圖3 陸地輸油管道泄漏事故多米諾效應(yīng)事件樹模型
表2 陸地輸油管道泄漏事故多米諾效應(yīng)事件樹因素表
對于陸地輸油管道泄漏產(chǎn)生的6個后果 (S1至S6),由于本文研究比較嚴重泄漏事故,所以財產(chǎn)損失、管道停輸是必定產(chǎn)生的,S6發(fā)生的概率為1。其余5個后果,準確確定其概率比較困難。
例如某輸油管道管徑508mm,輸送介質(zhì)為汽油,流量為450m3/h,由于第三方施工破壞發(fā)生泄漏,10min后發(fā)現(xiàn)并采取相應(yīng)措施。泄漏段左側(cè)200m有飲用水庫,右側(cè)200m有河流。
由于輸油管道周圍存在河流、水庫,根據(jù)圖3,泄漏輸油管道由于擴散作用,容易引起污染河流、影響居民飲用水的多米諾效應(yīng)。此處以污染河流為例,計算其發(fā)生的概率,計算步驟如下:
1)選擇專家組,獲得專家對污染河流的判斷。選取10名專家,評估意見為可能性大3票、中等3票、較大4票。
2)將自然語言轉(zhuǎn)化為模糊數(shù)。根據(jù)自然語言的模糊數(shù)[48]表達形式如圖4,得到3種判斷語言的隸屬函數(shù)式(3)~(5),由此可將專家判斷語言轉(zhuǎn)為模糊數(shù)。
圖4 失效可能性等級與模糊數(shù)對應(yīng)關(guān)系
其中下角標H、M、FH分別表示大、中等、較大。
綜合處理10位專家意見,選用模糊截集與平均算法組合專家的不同意見,設(shè)(3),(4),(5)式的 α截集分別為 Hα=[H1,H2],Mα=[m1,m2],F(xiàn)Hα=[h1,h2],其中 H1,H2,m1,m2,h1,h2分別為(3),(4),(5)的 α 截集的上限和下限。對fH(x)而言,令α=(x-0.1)/0.1,則H1=0.1α+0.7,H2=0.9-0.1α。 同理可得 m1=0.1α+0.4,m2=0.6-0.1α;h1=0.1α+0.5,h2=0.8-0.1α;。 在 α 的截集下,10位專家意見的總模糊數(shù)Z為:
3)將模糊數(shù)轉(zhuǎn)化為模糊可能性值。由于需要后果發(fā)生的概率大小,所以必須先將模糊數(shù)轉(zhuǎn)換為模糊可能性值(FPS),再將模糊可能性值轉(zhuǎn)換為模糊概率。根據(jù)Cheng和Hwang[49]提出的左右模糊排序法將模糊數(shù)轉(zhuǎn)化為模糊可能性值。該方法定義的最小模糊集和最大模糊集為:
則模糊數(shù)W的左右模糊可能性值分別為:
則可得到模糊數(shù)W的FFR=1.23×10-2,即污染河流發(fā)生概率為1.23×10-2。同理可得影響居民飲用水的概率為4.39×10-3。污染河流發(fā)生的概率大于影響居民用水,表明在發(fā)生泄漏時,河流污染可能會比污染水庫提前發(fā)生,因此需要調(diào)節(jié)搶險措施。同時,此結(jié)果與圖1事故后果比例圖相符,說明此方法具有一定適用性。對于泄漏的其他情況,同理可以將實際管道泄漏情況與圖2圖3結(jié)合,請求專家評判,再由模糊理論獲得事故發(fā)生的概率,由此可以指導(dǎo)類似輸油管道泄漏事故的處理工作。
1)輸油管道泄漏事故多米諾效應(yīng)有4種主要后果以及2種組合后果分別是:人員傷亡、污染河流海域、污染土地、影響居民用水、污染河流和人員傷亡、污染河流和影響居民用水。其中人員傷亡、污染河流、污染土地最容易發(fā)生,所以在選線與事故處理時,應(yīng)多關(guān)注管道周圍人類活動場所、河流、土地的情況。
2)模糊理論與專家判斷法處理輸油管道事故多米諾效應(yīng)后果有一定適用性,當輸油管道泄漏發(fā)生,可以運用此法計算多米諾后果發(fā)生概率,從而起到對事故處理的指導(dǎo)作用。
[1]鮑丹.中石化青島輸油管道爆炸事故原因查明認定責(zé)任事故[N].人民日報,2014-01-10.
[2]周檬.印尼石油管道爆炸起火致5死18傷[N].新華社,2012-10-13.
[3]OJEC.Council directive 82/501/EEC of 24 June 1982 on the control of major Accident Hazards of Certain Industrial Activities[M].Brussels:Official Journal of the European Communities,1982.
[4]OJEC.Council directive 96/82/EC of 9 December 1996 on the control of major-accident hazards involving dangerous substances[M].Brussels:Official Journal of the European Communities,1997:97-107.
[5]胡小榮,龔時雨,曾士勇.基于危險的事故場景及形式化描述[J].工業(yè)安全與環(huán)保,2004,30(5):41-43.
[6]張永強,劉茂.多米諾效應(yīng)的風(fēng)險分析方法[J].安全與環(huán)境學(xué)報,2008,8(6):152-158.
[7]王艷華,戴雪松,鮮魚小東.化學(xué)工業(yè)重大事故的多米諾效應(yīng)分[J].中國安全科學(xué)學(xué)報,2008,18(5):129-136.
[8]劉艷華.基于多米諾效應(yīng)的城市燃氣管網(wǎng)事故后果研究[D].成都:西南石油大學(xué),2009.
[9]張青松,楊一楠,鄭薇.基于貝葉斯概率的油庫池火災(zāi)多米諾效應(yīng)分析[J].安全與環(huán)境學(xué)報,2013,13(3):239-242.
[10]師立晨,劉驥,魏利軍等.重大危險源多米諾效應(yīng)的后果分析[J].中國安全生產(chǎn)科技,2007,3(6):45-48.
[11]王豐豐.尼輸油管道爆炸燒死數(shù)百人[N].京華日報,2006-12-27.
[12]曾多禮.格拉成品油管追火災(zāi)事故分析及預(yù)防[J].油氣儲運,2004,23(8):40-41.
[13]帥健.油氣管道事故原因分析及分類方法研究[J].中國安全科學(xué)學(xué)報,2013,23(7):109-114.
[14]王小華.廣州石化輸油管漏油起火[N].新快報,2004-07-24.
[15]莊浩濱.山東臨邑至濟南輸油管線遭盜油[N].新華社,2006-08-10.
[16]黃培昭.尼日利亞發(fā)生輸油管道爆炸事故[N].人民日報,2006-05-13.
[17]史先振.尼日利亞油管大爆炸:漏的是汽油[N].新華每日電訊,2006-12-18.
[18]智利海底輸油管破裂大量石油泄漏[N].新華社,2007-06-01.
[19]朝鮮汽油管道事故已造成100多人死亡[N].俄羅斯新聞網(wǎng) , 2007-06-20.http://rusnews.cn/guojiyaowen/guoji_shehui/20070620/41813670-print.html.
[20]趙丹丹.山東輸油管遭盜竊破壞 近30噸原油泄漏[N].新聞晨報,2007-08-02.
[21]殷亮.波蘭石油管道泄漏造成維斯瓦河污染[N].國際在線,2007-12-11.
[22]李慧穎.廈門海滄碼頭輸油管泄漏[N].新京報,2008-01-16.
[23]魏向楠.尼日利亞輸油管道爆炸上百人死亡[N].石油商報,2008-05-21.
[24]李彪.近百噸柴油外泄污染10公里河水[N].華商報,2008-07-27.
[25]吳松濤.蘭成渝輸油管四川綿陽段被鉆破致柴油泄漏[N].四川新聞,2009-03-12.
[26]梁旭.俄羅斯輸油管發(fā)生爆炸 200噸石油泄漏[N].環(huán)球時報,2009-12-03.
[27]鐘智維.陝西華縣中石油輸油管洩漏形成污染帶進入黃河[N].香港商報,2010-01-04.
[28]張芳玲.中石化魯寧線淮河大橋段發(fā)生輸油管道泄漏[N].新華社,2010-04-15.
[29]尚前名,王仁貴.七萬公里油氣管道誰來護衛(wèi)[J].瞭望,2010(19):12-14.
[30]美國埃索石油公司輸油管在英國海域發(fā)生泄漏[N].中國新聞,2010-06-22.
[31]謝萍.陜西吳起輸油管道泄漏原油3至5噸 下游未受污染[N].中國新聞,2010-08-13.
[32]戚燕凌.土耳其東南部油管爆炸致兩人死亡[N].新華社,2010-08-11.
[33]唐葵陽.遼中縣一條石油管線發(fā)生泄漏 玉米地里飄起黑油[N].沈陽晚報,2010-08-28.
[34]諶融.墨西哥輸油管道被竊賊割破 逾5萬升柴油泄漏[N].中國日報,2011-02-10.
[35]王軍.美國蒙大拿州輸油管泄漏河流受污染[N].新華社,2011-07-03.
[36]楊柳.美國通往加拿大管道泄露 80萬加侖石油泄入河流[N].中國新聞,2011-07-28.
[37]落落.肯尼亞輸油管線爆炸的慘痛教訓(xùn)[J].廣東安全生產(chǎn),2011,19:50.
[38]楊彥宇.哥倫比亞一輸油管道爆炸 致百余人死傷[N].中國新聞,2012-12-24.
[39]懷若谷.“秦京輸油管道泄漏”追蹤[N].京華時報,2012-03-30.
[40]張海鵬.銅川村民取土致輸油管破裂 17噸原油泄漏未入河道[N].中國日報,2012-04-15.
[41]王歡.美油管破裂再陷“漏油門”居民撤離未構(gòu)成即時威脅[N].中國新聞,2013-04-01.
[42]長慶采油一廠輸油管線破裂原油外泄 王瑤水庫水質(zhì)遭受污染[N].延安日報,2013-07-15.
[43]王海青.陜西子長山體滑坡致輸油管線破裂泄漏約10噸原油[N].中新社,2012-08-14.
[44]泰國國家石油全球化工公司羅勇海底輸油管洩漏 [N].泰國世界日報,2013-08-29.
[45]哥倫比亞最大輸油管受爆炸襲擊 發(fā)生原油泄漏[N].法新社,2013-07-30.
[46]Cozzani V,Antonioni G,Spadom G.Quantitative assessment of domino scenarios by a GIS-based software tool[J].Journal of Loss Prevention in the Process Industries,2006,19(5):463-477.
[47]劉雙躍.安全評價[M].北京:冶金工業(yè)出版社2010.
[48]余樹榮,楊慧來.基于梯形模糊數(shù)的長輸管道故障樹分析.蘭州理工大學(xué)學(xué)報[J].2009,35(5):62-65.
[49]Shu Jen Chen,Ching Lai Hwang,et al.Fuzzy multiple attribute decision making methods and applications[M].Germany SpringerVerlag Berlin Herdelberg,1990.