周 杰

（1.周口理工職業(yè)學(xué)院，河南 周口 466000；2.周口技師學(xué)院，河南 周口 466000）

由于目前鹽業(yè)化工企業(yè)主要采用鹽鹵作為生產(chǎn)原料，大多數(shù)鹽腔都處于地下，且鹽化廠的地理位置使得鹽鹵管道是鹽化廠收獲鹵水的主要方法，因此鹵水集中輸送成為提取鹽鹵的重要環(huán)節(jié)，在生產(chǎn)中常需要長距離運輸。然而隨著管道數(shù)增多、管道輸送距離拉長，以及基本管線處于露天狀態(tài)，導(dǎo)致管線遭遇破壞和自然環(huán)境的腐蝕影響，泄露的風(fēng)險和狀態(tài)與日增加。由于鹵水具有一定的腐蝕性和流動性，使得鹵水泄漏對周圍村莊和環(huán)境造成嚴重威脅，同時對企業(yè)造成巨大的經(jīng)濟損失[1]，因此鹽鹵管道的安全檢測至關(guān)重要。

我國鹽鹵管道檢測技術(shù)主要有人工檢測、半自動檢測、全自動化檢測3 種。其中，人工檢測速度較慢，周期較長，給企業(yè)帶來的損失大大增加；半自動檢測方法較為簡單，但魯棒性不強，采用物理量檢測極易受到周圍環(huán)境的干擾；全自動化檢測是未來的發(fā)展趨勢，主要分為流量、聲波、壓力、光纖、磁電傳感等[2]。

流量傳感型對流入流出量進行計量，但是受到管道形變和溫度影響，精確度不高，且無法準確定位到泄漏點；壓力傳感型對管道內(nèi)壓力信號進行處理，多采用均值濾波等進行消噪、提取特征值，與此同時也會大大降低泄漏信息的特征、影響定位精度；聲波傳感型以負壓波為主，通過管壁負壓波進行定位，由于長距離問題，負壓波在傳輸過程中衰減較為嚴重，對泄漏變化極為敏感，易受到干擾[3]。

目前負壓波檢測技術(shù)的精度受到管道、材料、尺寸、流量等影響，造成負壓波適用場景不能通用。本文提出一種基于CUSUM 算法的負壓波管道無損檢測方法，具有較強的自身抗噪性能，通過歷史數(shù)據(jù)預(yù)處理，解決了負壓波抗噪能力較差的問題，可靠性強，靈敏度高。該方法對管道泄漏檢測定位具有一定的工程實際意義。

1 負壓波法基本理論

1.1 負壓波檢測原理

當(dāng)鹵水管道出現(xiàn)裂紋、造成泄漏時，管道內(nèi)本來的穩(wěn)定流體將會沿著裂紋噴射向外擠壓，造成流體流失，使得管內(nèi)在泄漏點壓力降低，導(dǎo)致流體膨脹且密度降低。流體以水擊波的形式往兩端傳播、擴散、反射、衰減，由于負壓波屬于定頻信號，因此反射衰減較慢，沿管壁達到管道盡頭[4]。當(dāng)負壓波傳遞到首末兩端時，會在一段時間中產(chǎn)生壓力陡降，通過首末兩端采集的同步時鐘差，根據(jù)具體的材質(zhì)，計算出實際負壓波在管道中的傳播速度，計算出泄漏點的具體位置[5]。負壓波傳播速度一般在1 000 m/s 左右，負壓波法定位迅速，防止事故擴大，爭取了寶貴的時間，降低了損失[6]。負壓波法原理見圖1。

圖1 負壓波法原理

1.2 負壓波法采集信號與定位

管道總長度為L，故障距離PT1 為x，假設(shè)負壓波傳播速度為α，管內(nèi)流體速度為V，負壓波到達兩端時間差為Δt，則泄漏點到達PT1 距離的計算公式為

負壓波法采集的主要信號包括流量、壓力、溫度時域信號。根據(jù)負壓波的結(jié)構(gòu)特征與所采集的壓力、流量信號進行分析，將首末兩端所采集的壓力、流量信號的結(jié)構(gòu)特征與負壓波的結(jié)構(gòu)特征進行對比。當(dāng)滿足相應(yīng)結(jié)構(gòu)特征時，綜合管內(nèi)壓力、流速等，計算負壓波在管道中不同段的取值情況，通過時間差對泄漏點位置進行定位[7]。

2 CUSUM 算法基本理論

2.1 CUSUM 算法原理

累積和（Cumulative Sum，CUSUM） 是統(tǒng)計學(xué)中一種重要的方法，于1945 年在英國劍橋大學(xué)誕生[8]。其本質(zhì)就是凸顯觀測值與水平值偏移程度的控制圖。主要統(tǒng)計一個時序數(shù)列中過程平均值的微小變化和未來可變性，具有提前預(yù)測的能力。CUSUM算法基于序貫分析中的序貫概率比進行檢測，在抽樣時，不預(yù)先指定樣本總數(shù)，而且先抽取少量樣本，然后停止或者繼續(xù)抽樣，直到滿足條件，根據(jù)最終的結(jié)果狀態(tài)確定樣本數(shù)[9]。

2.2 基于CUSUM 算法的負壓波檢測

由于負壓波攜帶振動信號等泄露信息往兩端傳播，假設(shè)負壓波整個時間域上存在m個時間觀察點，可知m個時間觀察點中觀察到的負壓波波形數(shù)據(jù)分別為X1，X2，X3，…，Xm，相互獨立且服從分布N（μ1，δ2）。模型主要通過對時間域中的m進行檢測，觀察其均值分布，假設(shè)在時間觀察點n＜m后的負壓波數(shù)據(jù)發(fā)生突變，則負壓波數(shù)據(jù)量Xn+1，Xn+2，…，且服從分布N（μ2，δ2），要求μ1＜μ2。在發(fā)生突變時，假設(shè)檢驗公式為

式中：初始值S0=0；zi為當(dāng)前負壓波觀察值相對于負壓波歷史基線的一個偏移量；Xi為當(dāng)前PT1 的觀測值；μi為歷史基線均值；σi為歷史基線的標準差；Si為所需的負壓波預(yù)警統(tǒng)計量；Si-1為上一階段的負壓波狀態(tài)統(tǒng)計預(yù)警量；k為參考值，是CUSUM 負壓波模型中待定參數(shù)。

當(dāng)輸鹵管道出現(xiàn)泄漏時，通過以往歷史數(shù)據(jù)的長期檢測，當(dāng)出現(xiàn)變化時，立刻對變化較為明顯的負壓波做出精確預(yù)警，從而檢測負壓波到兩端PT1、PT2 的時刻，通過同步時間差，即可精準確定泄漏點位置。

3 負壓波檢測及識別

3.1 程序設(shè)計流程

基于CUSUM 算法的負壓波檢測具有良好的抗噪性，但在復(fù)雜工況下，由于現(xiàn)場大量噪聲夾雜時獲取原始信號較為困難，因此需要采用消噪方法對環(huán)境噪聲進行先一步處理。本程序采用小波去噪，結(jié)合固有工況，設(shè)置特定小波基db4，通過長期的正態(tài)分布記錄數(shù)據(jù)，保證正常工況下數(shù)據(jù)的穩(wěn)定性，通過MATLAB 程序設(shè)計，進一步驗證算法的準確性，MATLAB 程序中CUSUM 算法負壓波檢測流程見圖2。

圖2 MATLAB 程序中CUSUM 算法負壓波檢測流程圖

3.2 檢測判別測試及結(jié)果分析

根據(jù)MATLAB 程序中CUSUM 算法負壓波檢測流程，首先對采集的時域數(shù)據(jù)進行小波去噪，得到較為平滑的數(shù)據(jù)片段，即輸鹵管道實測負壓波波形數(shù)據(jù)片段，見圖3。

圖3 輸鹵管道實測負壓波波形數(shù)據(jù)片段

將實測負壓波波形數(shù)據(jù)通過CUSUM 算法檢測，判斷數(shù)據(jù)是否超出歷史數(shù)據(jù)基線裕度，當(dāng)達到泄漏條件時，即可采集雙端檢測的負壓波到達時間，通過距離公式進行計算，完成管道泄漏點的定位[10]。

通過輸鹵管道智能檢測系統(tǒng)數(shù)據(jù)采集，選擇STM32F767 為主控芯片，建立8 個數(shù)據(jù)采集點。在每段輸鹵管道上下游安裝終端，采集壓力信號和振動信號，并將信號進行存儲，采取圖2 中CUSUM算法負壓波檢測流程，通過調(diào)試，得到輸鹵管道檢測模型檢測判別正確率，見表1。

表1 輸鹵管道檢測模型檢測判別正確率

從測試結(jié)果來看，基于CUSUM 算法的負壓波檢測對于鹽化工輸鹵管道的泄漏檢測具有較高的準確性，該檢測模型能夠根據(jù)不同場景調(diào)整修改參數(shù)，從而滿足不同現(xiàn)場的要求，具有一定的工程實用價值。

4 結(jié)束語

本文從算法層面對負壓波檢測輸鹵管道工泄漏檢測提出一種基于CUSUM 算法的改進方法，對歷史數(shù)據(jù)進行預(yù)處理，通過負壓波進行泄漏點定位，具有良好的抗噪性能，同時為企業(yè)爭取了大量的時間。相對于傳統(tǒng)的負壓波檢測只能檢測到較大泄漏的缺點，該算法可以檢測到微小泄漏，對于不同場景工況下，通過設(shè)置閾值，正確率高于90%，滿足實際工況要求。