基于負(fù)壓波信號消除的管道泄漏定位算法

謝文 徐自強(qiáng) 袁宇飛

【摘 要】對于帶壓輸氣管道，泄漏會產(chǎn)生負(fù)壓波（NPW），并沿管道傳播。利用壓電陶瓷傳感器檢測管道表面的NPW，可以分析獲得泄漏位置。然而，由于NPW信號的反射和持續(xù)時間較長，傳統(tǒng)的定位算法難以獲得高分辨率的精確定位結(jié)果，如傳統(tǒng)的時間反轉(zhuǎn)（TR）定位算法。文章提出了一種適用于管道泄漏檢測的TR定位方法。在該方法中，為了提高定位精度和減小由反射NPW引起的定位誤差，提出了一種TR自適應(yīng)對消方法，并通過實驗研究該方法的性能。實驗中，在一條長55.8 m的PVC管道上安裝6個壓電陶瓷傳感器，其中有2個手動控制的泄漏點用于記錄NPW信號。實驗結(jié)果表明，該方法能夠高分辨率、高精度地識別管道泄漏位置，這是傳統(tǒng)定位方法無法實現(xiàn)的。

【關(guān)鍵詞】時間反轉(zhuǎn);負(fù)壓波;定位;壓電陶瓷傳感器;管道泄漏

【中圖分類號】TP13 【文獻(xiàn)標(biāo)識碼】A 【文章編號】1674-0688（2020）12-0053-04

0 前言

結(jié)構(gòu)健康監(jiān)測（SHM）是研究結(jié)構(gòu)健康和耐久性的跨學(xué)科工程領(lǐng)域。它集遙感、智能材料、信號處理等功能于一體。SHM尤其適用于橋梁和水壩等遠(yuǎn)程監(jiān)控大型基礎(chǔ)設(shè)施系統(tǒng)，以及高姿態(tài)機(jī)械系統(tǒng)，如飛機(jī)、航天器、船舶、海上結(jié)構(gòu)物和管道，這些系統(tǒng)的性能至關(guān)重要，但現(xiàn)場監(jiān)控很困難，甚至無法監(jiān)控[1]。

由于管道泄漏每年在世界范圍內(nèi)造成大量的災(zāi)難性事故，引起廣泛關(guān)注[2-4]。因此，管道泄漏監(jiān)測成為SHM的一個重要研究領(lǐng)域。當(dāng)發(fā)生泄漏時，負(fù)壓波（NPW）從泄漏點向管道兩側(cè)傳播[5]。通過安裝在管道表面的傳感器記錄NPW信號，采用信號處理算法對NPW信號進(jìn)行處理，以定位泄漏點。例如，Jia Z G等人[6]利用負(fù)壓波引起的環(huán)向應(yīng)變開發(fā)了定位泄漏方法。Zhu J等人[7]通過計算負(fù)壓波的到達(dá)時間對泄漏進(jìn)行定位，并對定位精度進(jìn)行分析。

時間反演是基于互易性的技術(shù)[8]。當(dāng)傳感器記錄的信號在時域（或在頻域中相位共軛）中被時間反轉(zhuǎn)并在傳感器位置發(fā)回時，在激勵位置自適應(yīng)的時空聚焦信號[9]。時間反轉(zhuǎn)技術(shù)有兩種應(yīng)用方式：一種是所謂的物理時間反轉(zhuǎn)。時間反轉(zhuǎn)信號被重新傳輸?shù)轿锢斫橘|(zhì)中，然后信號將圍繞源的物理位置重新聚焦。這種效果對于許多需要將波的能量物理地集中在期望位置的應(yīng)用極具吸引力。另一種是計算時間反轉(zhuǎn)，信號通過計算被重新輻射到感興趣的領(lǐng)域，而不是在真實的介質(zhì)中實現(xiàn)。在計算過程中，時間反演信號與頻域中的傳遞函數(shù)（通常使用格林函數(shù)作為傳遞函數(shù)）相乘[10]。由于頻域的乘法等于時域的卷積，同樣的過程可以通過將時間反轉(zhuǎn)信號與信道沖激響應(yīng)卷積來完成。近年來，時間反轉(zhuǎn)技術(shù)在無損檢測（NDT）和結(jié)構(gòu)健康監(jiān)測（SHM）中得到廣泛應(yīng)用[11-13]。

然而，為了監(jiān)測整個管道，傳感器必須安裝在管道末端。端部反射的NPW信號導(dǎo)致定位精度下降[7]。而且，當(dāng)接收信號的持續(xù)時間較長時，反射信號與原入射信號的疊加，將使基于信號最大值的TR定位方法準(zhǔn)確率進(jìn)一步下降。本文提出了一種新的定位方法，該方法通過將時間反轉(zhuǎn)信號與特定系數(shù)疊加，以確定泄漏區(qū)域。由于反時信號與系數(shù)相乘，基于該方法的定位函數(shù)值隨觀測點的移動而迅速減小，準(zhǔn)確率得到提高。此外，由于反射的NPW可以互相抵消，因此提高了定位精度，特別是對于末端泄漏的定位。為了研究該方法在被動檢測應(yīng)用中的性能，而進(jìn)行了管道泄漏檢測實驗。實測結(jié)果表明，該方法能準(zhǔn)確地識別55.8 m管道中的兩個泄漏位置。

1 定位方法

針對管道泄漏實驗，將利用管道模型描述所提出的方法。對于帶壓輸氣管道，泄漏會產(chǎn)生NPW，并沿管道傳播。壓電陶瓷傳感器可以檢測出管道表面的NPW。我們假設(shè)沿管道使用N個傳感器，第n個傳感器位于rn，如圖1所示。我們進(jìn)一步假設(shè)泄漏發(fā)生在rL。

為了方便分析，我們在時域中描述了所提出的方法。假設(shè)泄漏產(chǎn)生一個NPW信號，泄漏發(fā)生時間為t=T。所有傳感器都是同步的。傳感器記錄的泄漏信號可以建模如下：

時間反轉(zhuǎn)信號y（t，rn，rL），公式（1）的時間反轉(zhuǎn)版本可以表示如下：

時間反轉(zhuǎn)過程實際上相當(dāng)于脈沖響應(yīng)的相關(guān)計算。因此，考慮從第n個傳感器到監(jiān)測域中任一點rk的后向傳播，一般觀測點rk的時間反轉(zhuǎn)信號可以表達(dá)如下：

其中，gc（rn，rk，t）是從傳感器n到點rk的計算信道脈沖響應(yīng)。

本方案對時間反演信號進(jìn)行如下處理：首先，乘以一個系數(shù)（-1）n;其次，對所有乘以系數(shù)后的時間反轉(zhuǎn)信號求和。相應(yīng)的結(jié)果可以寫成如下：

由于互易原理，并假設(shè)計算信道響應(yīng)函數(shù)與測量數(shù)據(jù)完全匹配，即gc（rn，rL，t）=gy（rn，rL，t）。然后，所有時間反轉(zhuǎn)的信號 f（t，rn，rk）（n=1…N）將集中在泄漏位置（rk=rs），且具有相似的波形。因此，所有乘以系數(shù)后的時間反轉(zhuǎn)信號將在rs處互相抵消，獲得接近零的振幅。

由于時間反轉(zhuǎn)信號在泄漏位置相互抵消，泄漏區(qū)域的信號能量比其他區(qū)域的低。因此，定位函數(shù)設(shè)計如下：

在定位圖中，泄漏位置的定位函數(shù)值將大于其他位置的定位函數(shù)值。

2 實驗結(jié)果

利用本檢測方法進(jìn)行管道泄漏試驗。全模型管道由總長55.8 m的PVC管段組成，如圖2所示。管道共有6段9.1 m直管段，采用10個90°彎頭接頭和5個0.2 m管段連接。將6個尺寸為15 mm×10 mm的壓電傳感器粘接在管道表面，具體位置見表1。使用兩個手動控制閥產(chǎn)生泄漏。閥門位置見表2。數(shù)據(jù)采集系統(tǒng)為NI-PXI-5105，由傳感器1的電壓信號觸發(fā)，觸發(fā)電平為-0.02 V，采樣速率為100 KS/s，壓電材料為APC855。

空氣通過壓縮機(jī)泵入管道，通過打開管道上的閥門產(chǎn)生泄漏事件;然后用6個壓電陶瓷傳感器檢測泄漏產(chǎn)生的NPW信號。圖3顯示來自泄漏L2的NPW信號。由于管道內(nèi)容物從泄漏點逸出，泄漏點處的內(nèi)部壓力顯著下降。管道內(nèi)容物同時從上游和下游向泄漏點移動。這種物質(zhì)流產(chǎn)生一個負(fù)壓波，從泄漏點向管道兩端傳播。管道內(nèi)壓力的降低伴隨著管道周長的收縮，導(dǎo)致管道壁上的應(yīng)變變化。

由于應(yīng)變的變化，直接安裝在管道壁上的壓電傳感器會產(chǎn)生相應(yīng)的電信號。如圖3所示，在發(fā)生泄漏之前，傳感器捕捉到的信號為0 V，因為在正常工作條件下，內(nèi)部壓力保持恒定。然后，由于NPW到達(dá)傳感器，產(chǎn)生向下的類脈沖波形。脈沖的初始下緣與NPW導(dǎo)致的內(nèi)壓降低有關(guān)。脈沖的上邊緣表示在不同基線壓力下的內(nèi)部壓力沉降。波形的負(fù)峰值表示NPW通過壓電傳感器的位置。最后，當(dāng)內(nèi)部壓力穩(wěn)定下來時，低頻信號使得壓電陶瓷傳感器輸出振幅返回到0 V。

在信號處理階段，利用傳統(tǒng)的時間反演定位算法和本文提出的算法分別對泄漏產(chǎn)生的NPW信號進(jìn)行處理。在重新構(gòu)建泄漏區(qū)域的過程中，我們用公式（6）設(shè)計信道脈沖響應(yīng)gc（rn，rk，t）[14]：

其中，δ（t-distancen，k）是狄拉克函數(shù)，distancen，k是rn和rk之間的距離，an，k是信道從點rn到定位點rk的衰減系數(shù)。根據(jù)公式（6），如果不同信道的衰減系數(shù)不同，衰減系數(shù)將對定位圖的定位函數(shù)值產(chǎn)生有害影響。為了消除這種影響，我們使各信道的衰減系數(shù)相同，即an，k=1。參數(shù)vg表示結(jié)構(gòu)中NPW的群速度。在本實驗中，vg=300 m/s[7]。

由于系統(tǒng)在電壓觸發(fā)模式下記錄了NPW信號，泄漏發(fā)生在采集系統(tǒng)開始獲取信號之前。因此，在泄漏發(fā)生時間和接收到的泄漏信號開始時間之間存在時間延遲，即第2部分的參數(shù)“T”。對接收到的NPW信號進(jìn)行時間反轉(zhuǎn)和重新發(fā)送后，時間反轉(zhuǎn)信號聚焦在t=-T時刻。

為了得到準(zhǔn)確的泄漏位置，采用基于時間反轉(zhuǎn)信號最大值的傳統(tǒng)TR定位方法對數(shù)據(jù)進(jìn)行處理。基于最大值的傳統(tǒng)TR定位算法的定位函數(shù)表示如下：

為了研究分辨率，在基于時間反轉(zhuǎn)信號最大值的傳統(tǒng)TR定位方法得到的結(jié)果中，對定位函數(shù)值進(jìn)行歸一化處理。如圖4所示，基于時間反轉(zhuǎn)信號最大值的傳統(tǒng)TR定位方法可以揭示泄漏位置。然而，由于信號的持續(xù)時間很長，時間反轉(zhuǎn)信號仍然在泄漏點以外的點上相互疊加。因此，輸出信號的最大值隨觀測點的移動而緩慢衰減，很多點的定位函數(shù)值都非常接近泄漏點的值。傳統(tǒng)的基于時間反轉(zhuǎn)信號最大值的TR定位方法所揭示的泄漏區(qū)域覆蓋范圍很大。從圖4可以看出，L1覆蓋層的-3 dB區(qū)域從長度等于0 m到長度等于19 m，L2的-3 dB區(qū)域從長度等于23 m到長度等于27 m。

所有時間反轉(zhuǎn)信號都會以相似的波形聚焦在泄漏位置。這意味著時間反轉(zhuǎn)信號可以通過公式（4）完全相互抵消。在泄漏位置，通過公式（5）對疊加信號進(jìn)行積分后。泄漏位置的定位函數(shù)值將大于其他位置的定位函數(shù)值。由于采用了時域積分，所以不需要計算無法得到的時延。如圖4所示，基于該方法繪制的圖很好地估計了泄漏位置。此外，該方法的分辨率優(yōu)于傳統(tǒng)的基于最大信號值的TR定位算法。如圖4所示，由本方案獲得的兩個泄漏的-3 dB區(qū)域長度均約為2 m。顯然，基于該方法得到的結(jié)果-3 dB區(qū)域比基于最大信號值的傳統(tǒng)TR定位算法要小。

由于管道兩端反射的信號，采用傳統(tǒng)的TR定位方法，兩次泄漏的定位結(jié)果分別為8.2 m和25 m。然而，在基于該方法的定位結(jié)果中，漏點分別位于7.2 m和24.8 m處，明顯高于傳統(tǒng)的TR定位方法。

表3列出了該方法與傳統(tǒng)TR定位方法的多次定位結(jié)果。采用本文方法，基于55.8 m的模型管道總長度，定位結(jié)果的最大誤差為位置預(yù)測的1.9%，各漏點結(jié)果的方差穩(wěn)定在5.4%以內(nèi)。基于傳統(tǒng)TR定位方法的定位誤差比基于該方法的定位誤差大。表3顯示，傳統(tǒng)的TR定位方法的最大誤差為4.4%。在模型管道上，NPW在0.05 s內(nèi)傳播15 m，這是NPW脈沖的典型持續(xù)時間，如圖3所示。由于入射和反射NPW的疊加，當(dāng)壓電傳感器與管道末端的距離小于7.5 m時，脈沖的負(fù)峰值將被延遲，這是典型持續(xù)時間內(nèi)NPW傳播距離的一半。在模型管道中，傳感器1和傳感器6位于NPW反射范圍內(nèi)。因此，采用傳統(tǒng)的TR定位方法，對泄漏L1的估計距離誤差要大于本方法。然而，在本方法中，傳感器1和傳感器6的反射NPW相互抵消。因此，反射的NPW對定位結(jié)果的影響很小。該方法的定位精度高于傳統(tǒng)的TR定位方法。

3 結(jié)語

傳統(tǒng)的TR定位方法在應(yīng)用于管道泄漏監(jiān)測等低頻信號被動檢測時，存在定位結(jié)果不準(zhǔn)確、分辨率低等問題。本文提出了一種新的管道泄漏時間反演定位方法。通過多信號對消和時間積分，該定位方法能夠準(zhǔn)確、高分辨率地定位泄漏源。將該定位方法應(yīng)用于含壓電陶瓷傳感器的模型輸氣管道泄漏檢測系統(tǒng)中。結(jié)果表明，該方法能提供高分辨率、高精度的55.8 m PVC管道泄漏定位圖。除管道泄漏監(jiān)測外，該方法還可用于其他被動監(jiān)測系統(tǒng)，如沖擊監(jiān)測系統(tǒng)。

參 考 文 獻(xiàn)

[1]Song G，Wang C，Wang B.Structural Health Mo-

nitoring（SHM）of Civil Structures[J].Appl. Sci，2017

（8）：789.

[2]Papadakis G A.Assessment of requirements on s-

afety management systems in EU regulations for the control of major hazard pipelines[J].Journal of Hazardous Materials，2000，78（1-3）：63-89.

[3]Hu J，Zhang L，Liang W.Detection of small leakage from long transportation pipeline with complex noise[J].Journal of Loss Prevention in the Process Industries，2011，24（4）：449-457.

[4]張學(xué)軍.基于差分進(jìn)化的輸油管道泄漏定位方法[J].信息技術(shù)，2012（8）：178-180.

[5]劉明亮.基于BP神經(jīng)網(wǎng)絡(luò)輸油泄漏監(jiān)測識別技術(shù)研究[J].信息技術(shù)，2006，3：49-51.

[6]Jia Z G，Ren L，Li H N，Ho S C，et al.Experimental study of pipeline leak detection based on hoop strain measurement[J].Structural Control and Health Monitoring，2015，22（5）：799-812.

[7]Zhu J，Ren L，Ho S C，et al.Gas pipeline leaka-

ge detection based on PZT sensors[J].Smart Ma-

ter.Struct，2017，26（2）：025022.

[8]Fink M.Time-reversal of ultrasonic fields-part I：basic principles[J].IEEE Trans. Ultrason.Ferroelectr.Freq.Control，1992，39（5）：555-566.

[9]Fink M.Time-reversed acoustics[J].Sci.Am，1997，

281（3）：91-97.

[10]Liu D，Kang G，Li L，et al.Electromagnetic time-

reversal imaging of a target in a cluttered environment[J].IEEE Transactions on Antennas and Pro-

pagation，2005，53（9）：3058-3066.

[11]Ing R K，Quieffin N，Catheline S，et al.In solid localization of finger impacts using acoustic time-

reversal process[J].Appl.Phys.Lett，2005，87（20）：204104.

[12]Huo L，Wang B，Chen D，Song G.Monitoring of Pre-Load on Rock Bolt Using Piezoceramic-Transducer Enabled Time Reversal Method[J].Sensors.2017，17（11）：2467.

[13]Liang Y，F(xiàn)eng Q，Li D.Loosening Monitoring of the Threaded Pipe Connection Using Time Reversal Technique and Piezoceramic Transducers[J].Sensors，2018，18（7）：2280.

[14]Cai J，Shi L，Yuan S.High spatial resolution loc-

alization for structural health monitoring based on virtual reversal[J].Smart Materials and Structures，2011，20（5）：901-904.