工業(yè)泵網(wǎng) 國內(nèi)對自然氣需求不斷增長促進了自然氣管道行業(yè)的快速發(fā)展。按照相關(guān)規(guī)劃，到2015年，中國自然氣管道總長將達到6萬km。自然氣管道中，球閥是其重要組成部門，自然氣管道上的球閥要求具有很好的密封性，不能存在內(nèi)漏。因此如何了解自然氣管道閥門內(nèi)漏特征并開發(fā)出相應(yīng)的內(nèi)漏檢測技術(shù)，實時、快速、精確檢測球閥是否發(fā)生內(nèi)漏，最大程度地減少泄漏，避免誘發(fā)重大安全事故的發(fā)生，將對自然氣管道安全、經(jīng)濟運行具有重要意義。
國外從20世紀80年代起就開展了閥門內(nèi)漏聲場特征研究，閥門內(nèi)漏檢測比較難題，但是聲發(fā)射檢測技術(shù)作為一種新型的無損檢測技術(shù)被證實是一種有效的檢測手段。目前海內(nèi)在閥門內(nèi)漏聲場特性研究方面剛剛起步，值得一提的是我國學者在閥門內(nèi)漏實驗和仿真研究方面也取得了一定的成果，戴光等以閘閥為研究對象，以水為實驗介質(zhì)，對閘閥內(nèi)漏過程進行了理論分析和實驗研究。張穎等以Lighthill氣動力聲學方程作為控制方程，采用有限差分法，對不同內(nèi)漏程度和壓力差下球閥氣體內(nèi)漏噴流聲場進行數(shù)值模擬，因為聲發(fā)射檢測技術(shù)具有檢測過程便捷、檢測范圍廣、檢測精度高等長處，因此更合用于自然氣管道閥門的內(nèi)漏檢測。
筆者分別以自然氣管道中常用的球閥和截止閥為研究對象，運用聲發(fā)射理論和信號分析與處理方法探討不同閥門類型和泄漏率等工況下閥門內(nèi)漏聲音在頻域范圍內(nèi)和時域范圍內(nèi)聲場特征，針對閥門發(fā)生內(nèi)漏時其聲場特征比較復雜，影響聲發(fā)射信號的分析與處理等難點題目。采取小波包分析方法對聲發(fā)射信號進行降噪處理，表征真實管道閥門內(nèi)漏聲場特征參數(shù)，為自然氣輸送行業(yè)管道閥門內(nèi)漏檢測和自然氣管道泄漏檢測提供理論和技術(shù)支持。
1球閥內(nèi)漏聲發(fā)射檢測基本理論
聲發(fā)射現(xiàn)象是指材料中局部區(qū)域應(yīng)力集中，快速開釋并產(chǎn)生瞬態(tài)彈性波的現(xiàn)象（AcousticEmission，簡稱AE）。聲發(fā)射信號攜帶著屬性變性、裂紋擴展和斷裂等信息傳播到物體表面，通過安裝在物體表面的聲發(fā)射檢測儀拾取聲發(fā)射信號，就可以對材料應(yīng)力特征進行分析。自然氣管道閥門內(nèi)漏過程中產(chǎn)生的聲信號嚴格意義上說不是聲發(fā)射現(xiàn)象，閥門在泄漏過程中不會產(chǎn)生能量，而是作為聲發(fā)射信號傳播的媒介。然而閥門在內(nèi)漏過程中因為上下壓差的影響在泄漏口處，泄漏氣體會產(chǎn)生渦流現(xiàn)象，渦流不僅會擾動氣體活動而且會產(chǎn)生不同頻率范圍的應(yīng)力波，該應(yīng)力波攜帶著泄漏信息傳播到閥門表面。因此閥門內(nèi)漏可以稱之為非嚴格意義上的聲發(fā)射現(xiàn)象。
球閥泄漏產(chǎn)生的聲發(fā)射信號在時域上呈現(xiàn)出一種連續(xù)變化的類型，且伴隨大量的環(huán)境干擾噪聲，可以表征球閥內(nèi)漏特征的聲發(fā)射信號經(jīng)常被沉沒在大量的干擾噪聲當中。為了從檢測信號中提取可以表征球閥內(nèi)漏特征的信息，小波包變換作為有效的信號處理方法已經(jīng)廣泛應(yīng)用于平穩(wěn)信號和非平穩(wěn)信號分析和處理中。小波變換采用二進變換的方法對每一層分解后的高頻段進行再分解，彌補了小波變換中高頻段局部性分解差的局限。因為管道泄漏聲發(fā)射信號為一寬頻信號，因此采用小波包變換能夠更正確的提取泄漏源信息。
小波包變換定義為：設(shè)｛hn｝n∈z是正交標準函數(shù)μ0（t）＝Φ（t）對應(yīng)的正交低通實系數(shù)濾波器，｛gn｝n∈z是正交小波函數(shù)μ1（t）＝φ（t）對應(yīng)的高通濾波器，其中g(shù)n＝（-1）nh1-n，則他們知足以下二標準方程和小波方程：
通過μ0、μ1、h、g在固定尺度下可定義一組稱為小波包的函數(shù)。為推廣二尺度方程，定義下列遞推關(guān)系：
式中：當n＝0時，μ0（t）＝Φ（t），μ1（t）＝φ（t）。以上定義的函數(shù)集合｛μn（t）｝m∈z為由μ0（t）＝Φ（t）所確定的小波包，由此，小波包｛μn（t）｝n∈z是包括尺度函數(shù)μ0（t）和小波母函數(shù)μ1（t）在內(nèi)的一個具有一定聯(lián)系的函數(shù)集合。
聲發(fā)射傳感器檢測信號x（t）在子空間Unj+1上的投影（小波包）系數(shù)可以表示為：
記x（t）在子空間U2nj和U2n+1j上的投影（小波包）系數(shù)分別為d2nj和d2n+1j，則可得：
小波包分解公式為：
小波包重構(gòu)公式為：
閥門內(nèi)漏過程中產(chǎn)生的聲發(fā)射信號在時域范圍內(nèi)呈現(xiàn)出一種連續(xù)變化的類型，為了在時域范圍內(nèi)了解不同泄漏率下聲場特征，可以采用聲發(fā)射信號均方根（AERMS）作為評價標準，例如一個聲發(fā)射信號含有N個樣本，x［0］，x［1］，…，x［n-1］均方根的方程如下：
同時通過頻域范圍內(nèi)的聲發(fā)射信號參數(shù)測量也可以表征球閥內(nèi)漏特征，例如功率譜密度（PSD）、頻率峰值、主頻率的相位值以及能量所集中的主要部分功率譜密度是指聲發(fā)射信號功率在頻譜圖上的分布情況，功率譜密度（PSD）可以由下式獲得：
式中，P［k］是功率譜密度，X［k］是一個聲發(fā)射信號的離散傅里葉變換，x［n］和T是樣本和周期。在以前的研究中得出球閥內(nèi)漏信號的聲功率直接近似等于泄漏信號的平均能量，這個關(guān)系可以簡單的定義為：
聲功率的指數(shù)與內(nèi)漏率的指數(shù)成線性關(guān)系，即：
設(shè)置帶通濾波帶通頻率為12.5～75kHz，重復以上實驗過程，并提取不同泄漏率下聲發(fā)射信號值，球閥和截止閥內(nèi)漏聲發(fā)射檢測實驗結(jié)果如圖6所示。通過實驗結(jié)果可以看出：內(nèi)漏過程產(chǎn)生的噴流與閥壁相互作用而產(chǎn)生高頻應(yīng)力波，該應(yīng)力波強度隨氣體泄漏率增加而增大，但是因為球閥和截止閥的內(nèi)部結(jié)構(gòu)不同，聲發(fā)射信號從泄漏源傳播到球閥外壁檢測傳感器探頭過程中傳播路徑和間隔不同，信號的衰減程度也不同，因此在統(tǒng)一泄漏率處兩者的AERMS有所不同。將球閥和截止閥不同泄漏率下檢測聲發(fā)射信號AERMS值帶入到等式（9），計算求得閥門內(nèi)漏率聲發(fā)射信號AERMS和體積泄漏率（Q）的雙對數(shù)關(guān)系曲線longAERMS＝blogQ+c，利用Matlab對二者進行擬合（見圖7），求得球閥和截止閥泄漏過程中雙對數(shù)系數(shù)為b＝0.3443、c＝0.916和b＝0.3713、c＝1.293。從圖7可以看出，2種結(jié)構(gòu)形式的閥門在內(nèi)漏過程中聲發(fā)射信號隨泄漏率變化規(guī)律基本保持一致，其雙對數(shù)曲線斜率基本一致，因為2種閥門結(jié)構(gòu)不同，聲發(fā)射信號在傳播過程中衰減程度不同，導致在統(tǒng)一泄漏率下對數(shù)值不同。
5結(jié)論
利用聲發(fā)射檢測技術(shù)對閥門氣體體積泄漏率檢測分析表明，閥門內(nèi)漏過程中產(chǎn)生的聲發(fā)射信號為一寬頻范圍，球閥和截止閥在密封圈劃傷寬度0.5mm、深度0.5mm時，其主要能量分布的頻率范圍為12.5～75kHz；泄漏過程中氣體泄漏產(chǎn)生的聲發(fā)射信號均方根（AERMS）與氣體泄漏率有關(guān)，且隨泄漏率增大而增大；通過實驗研究確定DN50球閥和截止閥密封圈處發(fā)生內(nèi)漏時，聲發(fā)射信號值與內(nèi)漏率的雙對數(shù)關(guān)系logAERMS=blogQ+c中系數(shù)分別為：b＝0.3443、c＝0.916和b＝0.3713、c＝1.293，因此可以采用聲發(fā)射技術(shù)對閥門內(nèi)漏進行無損檢測。
在分析閥門內(nèi)漏過程中氣體體積泄漏率與聲發(fā)射信號特征參數(shù)均方根（AERMS）函數(shù)關(guān)系時，選取了有代表性的球閥和截止閥作為研究對象，因為不同種類閥門結(jié)構(gòu)不同，在一定程度上不能完全分析出其對應(yīng)關(guān)系，故在以后的工作中選取更多的閥門如閘閥、針閥等作為研究對象，以便使分析結(jié)果更為正確。