陳國華1,2 徐佳1,2 方順1,2 吳悅1,2 李明陽(yáng)3

1.華南理工大學(xué)安全科學(xué)與工程研究所；2.廣東省安全生產(chǎn)科技協(xié)同創(chuàng )新中心；3.廣東省特種設備檢測研究院珠海檢測院

摘要：隨著(zhù)城鎮天然氣需求量增加，燃氣管道并行鋪設現象逐漸增多。為避免低概率、高風(fēng)險事故發(fā)生，開(kāi)展并行燃氣管道失效概率研究，分析失效概率影響因素及參數敏感性。結果表明：并行管道失效概率在設計階段受管道間距影響最顯著(zhù)，在運行階段則受氣云半徑影響最為顯著(zhù)；由于其對源管道內壓敏感度更高，因此應盡量控制運行階段并行管道內壓波動(dòng)，避免峰值過(guò)大。

關(guān)鍵詞：并行燃氣管道；蒙特卡洛法；失效概率；參數敏感性

2012年，加拿大哥倫比亞省北部一條燃氣管道破裂爆炸，導致與其并行間距為3 m的一條燃氣管道發(fā)生破裂著(zhù)火，此次事件形成了一個(gè)長(cháng)17 m，寬7.6 m，深1.1 m的火坑[1]。近年來(lái)，由于國內城鎮燃氣需求量增加，燃氣管道并行鋪設現象逐漸增多，一旦其中某條管道因泄漏引發(fā)的蒸氣云爆炸，極易對其臨近鋪設的燃氣管道造成沖擊進(jìn)而也發(fā)生失效，產(chǎn)生多米諾效應并導致事故后果進(jìn)一步擴大[2]。因此，基于已構建的并行管道失效概率模型[3]，結合蒙特卡洛模擬法進(jìn)行計算，對并行管道失效概率進(jìn)行影響因素分析及參數敏感性分析。

1  并行管道失效概率影響因素

在并行管道設計階段，影響其事故發(fā)生概率的主要因素為管道間距、目標管道埋深、目標管道管徑、目標管道壁厚與土壤重度。管道運行參數中影響其失效概率的因素主要包括初始失效管道（以下簡(jiǎn)稱(chēng)源管道）泄漏氣云半徑、源管道內壓及目標管道內壓。

2  并行管道失效概率計算

管材為X70鋼的兩條燃氣管道并行鋪設，管材屈服強度為480 MPa。以源管道發(fā)生大孔泄漏為初始失效場(chǎng)景，該場(chǎng)景下并行管道失效概率計算模型相關(guān)隨機變量分布規律如表 1所示。采用蒙特卡洛抽樣方法對各隨機變量進(jìn)行107次隨機抽樣，隨后將各組隨機抽樣值代入并行管道失效概率計算模型進(jìn)行計算，得到該場(chǎng)景不同參數影響下的并行管道失效概率曲線(xiàn)，如圖 1、圖 2所示。

表 1 并行管道失效概率模型相關(guān)隨機變量統計分布[4-6]

圖 1 不同設計參數影響下的并行管道失效概率曲線(xiàn)

圖 2 不同運行參數影響下的并行管道失效概率曲線(xiàn)

3  并行管道失效概率影響因素分析

由圖 1（a）和（b）可知，增大管道間距和目標管道壁厚均可以使并行管道失效概率下降。增大管道間距使目標管道受爆炸載荷影響減少，進(jìn)而使其損傷幾率下降。增大管道壁厚，有利于提高目標管道抵抗外部荷載的能力。同時(shí)，壁厚增大也會(huì )減少其因管體受腐蝕、施工挖掘與其他環(huán)境因素而遭受損傷的可能，進(jìn)而有效降低并行管道失效概率。

由圖 1（c）和（d）可知，目標管道埋深與管徑對并行管道失效概率有著(zhù)相反的影響。對于管道并行鋪設場(chǎng)景，增加埋深可使其降低遭受施工挖掘、地面違規占壓等第三方損壞的概率。由于燃氣泄漏過(guò)程中土壤對泄漏燃氣的黏滯作用以及吸收作用，增加埋深也有助于減少燃氣泄漏量，降低目標管道失效概率。由圖 1（e）可知，土壤重度對并行管道失效概率影響較小。

在管道運行階段，泄漏氣云半徑及管道內壓為并行管道失效概率主要影響因素。如圖 2（a）所示，氣云半徑與并行管道失效概率呈正相關(guān)關(guān)系。而氣云半徑大小主要取決于源管道泄漏時(shí)間，為探究該參數與并行管道失效概率之間的關(guān)系，基于文獻[2]的管道燃氣泄漏擴散模型，獲得源管道泄漏時(shí)間與并行管道失效概率關(guān)系曲線(xiàn)。如圖 2（b）所示，并行管道失效概率隨著(zhù)源管道泄漏時(shí)間的增加而逐漸增大。在源管道泄漏時(shí)間較短時(shí)，泄漏氣云半徑相對較小，因此目標管道失效概率及并行管道失效概率均較小。而在源管道發(fā)生泄漏約2200 s后，由于其產(chǎn)生的爆炸載荷超出管道可承受的臨界值，導致目標管道必然發(fā)生失效，目標管道失效概率值趨近于1，此時(shí)，并行管道失效概率主要取決于源管道失效概率。由圖 2（c）和（d）可知，增高目標管道內壓促使管道抵抗外部載荷的能力增強，從而有利于使并行管道失效概率下降，而增大源管道內壓，使源管道失效泄漏時(shí)燃氣的泄放速率上升，增加了源管道燃氣泄漏量，目標管道失效概率增大，進(jìn)而使并行管道失效概率上升。

4  并行管道失效概率參數敏感性分析

為了提高管道故障預防和事故控制的有效性，對各參數進(jìn)行敏感性分析。敏感性大小用敏感性指數αi表示，其計算公式如下[7]：

式中，Ci表示上述影響參數中某參數的取值，ΔCi為一定范圍內該影響參數的變化量，Fp表示并行管道失效概率，ΔFp為由ΔCi引起Fp的變化量。

令取值分別在5%至25%共5個(gè)階梯范圍內波動(dòng)，得到如圖 3、圖 4所示敏感性分析結果。由圖 3可知，并行管道設計階段各參數敏感性大小排序依次為：管道水平間距、土壤埋深、管道壁厚、管道內徑、土壤重度。在規定波動(dòng)范圍內，管道間距和土壤埋深對并行管道失效概率影響顯著(zhù)。

圖 3 并行管道設計階段參數敏感性分析結果

圖 4 并行管道運行階段參數敏感性分析結果

由圖 4可以看出，在規定波動(dòng)范圍內，并行管道運行階段各參數敏感性大小排序依次為：氣云半徑、源管道內壓、目標管道內壓。并行管道失效概率與源管道內壓之間呈正相關(guān)關(guān)系，與目標管道內壓之間則呈負相關(guān)關(guān)系，但由敏感性分析結果對比可知，并行管道失效概率對源管道內壓敏感性更高。

5  結論

（1）并行管道設計階段，各參數對并行管道失效概率影響程度從大到小依次為管道水平間距、土壤埋深、管道壁厚、管道內徑、土壤重度。

（2）并行管道運行階段，影響程度大小排序為氣云半徑、源管道內壓、目標管道內壓。

（3）并行管道失效概率與源管道內壓呈正相關(guān)，與目標管道內壓則呈負相關(guān)，但對源管道內壓敏感性更高，因此在管道運行過(guò)程中應盡量控制運行階段并行管道的內壓波動(dòng)，避免峰值過(guò)大。

參考文獻：

[1]Ñanalysis based on pipeline crater models and historical accidents[J]. Journal of loss prevention in the process industries, 2016, 43：315-331.

[2]RAMIREZ C J G, PASTOR E, CASAL J, MUÑOZ G F. Analysis of domino effect in pipelines[J]. Journal of hazardous materials, 2015, 298：210-220.

[3]徐佳. 埋地并行交叉天然氣管道爆炸耦合風(fēng)險分析及評估研究[D].廣州：華南理工大學(xué), 2020.

[4]中華人民共和國住房和城鄉建設部, 中華人民共和國國家質(zhì)量監督檢驗檢疫總局.輸氣管道工程設計規范：GB 50251-2015 [S]. 北京：中國計劃出版社, 2015.

[5]中華人民共和國住房和城鄉建設部, 中華人民共和國國家質(zhì)量監督檢驗檢疫總局.給水排水管道工程施工及驗收規范： GB 50268-2008[S]. 北京：中國標準出版社，2008.

[6]中華人民共和國國家質(zhì)量監督檢驗檢疫總局, 中國國家標準化管理委員會(huì ). 壓力管道規范動(dòng)力管道：GB/T 32270-2015[S]. 北京：中國標準出版社, 2015.

[7]韓文海，周晶. 腐蝕海底管道可靠性分析[J].石油學(xué)報，2015，36(4)：516-520.

基金項目：2020年度珠海市科技計劃項目，ZH22036205200013PWC；2019年廣東省中央引導地方科技發(fā)展專(zhuān)項資金/廣東省省級科技計劃項目，2019B020208012。

作者簡(jiǎn)介：陳國華，1967年生，教授，1993年博士畢業(yè)于南京工業(yè)大學(xué)學(xué)�；�工過(guò)程機械專(zhuān)業(yè)，現主要從事過(guò)程裝備安全可靠性及風(fēng)險評價(jià)技術(shù)方向的研究工作。聯(lián)系方式：13660261735，mmghchen@scut.edu.cn。