作者：焦建瑛，張濤，王嵩梅， 何少平，任立坤，陳鵬第一作者單位：北京市燃?xì)饧瘓F(tuán)有限責(zé)任公司
摘自《煤氣與熱力》2021年3月刊
1 概述 隨著社會(huì)經(jīng)濟(jì)發(fā)展，天然氣在我國(guó)經(jīng)濟(jì)和能源生產(chǎn)中的地位越來(lái)越重要。作為城市公共基礎(chǔ)設(shè)施的主要組成部分，城鎮(zhèn)燃?xì)夤艿腊踩\(yùn)行不僅關(guān)系城市經(jīng)濟(jì)發(fā)展，更關(guān)系到居民的生命財(cái)產(chǎn)安全［1］。 我國(guó)城市大部分燃?xì)夤艿赖姆蹠r(shí)間已超出30 a，正逐步進(jìn)入事故多發(fā)階段［2］。因此，有效評(píng)估燃?xì)夤艿?/span>泄漏風(fēng)險(xiǎn)非常必要。為了更好地評(píng)估燃?xì)夤艿赖慕】禒顩r，相關(guān)部門采取了很多的檢測(cè)措施，在特定時(shí)期內(nèi)對(duì)燃?xì)夤艿栏g程度進(jìn)行檢測(cè)。但檢測(cè)需要耗費(fèi)大量的人力物力，且存在檢測(cè)頻率低等問(wèn)題［3］。本文以燃?xì)夤艿佬孤╋L(fēng)險(xiǎn)影響因素為出發(fā)點(diǎn)，建立燃?xì)夤艿佬孤╋L(fēng)險(xiǎn)評(píng)價(jià)體系。 2 體系架構(gòu) 燃?xì)夤艿佬孤╋L(fēng)險(xiǎn)評(píng)價(jià)體系包括風(fēng)險(xiǎn)預(yù)測(cè)層、風(fēng)險(xiǎn)篩查層、風(fēng)險(xiǎn)預(yù)警層，架構(gòu)見(jiàn)圖1。風(fēng)險(xiǎn)預(yù)測(cè)層：將泄漏風(fēng)險(xiǎn)影響因素歸納為10個(gè)風(fēng)險(xiǎn)強(qiáng)度指標(biāo)，對(duì)評(píng)價(jià)時(shí)間段內(nèi)各風(fēng)險(xiǎn)強(qiáng)度指標(biāo)進(jìn)行評(píng)分，量化評(píng)價(jià)風(fēng)險(xiǎn)強(qiáng)度指標(biāo)。風(fēng)險(xiǎn)篩查層：基于風(fēng)險(xiǎn)強(qiáng)度指標(biāo)評(píng)分，計(jì)算綜合風(fēng)險(xiǎn)得分，對(duì)高風(fēng)險(xiǎn)區(qū)域進(jìn)行篩查，并量化各風(fēng)險(xiǎn)強(qiáng)度指標(biāo)的風(fēng)險(xiǎn)強(qiáng)度貢獻(xiàn)度。風(fēng)險(xiǎn)預(yù)警層：在時(shí)序上實(shí)現(xiàn)快速迭代，對(duì)未來(lái)綜合風(fēng)險(xiǎn)得分做出預(yù)判，以實(shí)現(xiàn)風(fēng)險(xiǎn)預(yù)警。 
圖1 燃?xì)夤艿佬孤╋L(fēng)險(xiǎn)評(píng)價(jià)體系架構(gòu)
3 風(fēng)險(xiǎn)預(yù)測(cè)層 通過(guò)整合燃?xì)夤艿佬孤╋L(fēng)險(xiǎn)影響因素，并結(jié)合燃?xì)夤艿佬孤┦鹿试?/span>［4］，將燃?xì)夤艿佬孤╋L(fēng)險(xiǎn)影響因素歸納為10個(gè)風(fēng)險(xiǎn)強(qiáng)度指標(biāo)：雜散電流風(fēng)險(xiǎn)強(qiáng)度、土壤腐蝕性風(fēng)險(xiǎn)強(qiáng)度、防腐層風(fēng)險(xiǎn)強(qiáng)度、管道本體風(fēng)險(xiǎn)強(qiáng)度、庭院管道環(huán)境風(fēng)險(xiǎn)強(qiáng)度（僅適用于庭院管道）、地面沉降風(fēng)險(xiǎn)強(qiáng)度、地質(zhì)擾動(dòng)風(fēng)險(xiǎn)強(qiáng)度、天氣風(fēng)險(xiǎn)強(qiáng)度、施工風(fēng)險(xiǎn)強(qiáng)度、穿跨越風(fēng)險(xiǎn)強(qiáng)度。 3.1 評(píng)價(jià)單元 在實(shí)際工程中，由管網(wǎng)編號(hào)（由壓力級(jí)制編碼、測(cè)量編碼、順序碼、地區(qū)編碼等組成）描述的管段是燃?xì)夤艿赖淖钚卧?，由于管段長(zhǎng)度過(guò)短，將管段作為風(fēng)險(xiǎn)評(píng)價(jià)單元并不合適。因此，將由多條管段組成的評(píng)價(jià)單元作為燃?xì)夤艿佬孤╋L(fēng)險(xiǎn)評(píng)價(jià)體系的評(píng)價(jià)對(duì)象。 通常，評(píng)價(jià)單元可分為固定長(zhǎng)度單元、可變長(zhǎng)度單元［5］。本文采用可變長(zhǎng)度單元，并將管道分為庭院管道、庭院外管道：庭院管道以小區(qū)燃?xì)夤艿雷鳛樵u(píng)價(jià)單元，庭院外管道以若干分配管段作為評(píng)價(jià)單元。
3.2 風(fēng)險(xiǎn)強(qiáng)度指標(biāo)
3.3 評(píng)價(jià)單元數(shù)據(jù)融合 第k（k=1~10）個(gè)風(fēng)險(xiǎn)強(qiáng)度指標(biāo)的風(fēng)險(xiǎn)強(qiáng)度值Sk的計(jì)算式為：
4 風(fēng)險(xiǎn)篩查層 4.1 高風(fēng)險(xiǎn)區(qū)域篩查 基于風(fēng)險(xiǎn)強(qiáng)度指標(biāo)評(píng)估燃?xì)夤艿佬孤╋L(fēng)險(xiǎn)的問(wèn)題屬于有監(jiān)督學(xué)習(xí)，高風(fēng)險(xiǎn)區(qū)域篩查的一般步驟為：①統(tǒng)計(jì)各風(fēng)險(xiǎn)強(qiáng)度上的事件數(shù)。②標(biāo)注標(biāo)簽：將泄漏樣本標(biāo)記為1，未泄漏樣本標(biāo)記為0。③樣本對(duì)照組數(shù)據(jù)集整理：正樣本集是發(fā)生管道泄漏事件的樣本集，負(fù)樣本集是未發(fā)生管道泄漏事件的樣本集。④模型訓(xùn)練：確定訓(xùn)練集、驗(yàn)證集和測(cè)試集。本研究中的正負(fù)樣本分布不均衡，生成訓(xùn)練集時(shí)采用的是組合-集成方法，反復(fù)多次抽取與正樣本數(shù)相同的負(fù)樣本，采用隨機(jī)森林模型建立泄漏風(fēng)險(xiǎn)二分類模型［6］。⑤泄漏風(fēng)險(xiǎn)評(píng)價(jià)：將二分類預(yù)測(cè)時(shí)的輸出概率作為評(píng)價(jià)單元對(duì)應(yīng)的綜合風(fēng)險(xiǎn)得分。綜合風(fēng)險(xiǎn)得分越高，則表明該評(píng)價(jià)單元的泄漏風(fēng)險(xiǎn)可能性越大。 4.2 各風(fēng)險(xiǎn)強(qiáng)度指標(biāo)的風(fēng)險(xiǎn)強(qiáng)度貢獻(xiàn)度
5 風(fēng)險(xiǎn)預(yù)警層 風(fēng)險(xiǎn)預(yù)警層并不是對(duì)某一風(fēng)險(xiǎn)強(qiáng)度指標(biāo)進(jìn)行預(yù)警，而是對(duì)評(píng)價(jià)單元的綜合風(fēng)險(xiǎn)得分進(jìn)行預(yù)警。在時(shí)序上實(shí)現(xiàn)快速迭代，基于風(fēng)險(xiǎn)時(shí)序數(shù)據(jù)對(duì)未來(lái)風(fēng)險(xiǎn)得分做出預(yù)判，實(shí)現(xiàn)風(fēng)險(xiǎn)預(yù)警。 6 結(jié)語(yǔ) 伴隨著機(jī)器學(xué)習(xí)技術(shù)的發(fā)展，利用內(nèi)外部的時(shí)空大數(shù)據(jù)，由宏觀至微觀，由預(yù)測(cè)篩查至預(yù)警的城鎮(zhèn)燃?xì)夤艿佬孤╋L(fēng)險(xiǎn)評(píng)價(jià)得以實(shí)現(xiàn)。軟檢測(cè)技術(shù)對(duì)大量容易測(cè)量的過(guò)程數(shù)據(jù)和外部數(shù)據(jù)進(jìn)行有監(jiān)督學(xué)習(xí)，突破傳統(tǒng)風(fēng)險(xiǎn)預(yù)測(cè)技術(shù)的小樣本、貧信息等限制，有效預(yù)測(cè)復(fù)雜的多因素耦合泄漏風(fēng)險(xiǎn)的發(fā)生和發(fā)展趨勢(shì)。 以燃?xì)夤艿佬孤╋L(fēng)險(xiǎn)影響因素為出發(fā)，建立包括風(fēng)險(xiǎn)預(yù)測(cè)層、風(fēng)險(xiǎn)篩查層、風(fēng)險(xiǎn)預(yù)警層的燃?xì)夤艿佬孤╋L(fēng)險(xiǎn)評(píng)價(jià)體系?？蓪?shí)現(xiàn)對(duì)評(píng)價(jià)時(shí)間段內(nèi)各風(fēng)險(xiǎn)強(qiáng)度指標(biāo)進(jìn)行評(píng)分，量化評(píng)價(jià)風(fēng)險(xiǎn)強(qiáng)度指標(biāo)。基于風(fēng)險(xiǎn)強(qiáng)度指標(biāo)評(píng)分，計(jì)算綜合風(fēng)險(xiǎn)得分，對(duì)高風(fēng)險(xiǎn)區(qū)域進(jìn)行篩查，并量化各風(fēng)險(xiǎn)強(qiáng)度指標(biāo)的風(fēng)險(xiǎn)強(qiáng)度貢獻(xiàn)度。在時(shí)序上實(shí)現(xiàn)快速迭代，對(duì)未來(lái)綜合風(fēng)險(xiǎn)得分做出預(yù)判，以實(shí)現(xiàn)風(fēng)險(xiǎn)預(yù)警。 本研究由北京市燃?xì)饧瘓F(tuán)有限責(zé)任公司“燃?xì)夤艿佬孤╋L(fēng)險(xiǎn)耦合軟檢測(cè)模型研究與應(yīng)用示范”項(xiàng)目提供資助。
參考文獻(xiàn)：
［1］陳利民，胡勇.淺析合肥燃?xì)饩用駪魞?nèi)安全風(fēng)險(xiǎn)和安檢措施［J］.城市燃?xì)猓?/span>2017（9）：24-28.
［2］王春雪，呂淑然.我國(guó)燃?xì)獗ㄑ芯楷F(xiàn)狀及熱點(diǎn)前沿可視化分析［J］.消防科學(xué)與技術(shù)，2016（11）：1620-1623. ［3］劉新想，楊永超，范建全.燃?xì)夤艿罊z測(cè)技術(shù)［J］.化工管理，2016（27）：95. ［4］賈雁群，李連星，葉茂，等.燃?xì)夤芫W(wǎng)泄漏事故原因分析與預(yù)防措施［J］.煤氣與熱力，2005（7）：59-61. ［5］馬濤，於孝春，沈士明.埋地管道風(fēng)險(xiǎn)評(píng)價(jià)單元?jiǎng)澐址椒ǎ?/span>C］//第一次全國(guó)工程風(fēng)險(xiǎn)分析技術(shù)學(xué)術(shù)會(huì)議論文集.南京：中國(guó)石化股份有限公司，2012：197-202. ［6］ANTON S D D，SINHA S，SCHOTTEN H D. Anomaly-based intrusion detection in industrial data with SVM and Random Forests［J］. Cryptography and Security，2019（7）：1-6.
