基于可靠性的定量風(fēng)險評價(jià)技術(shù)及應用

——以某天然氣管道地區等級升級為例 

周亞薇  張振永  田姍姍 

中國石油管道局工程有限公司設計分公司

摘 要：隨著(zhù)我國社會(huì )、經(jīng)濟高速發(fā)展，一些在役天然氣管道通過(guò)的地區逐漸由人口稀少的一級、二級地區發(fā)展成為人口密集的三級、四級地區（高后果區），一旦管道發(fā)生泄漏事故，其造成的后果將更加嚴重。因此，有必要建立相應的定量風(fēng)險評價(jià)技術(shù)，確定管道實(shí)際風(fēng)險水平�；�于可靠性的定量風(fēng)險評價(jià)技術(shù)實(shí)質(zhì)是采用基于可靠性的極限狀態(tài)方法計算評價(jià)管段的失效概率，結合建立的國內天然氣管道失效后果模型，綜合計算管道失效風(fēng)險。之后，引入風(fēng)險可接受準則，判定地區等級升級后管道個(gè)體風(fēng)險和社會(huì )風(fēng)險的可接受性，并有針對性地制定風(fēng)險消減/防護措施。該技術(shù)不僅可以實(shí)現天然氣管道失效風(fēng)險的定量計算和評價(jià)，且能夠對不同風(fēng)險消減/防護措施進(jìn)行效果檢查及對比，在保障管道風(fēng)險可接受的前提下，實(shí)現技術(shù)性和和經(jīng)濟性的最優(yōu)化。

關(guān)鍵詞：天然氣管道；地區等級；可靠性；定量風(fēng)險評價(jià)；風(fēng)險可接受準則

 

通常情況下，天然氣管道敷設路由會(huì )盡量避開(kāi)人口密集、交通頻繁、地下設施多的區段。在管道設計階段，輸氣管道通過(guò)的地區將按照沿線(xiàn)居民戶(hù)數和/或建筑物的密集程度劃分符合標準要求的地區等級，并依據地區等級，結合地區的發(fā)展規劃等選擇相應的設計系數等管道設計方案。但是，隨著(zhù)我國社會(huì )、經(jīng)濟高速發(fā)展，許多管道沿線(xiàn)的建筑物和人口密度與建設初期相比均發(fā)生了較大變化，一些在役天然氣管道通過(guò)的地區逐漸由人口稀少的一級、二級地區發(fā)展成為人口密集的三級、四級地區，地區等級升級使得原始的管道設計方案與目前現狀不匹配^[1-3]。

由于天然氣管道的輸送介質(zhì)易燃、易爆，管道失效后因火災熱輻射造成的潛在影響半徑高達上百米至幾百米不等，管道周邊人口密度顯著(zhù)增加后，管道失效所帶來(lái)的失效后果將更加嚴重；另一方面，地區等級升級將導致沿線(xiàn)人員活動(dòng)、設備撞擊等帶來(lái)的管道損壞更加頻繁，管道的失效概率也隨之增大。因此，地區等級升級將導致管道失效風(fēng)險（失效概率與失效后果的乘積）大幅提升。為了準確量化管道風(fēng)險，并針對性地制定風(fēng)險管控措施，本文提出了一套基于可靠性的定量風(fēng)險評價(jià)流程，依據定量的失效概率與失效后果分析，對輸氣管道地區等級升級進(jìn)行定量風(fēng)險評價(jià)，通過(guò)個(gè)體風(fēng)險和社會(huì )風(fēng)險水平的判定，制訂風(fēng)險管控措施。

1　基于可靠性的定量風(fēng)險評價(jià)

管道風(fēng)險評價(jià)方法一般分為定性、半定量和定量風(fēng)險評價(jià)方法^[4]。定性評價(jià)和半定量評價(jià)方法簡(jiǎn)單易用，但不能定量計算管道的失效概率、后果及風(fēng)險，難以評判管道風(fēng)險是否符合國家行業(yè)規定的風(fēng)險可接受水平，也無(wú)法通過(guò)量化計算科學(xué)地確定風(fēng)險管控措施。定量風(fēng)險評價(jià)（Quantitative Risk Assessment，QRA）方法是管道風(fēng)險評價(jià)的高級階段，它將管道的失效概率和事故后果的進(jìn)行定量計算，實(shí)現了對管道風(fēng)險的精確描述^[5]。根據失效概率的計算方法，定量風(fēng)險評價(jià)可以分為基于失效統計的定量風(fēng)險評價(jià)技術(shù)和基于可靠性的定量風(fēng)險評價(jià)技術(shù)。

基于失效統計的定量風(fēng)險評價(jià)技術(shù)是通過(guò)對管道失效事故統計分析，建立基線(xiàn)失效概率，并根據待評價(jià)管道的實(shí)際情況分配不同的修正因子，以此來(lái)計算管道的失效概率^[6]�；�于可靠性的天然氣管道定量風(fēng)險評價(jià)技術(shù)實(shí)質(zhì)是采用基于可靠性的極限狀態(tài)方法，針對天然氣管道評價(jià)管段，通過(guò)分析管道沿線(xiàn)環(huán)境和荷載狀況，確定可能導致管道失效的主要極限狀態(tài)和狀態(tài)方程，采用應力-強度分布干涉理論計算管段失效概率；失效后果模型考慮了在一定的泄漏頻率、泄漏量、立即點(diǎn)燃情景下，熱輻射引起管道周?chē)�人員傷亡的程度，從而定量計算管道風(fēng)險。

基于可靠性的定量風(fēng)險評價(jià)方法工作流程見(jiàn)圖1，主要包括管段劃分、失效概率計算、失效后果計算、風(fēng)險計算、風(fēng)險評價(jià)和風(fēng)險決策等。

 

圖1 基于可靠性的定量風(fēng)險評價(jià)流程圖

1.1　失效概率計算

近年來(lái)，中石油集團公司開(kāi)展了天然氣管道基于可靠性的設計和評價(jià)方法研究，收集了國內近4萬(wàn)公里已建天然氣管道數據，構建了國內天然氣管道管材、施工和運行維護等變量參數數據庫，基于國內管道數據建立了系統了天然氣管道基于可靠性的設計流程，并形成了相應標準^[7-10]。

根據建立的極限狀態(tài)方程，利用蒙特卡洛方法，對管道的極端極限狀態(tài)（大孔泄漏和破裂）進(jìn)行模擬仿真計算，以此確定特定管道設計工況的失效概率。為了簡(jiǎn)化分析，僅計算腐蝕和第三方破壞失效概率。根據國內外數據統計，腐蝕和施工設備等撞擊造成的大泄漏和破裂約占所有失效的60%～76%。為了統籌考慮導致管道失效的其他因素，腐蝕與設備沖擊的失效概率應分別在計算結果的基礎上放大1.5倍^[11-13]。管道極端極限狀態(tài)下的總失效概率計算如下^[14]：

                            （1）

式中P_ULS表示極端極限狀態(tài)失效概率，/km·a；P_LL和P_RU分別為大泄漏失效概率和破裂失效概率，/km·a。

腐蝕缺陷的尺寸是隨時(shí)間增長(cháng)的，腐蝕失效概率與時(shí)間具有相關(guān)性；第三方設備撞擊造成的失效概率與時(shí)間無(wú)關(guān)，管道每公里每年的失效概率等于每公里每年的沖擊次數和每次沖擊造成的失效概率之積。管道腐蝕失效概率和第三方破壞失效概率分別使用加拿大C-FER公司PRISM^TM軟件中的時(shí)間相關(guān)模型和時(shí)間無(wú)關(guān)模型進(jìn)行年失效概率的計算。

1.2　失效后果計算

失效后果模型考慮了在一定的泄漏頻率、泄漏量、立即點(diǎn)燃情景下，熱輻射引起管道周?chē)�人員傷亡的程度。根據天然氣研究協(xié)會(huì )（Gas Research Institute， GRI）研究成果，將失效后果定義為死亡人數的函數，與天然氣著(zhù)火釋放的熱量有關(guān)。由于破裂的失效后果遠大于泄漏的后果，因此破裂是失效后果的主要控制因素，管道破裂情況下的預期死亡人數計算如下：

                                 （2）

                               （3）

式中，為預期死亡人數；為著(zhù)火概率；為失效事故影響范圍；為人口密度，人/公頃；為實(shí)際占有概率（失效事故發(fā)生時(shí)，公眾出現在事故影響范圍內的概率），D為管道的直徑，mm。

根據Stephen^[15]等人建立的模型，危害區域假定為圓形。圖2中兩個(gè)危害區域和相應的半徑定義了對應的熱強度的上下限門(mén)檻值，上限門(mén)檻值確定的范圍內，假定致死率是100%，在下限門(mén)檻值確定的范圍外，致死率是0。在這兩個(gè)門(mén)檻值之間，室外的致死率是50%，室內的致死率是25%，管道破裂條件下的死亡人數表達為式（4）：

        （4）

                                   （5）

根據對國內陸上天然氣管道沿線(xiàn)人口密度及其活動(dòng)頻率調查分析，P_in=80%；P_out=20%；對于室外情況，下限和上限分別是12.6 kw/m²和31.6 kw/m²；對于室內情況，下限和上限分別是15.8 kw/m²和31.6 kw/m²；=40%。由此得：

                            （6）

式中，為人口密度，人/公頃；為管道運行壓力，MPa；為管道直徑，mm。其中，人口密度根據評價(jià)初期人口數據調研收集和潛在影響范圍綜合確定。

 

圖2 估計預期死亡人數的危害區域

1.3　風(fēng)險可接受準則

輸氣管道風(fēng)險可接受準則工程上較多地采用最低合理可行（As Low As Reasonably Practically，ALARP）原則，ALA RP原則認為任何工業(yè)系統都存在風(fēng)險，不可能通過(guò)預防措施徹底消除，當系統的風(fēng)險水平超低時(shí)，要進(jìn)一步降低就很困難。為此所花費的成本往往呈上升趨勢，也可以表示為安全風(fēng)險改進(jìn)措施投資的邊際效益遞減，趨于零，最終為負值。因此，必須在工業(yè)系統的風(fēng)險水平和成本之間做出折衷^[16]。

分別用個(gè)體風(fēng)險與社會(huì )風(fēng)險來(lái)衡量輸氣管道地區等級升級后的管道風(fēng)險^[17]。

（1）個(gè)體風(fēng)險是指在評價(jià)位置長(cháng)期生活、工作的，并未采取任何防護措施的人員遭受特定危害而死亡的概率�？梢园聪率接嬎悖�

                                    （7）

式中P為失效概率；P_i為點(diǎn)燃概率；L_ir為相互作用長(cháng)度，該長(cháng)度定義為事故有可能影響所考慮位置的管段長(cháng)度，計算方法參考圖3；為占用概率。

 

a-相互作用長(cháng)度為2×R； b-相互作用長(cháng)度為2×

1-觀(guān)察點(diǎn)位置； 2-環(huán)向影響半徑； 3-管道； 4-相互作用長(cháng)度

圖3 相互作用長(cháng)度計算示意圖

基于ALARP原則，根據《油氣輸送管道風(fēng)險評價(jià)導則》（SY/T 6859-2012），個(gè)體風(fēng)險分別以1×10^-4和1×10^-6為界分為不可接受區、可接受區（即最低合理可行區ALARP）和廣泛接受區^[18]三個(gè)區域（圖4）。

 

圖4 個(gè)體風(fēng)險可接受標準推薦值

（2）社會(huì )風(fēng)險用于描述事故發(fā)生的可能性和災害導致人員傷亡數量之間的關(guān)系，或者解釋為每公里管道事故發(fā)生概率（每年）F和事故導致的死亡人數N之間的關(guān)系。F即為極端極限狀態(tài)失效概率P_ULS。

社會(huì )風(fēng)險的ALARP準則是建立在F-N曲線(xiàn)的基礎上，根據事故頻率F與死亡人數N分析輸氣管道失效的社會(huì )風(fēng)險。結合《油氣輸送管道風(fēng)險評價(jià)導則》（SY/T 6859-2012），社會(huì )風(fēng)險可接受標準（F-N曲線(xiàn)）如圖5所示。

 

圖5 社會(huì )風(fēng)險可接受標準（F-N）曲線(xiàn)

3　地區等級升級地區應用案例

某管線(xiàn)鋼管直徑為D610mm，直管段全部采用L415螺旋縫埋弧焊鋼管，設計壓力為6.3 MPa。在施工圖設計階段，該評價(jià)段管道所處地區等級為二級地區，鋼管設計壁厚為9.5 mm。評價(jià)段管道全部采用溝埋敷設，根據施工圖信息，評價(jià)段管道管頂平均埋深為1.5 m。目前，評價(jià)段管道兩側5米間距處均建有兩排多層民房（見(jiàn)圖6），人口密度與施工建設活動(dòng)增多，導致管段所處區域的地區等級升級至三級地區。

圖6  評價(jià)段管道現場(chǎng)情況

3.1  管段風(fēng)險評價(jià)

根據GB 32167-2015《油氣輸送管道完整性管理規范》規定的潛在影響半徑計算方法，將該評價(jià)段的潛在影響范圍識別為圖7的矩形區域，評價(jià)段長(cháng)度400 m，潛在影響半徑153 m，潛在影響范圍面積約12.4公頃。根據收集的潛在影響范圍內人口數據，計算得地區等級升級后人口密度為16.9人/公頃。

 

圖7 潛在影響區域

結合國內的管材、焊接、腐蝕和運行維護等數據的統計分析結果，建立參數的計算模型，選擇合適的設計參數及維修計劃等，利用軟件計算失效概率，計算模擬次數為1億次，假定管道設計壽命為30年，分別取30年模擬計算中外部腐蝕失效概率（圖8）和第三方破壞失效概率（圖9）的最大失效概率（大孔泄漏、破裂），計算得總失效概率為1.32E-06（表1）。采用已有失效后果模型計算評價(jià)范圍內的預期死亡人數N=10人。

 

圖8 外部腐蝕失效概率計算結果/km·a

 

圖9  第三方破壞失效概率計算結果/km·a

                                      表1 失效概率計算結果                             /km·a

采用已有模型，計算得出個(gè)體風(fēng)險為4.85E-05，處于個(gè)體風(fēng)險可接受區；將計算所得的死亡人數和事故發(fā)生概率通過(guò)F-N曲線(xiàn)進(jìn)行風(fēng)險評價(jià)。由圖10可見(jiàn)，評價(jià)管段的社會(huì )風(fēng)險位于F-N曲線(xiàn)的不可接受區，說(shuō)明該管段所處地區等級升級后，社會(huì )風(fēng)險水平略高于標準規定，需要采取相應的風(fēng)險減緩措施。

 

圖10  社會(huì )風(fēng)險評價(jià)結果

3.2　風(fēng)險管控措施

對于風(fēng)險水平高的管道，應采取風(fēng)險減緩措施，并對采取措施后的工況再次進(jìn)行風(fēng)險評價(jià)，并衡量不同措施之間的經(jīng)濟效益。常用的措施可為下列任意一種或幾種的組合形式：增加第三方防護措施，選擇混凝土蓋板防護；縮短內檢測的時(shí)間間隔；降低運行壓力；增加管道壁厚；加大管道埋深；改變管道路由等。

針對上述措施，分別對風(fēng)險重新評價(jià)，結果顯示選擇上述任意一種措施后，均能將社會(huì )風(fēng)險降低至可接受區范圍內；但是考慮到縮短內檢測間隔至5年/次將大幅增加管道運營(yíng)資金投入；降低運行壓力可能影響到下游用戶(hù)的用氣需求以及降低管道運營(yíng)盈利；增加壁厚和加大埋深都將導致大量的工程投入和實(shí)施困難等，最終建議對該地區等級升級段采取增加混凝土蓋板進(jìn)行第三方防護的措施，以此保證該評價(jià)段管道風(fēng)險水平符合標準的規定。

4　結論及建議

（1）在役天然氣管道的地區等級升級導致原始的管道設計方案與目前現狀不匹配，存在管道失效風(fēng)險隱患，通過(guò)定量風(fēng)險評價(jià)技術(shù)能夠準確地量化和評價(jià)管道風(fēng)險，并有針對性地制定風(fēng)險管控措施，可以明確管控重點(diǎn)，且降低管理成本。

（2）基于可靠性的設計方法通過(guò)對國內已建天然氣管道及其沿線(xiàn)相關(guān)數據的統計分析，結合國內實(shí)際科學(xué)地進(jìn)行失效概率的定量計算，避免了傳統失效可能性計算過(guò)程中主觀(guān)因素的影響，使量化結果更加適用于我國的天然氣管道的建設實(shí)際。

（3）風(fēng)險管控措施種類(lèi)較多，依據風(fēng)險可接受準則，對于風(fēng)險評價(jià)結果顯示需采取風(fēng)險減緩措施的地區，應分別針對特定風(fēng)險管控措施的定量分析確定其有效性的基礎上，結合經(jīng)濟性和可行性綜合確定風(fēng)險減緩措施。

（4）對于地區等級暫時(shí)未發(fā)生變化，但根據地區規劃等發(fā)展需求有可能導致地區等級升級的區域，可通過(guò)基于可靠性的定量風(fēng)險評價(jià)技術(shù)，確定其未來(lái)可能出現的風(fēng)險水平，并制定相應的風(fēng)險防控措施。

 

參考文獻：

[1] 董紹華，王東營(yíng)，費凡，等. 管道地區等級升級與公共安全風(fēng)險管控 [J]. 油氣儲運，2014，33(11)：

   1164-1170.
 [2] 姚安林，周立國，汪龍，等. 天然氣長(cháng)輸管道地區等級升級管理與風(fēng)險評價(jià) [J]. 天然氣工業(yè)，2017，37(1)： 

   124-130.

 [3] 中國石油集團工程設計有限責任公司西南分公司. 輸氣管道工程設計規范：GB50251-2015 [S]. 北京：中國計劃出版社，

   2015.

 [4] 王其磊，程五一，張麗麗，等. 管道量化風(fēng)險評價(jià)技術(shù)與應用實(shí)例 [J]. 油氣儲運，2011，30(7)：494-496.

 [5] 張圣柱，吳宗之. 油氣管道風(fēng)險評價(jià)與安全管理 [M]. 北京：化學(xué)工業(yè)出版社，2016.

 [6] 張華兵. 基于失效庫的在役天然氣長(cháng)輸管道定量風(fēng)險評價(jià)技術(shù)研究 [D]. 中國地質(zhì)大學(xué)（北京），2013.

 [7] ZHANG W W，ZHANG Z Y，ZHANG J Y，et al. Research on reliability-based design technique of China domestic onshore gas   pipeline [J]. Guangzhou：Proceedings of the 10th International Conference on Reliability, Maintainability and Safety，2014：IEEE  32999

[8] ZHANG J Y，ZHANG Z Y，YU Z F，et al. Building a target reliability adaptive to China onshore natural gas pipeline [C]. Calgary：

   Proceedings of the 10th International Pipeline Conference，2014：IPC 2014-33129.

[9] 溫凱，張文偉，宮敬，等. 天然氣管道可靠性的計算方法 [J]. 油氣儲運，2014，33（7）：729-733.

[10] ZHOU J，ROTHWELL B，NESSIM M et al. Development of reliability-based design and assessment standards for onshore

   natural gas transmission pipelines [C]. Calgary：Proceedings of the 6th International Pipeline Conference，2006：IPC 2006-10045.

[11] ZIMMERMAN T，NESSIM M，MCLAMB M，et al. Target reliability levels for onshore gas pipelines [C]. Calgary：Proceedings

   of the 9th International Pipeline Conference，2002：IPC 2002-27213.

[12] NESSIM M.，ZIMMERMAN T，GLOVER A，et al. Reliability-based limit states design for onshore pipelines [C]. Calgary：

   Proceedings of the 9th International Pipeline Conference，2012：IPC 2012-27125.

[13] LU J，WU W，ZHANG Z Y，et al. Probability calculation of equipment impact based on reliability method [C]. Calgary： 

   Proceedings of the 10th International Pipeline Conference，2014：IPC 2014-33147.

[14] 張振永，周亞薇，張金源. 國內天然氣管道強度設計系數的評估研究 [J]. 天然氣工業(yè)，2017，37(04)：116-122.

[15] Stephens, M.J., Leewis, K. and Moore, D.K. 2002. “A Model for Sizing High Consequence Areas Associated with Natural Gas

   Pipelines[C]. Proceedings, IPC02, September, Calgary, Alberta, Canada, Paper No. IPC02-27073.

[16] 馬志祥. 油氣長(cháng)輸管道的風(fēng)險管理 [J]. 油氣儲運，2005，24(02)：1-7+4-61.

 [17] 邢志祥.天然氣長(cháng)輸管道的定量風(fēng)險評價(jià)方法[J].石油機械，2008，36(4)：15-17.
 [18] 中國石油管集團石油管工程技術(shù)研究院. 油氣輸送管道風(fēng)險評價(jià)導則：SYT 6859-2012 [S]. 北京：石油工業(yè)出版社，2012.

   

作者：周亞薇，女，1986年生，工程師，2012年碩士畢業(yè)于中國地質(zhì)大學(xué)（北京）礦產(chǎn)普查與勘探專(zhuān)業(yè)，現主要從事油氣長(cháng)輸管道線(xiàn)路工程設計和管道材料、可靠性以及風(fēng)險評估相關(guān)研究。

《管道保護》2017年第6期（總第37期）