王文東 胡緒權(quán) 徐正玉 王浩文 付志紅

（1. 重慶大學(xué)電氣工程學(xué)院 重慶 400044 2. 重慶璀陸探測(cè)技術(shù)有限公司 重慶 402660 3. 長(zhǎng)江師范學(xué)院 重慶 408100）

0 引言

我國(guó)普遍采用碳鋼作為接地網(wǎng)材料，常年埋于地下容易發(fā)生腐蝕，致使其電氣性能惡化，嚴(yán)重時(shí)直接危及電網(wǎng)的穩(wěn)定運(yùn)行[1-6]。查找接地網(wǎng)的斷點(diǎn)及嚴(yán)重腐蝕段已成為電力部門一項(xiàng)重大的預(yù)防事故措施。國(guó)家能源局給出了DL/T 1532—2016《接地網(wǎng)腐蝕診斷技術(shù)導(dǎo)則》，用于指導(dǎo)敞開式變電站、換流站等電力工程接地網(wǎng)腐蝕狀態(tài)的評(píng)估[7]，導(dǎo)則給出的方法要求接地網(wǎng)的全部或者部分拓?fù)浣Y(jié)構(gòu)已知。

目前，接地網(wǎng)拓?fù)浣Y(jié)構(gòu)探測(cè)方法主要有探地雷達(dá)法[8]和電磁場(chǎng)分析法[9-12]。探地雷達(dá)探測(cè)效果容易受目標(biāo)體尺寸、深度和周圍介質(zhì)的影響，如何準(zhǔn)確地從回波數(shù)據(jù)和圖像中識(shí)別接地網(wǎng)結(jié)構(gòu)有待進(jìn)一步研究。文獻(xiàn)[9]提出了一種基于電磁場(chǎng)理論的方法，使用特種電源通過(guò)引下線向接地網(wǎng)注入和抽出數(shù)百赫茲異頻正弦電流，接收線圈采集地表磁場(chǎng)，利用模擬前置放大器、工頻陷波器和帶通濾波器從背景噪聲中分離出感應(yīng)磁場(chǎng)，根據(jù)地面磁感應(yīng)強(qiáng)度峰值位置定位導(dǎo)體。但實(shí)際測(cè)量得到的磁場(chǎng)數(shù)據(jù)受地面設(shè)備和地下管線等干擾會(huì)出現(xiàn)異常凸起和跌落點(diǎn)，磁場(chǎng)峰值定位會(huì)導(dǎo)致對(duì)導(dǎo)體的漏判和誤判，而且遠(yuǎn)離電流注入點(diǎn)的區(qū)域磁場(chǎng)峰值不明顯，難以識(shí)別，需要對(duì)接地網(wǎng)分區(qū)分塊測(cè)量。

基于電磁場(chǎng)分析理論，文獻(xiàn)[13-15]研究了基于磁場(chǎng)計(jì)算的網(wǎng)格定位方法。文獻(xiàn)[13]揭示了變電站接地網(wǎng)正逆問(wèn)題的內(nèi)在聯(lián)系，采用正則化方法計(jì)算接地網(wǎng)回路電流判斷腐蝕情況，但求解過(guò)程復(fù)雜。文獻(xiàn)[14]提出了一種微分法，利用磁場(chǎng)形函數(shù)水平分量偶數(shù)階導(dǎo)數(shù)和垂直分量奇數(shù)階導(dǎo)數(shù)的主峰特性，對(duì)注入電流后接地網(wǎng)上方磁場(chǎng)測(cè)量數(shù)據(jù)進(jìn)行濾波和數(shù)值微分處理，診斷得到精確的載流網(wǎng)格拓?fù)浣Y(jié)構(gòu)。但微分對(duì)噪聲敏感，會(huì)放大數(shù)據(jù)測(cè)量誤差，降低了抗干擾能力，因此對(duì)數(shù)據(jù)的精確性提出了很高的要求。文獻(xiàn)[15]結(jié)合磁梯度張量特征算子將磁場(chǎng)測(cè)面數(shù)據(jù)的邊緣表征為梯度矢量模極大值，引入小波邊緣檢測(cè)場(chǎng)位分離方法識(shí)別接地網(wǎng)孔。但計(jì)算梯度進(jìn)行邊緣檢測(cè)同樣依賴測(cè)量數(shù)據(jù)的精確性，而且只能處理磁場(chǎng)垂直分量。

數(shù)學(xué)形態(tài)學(xué)作為一種行之有效的非線性濾波工具，廣泛應(yīng)用于圖像處理[16]、電弧識(shí)別[17]、信號(hào)分析[18]、電力系統(tǒng)保護(hù)[19]和局放檢測(cè)[20]等領(lǐng)域，但目前尚無(wú)針對(duì)形態(tài)學(xué)接地網(wǎng)拓?fù)涑上穹椒ǖ难芯俊＝M合形態(tài)學(xué)濾波器能夠?yàn)V除白噪聲脈沖噪聲，在信號(hào)伴隨嚴(yán)重噪聲甚至發(fā)生畸變的情況下，依然可以識(shí)別基本幾何結(jié)構(gòu)，實(shí)現(xiàn)從復(fù)雜電磁環(huán)境下磁場(chǎng)數(shù)據(jù)中提取接地網(wǎng)骨架和重建拓?fù)洹?/p>

本文提出了一種基于形態(tài)濾波的接地網(wǎng)拓?fù)浣Y(jié)構(gòu)檢測(cè)方法，避免了復(fù)雜的磁場(chǎng)計(jì)算，解決了弱磁場(chǎng)峰值難以識(shí)別和強(qiáng)干擾下易誤判的問(wèn)題，實(shí)現(xiàn)了直接提取出精確、完整的接地網(wǎng)拓?fù)浣Y(jié)構(gòu)。

1 載流導(dǎo)體的空間磁場(chǎng)分布

1.1 載流單導(dǎo)體模型

接地導(dǎo)體的尺寸相對(duì)于地下空間而言很小，將單段接地導(dǎo)體等效為不計(jì)尺寸的有限長(zhǎng)載流直導(dǎo)線。在實(shí)際測(cè)量中，相對(duì)于測(cè)量系統(tǒng)所采用的坐標(biāo)系，接地網(wǎng)導(dǎo)體的走向可能是未知的。因此考慮一段與平面坐標(biāo)軸有任意夾角α的有限長(zhǎng)載流直導(dǎo)體，其模型如圖1所示，導(dǎo)體埋深為h，土壤為單層均勻土壤，磁導(dǎo)率為0μ。P是地表任意一點(diǎn)，與導(dǎo)體兩端點(diǎn)連線同導(dǎo)體夾角分別為θ1、θ2，與導(dǎo)體垂直距離r，磁感應(yīng)強(qiáng)度方向矢量e。

圖1 單根有限長(zhǎng)載流導(dǎo)體模型Fig.1 Single conductor current-carrying model with limited length

以導(dǎo)體軸向作為y′方向建立x′Oy′坐標(biāo)系，導(dǎo)體流過(guò)電流I，忽略土壤的散流效應(yīng)，由畢奧-薩法爾定律和坐標(biāo)變換方法計(jì)算得到任意坐標(biāo)系xOy下載流導(dǎo)體在空間P點(diǎn)的磁感應(yīng)強(qiáng)度為

其中

式中，(x,y)和(x′,y′)分別為xOy和x′Oy′坐標(biāo)系下的坐標(biāo)；(x0,y0)和(x1,y1)為直導(dǎo)體兩端點(diǎn)坐標(biāo)；ex、ey為方向矢量。

1.2 載流接地網(wǎng)模型

網(wǎng)格狀接地網(wǎng)由多段導(dǎo)體組成，地表面磁場(chǎng)分布可以描述為各段載流導(dǎo)體的疊加。由式（1）可得載流接地網(wǎng)在空間中產(chǎn)生的磁感應(yīng)強(qiáng)度水平分量為

式中，n為導(dǎo)體總段數(shù)；Bxi和Byi分別為第i段導(dǎo)體產(chǎn)生的空間磁感應(yīng)強(qiáng)度x、y分量。

當(dāng)導(dǎo)體與某一坐標(biāo)軸平行時(shí)，式（1）中該方向的感應(yīng)磁場(chǎng)為0，該導(dǎo)體的感應(yīng)磁場(chǎng)便不會(huì)在式（2）中體現(xiàn)，故只采用單一方向的磁感應(yīng)強(qiáng)度無(wú)法全面描述載流接地網(wǎng)所有導(dǎo)體的空間磁感應(yīng)強(qiáng)度分布。因此，本文采用x和y方向磁場(chǎng)共同描述。將兩方向磁感應(yīng)強(qiáng)度Bx、By求矢量和，使用磁通密度模值描述載流接地網(wǎng)的空間磁場(chǎng)分布，可以兼顧沿各個(gè)方向埋設(shè)的導(dǎo)體，磁通密度模MBf為

對(duì)于如圖2a所示的網(wǎng)格大小為5m的5×5的接地網(wǎng)模型，選取左下角A為電流注入點(diǎn)，右上角C為電流抽出點(diǎn)，埋深h=0.8m，計(jì)算出其空間磁場(chǎng)分布如圖2b所示。

圖2 載流接地網(wǎng)地表磁場(chǎng)分布Fig.2 Magnetic field induced by current-carrying grids

在接地網(wǎng)區(qū)域(x,y)∈R上，對(duì)于任意的常數(shù)x0、y0，當(dāng)磁通密度模梯度?MBf滿足則稱在該點(diǎn)取得局部極大值。

磁通密度模局部極大值在導(dǎo)體正上方取得，局部極大值點(diǎn)的集合構(gòu)成連通圖，即為載流導(dǎo)體構(gòu)成的接地網(wǎng)拓?fù)?。由于接地網(wǎng)電流分布不均，遠(yuǎn)離電流注入點(diǎn)的導(dǎo)體正上方的局部極大值很快衰減到全局極大值的1/10左右，甚至小于注入點(diǎn)附近的局部極小值，因此弱磁場(chǎng)區(qū)域拓?fù)潆y以被有效識(shí)別。

2 形態(tài)濾波接地網(wǎng)拓?fù)涮崛》椒?/h2>

數(shù)學(xué)形態(tài)學(xué)的基礎(chǔ)是集合論，具有完備的數(shù)學(xué)基礎(chǔ)，其本質(zhì)是非線性濾波器。形態(tài)濾波器可以簡(jiǎn)化圖像數(shù)據(jù)，保持基本形狀特性并去除不相干結(jié)構(gòu)[14-16]，因此形態(tài)濾波方法可以很好地排除異常數(shù)據(jù)干擾，獲得接地網(wǎng)拓?fù)浣Y(jié)構(gòu)。

2.1 形態(tài)濾波接地網(wǎng)拓?fù)涮崛≡?/h3>

一般而言，灰值數(shù)學(xué)形態(tài)學(xué)的處理對(duì)象是圖像信號(hào)波形的拓?fù)涮匦?，本文將磁通密度模值MBf函數(shù)視為灰度圖像，作為輸入信號(hào)。

對(duì)于輸入信號(hào)MBf(x,y)和結(jié)構(gòu)元素g，定義形態(tài)學(xué)運(yùn)算操作。

1）腐蝕

式中，（x,y）、（s,t）分別為形態(tài)學(xué)運(yùn)算前后的坐標(biāo)；Θ為數(shù)學(xué)形態(tài)學(xué)“腐蝕”運(yùn)算符；Df和Dg分別為信號(hào)定義域和結(jié)構(gòu)元素定義域。

從幾何上講，腐蝕運(yùn)算使用結(jié)構(gòu)元素在信號(hào)下部進(jìn)行滑動(dòng)，消除信號(hào)中向下的尖峰，起到濾波效果。

2）膨脹

式中，⊕為數(shù)學(xué)形態(tài)學(xué)“膨脹”運(yùn)算符。

膨脹運(yùn)算是腐蝕運(yùn)算的對(duì)偶運(yùn)算，當(dāng)結(jié)構(gòu)元素在信號(hào)定義域內(nèi)滑動(dòng)時(shí)，利用結(jié)構(gòu)元素的反射，在當(dāng)前結(jié)構(gòu)元素的定義域內(nèi)，上推結(jié)構(gòu)元素，使結(jié)構(gòu)元素超過(guò)信號(hào)的最小值，從上部對(duì)信號(hào)進(jìn)行濾波，消除向上的異常尖峰噪聲。

（1）先腐蝕再膨脹的迭代運(yùn)算為開運(yùn)算

式中，○為數(shù)學(xué)形態(tài)學(xué)“開”運(yùn)算符。

（2）先膨脹再腐蝕為閉運(yùn)算

式中，?為數(shù)學(xué)形態(tài)學(xué)“閉”運(yùn)算符。

接地網(wǎng)拓?fù)涫且粋€(gè)封閉圖形，由基爾霍夫電流定律和基爾霍夫電壓定律可知，電流的流通路徑也必然是閉合的，因此不存在孤立的支路和節(jié)點(diǎn)。地表磁場(chǎng)測(cè)量數(shù)據(jù)若出現(xiàn)孤立的極大值區(qū)域或不連續(xù)的支路，必然由噪聲或干擾導(dǎo)致。當(dāng)接地導(dǎo)體出現(xiàn)斷點(diǎn)或嚴(yán)重腐蝕，其上電流為零或近似為零，整段導(dǎo)體上方磁場(chǎng)降低為背景磁場(chǎng)，不存在段中局部磁場(chǎng)跌落的情況，如果出現(xiàn)，也必然為噪聲或干擾導(dǎo)致。

形態(tài)學(xué)開運(yùn)算作為一個(gè)基于幾何運(yùn)算的濾波器，消除比結(jié)構(gòu)元素小的孤立點(diǎn)、毛刺和小橋，而總的位置和形狀不變。開運(yùn)算只消除磁場(chǎng)數(shù)據(jù)中孤立的噪點(diǎn)而不改變?cè)型負(fù)浣Y(jié)構(gòu)。閉運(yùn)算通過(guò)填充凹角來(lái)濾波，能夠填平前景內(nèi)的小裂縫，當(dāng)載流導(dǎo)體正上方磁場(chǎng)出現(xiàn)異常跌落，閉運(yùn)算對(duì)其進(jìn)行填補(bǔ)，增強(qiáng)連接，保證拓?fù)涞耐暾浴?/p>

結(jié)構(gòu)元素的大小和形狀決定了濾波效果，選擇合理的結(jié)構(gòu)元素能夠有效濾除噪聲并保持原信號(hào)的結(jié)構(gòu)特征。

2.2 改進(jìn)的分水嶺算法

分水嶺算法是一種基于拓?fù)淅碚摰臄?shù)學(xué)形態(tài)學(xué)分割方法，將輸入信號(hào)各點(diǎn)的值視為該點(diǎn)的海拔，將輸入信號(hào)曲面視為地形圖，運(yùn)用地理學(xué)上的測(cè)線重構(gòu)思想找到邊界實(shí)現(xiàn)分水嶺變換[21]。

經(jīng)典的分水嶺算法是L. Vincent等提出的浸入模擬法[22]，局部極小值點(diǎn)作為集水盆的起始點(diǎn)，隨著水位不斷上升，集水盆的區(qū)域向外擴(kuò)散，到一定高度后兩個(gè)或多個(gè)集水盆發(fā)生匯合，匯合處即為分水嶺。

使用分水嶺算法處理接地網(wǎng)地表感應(yīng)磁場(chǎng)函數(shù)，可以獲得精確、封閉、連通的拓?fù)浣Y(jié)構(gòu)圖。但是由于噪聲影響和插值方法限制，實(shí)際數(shù)據(jù)中會(huì)出現(xiàn)大量偽極小值，進(jìn)而產(chǎn)生偽分水嶺，拓?fù)鋱D上表現(xiàn)為偽導(dǎo)體。

分水嶺算法示意圖如圖3所示。浸入式分水嶺算法要求浸水閾值h1、h2，地表磁通密度模局部極大值的變化范圍較大，不恰當(dāng)?shù)拈撝禃?huì)導(dǎo)致有用分水嶺丟失。

圖3 分水嶺算法示意圖（浸入模擬）Fig.3 Watershed algorithm(immersion simulation)

針對(duì)感應(yīng)磁場(chǎng)的數(shù)據(jù)特點(diǎn)，進(jìn)行分水嶺算法處理之前首先進(jìn)行背景補(bǔ)償和信號(hào)增強(qiáng)，實(shí)現(xiàn)數(shù)據(jù)規(guī)范化。對(duì)于磁通密度模函數(shù)MBf，使用濾波器W對(duì)對(duì)其進(jìn)行卷積。

濾波器wuv=1，U、V為卷積核大小，卷積核需略大于最大網(wǎng)格間距，以使其能夠同時(shí)覆蓋局部極大值和局部極小值，保證后續(xù)數(shù)據(jù)增強(qiáng)效果。對(duì)函數(shù)MBf等距抽樣選點(diǎn)得到一個(gè)二維矩陣，卷積獲得該點(diǎn)附近領(lǐng)域U×V點(diǎn)個(gè)值的累加和，實(shí)際上為避免峰值影響，應(yīng)選取其中值最小的20%的點(diǎn)求算數(shù)平均，定義求取到的點(diǎn)(x,y)的背景值為MBfh(x,y)。

背景值越小，說(shuō)明該點(diǎn)越遠(yuǎn)離導(dǎo)體正上方，越遠(yuǎn)離電流注入和抽出點(diǎn)。使用距離反比法進(jìn)行背景補(bǔ)償，點(diǎn)(x,y)處的距離d和補(bǔ)償函數(shù)分別定義為

MBf′(x,y)為補(bǔ)償后的磁場(chǎng)，λ、p、k1、k2取值應(yīng)使得分段函數(shù)連續(xù)，一般而言，k1可取3~8之間的的一個(gè)常數(shù)，k2可取0.8~1.5。全局閾值濾波器d0用于濾除數(shù)據(jù)中的微小噪點(diǎn)，分界點(diǎn)d1~d3為先驗(yàn)值，根據(jù)實(shí)驗(yàn)結(jié)果確定取值。

圖4 a給出了圖2b三維圖中y=2.5m處的磁通密度模，對(duì)其進(jìn)行背景補(bǔ)償和數(shù)據(jù)增強(qiáng)處理，補(bǔ)償倍數(shù)1~6.2倍，d1和d3分別取距離最大值dmax的80%和10%，結(jié)果如圖4b所示。

圖4 y=2.5m磁通密度模規(guī)范化結(jié)果Fig.4 Normalization results for magnetic flux at y=2.5m

2.3 接地網(wǎng)拓?fù)浣Y(jié)構(gòu)提取流程

接地網(wǎng)拓?fù)浣Y(jié)構(gòu)提取流程分為數(shù)據(jù)預(yù)處理和組合形態(tài)濾波獲得拓?fù)鋬蓚€(gè)部分。圖5給出了基于形態(tài)濾波的接地網(wǎng)拓?fù)涮崛×鞒獭?/p>

圖5 接地網(wǎng)拓?fù)涮崛×鞒蘁ig.5 Topology extraction process

采集到地表磁場(chǎng)數(shù)據(jù)后，首先對(duì)數(shù)據(jù)預(yù)處理。變電站測(cè)量到的磁場(chǎng)數(shù)據(jù)幅值較小，信號(hào)容易被以工頻噪聲為主的背景噪聲淹沒。使用模擬濾波器濾除信號(hào)中的干擾較強(qiáng)的工頻以及低次諧波，再通過(guò)高速數(shù)據(jù)采集卡對(duì)信號(hào)采樣后經(jīng)過(guò)數(shù)字帶通濾波器獲得與激勵(lì)源同頻磁場(chǎng)信號(hào)，每個(gè)測(cè)點(diǎn)保存一個(gè)磁場(chǎng)數(shù)據(jù)幅值用于計(jì)算和成像。對(duì)原始磁場(chǎng)數(shù)據(jù)二次曲面插值，近似得到整個(gè)測(cè)面的數(shù)據(jù)。再對(duì)數(shù)據(jù)進(jìn)行高斯平滑，濾除測(cè)面上的高斯噪聲。

然后構(gòu)建復(fù)合形態(tài)濾波器，其構(gòu)成包括形態(tài)開(Og1)、形態(tài)閉(Cg1)、腐蝕重建(ERg1)、膨脹重建(DRg1)、數(shù)據(jù)規(guī)范化(DN)、形態(tài)開-閉(OCg2)、形態(tài)學(xué)分水嶺邊界提取(WA)。

預(yù)處理后的數(shù)據(jù)歸一化到灰度值0~255，經(jīng)過(guò)形態(tài)開濾波器消除孤立凸起點(diǎn)，形態(tài)閉濾波器填充結(jié)構(gòu)裂縫，腐蝕重建和膨脹重建增強(qiáng)數(shù)據(jù)，數(shù)據(jù)規(guī)范化增強(qiáng)弱磁場(chǎng)區(qū)域峰值，小結(jié)構(gòu)元素g2形態(tài)開-閉濾波消除數(shù)據(jù)規(guī)范化產(chǎn)生的噪點(diǎn)，通過(guò)分水嶺算法提取出接地網(wǎng)拓?fù)洹?/p>

3 仿真算例

建立圖2a中接地網(wǎng)的多物理場(chǎng)仿真模型，仿真模型如圖6a所示，導(dǎo)體截面積6cm×5mm，網(wǎng)格大小5m×5m，埋深0.8m，10A電流從A點(diǎn)注入、C點(diǎn)抽出。單層土壤，厚度50m，土壤電阻率50Ω·m，土壤磁導(dǎo)率近似取真空磁導(dǎo)率0μ。圖6b給出了仿真結(jié)果，將磁通密度模轉(zhuǎn)換得到大小為351×351像素灰度圖。對(duì)灰度圖數(shù)據(jù)分別使用文中提出的形態(tài)學(xué)方法、canny邊緣檢測(cè)算子、直接使用分水嶺算法處理，結(jié)果如圖7所示。網(wǎng)格間距為50像素，結(jié)構(gòu)元素尺寸不能大于該值，否則將造成局部極大值之間的粘連；結(jié)構(gòu)元素尺寸同樣不宜小于網(wǎng)格間距的1/2，否則不能實(shí)現(xiàn)網(wǎng)格中心異常凸起的去除。綜合考慮，形態(tài)學(xué)處理時(shí)選擇結(jié)構(gòu)元素大小為31。

圖6 仿真結(jié)果Fig.6 Simulation results

圖7 拓?fù)浣Y(jié)構(gòu)提取結(jié)果Fig.7 Topology extraction results

可以看出，直接使用canny邊緣檢測(cè)算法處理，存在磁場(chǎng)較弱的區(qū)域的接地網(wǎng)結(jié)構(gòu)不能有效識(shí)別的情況，而直接使用分水嶺算法，同樣不能有效識(shí)別出弱磁場(chǎng)區(qū)域的接地導(dǎo)體，而且存在一定程度的過(guò)分割情況。定義識(shí)別到的導(dǎo)體段數(shù)與導(dǎo)體總段數(shù)之比為拓?fù)渫暾?，定義識(shí)別到“偽導(dǎo)體”段數(shù)導(dǎo)體總段數(shù)之比為誤判率，對(duì)比三種方法拓?fù)涮崛〗Y(jié)果。

仿真結(jié)果對(duì)比見表1，傳統(tǒng)的邊緣檢測(cè)算法canny算子提取的拓?fù)渫暾蕛H為不到90%，弱磁場(chǎng)區(qū)域的拓?fù)浣Y(jié)構(gòu)丟失，而本文的形態(tài)學(xué)方法達(dá)到了100%。相比于直接使用分水嶺算法，形態(tài)學(xué)處理后提取到的拓?fù)涓暾?，而且基本消除了過(guò)分割情況。仿真結(jié)果表明，數(shù)學(xué)形態(tài)學(xué)方法可以有效處理弱磁場(chǎng)區(qū)域數(shù)據(jù)，提取到更完整的接地網(wǎng)拓?fù)浣Y(jié)構(gòu)。

表1 仿真結(jié)果對(duì)比Tab.1 Comparison of simulation results

4 實(shí)驗(yàn)分析

在實(shí)驗(yàn)場(chǎng)開展模擬實(shí)驗(yàn)。整個(gè)實(shí)驗(yàn)接地網(wǎng)長(zhǎng)13.6m，寬10.2m，由4×3個(gè)網(wǎng)格構(gòu)成，網(wǎng)格間距3.2~3.4m，接地導(dǎo)體使用60mm×8mm鍍鋅扁鋼，埋設(shè)深度0.8m。

接地網(wǎng)模擬實(shí)驗(yàn)如圖8所示。從圖8中A、B兩點(diǎn)注入和抽出頻率9kHz、幅值9.9A的雙極性電流。采用直徑15cm，匝數(shù)600的接收線圈采集地表感應(yīng)磁場(chǎng)，選用數(shù)字多用表34401A測(cè)量和記錄接收線圈感應(yīng)電壓有效值。接收線圈感應(yīng)電壓有效值可以表征地表磁感應(yīng)強(qiáng)度，磁場(chǎng)測(cè)量區(qū)域應(yīng)遠(yuǎn)離電流注入點(diǎn)[4]，選擇網(wǎng)格EFGH進(jìn)行測(cè)量，測(cè)點(diǎn)間隔10cm。測(cè)量區(qū)域內(nèi)，各段導(dǎo)體的正上方均布置測(cè)點(diǎn)。

圖8 接地網(wǎng)模擬實(shí)驗(yàn)Fig.8 Simulation experiment of grounding grids

使用感應(yīng)電壓有效值成像，以1cm為間隔對(duì)原始數(shù)據(jù)進(jìn)行插值，測(cè)量結(jié)果如圖9a所示。形態(tài)學(xué)方法處理結(jié)果如圖9b所示，處理后提取到單像素寬的拓?fù)浣Y(jié)構(gòu)圖。F點(diǎn)靠近電流注入點(diǎn)，直線距離僅為3.6m，磁場(chǎng)受到強(qiáng)干擾導(dǎo)致定位出現(xiàn)較大偏移，其余導(dǎo)體定位較為準(zhǔn)確。實(shí)際值與測(cè)量值對(duì)比見表2。

表2 實(shí)驗(yàn)結(jié)果對(duì)比Tab.2 Comparison of simulation results

圖9 形態(tài)濾波處理結(jié)果Fig.9 Processing results based on morphological filtration

將提取得到的拓?fù)浣Y(jié)構(gòu)與實(shí)際施工圖紙比較，以網(wǎng)格大小3.4m×3.4m作為基準(zhǔn)值計(jì)算定位誤差。F點(diǎn)距離電流注入點(diǎn)較近，在強(qiáng)干擾的情況下，形態(tài)學(xué)方法仍然能夠定位端點(diǎn)，且誤差不超過(guò)網(wǎng)格大小的10%。遠(yuǎn)離電流注入點(diǎn)的區(qū)域，各導(dǎo)體和端點(diǎn)的定位誤差不超過(guò)3%。實(shí)驗(yàn)結(jié)果驗(yàn)證了使用形態(tài)濾波方法提取接地網(wǎng)拓?fù)浣Y(jié)構(gòu)的可行性。

5 結(jié)論

本文提出了一種基于形態(tài)濾波的接地網(wǎng)拓?fù)浣Y(jié)構(gòu)提取方法，通過(guò)構(gòu)建復(fù)合形態(tài)濾波器，直接從磁場(chǎng)數(shù)據(jù)中識(shí)別和定位接地網(wǎng)。在數(shù)據(jù)精度有限的情況下得到較為精確的結(jié)果。仿真和模擬地網(wǎng)實(shí)驗(yàn)結(jié)果驗(yàn)證了該方法的可行性，為應(yīng)用于實(shí)際變電站接地網(wǎng)拓?fù)錂z測(cè)奠定了基礎(chǔ)。