吳莉佳，胡 博，于潤(rùn)橋，郝延松

(無損檢測(cè)技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室(南昌航空大學(xué))，南昌 330063）

0 引言

埋地金屬管道在工業(yè)生產(chǎn)和工業(yè)運(yùn)輸中廣泛使用。由于其具有腐蝕泄漏甚至爆炸的潛在危險(xiǎn)性，對(duì)其腐蝕性檢測(cè)確保安全運(yùn)行意義重大。埋地金屬管道的腐蝕性檢測(cè)可分為管道內(nèi)腐蝕檢測(cè)和管道外腐蝕檢測(cè)。管道內(nèi)腐蝕檢測(cè)技術(shù)主要包括智能清管器超聲波法［1］和智能清管器漏磁法［2］。清管器超聲波法檢測(cè)準(zhǔn)確度高，但不適合在輸氣管線和含蠟高的輸油管線上進(jìn)行檢測(cè);清管器漏磁法檢測(cè)精度不如清管器超聲波法高，且對(duì)管道上的軸向裂縫檢測(cè)有一定難度。內(nèi)檢測(cè)方法均無法在非開挖不停輸?shù)臈l件下進(jìn)行檢測(cè)，且易出現(xiàn)內(nèi)置爬行器卡殼堵塞現(xiàn)象。管道外腐蝕檢測(cè)主要是通過對(duì)陰極保護(hù)系統(tǒng)進(jìn)行檢測(cè)［3-6］，而即使在良好的陰極保護(hù)下，管道有時(shí)也會(huì)發(fā)生腐蝕或泄漏。所以，僅通過對(duì)防腐層的檢測(cè)判斷管道的腐蝕狀況并不是最理想的方法。

瞬變電磁法是一種人工源的時(shí)間域電磁法，它是利用階躍波形電磁脈沖激發(fā)，利用不接地回線向地下發(fā)射一次場(chǎng)，在一次場(chǎng)的間歇期間，測(cè)量出地下介質(zhì)產(chǎn)生的感應(yīng)二次場(chǎng)隨時(shí)間的變化，來達(dá)到探測(cè)地下地質(zhì)體的目的［7］。該技術(shù)具有簡(jiǎn)單易行、信息豐富、精度較高、投資小、見效快等優(yōu)點(diǎn)。目前，瞬變電磁法在地質(zhì)礦產(chǎn)、石油、煤礦、考古、工程物探等行業(yè)中已被廣泛應(yīng)用［8-12］。瞬變電磁法檢測(cè)埋地金屬管道是一個(gè)比較新的應(yīng)用技術(shù):黃桂柏［13］在大慶油田開展了瞬變電磁法檢測(cè)管體腐蝕狀況的探索性試驗(yàn);高永才等［14］對(duì)瞬變電磁法金屬管道腐蝕檢測(cè)模型進(jìn)行了理論分析;石仁委［15］研究了基于瞬變磁場(chǎng)對(duì)理地管道腐蝕地面檢測(cè)技術(shù)，并在勝利油田進(jìn)行了工程檢測(cè)。但以上都是使用對(duì)電磁場(chǎng)的反演擬合方法來檢測(cè)埋地金屬管道，用瞬變電磁檢測(cè)埋地金屬管道的正演研究成果甚少。

在國(guó)外，R.U.B?rner等［16］提出了在磁偶極源和電偶極源下三維瞬變電磁場(chǎng)的數(shù)值計(jì)算方法;A.Swidinsky等［17］計(jì)算了埋在層狀地質(zhì)下電阻薄板的三維電磁場(chǎng)響應(yīng)。在國(guó)內(nèi)，閆述等［18］在時(shí)間域?qū)ω?fù)階躍脈沖激發(fā)的二維瞬態(tài)場(chǎng)進(jìn)行了數(shù)值分析;章飛亮［19］研究了考慮激電效應(yīng)的瞬變電磁場(chǎng)正演理論，并探討了不同模型下激電效應(yīng)對(duì)瞬變電磁響應(yīng)的影響規(guī)律;徐凌華等［20］用有限單元法計(jì)算了二維大地電磁響應(yīng)。由于三維正演問題本身的復(fù)雜性，以及計(jì)算量、計(jì)算速度、存儲(chǔ)量的限制，目前二維反演方法仍是實(shí)際資料定量解釋的主要手段。正演問題作為瞬變電磁檢測(cè)埋地金屬管道的理論基礎(chǔ)和研究瞬變電磁場(chǎng)響應(yīng)特征的基本手段，對(duì)其進(jìn)行深入研究有著重要意義。

本研究利用ANSYS有限元仿真軟件建立二維模型，對(duì)比計(jì)算地下埋有金屬管道后地面磁場(chǎng)分布，計(jì)算不同腐蝕程度的埋地金屬管道地面磁場(chǎng)分布;從理論上驗(yàn)證瞬變電磁法檢測(cè)埋地金屬管道腐蝕失效的可行性，分析瞬變電磁法對(duì)埋地金屬管道腐蝕的探測(cè)能力。

1 模型參數(shù)

埋地金屬管道模型瞬變電磁場(chǎng)仿真的實(shí)際問題是三維的，但三維計(jì)算復(fù)雜、計(jì)算量大、模擬難度非常大。鑒于本研究?jī)H考察某個(gè)檢測(cè)點(diǎn)剖面的磁場(chǎng)，設(shè)管道中心軸方向上的磁場(chǎng)不變，從而將三維問題簡(jiǎn)化為二維模型(圖1)。二維模型中，埋地金屬管道在某個(gè)檢測(cè)點(diǎn)的剖面簡(jiǎn)化為圓環(huán)，管壁缺陷簡(jiǎn)化為環(huán)狀缺陷［21］。分別模擬計(jì)算壁厚為10、8、4、2 mm的埋地金屬管道的瞬變電磁場(chǎng)。激勵(lì)線圈簡(jiǎn)化為等效磁偶極形式的環(huán)向電流。地表簡(jiǎn)化為水平地面。管道內(nèi)填充空氣，地表上方填充空氣，地表下方與管道之間填充土壤。

圖1 埋地金屬管道示意圖Fig.1 Diagram of buried metal pipeline

用ANSYS模擬埋地金屬管道模型的瞬變電磁場(chǎng)時(shí)，整個(gè)模型采用基于節(jié)點(diǎn)的磁場(chǎng)單元plane53建立二維模型，設(shè)置4種材料:空氣、管道、線圈、土壤。4種材料的屬性見表1。由于實(shí)際問題屬于開域問題，但有限元仿真只適用于封閉的邊界區(qū)域電磁場(chǎng)計(jì)算，所以仿真中需要對(duì)開域情況做一定處理。本模型對(duì)最外層所有空氣施加磁力線平行邊界條件，用來模擬無限遠(yuǎn)場(chǎng)狀態(tài)，來滿足電磁場(chǎng)的有限元計(jì)算的要求。模型計(jì)算范圍設(shè)定為2 m×2.3 m的長(zhǎng)方形。其中，上半部分2 m×0.3 m的區(qū)域?yàn)榭諝?，下半部? m×2 m的區(qū)域?yàn)槁裼薪饘俟艿赖耐寥?，土壤覆蓋的范圍達(dá)到管道外徑的3～5倍，可以認(rèn)為邊界的影響很小。模型原點(diǎn)設(shè)在管道橫截面的中心，設(shè)管道軸心為y軸，橫截面為x-z面。管道的埋深即管道橫截面的中心距地面1 m，管道內(nèi)徑為18 mm。激勵(lì)線圈水平放置于地表面，為1 m×1 m方形線圈，線圈橫截面積為10 mm×20 mm，共30匝，填充率為0.95。激勵(lì)信號(hào)加載在線圈上，激勵(lì)信號(hào)電壓有效值為12 V。當(dāng)t=0時(shí)激勵(lì)信號(hào)關(guān)斷，關(guān)斷時(shí)間為10 μs。線圈耦合自由度CURR(電流)，模型計(jì)算求解后可獲得其磁場(chǎng)強(qiáng)度及磁力線圖等結(jié)果。仿真模型網(wǎng)格劃分等級(jí)為2(可選1～10，1最細(xì)，10最粗)，且采用自由劃分方式。

表1 材料的參數(shù)表Table 1 Parameters of materials

2 結(jié)果與討論

2.1 有無管道的瞬變電磁場(chǎng)對(duì)比

利用ANSYS仿真軟件中POST1通用后處理器，可以對(duì)某一模型查看任意時(shí)刻的磁力線分布圖。圖2為信號(hào)關(guān)斷后t=0.18 s時(shí)地下沒有金屬管道的磁力線圖，圖中模型分為上、下2部分，上半部分為空氣，下半部分為土壤。交界處為地表面，激勵(lì)線圈水平放置在于地表面，在圖中為上下部分交界處的2個(gè)小矩形。圖3為信號(hào)關(guān)斷后t=0.18 s時(shí)金屬管道埋深1 m、壁厚8 mm的磁力線分布圖，圖右側(cè)A～I(xiàn)的數(shù)值從外至里依次表示圖中9條磁力線上磁場(chǎng)強(qiáng)度的大小。

圖2、圖3中最里層磁力線數(shù)值最大，最外層磁力線數(shù)值最小，說明土壤中磁場(chǎng)強(qiáng)度從里至外逐漸減小。從磁力線分布分析，圖2中地下沒有金屬管道，磁力線均勻分布在土壤中;圖3中金屬管道埋深1 m、壁厚8 mm的9條磁力線中有6條集中在管壁上。說明當(dāng)?shù)叵聸]有金屬管道時(shí)土壤中磁場(chǎng)強(qiáng)度依次均勻減弱，當(dāng)?shù)叵侣裼薪饘俟艿罆r(shí)管壁上磁場(chǎng)強(qiáng)度較強(qiáng)，土壤中磁場(chǎng)強(qiáng)度較弱。從數(shù)值上分析，地下埋有金屬管道的磁場(chǎng)強(qiáng)度比地下沒有金屬管道的磁場(chǎng)強(qiáng)度強(qiáng)約3倍。綜上得出:信號(hào)關(guān)斷后，地下沒有金屬管道時(shí)磁場(chǎng)較弱且由內(nèi)至外均勻減弱;地下埋有金屬管道時(shí)埋地金屬管道上磁場(chǎng)較強(qiáng)，土壤中磁場(chǎng)較弱。

圖2 地下無管道的磁力線Fig.2 Magnetic lines of no pipeline in the underground

2.2 不同時(shí)刻的瞬變電磁場(chǎng)

為了查看在激勵(lì)信號(hào)關(guān)斷后磁場(chǎng)變化情況，實(shí)驗(yàn)利用ANSYS仿真軟件中POST 1通用后處理器查看埋地金屬管道壁厚8 mm的模型中不同時(shí)刻磁力線分布圖。圖4為埋深1 m壁厚8 mm的管道在激勵(lì)信號(hào)關(guān)斷后 t分別為 0.02、0.03、0.04、0.06 s磁力線分布圖，圖中共 9 條磁力線，圖右側(cè)A～I(xiàn)的數(shù)值從外至里依次表示圖中9條磁力線上磁場(chǎng)強(qiáng)度的大小由圖可見，信號(hào)關(guān)斷后，隨時(shí)間線圈周圍的磁力線隨時(shí)間增加逐漸減少，埋地金屬管壁上磁力線隨時(shí)間增加逐漸增多。說明信號(hào)關(guān)斷后，線圈周圍的磁場(chǎng)強(qiáng)度隨時(shí)間逐漸減弱，管壁上的磁場(chǎng)強(qiáng)度隨時(shí)間逐漸增強(qiáng)。

圖3 管道埋深1 m壁厚8 mm的磁力線Fig.3 Magnetic lines of 8 mm thickness and 1 m in buried depth pipeline

2.3 不同腐蝕程度的瞬變電磁場(chǎng)

管道腐蝕后壁厚發(fā)生變化，腐蝕程度不同即壁厚不同。為比較激勵(lì)信號(hào)關(guān)斷后管道腐蝕程度不同的磁場(chǎng)強(qiáng)度，對(duì)壁厚為10、8、4、2 mm埋地金屬管道分別建立4個(gè)模型，并計(jì)算激勵(lì)信號(hào)關(guān)斷后管道中心z軸磁場(chǎng)強(qiáng)度(Hz)的大小。計(jì)算時(shí)間從信號(hào)關(guān)斷后 0.000 01 s到 0.180 01 s，歷時(shí)0.18 s，每 100 μs取一個(gè)值，共取1800 個(gè)值。

利用POST26時(shí)間歷程后處理器將各模型計(jì)算出的管道中心z軸磁場(chǎng)強(qiáng)度(Hz)的大小列表，再用MATLAB畫出圖像，見圖5。圖中橫坐標(biāo)為時(shí)間的自然對(duì)數(shù);縱坐標(biāo)為管道中心z軸磁場(chǎng)強(qiáng)度(Hz)幅值的自然對(duì)數(shù)。圖中拐點(diǎn)前半部分為負(fù)值取反后所得，拐點(diǎn)后半部分為正值直接得出。從圖中可以看出，4條曲線尾支明顯分離。說明管內(nèi)徑相同時(shí)，壁厚不同的脈沖瞬變響應(yīng)曲線是可分的。

圖4 磁力線分布圖Fig.4 Magnetic force line distribution

表2 不同壁厚的磁場(chǎng)強(qiáng)度(Hz)Table 2 Magnetic field intensity(Hz)of various metal pipeline wall thicknesses

圖5 管內(nèi)徑18 mm不同壁厚仿真瞬變響應(yīng)幅值對(duì)比Fig.5 Amplitude comparison of TEM field response with various metal pipeline wall thickness of 18 mm inner diameter

將t=0.18 s時(shí)各模型的磁場(chǎng)強(qiáng)度計(jì)算結(jié)果列表，見表2。由表可見，各模型磁場(chǎng)強(qiáng)度相差1E-3 A/m以上。說明當(dāng)線圈中心點(diǎn)磁場(chǎng)強(qiáng)度差值大于1E-3 A/m時(shí)，能在激勵(lì)信號(hào)關(guān)斷0.18 s后檢測(cè)出管道埋深在1 m以內(nèi)，壁厚相差大于2 mm的埋地金屬管道。

3 結(jié)論

1)地下沒有金屬管道時(shí)，土壤中磁場(chǎng)較弱且由內(nèi)至外均勻減弱;地下埋有金屬管道時(shí)，埋地金屬管道管壁上磁場(chǎng)較強(qiáng)。

2)激勵(lì)信號(hào)關(guān)斷后，隨時(shí)間增加，激勵(lì)線圈周圍磁場(chǎng)強(qiáng)度逐漸減弱，埋地金屬管道管壁上的磁場(chǎng)強(qiáng)度逐漸增強(qiáng)。

3)當(dāng)金屬管道埋深一定、腐蝕程度(即壁厚)不同時(shí)，磁場(chǎng)強(qiáng)度曲線尾支明顯分離，且管壁越厚磁場(chǎng)強(qiáng)度越大。

4)當(dāng)線圈中心點(diǎn)磁場(chǎng)強(qiáng)度大于1E-3 A/m時(shí)，能在激勵(lì)信號(hào)關(guān)斷0.18 s后檢測(cè)出管道埋深在1 m以內(nèi)，壁厚相差大于2 mm的埋地金屬管道。

［1］Hartmut G.Ultrasonic or MFL inspection:which technology is better for You?［J］.Pipeline ＆ Gas Journal，2003，30(10):342-371.

［2］Vijay B，Lynann C.Residual magnetic flux leakage:A possible tool for studying pipeline defects［J］.Journal of Nondestructive Evaluation，2003，22(4):117-125.

［3］Segall S M，Reid R G.Use of an Integrated CIPS/DCVG Survey in the ECDA Process［C］//NACE International.Corrosion 2006.San Diego Ca:NACE International，2006:12-16.

［4］鄭耿鉑，余越泉，王健健，等.PCM管道電流檢測(cè)系統(tǒng)介紹及應(yīng)用［J］.腐蝕與防護(hù)，2002，23(1):21-23.

［5］董寶山.變頻—選頻技術(shù)在防腐層檢測(cè)中的應(yīng)用［J］.石油規(guī)劃設(shè)計(jì)，2002，13(1):38-39.

［6］周琰，靳世久，孫墨，等.埋地管道防腐層缺陷DCVG檢測(cè)技術(shù)研究及應(yīng)用［J］.管道技術(shù)與設(shè)備，2001，9(5):38-40.

［7］李貅.瞬變電磁測(cè)深的理論與應(yīng)用［M］.西安:陜西科學(xué)技術(shù)出版社，2002:1-4.

［8］牛之璉.時(shí)間域電磁法原理［M］.長(zhǎng)沙:中南大學(xué)出版社，2007:1-6.

［9］吳克峰，樊建華.四維瞬變電磁法油藏監(jiān)測(cè)正演模擬［J］.內(nèi)蒙古石油化工，2006(5):179-181.

［10］薛國(guó)強(qiáng)，秦克章，李貅，等.西藏沙讓特大型鉬礦的發(fā)現(xiàn)與瞬變電磁法探查［J］.地球物理學(xué)進(jìn)展，2011，26(3):954-960.

［11］Zhdanov M S.Electromagnetic geophysics:Notes from the past and the road ahead［J］.Geophysics，2010，75(5):A49-A66.

［12］Xue G Q，Yan Y J，Li X.Control of the waveform dispersion effect and applications in a TEM imaging technique for identifying underground objects［J］.Journal of Geophysics and Engineering，2011，8(1):195-201.

［13］黃桂柏.瞬變電磁法在埋地鋼質(zhì)管道腐蝕檢測(cè)中的應(yīng)用［J］.石油化工腐蝕與防護(hù)，2001，18(4):39-42.

［14］高永才，李永年，王緒本，等.瞬變電磁法金屬管道腐蝕檢測(cè)理論初探［J］.物探化探計(jì)算技術(shù)，2005，27(1):29-33.

［15］石仁委.埋地管道壁厚的瞬變電磁檢測(cè)技術(shù)研究［J］.石油化工腐蝕與防護(hù)，2007，24(2):10-13.

［16］B?rner R U，Ernst O G，Spitzer K.Fast 3D simulation of transient electromagnetic fieldsby modelreduction in the frequencydomain using Krylov subspace projection［J］.Geophysical Journal International，2008，173(3):766-780.

［17］Swidinsky A，Edward R N.The transient electromagnetic response of a resistive sheet:An extension to three dimensions［J］.Geophysical Journal International，2010，182(2):663-674.

［18］閆述，陳明生，傅君眉.瞬變電磁場(chǎng)的直接時(shí)域數(shù)據(jù)分析［J］.地球物理學(xué)報(bào)，2002，45(2):275-284.

［19］章飛亮.考慮激電效應(yīng)的瞬變電磁二維正演研究［D］.北京:中國(guó)地質(zhì)大學(xué)，2012:198-263.

［20］徐凌華，童孝忠，柳建新，等.基于有限單元法的二維/三維大地電磁正演模擬策略［J］.物探化探計(jì)算技術(shù)，2009，31(5):421-425.

［21］張路根，胡智，葛建明，等.管道缺陷寬度對(duì)渦流檢測(cè)信號(hào)影響的仿真研究［J］.失效分析與預(yù)防，2011，6(4):232-235.