結合三維圖形技術的綜合法泄漏檢測系統在中低壓輸水管道上的應用研究

朱智偉

（河南省水利勘測設計研究有限公司，鄭州450000）

0 引 言

在輸水管道運行過程中，由于管道老化、地質條件變化以及第三方破壞等原因，管道泄漏事故經常發生，造成停工停產、資源浪費、財產損失［1］，目前我國長距離輸水管道泄漏檢測應用及技術水準均不高。因此，管道運營單位迫切需要行之有效的技術手段，實時監控管道運行狀態，并在泄漏發生時第一時間發現并定位。

對管道泄漏檢測技術的研究已有幾十年的歷時，考慮到檢測的可靠性、可行性和經濟性大致可分為兩大類，一類是直接探測法，包含直接觀察法［2］、電纜法、管內檢測法、光纜法［3］、聲波法［4］等；一類是水力參數檢測法，包含水擊波法［5］、瞬態模型法［6］、質量平衡法［7］和壓力法［8］等。直接探測法基于物理方法（非軟件算法），通過現場探測器或傳感器探測并發送泄漏報警信號至調度室，顯示并報警。此類方法消耗人力、物力大，泄漏報警時間長，靈敏度差。水力參數檢測法是基于傳感器、數據采集系統、仿真軟件為一體的檢測方法是漏損檢測及管控的發展方向。

然而大部分長距離輸水管道運行壓力一般不高，當利用傳統的參數法作為檢測方法時，當管道發生泄漏時所產生的負壓波波動不易被安裝在管道上的壓力傳感器所捕捉到，易造成負壓波法泄漏檢測系統漏報。結合三維圖形技術的綜合法泄漏檢測系統是一種結合了多種方法于一體的系統，通過三維角度對壓力波動分析，以及基于模擬仿真的流量算法相結合，對泄漏進行判斷和定位。

通過對結合三維圖形技術的綜合法泄漏檢測系統在低壓輸水管道上的運用研究，對比分析了多種檢測方法，特別是傳統負壓波法、基于三維圖形的負壓波法及流量平衡法的區別。

1 水力參數泄漏檢測方法

水力參數檢測法是基于內部的泄漏檢測系統利用現場傳感器數據監測管道內部參數。輸入到泄漏檢測系統中的數據如：壓力、流量、黏度、密度、溫度等通過模型計算可推斷出是否存在輸送介質的泄漏。

（1）管道（流量）平衡法。流量平衡法是以計算一段時間內儀表測量的差值，當該差值持續一段時間的情況下，并與所設定的門限值長時間存在差值［9］。這種不平衡主要由一系列的時間周期或窗口來監視，從而檢測不同大小的管輸物的泄漏。

（2）實時瞬態模型法。簡稱為RTTM，泄漏檢測時通過水力學模型，實時模擬管道內流量、壓力和溫度。該模型基于質量、動量以及能量守恒定律。模型需要依據管道參數（長度、直徑、壁厚、拓撲線路、管道粗糙度、泵、閥門、設備位置等）以及所輸送介質的屬性（體積模量、黏度、密度）進行配置。經過不斷的測試表明，當現場情況達到準確的情況下，RTTM 泄漏檢測系統可以檢測到小至1%的泄漏。但問題在于是如何保證模型計算的準確性，以使其運算結果與管道內部正在發生的情況保持一致。

（3）壓力/流量監控法。壓力/流量監控法是最基礎的內部檢測法，該方法主要通過分析壓力和流量的變化，并輔助以設定的閾值，來確定管道是否發生泄漏。當管道發生泄漏時，管道壓力和流量都將發生改變，即管道壓力下降，流量差變大。壓力/流量監控法就是通過檢測管道的壓力和流量的綜合變化，再通過和設定的閾值進行比較，來確定管道是否發生泄漏。

（4）聲波/負壓波法。聲波/負壓波技術是利用管道內液體外泄后，與泄漏點位置發生振動撞擊而形成稀疏波的原理。泄漏造成管道泄漏點壓力下降，因而產生負壓波或稀疏波，分別向管道兩端傳播。基于此原理的泄漏檢測系統，需配置有高頻、高分辨率的檢測儀表檢測泄漏所產生的負壓波，并需具有較強大分析能力的軟件，以從含有噪音的信號中，提取出負壓波對泄漏進行判定和定位。

（5）統計分析法。統計分析系統是將校正后的體積平衡與序貫概率比檢驗SPRT結合使用以提供可靠的泄漏檢測。SPRT是一種計算概率的假設測試方法，用于確定泄漏（H1）和無泄漏（H0）的概率。SPRT 計算發生泄漏概率與無泄漏概率之比，并確定校正后的體積平衡是否以預定的概率增加。

統計分析系統將SPRT 應用到檢測管道出口和入口整體流量和壓力的變化。雖然每條管道的控制以及操作不盡相同，但當管道發生泄漏后，管道壓力和流量之間的關系將始終發生改變。例如，泄漏通常會導致管道壓力降低，并導致管道入口和出口的流量差變大。統計分析法泄漏檢測系統旨在識別這樣的變化模式，確定泄漏是基于對固定周期采樣的概率計算。雖然操作改變會導致管道中的流量和壓力發生波動，但從總體上來說，管網入口和出口的總流量差異可通過管網內部的管存變化修正后達到平衡狀態。如管網中發生泄漏，且SPRT 檢測到該泄漏引起的偏差，則系統無法保持平衡，SPRT 與模式識別相結合的統計分析法泄漏檢測系統可提供較高的系統可靠性，即最低的誤報警數量。

2 三維圖形技術的綜合法泄漏檢測

自20 世紀80年代開始，負壓波在管道泄漏檢測行業內已有應用，其原理是對管道輸送過程中，因泄漏所引起的壓力波動變化進行分析，而對泄漏以及泄漏位置做出判斷和定位；隨著近年來對系統所涉及的軟件分析能力以及硬件性能的不斷提高，負壓波法在使用過程中，也在不斷地更新和發展。

常用的水力參數檢測法各有優缺點，沒有一種技術適用于所有管道，使用一種以上技術的泄漏檢測系統能夠更廣泛地適用于管道的各種運行工況。［10］

2.1 傳統負壓波法

傳統負壓波法泄漏檢測系統，通過傳感器或自動控制系統提取管道入口和出口的壓力信號，通過檢測泄漏所產生的壓降波，并通過壓降波分別到達管道兩端的時間差以及波速，可計算出泄漏發生的位置。

負壓波檢測法在管道發生泄漏時，泄漏點的壓力將會下降并產生負壓波，該負壓波將沿管道向兩端傳播并最終由安裝在管道兩端的壓力傳感器檢測到該壓降波［11］。通過負壓波到達管道兩端的時間，可對發生的泄漏進行定位［12］。

傳統負壓波法泄漏檢測系統對發生迅速的即時泄漏有著較好的檢測效果；尤其是在中高壓的液體管道上；但針對低壓管道，由于管道整體運行壓力較小，泄漏所造成的瞬間壓降不明顯，很大可能性將掩蓋于管道運行的各種噪聲中；系統若無法在第一時間捕捉到泄漏特征，則將無法對該泄漏進行報警和定位。

由于判斷機制的原因，傳統負壓波在判斷管道是否發生泄漏時，通常使用的是對系統所檢測靈敏度的值進行手動設置，當管道內所發生的波動值大于所設的門限值時，系統即認定管道發生了泄漏并提供泄漏信息；即便是當管道處于瞬態情況下而帶來的壓力波動，當其波動值大于門限值時，傳統負壓波法依然會提供泄漏報警。

傳統負壓波法是通過負壓波的傳播速度、波傳播至管道兩端的時間差以及管長之間的關系，對所發生的泄漏進行定位；但基于傳統負壓波法對波動強度分析的不足，在較為復雜的低壓管道上，無法準確從包含大量噪音的原始信號中，分離出泄漏信號，將會導致系統產生大量誤報警的同時，也會對泄漏定位帶來較為明顯的誤差。

綜上所述，傳統負壓波法泄漏檢測系統存在如下幾個問題：①無法判斷管道運行狀態而導致系統在管道處于瞬態狀態下無法正常工作；②由于噪音以及信號的衰減，系統漏報可能性較大；③管道壓力較低，影響泄漏報警及定位。

2.2 基于三維圖形的負壓波法

結合三維圖形處理技術的綜合法泄漏檢測系統是從壓力、流量以及水力學變化多角度對管道運行狀態進行分析和判斷，是結合了多種方法于一體的泄漏檢測系統，可有效的適用于管道上的各種工況并檢測發生在其中的泄漏。

2.2.1 瞬態檢測

傳統負壓波法無法查看和分析整體動態壓力變化，忽略了管道全局的重要變化信息；綜合法泄漏檢測系統將實時監測管道壓力、流量的變化，監測壓力波在管道上傳播的整個過程以及管道上每一點的壓力變化情況，經過3 種綜合算法過濾掉數據中的噪音，生成以壓力強度、管長和時間為軸的三維圖形。通過對管道整體動態壓力變化進行分析，同時綜合泄漏點壓力變化的分析，全局化分析因泄漏引起的壓力波動，使得三維圖形處理技術可有效清晰的捕捉到發生在低壓管道上因泄漏所引起的微小壓力波動，從而對泄漏進行報警和定位。

圖1所示為基于三維圖形的負壓波法對泄漏判斷的運算過程。當管道處于瞬態下，壓力波動較為明顯且劇烈，綜合法泄漏檢測系統通過三維圖形法對比管道每一點的壓力波動與整條管道壓力波動的比率；當某一點的壓力波動比例大于整條管道壓力波動比例一定值時，系統即判斷管道發生變化。

圖1 三維圖形法運算過程Fig.1 Operation process of 3D graphics method

因此，系統將在管道所有狀態下工作，針對壓力波動進行分析，以此保證系統性能不受管道狀態的影響。圖2 所示為管道處于瞬態時系統檢測到泄漏所帶來的壓力波動原始數據。

圖2 管道處于瞬態系統檢測到泄漏（原始數據）Fig.2 The pipeline is in transient state and the system detects leakage（raw data）

圖3 為經過降噪濾波后檢測到泄漏，該圖標中分別展示了三類數據圖形，第一類為原始數據圖形，X軸為時間，Y軸為壓力強度；第二類為經過降噪濾波以及數學計算后的數據圖形，在該圖形中，X軸為時間，Y軸為對壓力波動進行二次求導處理后的波動值；第三類數據圖形，為系統依據所有管道信息生成的三維圖形，X軸為管長，Y軸為時間，Z軸為壓力波動強度，亮度越高壓力波動越強。

圖3 系統檢測到泄漏Fig.3 The system has detected leaks

2.2.2 降噪濾波保證數據的完整性

只有準確地捕捉到泄漏引發的壓力突降特征點，才能準確的檢測到泄漏并精確地判斷出壓力波傳播到上下游的時間差，從而提高泄漏檢測定位的精度。但由于在輸送過程中，例如電磁干擾以及泵所帶來的噪音等因素，使得系統所采集到的數據包含在大量的噪音之中，因此如何從噪音中進行有效準確的降噪濾波，對系統提供準確定位起著至關重要的作用。

對所采集數據的降噪濾波將直接影響系統的檢測性能以及使用感受。因此使用小波變換原理對數據進行降噪濾波。

（1）小波降噪原理［13］。在實際運行過程中，泄漏信號通常表現為低頻信號，而因瞬態所引起的信號或噪音信號，通常表現為高頻信號，小波降噪濾波的基本意義為：利用噪音與泄漏信號在各個坐標軸上的小波（變換模量極大值）譜具有不同的表現這一特征，將噪音波譜分量去除，然后利用小波變換重構算法，重構出原信號。如圖4 所示分別為對泄漏信號進行小波降噪濾波前后的對比趨勢圖，其數據坐標軸分別為X軸表示波動頻率，Y軸表示不同頻段下的壓力強度。

圖4 出入口壓力頻率分布圖Fig.4 The diagram of distribution between pressure and frequency of inlet and outlet

（2）小波變換原理［14］。小波變換由小波變換和小波級數兩部分組成。其中，小波變換的定義為：

式中：f(t)為平方可積的信號；ψ(t)為震蕩衰減且具有緊支集的函數稱為基本小波；基本小波中參數b起著平移作用，而參數a使得窗口函數的大小發生變化。

對于一個模擬信號f(t)，更主要的是把它展示成為小波級數，即

式中：Cj，k為小波變換系數；ψj，k(t)為由基本小波函數經平移和收縮得到的函數。

小波變換實質上是把函數f(t)展開為滿足一定條件的基本小波函數的線形組合。

2.3 流量平衡法泄漏檢測

流量平衡法是目前泄漏檢測方法中使用歷史最悠久，穩定性最高的方法；該方法對硬件的要求較為嚴格，傳統檢測方法的性能受硬件條件的制約，不能達到較好的效果；通過對管道輸送過程的了解和研究，針對性地提出多角度解決方案可有效地提高其可靠性和靈敏性。

2.3.1 傳統流量平衡法

依據泄漏發生時，管道流量差將變大的原理；傳統流量平衡法在于計算管道入口和出口的流量差，并對系統進行檢測閾值的設置；當流量差值大于所設閾值時，系統即提供泄漏報警以及相關信息。

傳統流量平衡法存在以下問題：

（1）誤報警。由于傳統流量平衡在判斷泄漏時，所采用的傳統的計算模式，因此所有數值皆為“絕對值”，無法將管道因瞬態或工況改變所帶來的異常變化考慮在內，造成誤報警。

（2）非獨立引擎。傳統流量平衡法需與壓力法配合使用，才可以檢測并對泄漏進行定位；無法作為一個獨立的引擎獨立對管道泄漏進行檢測。

（3）硬件制約。傳統流量平衡法在計算過程中，受硬件制約較大；即現場所提供信號的精度將大大影響系統的所檢測泄漏的大小。

2.3.2 智能流量平衡法

泄漏檢測系統中所使用的流量平衡法，屬于算法流量平衡，而非計算絕對值的流量平衡。在檢測過程中，系統對所采集的數據和管道狀態進行分析和判斷，并依據變化進行準確的泄漏檢測和定位［15］。

智能流量平衡法在計算管道流量差的同時，引入了SPRT運算技術和模式識別技術；系統根據所檢測管道提取現場數據，并實時計算管道發生泄漏的概率，SPRT計算公式為：

式中：τ(t)為泄漏檢測量；σi(t)為出口流量；σ0(t)為入口流量；Δσp(t)管存變化量。

通過假設檢驗算法判斷管道流量差是否增加，如下所示：

圖5 中如果大部分數據落在陰影區，H0大概率為真（無泄漏），如果平均值為1，假設檢驗無法判斷哪一個假定為真，因為H0的假定和H1的假定相等，如果平均值等于2，H1大概率為真（泄漏）。以時間間隔為t，概率計算還可表達為：

圖5 SPRT概率計算Fig.5 SPRT probability calculation

式中：λ(t)為管道在t秒時泄漏概率值；p1(t)為管道未發生泄漏的概率；p0(t)管道發生泄漏的概率。

式（4）也可表達為：

式中：λ(t- 1)為管道在t-1 秒時泄漏概率值；為管道常數；τ(t)為t秒時管道的流量差；m為管存修正量；Δm為系統可檢測的最小泄漏量。

通過系統中的概率統計分析法和模擬仿真法實時對管道管存的模擬計算，結合公式實時對所采集的數據進行計算，可準確計算出掃描周期內，管道是否發生泄漏。

3 綜合法泄漏檢測系統提高定位精度

3.1 提高定位精度

傳統泄漏檢測定位方法，即計算壓力波傳播至管道兩端的時間差與波速之間的關系，即可對泄漏位置進行定位，其公式表示為：

即：

式中：L1為泄漏點距離管道入口距離；L2為泄漏點距離管道出口距離；c為波速；t1為壓力波從泄漏點傳至管道入口所需的時間；t2為壓力波從泄漏點傳至管道出口所需的時間。

管道總長L=L1+L2，則泄漏定位計算公式為：

以上計算公式為傳統計算方法，在計算過程中，不可避免地存在著一定的盲目性和較大的計算量。特別是波速c事實上并非常數，而是隨外界條件發生變化，會加大計算過程中的誤差值而導致系統對定位誤差偏大。

如將壓力波在水中傳播的速度看作位置的函數c(x)，則壓力波從泄漏點傳播至管道兩端計量點的時間分別為：

由上式可得出：

式中：x0為泄漏點距離上游壓力傳感器的管道長度。

在式（11）中只有Δt是已知量，而x0和c(x)均為未知量，因此x是不可解得的。可通過Gauss Legendre公式［16］進行計算：

式中：n為節點數；Ak、xk為系數。

使用Gauss Legendre公式對式（11）進行處理，可得出：

通過以上公式的計算，可很大程度上減小計算量和盲目性，以此從數學的角度提高了計算的精確性。

3.2 減小誤報警

誤報警通常指不是由于直接發生的管輸物漏失、其他緊急工況或異常運行工況所引起的報警。為保障系統在正常運行過程中，不受操作影響而帶來更多的誤報警，使用負壓波方向判斷法，可區分出壓力波的傳播方向，如圖6所示。

圖6 方向判斷法Fig.6 Direction judgment method

通過系統算法以及數據分析可在低壓管道運行時提供較好檢測靈敏度；綜合法泄漏檢測系統中的方向判斷算法，可對管道內部的壓力波動進行甄別，區分出非泄漏引起的壓力波動，使得系統具有較好的可靠性。

通過方向判斷法，對不同原因造成的壓力波動進行判斷，區分出非泄漏造成的壓力波動，使系統在管道瞬態（啟停泵，閥門開關）過程中，不會出現誤報警且系統正常對發生在該瞬態下的泄漏進行有效的判斷。

4 結 論

水力參數泄漏檢測方法是輸水管道泄漏檢測發展的方向，結合三維圖形處理技術的綜合法泄漏檢測系統結合多種方法和算法，對于長距離中低壓輸水管道能夠有效、快速檢出泄漏并對泄漏點精確定位，是未來泄漏檢測的發展方向。

（1）由于中低壓管道，泄漏所造成的負壓波壓力較小以及衰減和噪音的影響，不易被傳感器所捕捉到，而存在漏報的可能性，因此傳統負壓波法檢測泄漏對于中低壓輸水管道適用性不強。

（2）經過本工程輸水管道驗證，結合了三維圖形的綜合法泄漏檢測系統可對中低壓管道進行有效的泄漏檢測并可提供穩定有效的泄漏報警。

（3）智能流量平衡法可修正彌補了傳統流量平衡無法定位，運算簡單的弊端，使得泄漏檢測系統具有靈敏度高、可靠性高的特征。

（4）基于三維圖形的負壓波法與智能流量平衡法相結合的綜合法泄漏檢測系統，可有效對管道上發生的大、小泄漏進行檢測，并提供準確的定位，是未來泄漏檢測主要選擇方向。 □