摘　要：客觀有效的進(jìn)行城市燃?xì)夤艿里L(fēng)險評估對于燃?xì)馄髽I(yè)的安全運(yùn)行和科學(xué)管理具有重要的意義，采用定量風(fēng)險評價方法，利用收集的歷史數(shù)據(jù)分析燃?xì)夤艿朗Ц怕?、失效后造成的?jīng)濟(jì)損失，使用雙曲型理論及回歸分析法建立城市燃?xì)夤艿离p曲型風(fēng)險決策模型來判斷燃?xì)夤艿赖娘L(fēng)險可接受性，并根據(jù)GIS4支術(shù)，以SUPERMAP為開發(fā)平臺，將燃?xì)夤艿里L(fēng)險水平進(jìn)行可視化表達(dá)。在實例中應(yīng)用雙曲型風(fēng)險決策模型和GIS技術(shù)，結(jié)果表明：在管道風(fēng)險評估中，該模型與GIS的結(jié)合能夠客觀的反應(yīng)燃?xì)夤艿赖娘L(fēng)險水平，對城市燃?xì)夤艿腊踩\(yùn)行提供科學(xué)決策依據(jù)。

關(guān)鍵詞：燃?xì)夤艿馈』貧w分析　雙曲型風(fēng)險決策模型　地理信息系統(tǒng)　可視化表達(dá)

City Gas Pipeline Risk Assessment Research Based on Hyperbolic Risk Decision Model and GIS

AbStract：To achieve the safe operation and scientific m8nagement，it is very important for the gas enterprise to do the objective risk assessment of city gas pipeline．The quantitative risk assessment method is adopted．Frirstly．analyse the gas pipeline failure probability and the economic loss by using the historical data；Next，establish the city gas pipeline hyperbolic risk decision model to judge the gas pipeline risk acceptability by using the hyperbolic theory and regression analysis method；Lastly，the gas pipeline risk level is visually expressed by using GIS technology and taking SUPERMAP as development platform．Examples show that the hyperbolic risk decision model and GIS technology applied to pipeline risk assessment can objectively reflect the gas pipeline risk and provide scientific decision for the safe operation of city gas pipeline．

Key words：gas pipeline  regression analysis  city gas pipeline hyperbolic risk decision model  GIS visualization

1　引言

隨著城市化進(jìn)程的不斷加快，燃?xì)夤?yīng)成為城市發(fā)展中非常重要的組成部分，由于燃?xì)夤艿赖拇罅柯裨O(shè)和長期使用，燃?xì)夤艿劳话l(fā)事故頻繁發(fā)生給人們的生活帶來極大的風(fēng)險。城市燃?xì)夤艿里L(fēng)險評價是燃?xì)馄髽I(yè)進(jìn)行安全生產(chǎn)和安全管理的一種新型工具。由于我國燃?xì)夤艿涝谠O(shè)備和管理水平上同圍外還存在一定的差距，所以不能把國外成熟的管道風(fēng)險管理技術(shù)直接應(yīng)用于我國燃?xì)夤艿赖娘L(fēng)險管理中。國內(nèi)學(xué)者已經(jīng)開發(fā)了多種風(fēng)險評價方法^[1]，如黃小美利用故障樹和事件樹相結(jié)合對城市燃?xì)夤艿老到y(tǒng)失效概率進(jìn)行分析，并在其碩士畢業(yè)論文中對城市燃?xì)夤艿里L(fēng)險進(jìn)行了評價研究^[2]，隋楠使用對KENT評分法的指標(biāo)體系加以修正的方法對城市燃?xì)夤艿肋M(jìn)行風(fēng)險評價研究^[3]。以上學(xué)者主要采用定性與定量相結(jié)合的風(fēng)險評價方法，主觀因素影響較大。目前依據(jù)已有歷史數(shù)據(jù)進(jìn)行定量風(fēng)險評價的方法較少。基于此點，本文將雙曲型風(fēng)險決策模型應(yīng)用于城市燃?xì)夤艿里L(fēng)險評估中，并將風(fēng)險評估結(jié)果在SUPERMAP中可視化表達(dá)，較好的克服了主觀因素影響風(fēng)險評價結(jié)果準(zhǔn)確性的問題，從客觀上真實的反應(yīng)燃?xì)夤艿赖娘L(fēng)險水平，提高了燃?xì)夤艿赖娘L(fēng)險管理效率。

2　雙曲型風(fēng)險決策模型

對于燃?xì)夤艿朗鹿实娘L(fēng)險評價方法有多種，主要有定性風(fēng)險評價法、半定量風(fēng)險評價法和定量風(fēng)險評價法。定量風(fēng)險評價方法是在將與風(fēng)險相關(guān)的參數(shù)進(jìn)行量化的基礎(chǔ)上進(jìn)行的風(fēng)險評價，能夠?qū)ξｋU源或系統(tǒng)的風(fēng)險程度進(jìn)行定量描述，與定性風(fēng)險評價相比，能夠?qū)崿F(xiàn)多層次描述危險源或系統(tǒng)的風(fēng)險程度，以支持安全管理決策^[4]。雙曲型風(fēng)險決策模型就是定量風(fēng)險評價方法中的一種。

風(fēng)險就是危險發(fā)生的概率以及危險一旦發(fā)生所造成的后果的綜合，因此風(fēng)險是由兩部分組成：一是危險事件發(fā)生的概率。二是一旦出現(xiàn)危險，其所造成的后果嚴(yán)重程度和損失的大小。兩者的綜合結(jié)果，即風(fēng)險函數(shù)：

R=f(P．L) 　　　　　　　　   (1)

式中：R表示風(fēng)險值，P表示事故發(fā)生的概率；L表示事故后果嚴(yán)重程度。該函數(shù)關(guān)系一般取為概率與后果的乘積：R=P×L。

根據(jù)風(fēng)險函數(shù)，雙曲型風(fēng)險決策模型的建立分為以下幾個步驟。

2．1　燃?xì)夤艿朗鹿拾l(fā)生的概率分析

按照W．KentMuhlbauer的分類方法，造成管道失效和事故的原因共4大類：①第三方破壞。第三方破壞原因是指燃?xì)夤艿酪蜃钚÷裆?、地面上的活動狀況(水平)、當(dāng)?shù)鼐用竦乃刭|(zhì)(公眾教育)、管道地上設(shè)備安全、線路狀況、巡線頻率等因素，這些因素以外力擠壓或人為破壞的形式使燃?xì)夤艿朗艿綋p壞。②管道腐蝕。腐蝕是導(dǎo)致燃?xì)夤艿?span lang="EN-US">(主要是鋼管)穿孔、破裂的重要破壞因素。對于埋地管道而言，腐蝕來自兩個方面：內(nèi)腐蝕和外腐蝕。③人工誤操作。人工誤操作包括設(shè)計誤操作、施工誤操作、技術(shù)誤操作、管理與監(jiān)督誤操作、運(yùn)行維護(hù)與檢測誤操作等，這些誤操作會導(dǎo)致燃?xì)夤芫W(wǎng)運(yùn)行的安全性降低，風(fēng)險增加。④管網(wǎng)設(shè)計。設(shè)計需考慮的因素有管道最大承受壓力、管內(nèi)水擊因素、設(shè)計厚度和選材選型、士層移動、線路選擇、施工方案等，考慮不周會導(dǎo)致燃?xì)馐鹿实陌l(fā)生。在實際風(fēng)險管理中，用普通的概率計算方法計算其風(fēng)險概率相當(dāng)困難，有些還不可能實現(xiàn)。因此，采用美國核管理委員會(NUREG)推薦的方法，通過計算失效速率得出風(fēng)險概率。由分析統(tǒng)計的相關(guān)資料可知，管道發(fā)生事故的風(fēng)險概率一般符合U型曲線(見圖l)，即初期處于磨合階段，管道工作并不穩(wěn)定，后期接近設(shè)計壽命，事故發(fā)生也較頻繁。若只考慮中間的穩(wěn)定期，則其風(fēng)險速率v可視為常數(shù)。

管道事故的風(fēng)險概率按下式計算：

在上式中，P表示燃?xì)夤艿朗鹿拾l(fā)生的概率，x表示長L的管道上在使用時間t內(nèi)發(fā)生的事故總次數(shù)，v表示風(fēng)險速率，單位為次／(km*a)，t表示使用時間，單位為a，L表示管道長度，單位為km。根據(jù)上述方法可以先統(tǒng)計各因素在某段時間(如一年內(nèi))某段長度的管道上所引起的事故總次數(shù)，然后計算出風(fēng)險速率，最后求出燃?xì)夤艿朗鹿拾l(fā)生的概率^[5]。這樣從客觀的角度上反應(yīng)出燃?xì)夤艿朗鹿拾l(fā)生的可能性。

2．2　燃?xì)夤艿朗鹿拾l(fā)生的后果分析

根據(jù)風(fēng)險函數(shù)可知，風(fēng)險值除了要求計算事故發(fā)生的概率，還要對燃?xì)夤艿朗鹿拾l(fā)生所造成的后果進(jìn)行分析。所謂后果分析就是定量評價管道失效所導(dǎo)致的危害事件組合的嚴(yán)重程度，主要用貨幣化的直接經(jīng)濟(jì)損失表達(dá)這種嚴(yán)重程度嘲。燃?xì)夤艿朗鹿试斐傻氖Ш蠊话憧梢詮膬煞矫孢M(jìn)行計算：①直接經(jīng)濟(jì)損失，主要包括人員傷亡、財物損毀所造成的經(jīng)濟(jì)損失。②間接經(jīng)濟(jì)損失，主要包括漏氣經(jīng)濟(jì)損失以及系統(tǒng)恢復(fù)的費(fèi)用。根據(jù)下列公式計算每次已發(fā)生的燃?xì)夤艿朗鹿仕斐傻慕?jīng)濟(jì)損失總額。

式中L表示經(jīng)濟(jì)損失總額，Lda表示卣接經(jīng)濟(jì)損失，Lib表示間接經(jīng)濟(jì)損失。在式4中人員傷亡等經(jīng)濟(jì)損失要根據(jù)不同公司、不同地區(qū)、不同情況而定，應(yīng)具體問題具體分析。

2．3　建立雙曲型風(fēng)險決策模型

通過上述步驟可得知，燃?xì)夤艿涝诿恳灰?guī)定時間段內(nèi)(如1年)發(fā)生事故的概率以及事故發(fā)生后所造成的經(jīng)濟(jì)損失，使用回歸分析法對數(shù)據(jù)進(jìn)行擬合。

回歸分析就是從某一組數(shù)據(jù)出發(fā)，判斷變量之間是否存在相關(guān)關(guān)系，確定它們之間的數(shù)學(xué)表達(dá)式(即回歸方程)，并對其可信程度作統(tǒng)計檢驗，然后利用確定的回歸方程進(jìn)行預(yù)測^[7-9]。

由風(fēng)險計算公式R=pxl可知燃?xì)夤艿朗鹿拾l(fā)生的概率和可能造成的經(jīng)濟(jì)損失呈倒冪數(shù)關(guān)系，因此回歸方程即為R=pxl。用回歸分析法進(jìn)行相關(guān)分析，然后計算預(yù)測誤差，最后確定燃?xì)夤艿赖娘L(fēng)險值，即可得到該城市的風(fēng)險曲線，每段燃?xì)夤艿赖娘L(fēng)險值與此比較，判斷其風(fēng)險可接受性，即可建立雙曲型風(fēng)險決策模型如下圖2所示，其中P表示燃?xì)夤艿朗鹿拾l(fā)生的概率，L表示燃?xì)夤艿朗鹿拾l(fā)生造成的損失。

通過圖2可知，在同一條曲線上的風(fēng)險值是相同的，曲線越是靠外表示風(fēng)險值越大，而且R1>R2>R3。圖2中燃?xì)夤艿赖娘L(fēng)險預(yù)測值R2表示可接受的風(fēng)險水平閾值，若風(fēng)險預(yù)測值為R3，則表示風(fēng)險可接受；若風(fēng)險預(yù)測值為R1則表示不可接受，需要有關(guān)部門采取相關(guān)措施對風(fēng)險進(jìn)行有效適當(dāng)?shù)目刂?。通過雙曲型風(fēng)險決策模型的相關(guān)分析可以為燃?xì)夤艿里L(fēng)險決策提供有效的依據(jù)。

3　GIS住燃?xì)夤艿里L(fēng)險評估中的應(yīng)用

G1S的數(shù)據(jù)庫是以圖形為基礎(chǔ)，根據(jù)不同的要求有不同的數(shù)據(jù)結(jié)構(gòu)^[10]。本文采用supermap組建來實現(xiàn)燃?xì)夤艿里L(fēng)險水平的可視化。supermap是一個標(biāo)準(zhǔn)的控件GIS產(chǎn)品，可以在多種語言上使用。Supermap中deskpro是一套運(yùn)行在桌面端的專業(yè)GIS軟件，提供了地圖編輯、屬性數(shù)據(jù)管理、分析與決策輔助事務(wù)處理、地圖輸出、報表打印、三維建模等方面的功能。supermap is.net是新一代網(wǎng)絡(luò)地理信息系統(tǒng)開發(fā)平臺，supemlap is.net由客戶端用戶界面表現(xiàn)組件、web服務(wù)器擴(kuò)展、GIS服務(wù)器、數(shù)據(jù)服務(wù)器以及遠(yuǎn)程管理器等多個組件組成，如圖3所示^[11-13]。

燃?xì)夤艿赖娘L(fēng)險評估結(jié)果主要使用supermap is.net和supermap deskpro來實現(xiàn)。先將燃?xì)夤艿赖牡貓D信息及風(fēng)險評估結(jié)果存入deskpro。再使用supermap is.net將相關(guān)數(shù)據(jù)顯示出來，提高燃?xì)夤艿赖娘L(fēng)險管理效率。

4　實例分析

以某市燃?xì)夤艿罏槔撌腥細(xì)夤艿赖目傞L度為823.5km。

(1)統(tǒng)計該市從2001年到2011年每年燃?xì)夤艿朗鹿拾l(fā)生的次數(shù)以及事故發(fā)生后造成的經(jīng)濟(jì)損失，按照上述步驟計算燃?xì)夤艿赖娘L(fēng)險概率以及事故經(jīng)濟(jì)損失，其結(jié)果如圖4所示。

(2)對上述數(shù)據(jù)用回歸分析法進(jìn)行擬合，可得該城市的等風(fēng)險曲線，如表1和圖5所示，運(yùn)算結(jié)果為R=P*L=24.765，即該城市燃?xì)夤艿里L(fēng)險值為24.765。

(3)按照公式(1)—(4)計箅每段燃?xì)夤艿赖娘L(fēng)險值，運(yùn)用雙曲型風(fēng)險決策模型，將每段燃?xì)夤艿赖娘L(fēng)險值與該市燃?xì)夤艿赖娘L(fēng)險值相比較，判斷其風(fēng)險可接受性。若某段燃?xì)夤艿赖娘L(fēng)險值小于24.765，則該段燃?xì)夤艿赖娘L(fēng)險可接受，若大于24.765，則該段燃?xì)夤艿赖娘L(fēng)險不可接受，應(yīng)及時采取相關(guān)措施降低燃?xì)夤艿赖娘L(fēng)險水平。

(4)實現(xiàn)燃?xì)夤艿里L(fēng)險水平的可視化。首先得到該市的燃?xì)夤艿赖臇鸥駡D像，經(jīng)轉(zhuǎn)換成矢量數(shù)據(jù)并作適當(dāng)處理后存入supermap deskpro。將燃?xì)夤艿赖幕拘畔⒓帮L(fēng)險水平等相關(guān)數(shù)據(jù)導(dǎo)入supermap deskpro，再使用supermap is.net將相關(guān)數(shù)據(jù)顯示出來，可以使用查詢條件查詢?nèi)細(xì)夤艿赖幕拘畔⒓捌滹L(fēng)險水平，這樣可以通過圖得知各段管道的風(fēng)險情況并及時采取相應(yīng)的措施。查詢結(jié)果如圖6所示。

5　結(jié)論

(1)采用定量風(fēng)險評價方法對燃?xì)夤艿肋M(jìn)行風(fēng)險評估，排除主觀因素影響，從客觀的角度上反映燃?xì)夤艿赖娘L(fēng)險水平；

(2)基于雙曲型理論和回歸分析法建立雙曲型風(fēng)險決策模型以判斷燃?xì)夤艿赖娘L(fēng)險可接受性，為燃?xì)夤艿赖娘L(fēng)險管理提供科學(xué)決策依據(jù)；

(3)將雙曲型風(fēng)險決策模型分析的結(jié)果運(yùn)用supermap deskpro以及supermap is.net實現(xiàn)可視化表達(dá)，記錄每段燃?xì)夤艿赖幕拘畔⒓帮L(fēng)險可接受性，及時、方便的獲取燃?xì)夤艿赖娘L(fēng)險水平，對燃?xì)馄髽I(yè)的安全運(yùn)行與科學(xué)管理有重要意義。

參考文獻(xiàn)

1李志鵬，李艷紅，胡國新等.燃?xì)夤艿赖娘L(fēng)險評估方法[J].煤氣與熱力，2004；4：228-230

2黃小美.城市燃?xì)夤艿老到y(tǒng)定量風(fēng)險評價及其應(yīng)用研究[D].重慶：重慶大學(xué)，2008

3隋楠，魏立新，隋溪等.城市燃?xì)夤芫W(wǎng)風(fēng)險評價體系研究[J].油氣儲運(yùn)，2010；29：（25301）12-14+89

4金星，洪延姬等．工程系統(tǒng)可靠性數(shù)值分析方法[M]．北京：國防工業(yè)出版社，2002

5郭章林，谷洪雁，劉海新等．GIS在埋地燃?xì)夤艿里L(fēng)險分析的應(yīng)用[J]．煤氣與熱力，2005；7：33-35

6黃小美，彭世尼，張才華等．燃?xì)夤艿蓝匡L(fēng)險評價及其軟件研制[J]．煤氣與熱力，2010；30(1)：B24-B29

7劉燕，張應(yīng)輝，杜學(xué)平．基于回歸分析的北京市天然氣供應(yīng)規(guī)模預(yù)測[J]．煤氣與熱力，2010；30(6)：A26-A29

8何燦芝，羅漢．應(yīng)用數(shù)理統(tǒng)計[M]．長沙：湖南大學(xué)出版社，2004：95-ll6

9王巖，隋思連，王愛青．?dāng)?shù)理統(tǒng)計與MATLAB工程數(shù)據(jù)分析[M]．北京：清華大學(xué)出版社，2006：126-160

10吳信才地理信息系統(tǒng)原理與方法[M]．北京：電子工業(yè)出版社，2002

11陳樺，李小兵，徐光輝．基于SuperMap GIS的地理信息服務(wù)系統(tǒng)的設(shè)計與實現(xiàn)[J].計算機(jī)工程與設(shè)計，2009；8：2030-2033

12鐘廣銳．基于Ajax和supermap IS．NET 2008的WebGIS開發(fā)[J]．地理空間信息，2009；2：12-14

13程林，王美玲，張毅．一種基于SuperMap GIS的改進(jìn)Dijkstra算法[J]．地球信息科學(xué)學(xué)報，2010；5：649-654

本文作者：郭章林　劉彥香　宮亮

作者單位：華北科技學(xué)院建筑工程學(xué)院

  河北工程大學(xué)經(jīng)濟(jì)管理學(xué)院

  山東理工大學(xué)農(nóng)業(yè)工程與食品科學(xué)學(xué)院