孫中旺，閆科樂，儲佳偉，劉衛(wèi)國，楊 磊，趙佳飛

(1.大連理工大學(xué) 能源與動力學(xué)院，遼寧 大連 116024；2.中國石油化工股份有限公司青島安全工程研究院，山東 青島 266071)

0 引言

隨著我國能源供給結(jié)構(gòu)的轉(zhuǎn)變以及對清潔能源的需求增加，天然氣在當(dāng)今化石燃料市場上占據(jù)優(yōu)勢，運輸天然氣最快捷、最可行的方式是通過管道[1-2]。從海底開采出的天然氣攜帶有沙礫、水及凝析油等雜質(zhì)，輸運過程中，這些雜質(zhì)的不斷堆積會形成固體沉積，隨著輸運里程的增加，沉積物的累計會加劇管道堵塞[3]。大量堵塞可能導(dǎo)致深水開發(fā)中管道破裂等災(zāi)難性事件，造成重大安全事故，威脅人身安全[4]。因此，當(dāng)管道形成堵塞時，需要快速精準(zhǔn)定位堵塞位置，采取合適的解堵措施，保障管道安全運行[5-6]。

目前，國內(nèi)外研究者已經(jīng)提出了一些不同的管道堵塞檢測方法和技術(shù)。Watanabe等提出了利用聲波技術(shù)檢測管道泄漏位點，通過在管道內(nèi)2個終端檢測到的信號的聲波脈沖響應(yīng)來估計管道的泄漏點[7-8]；隨后，Koyama等將該方法擴展到管道堵塞和堵塞程度的檢測上[9]；接著，Sharp等研究了聲波法檢測管道堵塞，這一方法理論上可行，但是聲波信號易受環(huán)境噪音等的干擾，如何提高抗干擾能力還是急需解決的問題[10]。Scott等提出的利用背壓法檢測天然氣管道堵塞，可以確定堵塞面積，但不能得到堵塞位置及長度[11]。為解決這一問題，Liu和Scott等將背壓法與特性曲線法結(jié)合，可對非完全堵塞位置進行估算[12]。

近年來，基于瞬態(tài)分析的壓力波脈沖法也被廣泛研究[13-14]。Gudmundsson等在輸水管道中利用快速關(guān)閉閥門產(chǎn)生水擊脈沖壓力波，通過壓力波反射信號可以檢測到管路內(nèi)由水合物沉積堵塞造成的管徑變化，并分析管內(nèi)流動狀態(tài)[15-16]。Sebastian開發(fā)的用于預(yù)測氣液兩相管道流動的壓力波傳播特性的模型表明，壓力波檢測堵塞的能力受到摩擦系數(shù)、管徑、管內(nèi)流速以及聲速的影響[17]。Adeleke等開發(fā)出一種快速定位堵塞位置以及準(zhǔn)確判斷堵塞程度的理論模型，其發(fā)現(xiàn)黏性阻力是影響準(zhǔn)確檢測堵塞程度的關(guān)鍵因素，但是對堵塞段定位及堵塞長度檢測影響不大[18-19]。為驗證壓力脈沖波法檢測管道堵塞的有效性，Chu等搭建了長為220 m的壓力波堵塞檢測實驗管路平臺，并模擬管道堵塞進行檢測，實驗結(jié)果證明該方法可有效地檢測長距離輸運管路堵塞位置、堵塞長度及堵塞面積[20-21]。

目前，對壓力脈沖波在復(fù)雜管道中的傳播特性及堵塞檢測研究很少，但這又是該技術(shù)現(xiàn)場化應(yīng)用不可避免的問題。因此為了解決這一問題，本文改造搭建一條長18.1 m的壓力波法復(fù)雜管道堵塞檢測實驗管路，利用閘板閥模擬堵塞，結(jié)合FFT諧波分析技術(shù)，對堵塞定位檢測，并檢驗該技術(shù)是否適用于低含水氣液混輸管道。以驗證利用壓力脈沖波法檢測復(fù)雜天然氣管道堵塞的可行性，為該技術(shù)在現(xiàn)場應(yīng)用提供數(shù)據(jù)依據(jù)和工程支撐。

1 實驗系統(tǒng)及實驗方法

1.1 實驗系統(tǒng)

圖1所示為復(fù)雜管路壓力波法堵塞檢測實驗裝置和實驗流程圖。該裝置主要由注氣系統(tǒng)、注液系統(tǒng)、壓力脈沖波產(chǎn)生系統(tǒng)、主體管路系統(tǒng)、單點堵塞系統(tǒng)、數(shù)據(jù)測量及采集系統(tǒng)6部分組成。

圖1 壓力波法堵塞檢測實驗裝置照片及流程Fig.1 Pictures and procedure of blockage detection experimental device with pressure wave method

注氣系統(tǒng)由空氣壓縮機及儲氣罐組成，利用空氣作為系統(tǒng)氣源。注液系統(tǒng)由儲液罐與蓄水池組成，作為系統(tǒng)水源。壓力脈沖波產(chǎn)生系統(tǒng)的核心部件是1個快速閉合電磁閥，當(dāng)其迅速工作時，在短時間內(nèi)將管道中的部分流體釋放到大氣中，并在管道入口產(chǎn)生1段壓力脈沖波，隨后，脈沖波向主體管道傳播。管路系統(tǒng)采用不銹鋼材料，管道直徑為100 mm，包含17個三通結(jié)構(gòu)，其中盲端1與電磁閥之間的管道長為8.25 m；盲端2與電磁閥之間的管道長為18.58 m。單點堵塞系統(tǒng)由距電磁閥8.89 m的閘板閥組成，可根據(jù)實驗要求通過調(diào)節(jié)閘板閥開度以模擬不同堵塞率。數(shù)據(jù)測量及采集系統(tǒng)由1個動態(tài)壓力傳感器D1、1個靜態(tài)壓力傳感器、1個高速數(shù)據(jù)采集卡及1個工業(yè)電腦組成，靜壓傳感器用于記錄管道內(nèi)的靜態(tài)壓力；動態(tài)壓力傳感器用于記錄由電磁閥快速閉合產(chǎn)生的瞬態(tài)動態(tài)壓力；利用高速數(shù)據(jù)采集卡獲取這些傳感器的實驗數(shù)據(jù)，并利用Labview數(shù)據(jù)采集程序進行同步記錄。

1.2 實驗步驟及方法

1) 檢查氣密性。閘板閥全開，保持管路通暢，關(guān)閉電磁閥。利用注氣系統(tǒng)向管路注氣并加壓到0.5 MPa，關(guān)閉單向閥，檢查管路氣密性。

2)加壓到實驗壓力。確認管路氣密性良好后，打開單向閥繼續(xù)注氣到實驗壓力1.5 MPa，隨后關(guān)閉注氣系統(tǒng)。

3)數(shù)據(jù)測量及采集。打開數(shù)據(jù)測量及采集系統(tǒng)，控制電磁閥快速開閉，產(chǎn)生壓力脈沖波，隨后打開數(shù)據(jù)測量，采集壓力波信號。當(dāng)壓力脈沖波沿管道傳播直至振幅衰減為0時，完成實驗。

4)氣液混輸管道堵塞檢測。先利用注液系統(tǒng)向管路注入實驗所需含水量，然后繼續(xù)向管路中注氣并加壓到0.5 MPa檢查氣密性，隨后重復(fù)上述實驗步驟。

2 實驗結(jié)果與分析

2.1 閘板閥全關(guān)(100%堵塞率)

首先將閘板閥全關(guān)，探究壓力脈沖波在復(fù)雜管道的傳播特性，實驗結(jié)果如圖2所示。圖2中為數(shù)據(jù)采集卡采集到的原始信號，橫坐標(biāo)為壓力波在管道內(nèi)的傳播時間，縱坐標(biāo)為動態(tài)壓力。由圖2可知，在復(fù)雜管道中采集到的是無規(guī)律的高低頻混疊壓力信號，隨著壓力波在管道內(nèi)的傳播，高頻信號逐漸衰減直至消失，且衰減速度高于低頻信號。但對于長距離輸運管道堵塞檢測而言，所得到的壓力波信號僅為單峰低頻信號[21]。

圖2 100%堵塞率時堵塞檢測原始信號Fig.2 Original signals of blockage detection at 100% blockage rate

本文實驗中，電磁閥發(fā)射的低頻負壓波沿著管道連續(xù)傳播，當(dāng)遇到三通結(jié)構(gòu)時會發(fā)生反射和透射。由于三通截面積由小變大，引起的反射波為正壓波，反射波的頻率發(fā)生變化，形成高頻波；引起的反射波逆向傳播到管道入口，再次反射。而三通引起的透射波仍為負壓波，沿管道繼續(xù)傳播，遇到盲端/堵塞時會引起反射負壓波。這些反射波逆向傳播時仍會經(jīng)過三通結(jié)構(gòu)并引起透射負壓波和反射正壓波，隨后繼續(xù)傳播到管道入口處再次反射。因此，壓力波在管道入口和堵塞點之間來回傳播反射，直至衰減為0。在管道入口至堵塞點之間分布著9個不同位置的三通結(jié)構(gòu)，靠近電磁閥的三通優(yōu)先產(chǎn)生反射波，隨著透射波的繼續(xù)傳播，后面的三通陸續(xù)產(chǎn)生反射波。

不同位置的三通產(chǎn)生的反射高頻正壓波由于時間差會相互疊加，增強振幅；盲端及堵塞點產(chǎn)生的反射低頻負壓波也會相互疊加，導(dǎo)致反射波振幅大于入射波；同時低頻壓力波和高頻壓力波的疊加，產(chǎn)生高低頻混疊信號。

壓力波在管道內(nèi)的多次反射和透射，使高頻信號迅速衰減；同時頻率越高，衰減系數(shù)越高，衰減也越快?？紤]到入射波振幅較小，隨著傳播過程的衰減振幅逐漸降低，壓力波能量下降，導(dǎo)致采集后期只有明顯的低頻壓力波信號。

由原始信號可知，信號采集受噪聲影響較小，主要干擾信號由管道內(nèi)三通結(jié)構(gòu)引起。為確定三通所產(chǎn)生的高頻信號的頻率范圍，采用FFT諧波變換法對原始信號進行不同頻率的低通濾波分析[22-23]，所得信號頻譜圖如圖3所示，圖3(a)，(b)，(c)分別代表100 Hz、150 Hz及200 Hz低通濾波變換后的信號頻譜。由圖3(a)可知，100 Hz濾波后的動壓信號曲線較為平滑，高頻信號均被過濾，因此所發(fā)射的入射波頻率低于100 Hz。對比3種不同的濾波信號，發(fā)現(xiàn)在100 Hz與150 Hz低通濾波后的信號曲線相差無幾，而200 Hz濾波變換后仍保留大量高頻信號，因此三通產(chǎn)生的反射高頻壓力波頻率集中在150～200 Hz之間。

圖3 100%堵塞率時FFT不同頻率低通濾波信號Fig.3 FFT low-pass filtered signals with different frequencies at 100% blocking rate

為定位堵塞點，對FFT 150 Hz低通濾波信號頻譜圖進行放大分析，如圖4所示。P1為動壓傳感器D1測得的入射波；P2為D1測得的盲端1反射波，由于盲端1與電磁閥之間的三通產(chǎn)生的反射正壓波與盲端1產(chǎn)生的反射負壓波相互疊加，導(dǎo)致P2振幅遠小于P1振幅；而P3為D1測得的堵塞點反射波，由于堵塞點產(chǎn)生的反射波與三通盲端產(chǎn)生的反射波疊加，使得P3振幅大于P1振幅。

圖4 100%堵塞率時FFT 150 Hz低通濾波信號放大圖Fig.4 Magnification of FFT 150 Hz low-pass filtered signals at 100% blocking rate

堵塞點位置計算公式如式(1)所示：

(1)

式中：L為堵塞點至D1間的計算距離,m；C為當(dāng)?shù)芈曀?，本文均?40 m/s；T1表示入射波時間,s；T′表示所計算的堵塞點反射波時間,s。

由式(1)得，盲端1距D1計算距離為7.66 m，而實際距離為7.75 m，因此，盲端1的定位誤差為1.11%。同理，根據(jù)P1與P3的時間差T1和T3，得到堵塞點至D1之間的計算距離為8.48 m，而實際距離為8.39 m，因此,堵塞點定位誤差為1.07%，誤差均在允許范圍內(nèi)。

由上述實驗結(jié)果可知，利用壓力波法定位復(fù)雜管道堵塞，外界干擾信號較小時，采用FFT諧波變換方法對采集到的高低頻混疊原始信號進行處理，同時在150 Hz濾波頻率下得到的頻譜信號更利于對堵塞進行定位。

2.2 閘板閥部分打開(不同堵塞率)

隨后進行堵塞率分別為93%，85%，70%和50%的4組堵塞檢測實驗，以驗證不同堵塞率下的堵塞定位效果及差異，實驗結(jié)果如圖5所示。實驗過程中，入射波遇到堵塞之后，會產(chǎn)生反射波和透射波，均為負壓波；透射波穿過堵塞段，沿管道繼續(xù)傳播到盲端2處，并反射。

圖5 不同堵塞率時堵塞檢測原始信號Fig.5 Original signals of blockage detection at different blockage rates

采用FFT 150 Hz低通濾波對不同堵塞率原始信號進行分析，結(jié)果如圖6所示。從圖6中可以看出在不同堵塞率時壓力波振幅衰減速率有差異，在堵塞率85%及93%時衰減速率較快，其中85%堵塞率衰減最快，而當(dāng)堵塞率70%及以下時，振幅衰減速率較慢。這是因為在堵塞百分比很高時，堵塞產(chǎn)生的入射波的反射率高，透射率低，此時反射波和透射波在管道內(nèi)會來回反射，高堵塞率增加反射次數(shù)，衰減快；同時閘板閥結(jié)構(gòu)會改變壓力波相位，導(dǎo)致入射波和反射波的相互疊加改變振幅，加快信號的衰減。

圖6 不同堵塞率時FFT 150 Hz低通濾波信號Fig.6 FFT 150 Hz low-pass filtered signals at different blocking rates

圖7為不同堵塞率時150 Hz低通濾波信號放大圖。圖7中，P1為入射波，P2為堵塞產(chǎn)生的反射波，P3為盲端2產(chǎn)生的反射波。堵塞及盲端2位置計算結(jié)果如表1所示。由表1可知，不同堵塞率下堵塞位置預(yù)測均較為準(zhǔn)確，最大誤差1.93%。盲端2位置預(yù)測在堵塞率85%時預(yù)測不夠準(zhǔn)確，誤差達到3.53%，可能與反射波相互疊加導(dǎo)致的振幅快速衰減有關(guān)，其他堵塞率時盲端2定位均較為準(zhǔn)確。

圖7 不同堵塞率時FFT 150 Hz低通濾波信號放大圖Fig.7 Magnification of FFT 150 Hz low-pass filtered signals at different blocking rates

表1 不同堵塞率時堵塞及盲端2定位結(jié)果及誤差Table 1 Location results and errors of blockage and blind flange 2 at different blocking rates

2.3 氣液混輸管道堵塞檢測

實際工程中，天然氣輸運管道內(nèi)含有少量的水，形成氣液混輸工況。為了檢驗壓力波堵塞檢測技術(shù)是否適用于氣液混輸管道，本文利用實驗管道進行1.5 MPa下，管道含水量8%，堵塞率100%的堵塞定位實驗，結(jié)果如圖8所示，其中，圖8(a)為原始信號，圖8(b)為FFT 150 Hz低通濾波信號。

由圖8可得，氣液混輸管道內(nèi)壓力波傳播特性與輸氣管道差異明顯。對于氣液混輸管道，入射波快速衰減，并且管內(nèi)長期處于負壓狀態(tài)，產(chǎn)生這種現(xiàn)象的原因主要有以下3方面：

圖8 100%堵塞率時氣液混輸管道堵塞檢測信號Fig.8 Detection signals of gas-liquid mixed transportation pipeline at 100% blockage rate

1)對于輸氣管道，假設(shè)管內(nèi)的任一界面均為理想界面，壓力波在管道內(nèi)傳播遇到理想界面后為全反射，不存在折射。而汽液混輸管道底部有一層液體，入射波到達汽水界面后不會完全反射，有一部分壓力波折射進水相，然后又會從水相折射回氣相，與原來的壓力波相互疊加，消耗了能量。

2)底部水層受到劇烈壓力波動后，產(chǎn)生強烈震蕩，導(dǎo)致部分壓力波能量轉(zhuǎn)化為水的動能，進一步消耗能量，加快振幅衰減。

3)堵塞產(chǎn)生的反射負壓波峰值高，與三通產(chǎn)生的正壓波疊加后仍然為負壓波。同時，在壓力波動瞬間，管道內(nèi)外強烈的壓差促進汽水界面的水層的激濺，分裂成為多個液滴，并且液滴迅速霧化，形成的小水珠分布在管內(nèi)各層空間中，由于水珠的拖曳力延緩管道壓力的恢復(fù)，造成管內(nèi)長期處于負壓狀態(tài)，隨著時間的推移，壓力恢復(fù)速率變慢，恢復(fù)時間增加。

對圖8(b)的信號進行放大，如圖9所示。其中，P1為入射波，P2為盲端1產(chǎn)生的反射波，P3為堵塞產(chǎn)生的反射波。根據(jù)堵塞計算理論得到堵塞段定位誤差為0.48%，盲端1的定位誤差為1.94%。說明壓力波堵塞檢測技術(shù)適用于氣液混輸管道。

圖9 100%堵塞率下氣液混輸管道FFT 150 Hz低通濾波信號放大圖Fig.9 Magnification of FFT 150 Hz low-pass filtered signals at 100% blocking rates of gas-liquid mixed transportation pipeline

3 結(jié)論

1)基于壓力波法定位復(fù)雜天然氣管道堵塞時，對采集到的無規(guī)律高低頻混疊動壓信號進行FFT諧波變換分析，發(fā)現(xiàn)三通導(dǎo)致的高頻信號的頻率集中在150～200 Hz之間，同時，外界干擾信號較小時，采用150 Hz低通濾波頻率所得到的頻譜信號最好，更利于對堵塞定位分析。

2)實驗結(jié)果表明，100%堵塞率下，堵塞定位誤差為1.07%；在堵塞率50%時，堵塞定位誤差最大達到1.93%，而在堵塞率85%時盲端2定位誤差最大，達到3.53%。另外對于含液率8%的氣液混輸管道，由于管內(nèi)液滴霧化的影響，所采集到的壓力波信號與輸氣管道有較大差異，此時，100%堵塞率時定位誤差為0.48%，亦可實現(xiàn)精準(zhǔn)定位。

3)實驗結(jié)果表明，壓力脈沖波法可以快速、精準(zhǔn)地檢測復(fù)雜管道堵塞，也適用于氣液混輸管道，具有較好的可靠性與實用性。但在工程現(xiàn)場，電磁環(huán)境更加復(fù)雜，信號干擾更明顯，下一步將提高傳感器抗干擾能力和信號去噪能力；同時，不同位置及不同種類的堵塞定位會有所差異，因此，爬/降坡管道以及流動管道堵塞定位，也是下一步研究重點。