基于有限個SCADA監(jiān)測點的天然氣管道泄漏檢測

摘 要

摘要：探討了實時模型法、負壓波法、模式識別法3種天然氣管道泄漏檢測技術(shù)的原理、特點、適用情況及應(yīng)用案例。關(guān)鍵詞：天然氣管道；泄漏檢測；SCADA監(jiān)測點；實時模型法；負壓波法；模式識
摘要：探討了實時模型法、負壓波法、模式識別法3種天然氣管道泄漏檢測技術(shù)的原理、特點、適用情況及應(yīng)用案例。
關(guān)鍵詞：天然氣管道；泄漏檢測；SCADA監(jiān)測點；實時模型法；負壓波法；模式識別法
Detection of Natural Gas Pipeline Leakage Based on Limited SCADA Monitoring Points
WANG Shuai，YU Jianjun，YAN Mingqing
Abstract：The principle，characteristics，applicability and application case of three kinds of detection technologies of natural gas pipeline leakage，including real-time model method，negative pressure wave method and pattern recognition method are introduced.
Key words：natural gas pipeline；leakage detection；SCADA monitoring point；real-time model method；negative pressure wave method；pattern recognition method
1 概述
    近些年來，我國長輸天然氣管道和城市天然氣管網(wǎng)相繼設(shè)置了SCADA系統(tǒng)。它可對天然氣管網(wǎng)一些節(jié)點的壓力和流量進行自動監(jiān)測，并實時將監(jiān)測點的壓力和流量信號傳回控制中心，以此來監(jiān)控整個天然氣輸配系統(tǒng)的運行狀況[1]。但SCADA系統(tǒng)的監(jiān)測裝置價格很高，在管網(wǎng)中設(shè)置的監(jiān)測點數(shù)量往往有限(約每5km設(shè)置1個監(jiān)測點)，而且目前缺乏有效的故障診斷方法，在遇到管道突發(fā)泄漏或小泄漏量事故時，難以及時找出事故發(fā)生的具體位置。因此利用有限個SCADA監(jiān)測點來識別天然氣輸配系統(tǒng)的泄漏事故有著重要的實際應(yīng)用意義。
2 利用2個監(jiān)測點檢測長輸管道泄漏
長輸管道具有分支少、流量大、管網(wǎng)運行壓力高等特點。因此，可以僅在管道首末兩端分別設(shè)置1個SCADA遠傳裝置，傳送首末端固定傳感器采集到的管道實時數(shù)據(jù)，利用控制中心的計算機進行數(shù)據(jù)分析，可以判定管道是否存在泄漏并定位泄漏點。固定傳感器采集到的數(shù)據(jù)包括管道監(jiān)測點的壓力、流量，管道中天然氣的溫度和相關(guān)聲波參數(shù)(主要應(yīng)用于負壓波法)等。利用兩個SCADA監(jiān)測點來監(jiān)測復(fù)雜的管道運行，可以減少人力和財力支出，在管道泄漏檢測中有較大的發(fā)展空間。目前在該領(lǐng)域研究較深的方法有實時模型法和負壓波法。
2.1 實時模型法
    ① 實時模型法的泄漏檢測及定位原理
    用模型在線檢測管道的壓力和流量，并與壓力、流量的實測值比較來診斷泄漏事故，這是實時模型法的基本思想。實時模型法以管網(wǎng)SCADA系統(tǒng)的實測參數(shù)為邊界條件，以管道沿程的熱力、動力平衡微分方程和氣體連續(xù)性方程建立管道的實時模型。在沒有分支的管道首末端各設(shè)置1個SCADA監(jiān)測點，采集實時數(shù)據(jù)進行模擬后加以分析，可以檢測出管道是否存在泄漏。
要建立管道的實時模型，需要把管道劃分成無數(shù)個管長為出的微元管段。由熱力學方程、動力學方程、連續(xù)性方程和狀態(tài)方程[2～3]，同時考慮氣體壓縮因子方程Z=Z(p，T)，得到方程組，見式(1)。
 
式中cp——天然氣的比定壓熱容，J/(kg·K)
    Dj——焦耳-湯姆遜系數(shù)，K/Pa
    p——天然氣壓力，Pa
    dx——微元管道的長度，m
    T——天然氣溫度，K
    K——管道的傳熱系數(shù)，W/(m2·K)
    D——管道內(nèi)徑，m
    T0——管道埋深處的土壤溫度，K
    qm——天然氣的質(zhì)量流量，kg/s
    ρ——天然氣密度，kg/m3
    g——重力加速度，m/s2
    λ——摩擦阻力系數(shù)
    v——天然氣流速，m/s
    Z——壓縮因子
    根據(jù)式(1)建立管道模型，利用SCADA系統(tǒng)同一時刻對管道兩端實時采集到的數(shù)據(jù)，分別進行兩組管道數(shù)據(jù)的仿真：由起點到終點的模型仿真和由終點到起點的模型仿真。如果管道運行正常，則以上兩組仿真數(shù)據(jù)與實測數(shù)據(jù)偏差較小(允許相對誤差為5%[4])。但當管道發(fā)生泄漏后，仿真數(shù)據(jù)將不能與實測數(shù)據(jù)中的1個或幾個吻合，由此可以判定管道發(fā)生了泄漏。當天然氣在管道泄漏點處有相同的邊界條件時，由流體的連續(xù)性可知，泄漏點處的壓力和溫度也相同。因此當管道泄漏時，根據(jù)管道起點邊界條件繪制的壓力曲線1和根據(jù)終點邊界條件繪制的壓力曲線2將相交于L點，見圖1，此時可判定L點對應(yīng)的管道位置x0就是泄漏點位置。泄漏量為起點到終點模型仿真在泄漏點處的流量與終點到起點模型仿真在泄漏點處的流量之差。
 

    ② 實時模型法的考慮因素及適用情況
    a. 實時模型法中決定精度的根本是模型的仿真程度，模型的敏感參數(shù)為壓力，因此在建立模型時需要充分考慮管道的熱環(huán)境和沿程阻力及局部阻力對壓力的影響，以保證模型的準確性。文獻中的模型缺少對沿程阻力變化和管道熱環(huán)境變化的考慮[5]。
    b. 實時模型法只適用于高壓干管且在同一時刻只有1個泄漏點的情況。對同一時刻發(fā)生多點泄漏或在某一點發(fā)生泄漏后其擾動尚未傳到管道兩端時又有新泄漏發(fā)生，以及中壓管網(wǎng)的情況，實時模型法并不適用。
2.2 負壓波法
    ① 負壓波法的泄漏檢測及定位原理
    當流體輸送管道發(fā)生泄漏時，其泄漏部位立即產(chǎn)生物質(zhì)損失，這將引起泄漏區(qū)的流體密度減小，壓力降低。由于流體具有連續(xù)性，管道中的流體不會立即改變速度，泄漏點和相鄰兩側(cè)區(qū)域之間的壓力差導致流體從上、下游區(qū)域向泄漏區(qū)填充，從而引起泄漏區(qū)相鄰區(qū)域的密度減小、壓力降低。這種現(xiàn)象依次向泄漏區(qū)上、下游擴散，在水力學上稱為負壓波(或減壓波)，其傳播速度就是聲波在管道流體中的傳播速度[6～7]。
    在管道首末兩端安裝壓力傳感器，并由SCADA系統(tǒng)記錄傳感器傳送的數(shù)值。當泄漏發(fā)生時，會在管道內(nèi)產(chǎn)生負壓波動，從泄漏點向上、下游傳播，并以指數(shù)規(guī)律衰減，這種壓力波動和正常壓力波動大不相同，具有幾乎垂直的前緣。首末兩端的壓力傳感器分別在不同時刻捕捉到負壓波信號，從而判斷管道發(fā)生了泄漏，根據(jù)負壓波傳到管道兩端的時間差和負壓波波速可以進行泄漏點定位[8]。
   負壓波法進行泄漏點定位見圖2。負壓波從泄漏點到達管道首、末端的時間分別為t1、t2。令△t=t1-t2，則存在如下關(guān)系：
 

 
式中x——泄漏點距離管道首端的距離，m
    a——負壓波波速，m/s
    L——被測管道長度，m
    v——流體速度，m/s
    △t——負壓波從泄漏點到達管道首、末端的時間差，s
式(2)中，a一般大于1000m/s，v為1.5～3.0m/s，因此v通常忽略不計，式(2)可以簡化為式(3)。在式(3)的3個變量中，L已知，a與△t待求。
 
   影響負壓波法對泄漏點進行準確定位的兩個關(guān)鍵因素[8～9]為：負壓波傳播到首、末兩端傳感器的時間差的準確確定和管道內(nèi)負壓波波速的確定，即負壓波法定位的誤差主要來自a和△t。在傳統(tǒng)的常波速泄漏定位法中，負壓波在管道中的傳播速度被視為定值，通常為1000～1200m/s。實際上負壓波的傳播速度計算見式(4)，準確的a值有助于泄漏點的準確定位。
 
式中W——流體的體積彈性系數(shù)，Pa
    E——管材的彈性模量，Pa
    δ——管壁厚度，m
    C——與管道約束條件有關(guān)的修正系數(shù)
    要得出較準確的△t，需要管道首、末端的系統(tǒng)時間必須同步?；谪搲翰ㄟM行泄漏檢測和定位的后續(xù)方法主要有相關(guān)分析法、時間序列分析法和小波變換法等[9～11]。
    ② 負壓波法的局限性及適用情況
    a. 負壓波法的局限性在于要求泄漏的發(fā)生是快速和突發(fā)性的。該方法的抗干擾能力差，易誤報警。
    b. 要求首、末端傳感器的時間精度始終保持相同，否則定位會失真。負壓波法要求泄漏點產(chǎn)生突然壓降，對快速突發(fā)的泄漏和量較大的泄漏敏感，對于緩慢發(fā)生的泄漏或已經(jīng)發(fā)生的泄漏則失效。
2.3 應(yīng)用案例
    實時模型法和負壓波法只用兩個SCADA監(jiān)測點即可檢測長輸管道泄漏，目前已經(jīng)取得實驗的成功并進行了實際應(yīng)用。
    ① 四川省某輸氣管道長23.5km，管道外徑為426mm，壁厚為7mm，運行壓力為1.5MPa，起點溫度為305.3K，終點溫度為296.5K，在距管道起點5.6km處人為制造了1個泄漏點。采用實時模型法進行泄漏檢測和定位。監(jiān)控軟件采用BWRS狀態(tài)方程，18s后發(fā)出了泄漏報警信號，47s后將泄漏點定位于距管道起點5.38km處。
    ② 采用負壓波法對英國某天然氣管道進行泄漏檢測和定位。管道長220km，內(nèi)徑為560mm，每隔10km安裝1個截止閥，兩個截止閥之間設(shè)置5個壓力測點。該系統(tǒng)能在60s內(nèi)檢測出直徑為6.4mm的泄漏點，泄漏點定位精度為152m。
    ③ 清華大學與勝利油田油氣集輸公司聯(lián)合研制的長輸管道泄漏監(jiān)測系統(tǒng)采用負壓波法進行泄漏檢測和定位，利用小波變換法和相關(guān)分析法同時進行檢測。該泄漏監(jiān)測系統(tǒng)對32km管道進行了在線檢測，檢測結(jié)果顯示：最小泄漏量為5m3/h，泄漏點的定位精度約為管道長度的2%。
3 城市天然氣管網(wǎng)泄漏檢測的模式識別法
    ① 模式識別法的原理
城市天然氣管網(wǎng)系統(tǒng)具有壓力級制偏低，分支多，管網(wǎng)復(fù)雜，多成環(huán)狀等特點。從理論方面來看，無論一個管網(wǎng)有多大，2個監(jiān)測點均無法唯一確定泄漏點的具體位置，泄漏點的定位至少需要3個壓力監(jiān)測點。實時模型法和負壓波法均不能解決城市天然氣管網(wǎng)的泄漏檢測問題，筆者認為基于計算機人工神經(jīng)網(wǎng)絡(luò)的模式識別法能夠解決這個難題。
    模式識別法可以根據(jù)管網(wǎng)中設(shè)置的3個監(jiān)測點的壓力變化來診斷故障。SCADA監(jiān)測點在管網(wǎng)中的布置方式會影響診斷效果，因此盡量使監(jiān)測點在管網(wǎng)中均勻布置。首先要對天然氣管網(wǎng)各種代表陛的故障情況與3個監(jiān)測點壓力變化之間的關(guān)系進行學習，建立管網(wǎng)壓力、流量的特征矢量，形成故障模式的特征空間(模式庫)并儲存。一旦天然氣管網(wǎng)出現(xiàn)故障，將3個監(jiān)測點的數(shù)據(jù)與特征空間進行對比，分析出故障發(fā)生的位置和類別。充分利用模式識別技術(shù)在非線性映射方面的優(yōu)勢，彌補其他數(shù)學模型如多元回歸分析等在處理強非線性問題上的不足，提高故障診斷的準確性[12]。
    特征值的提取與選擇、訓練學習和分類識別是任何模式識別系統(tǒng)的3個核心問題[13]。一個功能較完善的識別系統(tǒng)在進行識別前，首先要進行學習，天然氣管網(wǎng)故障診斷模式識別系統(tǒng)的原理見圖3。
 

    一個測點的測量參數(shù)可構(gòu)成表征該測點所在管道運行狀態(tài)的特征矢量，天然氣管網(wǎng)不同部位發(fā)生故障，同—測點的特征矢量會有所不同?？梢酝ㄟ^試驗(廣義的試驗，如數(shù)值模擬等)得出管網(wǎng)發(fā)生故障時特征矢量的數(shù)值特性知識(模式分類)。在日常運行中將實時監(jiān)測得到的特征矢量值與特征空間進行對比，判斷監(jiān)測得到的特征矢量值屬于何種管網(wǎng)故障的特征矢量值范圍，從而做出對故障管道的診斷。特征空間和特征矢量的建立決定了模式識別的準確性和精確性，因此在建立特征矢量初期，要考慮管網(wǎng)壓力損失對特征矢量的影響，為了增加模式識別的可靠性，應(yīng)該選取適當?shù)奶卣魇噶烤S度和分量類別[14]。
    目前針對各種特征矢量提出的模式匹配方法的研究越來越深入。針對天然氣管網(wǎng)的特殊性，判別域代數(shù)界面方程法中的位勢函數(shù)法、統(tǒng)計判決法中的正態(tài)模式參數(shù)判別法以及人工神經(jīng)網(wǎng)絡(luò)法這3種方法比較適合天然氣管網(wǎng)故障診斷的模式識別[12]。采集分布較好、分布范圍較廣的實際數(shù)據(jù)非常困難，可以模擬部分數(shù)據(jù)以補充實際數(shù)據(jù)的不足，對于研究和實際應(yīng)用都非常必要。
   ② 應(yīng)用案例
   文獻[15]將模式識別法應(yīng)用于城市給水系統(tǒng)，并且全部采用水力分析來模擬所需數(shù)據(jù)，即采用管網(wǎng)局部破壞狀態(tài)下的水力分析來得到所需數(shù)據(jù)。雖然給水系統(tǒng)和天然氣系統(tǒng)在水力學工況上有很大差異，但是利用模式識別的方法和思路可以應(yīng)用到天然氣管網(wǎng)中。
    ③ 小結(jié)
    a. 模式識別的算法有很多種，算法復(fù)雜度的度量應(yīng)獨立于程序語言種類，這樣能可靠地進行復(fù)雜度比較，可提供數(shù)據(jù)內(nèi)在的固有的信息量[13]。
    b. 天然氣管網(wǎng)用氣的不均勻性要求對管網(wǎng)訓練時區(qū)分區(qū)域和時間段進行設(shè)置，考慮用氣量的變化得到管網(wǎng)各種代表性狀態(tài)的數(shù)據(jù)，建立相應(yīng)的特征空間。比如可將一天分為幾個時段，區(qū)分周末和工作日、冬季和夏季，有助于特征空間的訓練學習及故障診斷，保證數(shù)據(jù)庫的詳實，避免系統(tǒng)誤報。
    c. 本文僅考慮設(shè)置3個監(jiān)測點的情況，若增至4個監(jiān)測點，其中的1個作為備用，可大大提高故障診斷系統(tǒng)的可靠性。即當其中1個監(jiān)測點出現(xiàn)問題或維修時，其余3個仍可以進行故障診斷。
4 結(jié)語
    目前城市天然氣管網(wǎng)利用有限個SCADA測點進行檢測的方法仍處于試驗階段，該方法有廣闊的應(yīng)用前景，研究人員可以加大科研投入。在天然氣管道泄漏檢測中，多種方法相結(jié)合能保證運行人員安全和減少不必要的經(jīng)濟損失。為保證定位精確和搶修及時，在管道泄漏點定位時還應(yīng)與其他方法相結(jié)合，如利用手持式激光探測儀、激光雷達等[16]。
參考文獻：
[1] 肖文暉，劉亞斌，王思存.城市燃氣管網(wǎng)SCADA系統(tǒng)的開發(fā)[J].煤氣與熱力，2003，23(1)：20-22.
[2] 胡憶為.天然氣管道泄漏分析及動態(tài)處理技術(shù)[J].天然氣工業(yè)，2001，21(4)：98-100.
[3] 張紅兵，李長俊，彭善碧.輸氣管道故障診斷中的實時模型法[J].天然氣工業(yè)，2005，25(10)：103-105.
[4] 肖建蘭，呂保和，王明賢，等.氣體管道泄漏模型的研究進展[J].煤氣與熱力，2006，26(2)：7-9.
[5] 李軍.基于瞬態(tài)模擬的燃氣管道泄漏檢測與定位技術(shù)研究(碩士學位論文)[D].天津：天津城市建設(shè)學院，2008：25-29.
[6] 陳春剛，王毅，楊振坤.長輸油氣管道泄漏檢測技術(shù)綜述[J].石油與天然氣化工，2002，31(1)：52-54.
[7] 高松巍，劉博，楊理踐.城市燃氣管道泄漏檢測及定位技術(shù)的研究[J].城市燃氣，2008(9)：14-18.
[8] 李宏，王蘭蘭，張冬生，等.負壓波輸油管道泄漏檢測技術(shù)研究[J].科學技術(shù)與工程，2011，11(2)：340-342.
[9] 鄧鴻英，王毅.負壓波管道泄漏檢測與定位技術(shù)[J].油氣儲運，2003，22(7)：30-33.
[10] 李長俊，劉恩斌，梁黨國，等.基于負壓波原理的輸油管道泄漏檢測技術(shù)研究及應(yīng)用[C]∥.中國國際石油天然氣安全技術(shù)管理高層研討會論文集.重慶：科學技術(shù)與工程雜志社，2005：67-75.
[11] 李新建，鄧雄.基于負壓波和流量平衡的管道泄漏監(jiān)測系統(tǒng)研究[J].石油和化工設(shè)備，2009(1)：167-175.
[12] 嚴銘卿.燃氣輸配工程分析[M].北京：石油工業(yè)出版社，2007：245-255.
[13] 孫即祥.現(xiàn)代模式識別[M].2版.北京：高等教育出版社，2008：1-6.
[14] 孫即祥.模式識別[M].北京：國防工業(yè)出版社，2009：2.
[15] 梁建文，肖笛，趙新華，等.給水管網(wǎng)故障實時診斷方法[J].水利學報，2001(12)：40-46.
[16] 李軍，徐永生，玉建軍.燃氣管道泄漏檢測新技術(shù)[J].煤氣與熱力，2007，27(7)：56-59.
 
(本文作者：王帥1 玉建軍1 嚴銘卿2 1.天津城市建設(shè)學院 天津 300384；2.中國市政工程華北設(shè)計研究總院 天津 300074)