摘要：為解決天然氣管道泄漏的不確定性問題以及確定參數(shù)的敏感性，考慮到管道泄漏事故的突發(fā)和隨機性，本研究基于小孔、管道和大孔三類泄漏模型進(jìn)行泄漏后果的不確定性分析，采用蒙特卡洛分析方法分析泄漏孔面積、管道運行壓力、氣體溫度、泄漏處距離初始站點的距離四個不確定性因素對泄漏速率的影響。確定性分析結(jié)果表明，泄漏孔徑50mm的小孔泄漏速率為27.0996kg/s，管道破裂的瞬時泄漏速率高達(dá)54204.6289kg/s。在不同泄漏模型中，各參數(shù)的敏感度也不同，但泄漏孔面積始終為主導(dǎo)參數(shù)。當(dāng)泄漏孔達(dá)到管道內(nèi)徑的一半時，由管道模型計算的初始泄漏速率具有兩個數(shù)量級(10)的優(yōu)勢。因此，在天然氣管道運輸工程中，要及時對管道中出現(xiàn)的孔、隙等缺陷進(jìn)行檢測并采取補救措施，防止其進(jìn)一步發(fā)展成為管道泄漏。

　　關(guān)鍵詞：安全工程;管道泄漏;蒙特卡洛;不確定性;泄漏模型

　　0引言2017年，馬丁等[1]前瞻性地對全球能源用量進(jìn)行了梳理，指出未來15年的能源需求量會持續(xù)增加，其中，天然氣仍然保持著較高水平的消耗量。為了滿足天然氣的高消耗量，天然氣運輸管道數(shù)量不斷增加[2]，管道泄漏事故也相應(yīng)地頻繁發(fā)生，嚴(yán)重威脅人們的安全、健康與經(jīng)濟(jì)的發(fā)展。因此，科學(xué)地評估管道泄漏事故的后果以及準(zhǔn)確地預(yù)測事故的發(fā)生對于預(yù)防災(zāi)害、開展緊急救援等具有重要意義。

　　天然氣論文范例： 天然氣儲運技術(shù)及其應(yīng)用發(fā)展前景

　　目前普遍認(rèn)為，天然氣管道失效后果的研究應(yīng)充分考慮不確定參數(shù)及其隨機性，進(jìn)行不確定性的后果分析[3~8]。不考慮泄漏事故的突發(fā)與隨機性，僅針對確定性場景或單一參數(shù)進(jìn)行后果分析存在片面性。現(xiàn)有文獻(xiàn)中，不乏有針對泄漏后果分析的不確定性研究，如2002年和2018年，沈斐敏等[4]和熊鴻斌等[5]在對氣體泄漏問題的蒙特卡洛分析中，先后建立了以泄漏介質(zhì)、泄漏模式、泄漏孔徑、燃?xì)鈮毫Α�環(huán)境溫度、風(fēng)速、泄漏處距起始點距離等作為不確定性參數(shù)的體系。

　　2014年，AlzbutasR等提出了參數(shù)敏感性是不確定性分析的一個重要內(nèi)容的思想。2018年，紀(jì)虹等[7]指出EGIG(歐洲輸氣管道事故數(shù)據(jù)組織)建立了以泄漏孔徑尺寸為劃分依據(jù)的小孔泄漏、大孔泄漏、管道斷裂三種天然氣管道泄漏模型。2015年，姜璐等引入以管道壓力、泄漏孔徑、環(huán)境風(fēng)速作為噴射火熱輻射傷害模型的不確定性參數(shù)，通過蒙特卡洛模擬及拉丁超立方抽樣方法計算死亡半徑，得到距離管道泄漏處不同位置的事故死亡半徑的概率分布的新思想。2021年，PetroD.Ndalila等[9]建立了天然氣管道泄漏口動壓測量的模型，發(fā)現(xiàn)泄漏率會隨著泄漏口尺寸的變化而改變。

　　2021年，魯寨軍等[10]建立了非平衡相變的CO泄漏模型，通過數(shù)值模擬分析了CO管道泄漏的瞬態(tài)行為。目前，國內(nèi)外對參數(shù)進(jìn)行敏感度分析的研究較少，風(fēng)險分析中的不確定性參數(shù)的樣本數(shù)少，對各參數(shù)概率分布的描述缺乏說服力。同時，考慮多個不確定性參數(shù)時，沒有分析各參數(shù)對結(jié)果的影響程度的差異性。鑒于此，本文針對泄漏后果系統(tǒng)地進(jìn)行了不確定性分析，考慮泄漏孔面積(or)、管道運行壓力(in)、氣體溫度、泄漏處距離初始站點的距離四個因素為管道泄漏的不確定性參數(shù)[1114]，泄漏速率作為泄漏后果參數(shù)。

　　首先，通過拉丁超立方抽樣方法，獲得大量的參數(shù)樣本集，由管道泄漏的物理模型計算各參數(shù)樣本集對應(yīng)的泄漏速率，定義多參數(shù)與泄漏速率之間的秩相關(guān)系數(shù)，以表征各不確定性參數(shù)對管道泄漏速率的影響程度;在此基礎(chǔ)上，采用概率分布函數(shù)描述各參數(shù)的不確定性，通過蒙特卡洛抽樣，獲得大量的參數(shù)樣本集，由管道泄漏的物理模型計算各參數(shù)樣本集對應(yīng)的泄漏速率，并統(tǒng)計分析泄漏速率的不確定性。

　　1不確定性分析方法

　　1.1泄漏后果的物理模型管道內(nèi)的氣體泄漏，可歸結(jié)為連續(xù)氣體的泄漏問題。氣體在壓力條件下從設(shè)備的裂口泄漏時，需要通過氣體流動標(biāo)準(zhǔn)方程計算。由于氣體發(fā)生亞音速流動、音速流動時滿足的基本假設(shè)不同，其標(biāo)準(zhǔn)方程的形式也有差異[3]。

　　根據(jù)孔徑尺寸變化引起的泄漏氣體流動方式的改變，將管道泄漏模型分為小孔、管道和大孔模型[7]三類，適用范圍為：1)小孔模型適用于泄漏孔徑小的管道，管道內(nèi)壓不影響泄漏速率，氣體膨脹過程等熵;2)管道模型適用于完全破裂的管道，氣體無等熵膨脹過程;3)大孔模型適用于各種泄漏孔徑的計算。

　　目前，最常用的不確定性分析方法是概率分析法和模糊集理論分析法，前者是定量分析方法，后者是定性分析方法。概率分析法通過概率分布量化變量的不確定性，能比較逼真地描述具有隨機性質(zhì)的事物的特點及物理實驗過程;而模糊集理論分析法研究概念模糊的事物，由于當(dāng)前研究尚不成熟，在實踐中應(yīng)用得較少。

　　概率分析法中最常用的是蒙特卡洛分析(Monte CarloAnalysis)，即通過對概率模型進(jìn)行抽樣實驗來計算所求參數(shù)的統(tǒng)計特性，得出該事件發(fā)生的頻率。如姜璐等通過蒙特卡洛模擬計算死亡半徑，得到距離管道泄漏處不同位置的事故死亡半徑的概率分布。MunirAliElfarra等[1通過蒙特卡洛模擬對風(fēng)速進(jìn)行隨機抽樣，計算出年風(fēng)能產(chǎn)量分布，呈現(xiàn)出較好的置信水平。S.Triambak等[1采用一個二維蒙特卡洛模型來研究COVID19類型的感染在一個可控人群中的傳播，并獲得了良好的增長指數(shù)。

　　2計算結(jié)果與討論

　　案例：某天然氣管道工程于2011年投產(chǎn)，管道全長288km，管道外徑為1016mm，設(shè)計壓力10MPa，運行壓力為～8.5，管道材質(zhì)L485螺旋焊縫管，設(shè)計輸量9×10/a。管道從某村東側(cè)穿過，房屋距離管道最近距離為25米，村落沿管道呈條帶狀，長年居住有155戶465人左右，小鎮(zhèn)住宅區(qū)，均為四層以上樓房。研究管段氣候?qū)贉貛О霛駶櫦撅L(fēng)型大陸氣候四季分明，雨熱同季，光熱條件好。月平均氣溫月份最低，平均7.4℃，月份最高，平均24℃，歷年平均氣溫9.4℃。

　　3結(jié)論

　　本文將泄漏孔面積(or)、管道運行壓力(in)、氣體溫度、泄漏處距離初始站點的距離作為天然氣管道泄漏事故的不確定性參數(shù)，考慮其概率分布函數(shù)，并結(jié)合一個工程實例進(jìn)行分析。確定性分析結(jié)果顯示，發(fā)生泄漏孔徑為50mm的小孔泄漏時，泄漏速率為27.0996kg/s;發(fā)生管道破裂時，瞬時泄漏速率高達(dá)54204.6289kg/s。

　　對參數(shù)進(jìn)行拉丁超立法抽樣，并計算Spearman秩相關(guān)系數(shù)，得出敏感度最大的參數(shù)為泄漏孔面積。蒙特卡洛分析結(jié)果表明，在實例的計算參數(shù)條件下，泄漏速率的累積概率分布函數(shù)在83%處有一個轉(zhuǎn)折點，說明多數(shù)情況下的泄漏速率Q<2000kg·s。為了減小管道附近區(qū)域燃燒、爆炸的可能性，要及時對天然氣管道中出現(xiàn)的孔、隙等缺陷進(jìn)行檢測，并對缺陷處采取補救措施，防止其進(jìn)一步發(fā)展成為管道泄漏。

　　參考文獻(xiàn)

　　[1]馬丁，單葆國.2030年世界能源展望——基于全球能源展望報告的對比研究[J].中國能源，2017，39(2)：21-24.MAD，SHANBG.WorldEnergyOutlook2030Brexit’sPotentialImpactsontheSinoUKEnergy[J].EnergyofChina，2017，39(2)：21-24.

　　[2]周祥君.天然氣管道工程施工建設(shè)質(zhì)量管理探析[J].全面腐蝕控制，2021，35(3)：84-86.ZHOUX.AnalysisonConstructionQualityManagementofNaturalGasPipelineEngineering[J].TotalCorrosionControl，201，35(3)：84-86.

　　[3]中國石油化工股份有限公司青島安全工程研究院.石化裝置定量風(fēng)險評估指南[M].北京：中國石化出版社，2007.QingdaoSafetyEngineeringResearchInstituteofChinaPetrohemicalCo.Ltd.Aguidanceforquantitativeriskassessmentinthepetrochemicalplant[M].Beijing：ChinaPetrochemicalPress，2007.

　　作者：陳華燕，劉演錕，敬思岑，雷鳴浩