一種接觸網(wǎng)絕緣子污閃預測方法的研究

王思華， 景 弘(蘭州交通大學 自動化與電氣工程學院， 甘肅 蘭州 730070)

隨著大氣污染的日益嚴重，接觸網(wǎng)絕緣子污穢程度亦隨之加重，其絕緣性能承受著雨雪霧等惡劣天氣及周邊廠礦企業(yè)工業(yè)粉塵、煤煙塵的實時考驗。絕緣子臟污潮濕導致的放電閃絡現(xiàn)象，將造成接觸網(wǎng)設備大面積停電，嚴重威脅鐵路運輸安全生產(chǎn)。污閃事故雖然次于占電力系統(tǒng)事故第一位的雷害事故，但是其造成的損失卻是雷害的10倍[1]。故研究相應的污閃預測方法對保證接觸網(wǎng)供電可靠性具有非常重要的意義。

目前污閃預測主要集中在以監(jiān)測絕緣子泄漏電流、光波導為核心的在線監(jiān)測方面[2-5]。存在監(jiān)測設備昂貴，實際運行中受周圍環(huán)境因素影響較大等缺陷。從20世紀90年代開始，國內(nèi)外許多學者開始研究運用數(shù)學方法對污閃進行預測，并取得了大量的研究成果。當前對絕緣子污閃預測的模型可歸納為基于模糊原理、神經(jīng)網(wǎng)絡和最小二乘支持向量機等3類模型。文獻[6-8]根據(jù)流經(jīng)絕緣子串的泄漏電流及電暈電流，利用模糊推理的分析處理方法評定絕緣子的劣化狀態(tài)；文獻[9-12]利用神經(jīng)網(wǎng)絡建立復雜環(huán)境信息與絕緣子閃絡電壓的神經(jīng)網(wǎng)絡模型，預測絕緣子的閃絡電壓；文獻[13-17]利用最小二乘支持向量機具有較傳統(tǒng)神經(jīng)網(wǎng)絡更好的處理小樣本的能力，對絕緣子的等值鹽密進行實時預測。以上研究多從定性的角度提出污閃預測方法，對影響污閃發(fā)生的具體因素及污閃預測方法進行了深入研究，得出污閃的發(fā)生受絕緣子污穢程度和氣象因素的影響。但受實驗條件、實際環(huán)境差異和污閃復雜性的限制，目前的研究成果尚無法絕對準確預測污閃發(fā)生，因此還需要繼續(xù)研究絕緣子污閃預測方法。

本文將絕緣子污穢程度與絕緣子污閃概率結(jié)合起來，首先探討人工魚群算法優(yōu)化BP神經(jīng)網(wǎng)絡的可行性，建立基于人工魚群BP神經(jīng)網(wǎng)絡算法的絕緣子污穢等級預測模型。選取月降水量、月均風速、月均相對濕度、月降塵量等4個最具代表性的氣象因子和絕緣子運行時間(月)作為模型輸入特征量。輸入上個月5個特征量的數(shù)據(jù)，預測當月絕緣子的污穢等級。然后建立基于絕緣子污穢等級的污閃概率預測模型，預測當月接觸網(wǎng)絕緣子污閃概率，提前進行預警并建立接觸網(wǎng)線路污閃短期預報響應機制。對評估接觸網(wǎng)線路可靠性和防污閃工作起到輔助參考作用，在一定程度上避免人力、物力、財力的浪費。

1 影響絕緣子污穢程度的氣象因素分析

由于接觸網(wǎng)絕緣設備多暴露在露天環(huán)境下且安裝高度相對較低，故設備表面積污易受外部氣象環(huán)境的影響。理論上考慮的影響因素越多越好，但根據(jù)可行性、實用性和直觀性原則，本文對影響絕緣子污穢程度的氣象因素分析如下。

(1) 相對濕度。濕度對絕緣子污穢程度影響很大，當空氣中相對濕度較低時絕緣子表面污穢物易形成干灰，很難凝結(jié)成污穢物；由于污穢物顆粒在絕緣子表面的碰撞和吸附過程可分為入射、碰撞、出射3個階段。當相對濕度逐漸增加時，絕緣子表面的水膜也逐漸變厚，污穢物顆粒出射時動能損失隨之變大，使污穢物顆粒出射速度變小，污穢物容易沉積在絕緣子表面。

(2) 降水。降水對絕緣子污穢程度的影響最為明顯，當進入夏秋季節(jié)(5～9月)時，大部分地區(qū)降雨較多，24 h內(nèi)降水量大于10 mm時，絕緣子表面ESDD受到降水的機械沖刷和溶解沖刷作用，對污穢絕緣子具有良好的沖洗效果，絕緣子污穢程度大幅度降低。但是降水沖刷污穢絕緣子形成的污流可能使絕緣子傘裙間氣隙橋接，引發(fā)閃絡事故。

(3) 風速。風速對絕緣子污穢程度的影響主要體現(xiàn)在風速大小上，當風速較小時，絕緣子表面的ESDD相對較低，隨著風速的不斷增大，污穢物顆粒的入射速度隨之增大，但污穢物顆粒的碰撞恢復系數(shù)將逐漸降低。致使污穢物顆粒出射速度增加，污穢顆粒難以吸附在絕緣子表面。

(4) 降塵量。降塵量是影響絕緣子污穢程度的重要因素之一。由于接觸網(wǎng)絕緣子一般架設在鐵路線路上方5～7 m，而裝載貨物的列車最高可達4.8 m。在列車運行時會揚起列車貨物中的煤塵、粉塵、輪軌摩擦產(chǎn)生的金屬塵埃、電弓的碳粉，這些污穢在風力和自身重力作用下沉降到絕緣子表面造成污穢程度加重。特別是接觸網(wǎng)周圍有廠礦企業(yè)時，工業(yè)污穢源造成絕緣子積污更為嚴重。

圖1為2006—2008年某供電段接觸網(wǎng)XP-70型懸式瓷絕緣子表面鹽密統(tǒng)計與該地區(qū)氣象因素變化情況的對比，時間節(jié)點按春夏秋冬4個季節(jié)排列。圖中列出了降水量、風速、相對濕度、降塵量和絕緣子表面鹽密的變化情況。從圖中可以看出，降水對絕緣子表面鹽密的影響最為明顯，3年的降水量高峰期基本對應當年絕緣子表面鹽密低谷期；絕緣子的鹽密高峰期與降塵量的增長期具有高度的一致性，均出現(xiàn)在每年的冬季和次年的初春。此外從圖中也可以看出，絕緣子的表面鹽密與相對濕度和風速呈現(xiàn)出一定的相關性。

2 絕緣子污穢等級預測模型的建立

2.1 建模思想

BP神經(jīng)網(wǎng)絡是一種有3層或3層以上的梯度下降型神經(jīng)網(wǎng)絡，上下層之間各神經(jīng)元實現(xiàn)了全連接。具有強大的逼近能力，可以實現(xiàn)輸入和輸出間的任意非線性映射，解決任何非線性函數(shù)的逼近問題。廣泛應用于模式識別、自動化系統(tǒng)、機器學習和圖像處理等領域。主要結(jié)構(gòu)包括輸入層、隱含層及輸出層，每層的神經(jīng)元稱為節(jié)點或單元。由于一個3層BP神經(jīng)網(wǎng)絡可以完成任意M維輸入到N維輸出的映射，故常用的是3層型BP神經(jīng)網(wǎng)絡，網(wǎng)絡模型見圖2。xm(i)為輸入層節(jié)點;wij為輸入層神經(jīng)元i與隱含層神經(jīng)元j之間的連接權(quán)值;vjk為隱含層神經(jīng)元j與輸出層神經(jīng)元k之間的連接權(quán)值;ym(k)為輸出層節(jié)點。各層之間的連接權(quán)值在樣本數(shù)據(jù)的學習過程中不斷修正，隱含層的個數(shù)根據(jù)實際需要可以是一層或多層。

BP神經(jīng)網(wǎng)絡盡管具有廣泛的實用性，但其算法存在著一些缺陷限制了它的使用范圍。比如BP神經(jīng)網(wǎng)絡易陷入局部最優(yōu)是其致命弱點，導致在求解最優(yōu)化非線性函數(shù)問題時函數(shù)的全局最優(yōu)解無法找到。同時算法還存在收斂速度慢，網(wǎng)絡學習和記憶性差等缺點。

2002年，李曉磊提出人工魚群算法AFSA[18]，該算法通過模仿魚群的覓食，聚群及追尾行為來尋求全局最優(yōu)解。該算法收斂速度快、全局尋優(yōu)能力強，具有良好的克服局部極值的能力，而且算法不受問題峰數(shù)增加和問題維數(shù)的影響，可以很好地彌補BP神經(jīng)網(wǎng)絡的缺陷。故可把BP神經(jīng)網(wǎng)絡和人工魚群算法有機的結(jié)合起來，用人工魚個體代表神經(jīng)網(wǎng)絡狀態(tài)，用食物濃度代表實際輸出值與期望值均方誤差的倒數(shù)。人工魚群按照魚群算法的步驟反復優(yōu)化BP神經(jīng)網(wǎng)絡的連接權(quán)值和閾值，以誤差函數(shù)E為最小尋優(yōu)對象，通過尋找最大食物濃度使均方誤差最小，從而達到獲取最優(yōu)神經(jīng)網(wǎng)絡結(jié)構(gòu)的目地。根據(jù)對影響絕緣子污穢程度的因素分析，選定BP神經(jīng)網(wǎng)絡的輸入、輸出節(jié)點，由此建立基于人工魚群-BP神經(jīng)網(wǎng)絡的絕緣子污穢等級預測模型。

2.2 算法結(jié)構(gòu)

本文選取3層型BP神經(jīng)網(wǎng)絡，設定輸入層神經(jīng)元個數(shù)為n，輸出層神經(jīng)元個數(shù)為m，隱含層神經(jīng)元個數(shù)為n1，期望輸出與實際輸出之間的誤差平方和作為誤差函數(shù)E。其中隱含層神經(jīng)元個數(shù)n1和誤差函數(shù)E分別為

( 1 )

( 2 )

利用人工魚群算法優(yōu)化BP神經(jīng)網(wǎng)絡的關鍵在于構(gòu)造人工魚個體模型。定義人工魚群的規(guī)模為S，魚群采用隨機分布，每條魚都代表一個BP神經(jīng)網(wǎng)絡，任意兩條人工魚p、q之差與之和仍代表一個BP神經(jīng)網(wǎng)絡[18-19]。兩條人工魚之間的距離為

( 3 )

式中：wj0為隱含層神經(jīng)元j的閾值；vk0為輸出層神經(jīng)元k的閾值。

2.3 算法流程

利用人工魚群算法優(yōu)化BP神經(jīng)網(wǎng)絡的流程為：

Step1設置初始化人工魚群算法的參數(shù)。即設定初始迭代次數(shù)Gen，人工魚移動的最大步長T，人工魚的視野Visual，擁擠度因子δ，誤差函數(shù)E等。

Step2計算初始魚群中各條人工魚當前位置的食物濃度FC，比較大小后將最大值送入到公告板中。其中FC=1/E。

Step3各條人工魚分別模擬追尾行為，聚群行為。選擇行動后FC值較大的行為來實際執(zhí)行。

Step4各條人工魚每迭代一次，比較自身FC值與公告板中FC值。若有人工魚的FC值優(yōu)于公告板中的FC值，則取而代之。最后將人工魚群算法找到的最優(yōu)解即公告板中人工魚狀態(tài)轉(zhuǎn)換成BP神經(jīng)網(wǎng)絡中對應參數(shù)即為訓練結(jié)果。

Step5算法的中止條件：判斷Gen是否達到最大迭代次數(shù)maxGen或已滿足解的誤差精度(通常設定誤差精度小于0.001)，只要二者中有一個條件符合就輸出公告板中的FC值，否則Gen+1后返回到Step3。具體流程圖見圖3。

2.4 行為描述

(1) 覓食行為

設人工魚p當前狀態(tài)為Fp，在其視野范圍內(nèi)隨機選擇一個狀態(tài)Fq。若人工魚Fq的食物濃度FCq大于Fp的食物濃度FCp，則意味著人工魚向著Fq的方向前進一步，否則人工魚在其視野里隨機選擇一個狀態(tài)向著該方向移動。人工魚優(yōu)化wij(p)的過程為

( 4 )

式中：T為人工魚移動的最大步長值；rand(T)為一個介于0和T之間的隨機數(shù)。同理可得出人工魚優(yōu)化wjk、vjk、wj0、vk0的過程。

(2) 聚群行為

人工魚p探索當前可視域內(nèi)的其他人工魚伙伴數(shù)目nf及中心位置F(采用歐式距離)，并計算該中心位置食物濃度值FCc。若FCc/nf>δFC，則表明中心位置的食物濃度高、魚群不擁擠，此時伙伴中心位置將會是人工魚的前進方向，人工魚朝著此方向游近一步，否則人工魚執(zhí)行覓食行為。人工魚wij(p)的變化過程為

wij(p+1)=

( 5 )

式中：dp,c為Fp距離中心位置的長度。

(3) 追尾行為

定義人工魚可視域內(nèi)所有伙伴中食物濃度最大的人工魚狀態(tài)為Fmax。若其滿足式( 6 )，則表明人工魚Fmax的食物濃度高且周圍不太擁擠。此時Fmax處將會是人工魚的前進方向，人工魚朝著此方向游近一步，否則人工魚執(zhí)行覓食行為。人工魚參數(shù)wij(p)的優(yōu)化過程為

FCmax>δFCp

( 6 )

wij(p+1)=wij(p)+

( 7 )

式中：wij(max)為最大人工魚狀態(tài)Fmax的參數(shù)；dp,max為人工魚當前狀態(tài)Fp與最大人工魚狀態(tài)Fmax之間的距離。

3 絕緣子污穢等級預測模型的應用

3.1 特征量的選取

文獻[20-24]研究表明絕緣子污穢程度不存在積累效應，污穢程度是一個隨機概率值，一年內(nèi)做周期變化。其不僅受污穢環(huán)境、氣象因素的影響，還與絕緣子運行時間有關。綜上，本文選取上個月的月降水量、月均風速、月均相對濕度、月降塵量和絕緣子運行時間(月)等5個指標作為網(wǎng)絡輸入特征量，網(wǎng)絡輸出當月絕緣子污穢等級。

3.2 使用的數(shù)據(jù)

本文收集了某供電段2008—2014年5個輸入特征量的全部數(shù)據(jù)。氣象數(shù)據(jù)來源于該供電段所在城市的氣象局，絕緣子數(shù)據(jù)來自該供電段接觸網(wǎng)絕緣子運維數(shù)據(jù)庫。將2008—2012年60個月的數(shù)據(jù)作為訓練樣本，2013—2014年24個月的數(shù)據(jù)作為檢測樣本。為了提高BP神經(jīng)網(wǎng)絡的計算精度和泛化能力，對數(shù)據(jù)做歸一化處理，處理方法[25-26]為

( 8 )

式中：xmax、xmin分別為指標x的最大和最小值；x*為指標x歸一化后的值，x*∈ [0，1]。

3.3 絕緣子污穢等級的預測

本文BP神經(jīng)網(wǎng)絡隱含層神經(jīng)元傳遞函數(shù)采用tansig函數(shù)，輸出層采用purelin函數(shù)。輸入層神經(jīng)元個數(shù)為5，輸出層神經(jīng)元個數(shù)為1，隱含層神經(jīng)元個數(shù)由式( 1 )初步確定為7。使用循環(huán)函數(shù)進行Matlab編程，當選取隱含層神經(jīng)元個數(shù)為5～9，誤差函數(shù)E為0.001時發(fā)現(xiàn)都能順利收斂，所用步數(shù)見表1。其中MSE為均方誤差，EPOCHS為訓練步數(shù)。從表1中可以看出，當隱含層神經(jīng)元個數(shù)為7時，此時收斂的最快，故確定隱含層神經(jīng)元個數(shù)為7。

表1 不同隱含層神經(jīng)元的均方誤差和訓練步數(shù)

使用人工魚群算法訓練BP神經(jīng)網(wǎng)絡，設定魚群規(guī)模S=100，人工魚的可視域范圍Visual為[0，2.5],擁擠度因子δ范圍為[0.95，2.98]，人工魚最大移動步長T=0.3，最大迭代次數(shù)kmax=200，網(wǎng)絡的誤差函數(shù)E=0.001。網(wǎng)絡輸入特征量為5個指標數(shù)據(jù)，期望輸出為絕緣子污穢等級，當訓練誤差趨向0時停止訓練。分別使用BP神經(jīng)網(wǎng)絡和人工魚群算法優(yōu)化BP神經(jīng)網(wǎng)絡進行訓練，結(jié)果見圖4、圖5。

從圖4、圖5可得：BP神經(jīng)網(wǎng)絡經(jīng)過166次訓練后才達到目標，而人工魚群算法優(yōu)化BP神經(jīng)網(wǎng)絡只需88次訓練即可。后者收斂速度快，全局搜索能力強，能更快的找到全局最優(yōu)點。從圖中還可以看出2種算法的收斂過程也有很大區(qū)別，BP神經(jīng)網(wǎng)絡在收斂過程中長時間停留在10-1左右，且停留了接近130次。而人工魚群算法優(yōu)化BP神經(jīng)網(wǎng)絡是以較大斜率收斂，很快收斂于目標誤差。雖然在局部最優(yōu)處也有停留，但是僅僅短暫停留了10次左右，相比于前者來說停留次數(shù)減少很多。綜上可以看出，人工魚群算法與BP神經(jīng)網(wǎng)絡的結(jié)合不僅能夠充分發(fā)揮魚群算法快速獲得全局最優(yōu)解的特性，同時克服了BP神經(jīng)網(wǎng)絡收斂速度慢、易陷入局部最優(yōu)解的缺陷，在相同的訓練次數(shù)下，得到的結(jié)果更精確。以XP-70型絕緣子為例，使用本文所提方法對絕緣子污穢等級進行預測，驗證本文所提方法的正確性和有效性。部分檢測樣本數(shù)據(jù)和絕緣子污穢等級預測結(jié)果見表2。

表2 部分檢測樣本數(shù)據(jù)和絕緣子污穢等級預測結(jié)果

從表2的預測結(jié)果可知，本文所提方法可以準確的預測出絕緣子污穢等級，證明本文所提方法的優(yōu)越性和可行性。其中絕緣子污穢等級參照GB/T 16434-1996《高壓架空線路和發(fā)電廠、變電所環(huán)境污區(qū)分級及外絕緣選擇標準》，分為0～Ⅳ級共5級，各級所對應的鹽密值分別為：0～0.03、0.03～0.06、0.06～0.10、0.10～0.25、0.25～0.35 mg/cm2。

4 接觸網(wǎng)絕緣子污閃概率的計算

4.1 污閃概率模型

絕緣子在工作電壓Uo下發(fā)生污閃主要由2個因素決定:大氣污染造成的絕緣子表面積污和能使表面污穢物質(zhì)充分受潮的氣象條件[27]。根據(jù)相關研究表明，絕緣子污穢程度是一個隨機概率值。由于接觸網(wǎng)線路上多種絕緣子并存，設絕緣子型號共有k類，絕緣子表面污穢程度(只考慮鹽密值)的概率密度函數(shù)為g(ρ)，第i類(i=1,2, …,k)絕緣子在某污穢程度下發(fā)生閃絡的概率為Pi(ρ)。則接觸網(wǎng)第i類單串絕緣子在等值附鹽密度處于大于a小于b區(qū)間時污閃概率為

( 9 )

式中：

(10)

(11)

(12)

(13)

(14)

(15)

式中：N為某地區(qū)未來一段時間惡劣天氣(霧、小雨、降雪等)的天數(shù)，本文中惡劣天氣的天數(shù)指24 h內(nèi)降水量小于10 mm的天數(shù)與霧日天數(shù)之和。在氣象學中定義小雨為24 h內(nèi)降水量小于10 mm的降水過程，等價到降雪天氣來說，24 h內(nèi)降水量(折合為融化后的雨水量)小于10 mm的天氣現(xiàn)象為小雪、中雪、大雪。定義霧為因大氣中懸浮的水汽凝結(jié)，能見度低于1 km的天氣現(xiàn)象。一天中只要出現(xiàn)霧，不論持續(xù)時間多長都算一個霧日，且氣象上是以北京時間20 h為日界，如果霧天氣從20 h前持續(xù)到20 h后則算兩個霧日。其中霧為水平降水，降雨、降雪為垂直降水，用降水量來區(qū)分降水的強度，單純的霧不作為降水量處理[28-30]；a、b為絕緣子污穢等級對應的鹽密值；ρ為等值附鹽密度；g(ρ)為Sforzini和Naito[31-32]經(jīng)過長期現(xiàn)場觀測發(fā)現(xiàn)絕緣子鹽密值服從對數(shù)正態(tài)分布規(guī)律，從而得出的概率密度函數(shù)[33]；μ為ρ的對數(shù)平均值；σ為ρ的對數(shù)標準差；z為絕緣子的標準方差與其50%閃絡電壓U50的比值；α為鹽密影響污閃電壓的特征指數(shù)；ρ0為切斷鹽密[34]，表示恒定電壓下運行的絕緣子當其表面等值鹽密小于ρ0時將不會發(fā)生閃絡；n為切斷參數(shù)，IEC推薦n值取4，但有關實驗數(shù)據(jù)表明n取4時對于污穢絕緣子來說偏高。Ivanov通過采集2 800個人工污穢實驗的結(jié)果進行比較得出對于針式絕緣子n取2.5最為適宜[35]；Houlgate和Risk對棒式和懸式絕緣子進行自然積污實驗得出n=2.1[36-37]；ρ50為50%閃絡概率對應的等值附鹽密度；A為與絕緣子形狀和污穢程度有關的系數(shù)；L為絕緣子爬電距離；β和h為威布爾參數(shù)。

由于接觸網(wǎng)不同型號的絕緣子在材質(zhì)和結(jié)構(gòu)上存在差異，絕緣子污穢特性不僅受外界環(huán)境的影響，還與絕緣子自身因素密切相關。在相同環(huán)境下，絕緣子表面的積污情況主要取決于絕緣子的材質(zhì)和結(jié)構(gòu)[38]，故接觸網(wǎng)不同型號的絕緣子的A值也不同。對于A的取值，Topalis通過對多種型式絕緣子進行人工污穢實驗得出：A是由絕緣子結(jié)構(gòu)決定的函數(shù)，且A與絕緣子的爬電距離之間存在著簡單的函數(shù)關系[39]，即

A=0.130L+1.947

(16)

表3列出了電氣化鐵路接觸網(wǎng)部分常用絕緣子的爬電距離。

表3 部分常用絕緣子的爬電距離

4.2 污閃概率計算

因為鐵路運維部門對接觸網(wǎng)絕緣部件的清掃提前制訂了專項清掃計劃，且平時天窗點緊張，故定義N取當月未來30 d惡劣天氣的天數(shù)較為合理。根據(jù)已建立的絕緣子污穢等級預測模型，在當月1日輸入上個月5個特征量的相關數(shù)據(jù)來預測當月絕緣子的污穢等級。在GB/T 16434—1996標準中找出當月污穢等級對應的鹽密值上下限，由此得到式( 9 )中參數(shù)a、b。查詢當?shù)貧庀缶株P于當月未來30天惡劣天氣預報，統(tǒng)計惡劣天氣天數(shù)，確定參數(shù)N。取μ=-2.617，σ=0.806，α=0.25，β=0.2，z=0.088[35]，n=2.1。由式(16)計算出接觸網(wǎng)相關型號絕緣子對應的系數(shù)A。將以上參數(shù)代入式( 9 )中計算在當月未來30 d接觸網(wǎng)相關型號絕緣子的污閃概率。

由于接觸網(wǎng)發(fā)生污閃時往往是多個絕緣子發(fā)生污閃，故預測接觸網(wǎng)線路發(fā)生污閃的概率更有實用意義。設一段距離(如10 km)的接觸網(wǎng)線路上，絕緣子型號共有y類，第i類(i=1,2, …,y)絕緣子共有mi串，相同型號的絕緣子串閃絡概率均相同。接觸網(wǎng)線路在某污穢程度下發(fā)生閃絡的概率為

Pline(ρ)=1-[1-P1(ρ)]m1[1-P2(ρ)]m2…

[1-Py(ρ)]my

(17)

將式(11)代入到式(17)中可得

Pline(ρ)=1-

(18)

則一段距離的接觸網(wǎng)線路在當月未來30 d發(fā)生污閃的概率為

(19)

式中：s、t為接觸網(wǎng)線路所在區(qū)域的污穢等級對應的鹽密值，線路污穢等級通過查詢污區(qū)圖即可得到。

4.3 污閃預警

污閃是一種對電氣化鐵路供電可靠性危害很大的頻繁性事故，為提高供電可靠性制定合理的防污閃措施是預防污閃發(fā)生的一個重要手段。為此可將預測的污閃概率按照發(fā)生的可能性分為5級，每級制定不同的防污閃措施。根據(jù)每月預測的污閃概率大小采取不同的響應措施，由此建立接觸網(wǎng)線路污閃短期預報響應機制，見表4。

表4 污閃短期預報響應機制

5 算例分析

以2015年2月多次遭遇惡劣天氣的某供電段一段長10 km的接觸網(wǎng)線路為例，計算絕緣子和接觸網(wǎng)線路的污閃概率。根據(jù)該供電段對此段線路接觸網(wǎng)供電系統(tǒng)事故調(diào)查材料，2月份頻發(fā)的惡劣天氣造成接觸網(wǎng)正饋線型號為PXN7-70EQ棒形懸式絕緣子和型號為QBSJ-25/16腕臂絕緣子發(fā)生污閃，累計中斷供電近3 h，為鐵路交通一般D21事故。此段線路接觸網(wǎng)運行有4種型號絕緣子，絕緣部件技術(shù)資料見表5。

表5 絕緣部件技術(shù)資料

調(diào)取4種類型絕緣子在2015年1月份5個輸入特征量的全部數(shù)據(jù)，見表6。查詢當?shù)貧庀缶?月份的天氣預報，確定參數(shù)N=5。根據(jù)以上數(shù)據(jù)來計算在2月份未來30 d 4種型號絕緣子的污閃概率和接觸網(wǎng)線路的污閃概率。

表6 輸入特征量的數(shù)據(jù)

5.1 絕緣子污閃概率計算

設腕臂絕緣子的污閃概率為P1，正饋線絕緣子的污閃概率為P2，下錨絕緣子的污閃概率為P3，懸式絕緣子的污閃概率為P4。把2015年1月份5個輸入特征量的數(shù)據(jù)輸入到絕緣子污穢等級預測模型中，預測結(jié)果見表7。由預測的污穢等級上下限對應得參數(shù)a、b的值。

表7 絕緣子污穢等級預測結(jié)果

將表5中4種類型絕緣子的爬電距離代入到式(16),求得系數(shù)A，計算結(jié)果見表8。

表8 4種型號絕緣子的系數(shù)A

將以上數(shù)據(jù)代入到式( 9 )～式(15)中，分別計算出4種類型絕緣子在2月份未來30 d的污閃概率，計算結(jié)果見表9。

表9 絕緣子的污閃概率

5.2 接觸網(wǎng)線路污閃概率計算

從該段的絕緣子運維數(shù)據(jù)庫中調(diào)取4種類型絕緣子的參數(shù)mi(m1=200，m2=210，m3=32，m4=1 010)。查詢該地區(qū)污區(qū)圖可知線路所在區(qū)域為Ⅲ級污區(qū)(0.05～0.10 mg/cm2)，即s=0.05，t=0.10。取ρ=0.08 mg/cm2，將前文算出的相關數(shù)據(jù)代入到式(17)～式(19)中，求得該段距離的接觸網(wǎng)線路在2月份未來30 d發(fā)生污閃的概率為P=0.821 7。對應表4可知在2月份此段接觸網(wǎng)線路容易發(fā)生污閃。應加強巡視，及時做好接觸網(wǎng)絕緣子表面污穢的清掃處理工作，避免接觸網(wǎng)線路污閃概率的增加，降低其給鐵路供電與行車安全帶來的風險。

5.3 分析對比

該供電段對2月份發(fā)生的污閃事故進行了分析總結(jié)，得出造成此段接觸網(wǎng)線路污閃的客觀原因為：多次出現(xiàn)大霧天氣，空氣非常潮濕(濕度：85%；溫度：-10～-5 ℃；風力：微風)。再加上進入冬季以來，天氣干燥久未降水，在該標區(qū)段附近K1931+550處(距離鐵路不到900 m)是個中型規(guī)模的水泥廠，日常排放的臟污顆粒被風吹到絕緣子上致使絕緣部件表面污穢程度較高，在表面污穢層受潮后絕緣電阻和絕緣性能大幅下降，產(chǎn)生絕緣子沿面氣體放電，導致絕緣子表層空氣絕緣擊穿，誘發(fā)接觸網(wǎng)大面積絕緣部件表面閃絡。現(xiàn)場巡視發(fā)現(xiàn)接觸網(wǎng)下行線支柱腕臂瓶底座、正饋線絕緣子與角鋼連接處有明顯的污閃放電痕跡。這與預測結(jié)果基本符合，由此證明本文所提方法的可行性和有效性。

6 結(jié)論

(1) 本文將接觸網(wǎng)絕緣子污穢等級預測方法與污閃概率進行融合，建立接觸網(wǎng)線路污閃短期預報響應機制。對鐵路運行和檢修起到了輔助指導作用，在一定程度上避免了人力、物力、財力的浪費。

(2) 絕緣子污穢程度雖然與降水量、風速、相對濕度、降塵量等氣象因素和自身因素有著密切關系，然而污穢程度也與相關空氣污染參數(shù)及絕緣子外絕緣參數(shù)有著重要聯(lián)系。因此通過監(jiān)測氣象和絕緣子自身因素實現(xiàn)對絕緣子污穢等級進行預測還需要進一步完善，后續(xù)的研究主要集中在綜合考慮氣象和絕緣子自身因素、空氣污染指數(shù)、外絕緣參數(shù)對污穢顆粒受力沉積的影響，實現(xiàn)基于這3類參數(shù)的接觸網(wǎng)絕緣子污穢等級預測。

參考文獻：

[1] 李群山.華中電網(wǎng)典型事故統(tǒng)計分析[J].華中電力，2002,16(5)：33-35．

LI Qunshan. Statistical Analysis of Typical Accidents in Central China Power Grid[J].Central China Electric Power，2002，16(5):33-35.

[2] 張銳,吳光亞,劉亞新,等.光技術(shù)在線監(jiān)測絕緣子鹽密和灰密的實現(xiàn)和應用[J]．高電壓技術(shù)，2010，36(6):1513-1519．

ZHANG Rui，WU Guangya，LIU Yaxin，et al.The Realization and Application of the Optical Technology on Line Monitoring of the Salt and Gray Density of Insulator[J].High Voltage Engineering，2010,36(6):1513-1519.

[3] 王燦林,王黎明,王珂，等.采用光放大器的絕緣子微放電紫外光遙測系統(tǒng)[J].高電壓技術(shù),2007,33(6):137-142．

WANG Canlin，WANG Liming，WANG Ke，et al. An Optical Amplifier for the Measurement of the Micro Discharge of the Insulator Using the Micro Discharge [J].High Voltage Engineering，2007，33(6):137-142.

[4] 李光，王成江，李如峰，等．基于聲發(fā)射波關聯(lián)維數(shù)的絕緣子放電識別研究[J].高壓電器，2013，49(5):79-85．

LI Guang，WANG Chengjiang，LI Rufeng，et al. Study on the Identification of Insulator Discharge Based on Correlation Dimension of Acoustic Emission Wave [J].High Voltage Electrical Apparatus，2013：49(5):79-85.

[5] 劉剛,胡倩楠,陳錫陽,等．由表面電導實現(xiàn)絕緣子污穢程度監(jiān)測的可行性分析[J].高電壓技術(shù)，2012，38(6):1321-1326．

LIU Gang，HU Qiannan，CHEN Xiyang，et al. Feasibility Analysis of Monitoring Insulator Contamination Level by Surface Conductance [J].High Voltage Engineering，2012，38(6):1321-1326.

[6] 熊蘭,劉鈺,林蔭宇,等．模糊變權(quán)法在絕緣子狀態(tài)綜合評判中的應用[J].電力系統(tǒng)及其自動化學報,2010,22(1):96-100．

XIONG Lan，LIU Yu，LIN Yinyu，et al. Application of Fuzzy Variable Weight Method in Comprehensive Evaluation of Insulator State[J].Electric Power System and Its Automation,2010，22(1):96-100.

[7] 張桂南,劉志剛,韓燁,等．接觸網(wǎng)棒式絕緣子故障檢測的快速模糊匹配方法[J].鐵道學報,2013,35(5):28-33．

ZHANG Guinan，LIU Zhigang，HAN Ye，et al. A Fast Fuzzy Matching Method for the Fault Detection of the Contact Net Rod Insulator[J].Journal of Railway，2013，35(5):28-33.

[8] ERENTURK K. Adaptive-Network-Based Fuzzy Inference System Application to Estimate the Flashover Voltage on Insulator[J]. Instrumentation Science and Technology，2013，37(4):446-461.

[9] 何彥諄．覆冰絕緣子閃絡預警特征參量及方法研究[D].重慶：重慶大學，2014:39-51.

[10] 趙世華,蔣興良,張志勁,等.染污玻璃絕緣子泄漏電流特性及其閃絡電壓預測[J].電網(wǎng)技術(shù)，2014，38(2):440-447．

ZHAO Shihua，JIANG Xingliang，ZHANG Zhijin，et al. Prediction of Stained Glass Insulator Leakage Current Characteristic and Flashover Voltage[J].Power System Technology，2014，38(2):440-447.

[11] 石巖,蔣興良,苑吉和．基于RBF網(wǎng)絡的覆冰絕緣子閃絡電壓預測模型[J].高電壓技術(shù)，2009，35(3):591-596．

SHI Yan，JIANG Xingliang，YUAN Jihe. Prediction Model of Icing Insulator Flashover Voltage Based on RBF Network [J]. High Voltage Engineering，2009，35(3):591-596.

[12] NASRAT L S，ALY S. Evaluation of Flashover Voltage Mechanisms on Hydrophobic Polymer Insulators with Artificial Neural Network Approach[J].International Journal of Electrical and Computer Engineering，2012，2(4):487-494.

[13] 吳軍,劉飛．基于最小二乘支持向量機的復雜氣象條件下絕緣子等值附鹽密度的預測[J].陜西電力，2011，39(2):9-12．

WU Jun，LIU Fei. Based on Least Squares Support Vector Machines for the Prediction of the Equivalent Salt Density of Insulators in Complex Meteorological Conditions[J].Shanxi Electric Power，2011，39(2):9-12.

[14] 靳軍．考慮氣候災害的架空線路停運模型研究[D]．北京：華北電力大學，2014:9-18.

[15] 夏青．線路絕緣子交流污穢預測的新特征量研究[D]．重慶：重慶大學，2011:16-49.

[16] 王勝輝．基于紫外成像的污穢懸式絕緣子放電檢測與評估[D]．北京：華北電力大學，2011:68-80.

[17] MAHDJOUBI A，ZEGNINI B，BELKHEIRI M. A LS-SVM(Least Squares Support Vector Machines)Approach for Predicting Critical Flashover Voltage of Polluted Insulators[J].Energy and Power Engineering ，2013，19(2):355-360.

[18] 李曉磊,錢積新．基于分解協(xié)調(diào)的人工魚群優(yōu)化算法研究[J].電路與系統(tǒng)學報，2003，8(1):1-6．

LI Xiaolei,QIAN Jixin. Research on Artificial Fish Swarm Optimization Algorithm Based on Decomposition and Coordination[J].Journal of Circuits and Systems，2003，8(1):1-6.

[19] 馬建偉,張國立,謝宏,等．利用人工魚群算法優(yōu)化前向神經(jīng)網(wǎng)絡[J].計算機應用，2004，24(10):21-23．

MA Jianwei，ZHANG Guoli，XIE Hong, et al.Application of Artificial Fish Swarm Optimization to Neural Network[J].Computer Applications, 2004，24(10):21-23.

[20] 王少華,胡文堂,龔堅鋼,等．浙江電網(wǎng)架空輸電線路絕緣子串的自然積污特性[J].高電壓技術(shù)，2014，40(4):1002-1009．

WANG Shaohua,HU Wentang,GONG Jiangang,et al.Natural Pollution Characteristics of Insulator Strings of Overhead Transmission Lines in Zhejiang Power Grid[J].High Voltage Engineering,2014,40(4): 1002-1009.

[21] 李清,胡文,李黎,等.絕緣子自然積污飽和趨勢及氣象因素影響分析[J].高電壓技術(shù),2014,40(11):3336-3342.

LI Qing，HU Wen，LI Li，et al. Analysis of the Natural Pollution Saturation Trend of Insulators and the Influence of Meteorological Factors [J].High Voltage Engineering,2014,40(11):3336-3342.

[22] 詹清華,阮江軍,鐘連宏．珠三角地區(qū)架空線路絕緣子的自然積污特性[J].高電壓技術(shù)，2010，36(11):2662-2667．

ZHAN Qinghua, RUAN Jiangjun，ZHONG Lianhong.Natural Pollution Characteristics of Overhead Line Insulators in Pearl River Delta Area [J].High Voltage Engineering，2010，36(11):2662-2667.

[23] 盧明,李瑤琴,李黎,等．基于歸一比較法的絕緣子積污特性研究[J].高電壓技術(shù)，2014，40(11):3374-3380．

LU Ming，LI Yaoqin，LI Li，et al. Study on the Fouling Characteristics of Insulators Based on the Normalized Comparison Method[J]. High Voltage Engineering，2014，40(11):3374-3380.

[24] 黃道春,李曉彬,徐濤,等．蒸汽霧中懸式絕緣子表面污層飽和附水量研究[J].高電壓技術(shù)，2014，40(11):3349-3356．

HUANG Daochun, LI Xiaobin，XU Tao，et al. Study on Saturated Water Content of Surface Layer of Suspended Insulator in Steam Fog [J].High Voltage Engineering，2014，40(11):3349-3356.

[25] 郭亞軍．綜合評價理論方法和應用[M].北京：科學出版社，2006:56-72．

[26] 鄭宏宇,鄧銀燕,賀瑞纏．綜合評價中數(shù)據(jù)變換方法的選擇[J].純粹數(shù)學與應用數(shù)學，2010，26(2):319-324．

ZHENG Hongyu, DENG Yinyan, HE Ruichan. Selection of Data Transformation Methods in Comprehensive Evaluation [J].Pure Mathematics and Applied Mathematics，2010，26(2):319-324.

[27] 關志成,王邵武,梁曦東,等．我國電力系統(tǒng)絕緣子污閃事故及其對策[J].高電壓技術(shù)，2000，26(6):37-39．

GUAN Zhicheng，WANG Shaowu，LIANG Xidong，et al. Insulator Pollution Flashover Accident of Power System in China and Its Countermeasures [J].High Voltage Engineering，2000，26(6):37-39.

[28] 高繼卿,楊曉光,董朝陽,等.氣候變化背景下中國北方干濕區(qū)降水資源變化特征分析[J].農(nóng)業(yè)工程學報，2015，31(12):99-110.

GAO Jiqing，YANG Xiaoguang，DONG Chaoyang ，et al. Analysis on the Change of Precipitation Resources in Dry and Wet Areas in Northern China under the Background of Climate Change[J].Journal of Agricultural Engineering，2015，31(12):99-110.

[29] 寧亮,錢永甫.中國年和季各等級日降水量的變化趨勢分析[J].高原氣象，2008，27(5):1010-1020.

NING Liang, QIAN Yongfu. Analyses on Trends of Annual and Seasonal Four Kinds of Daily Precipitation in China[J].Plateau Meteorology，2008，27(5): 1010-1020.

[30] 王靜，楊曉光，李勇，等.氣候變化背景下中國農(nóng)業(yè)氣候資源變化Ⅵ. 黑龍江省三江平原地區(qū)降水資源變化特征及其對春玉米生產(chǎn)的可能影響[J].應用生態(tài)學報，2011，22(6):1511-1522.

WANG Jing，YANG Xiaoguang，LI Yong，et al. Changes of China Agricultural Climate Resources Under the Background of Global Climate Change Ⅵ. Change Characteristics of Precipitation Resource and Its Possible Effect on Maize Production in Sanjiang Plain of Heilongjinag Province[J].Chinese Journal of Applied Ecology,2011，22(6):1511-1522.

[31] SFORZINI M,CORTINA G,MARRONE G.A Statistical Approach for Insulator Design in Polluted Areas[J].IEEE Transactions on Power Apparatus and Systems，1983，102(9):3157-3166.

[32] NAITO K,MIZUNO Y,NAGANAWA W.A Study on Proabilistic Assessment of Contamination Flashover of High-Voltage Insulator[J]. IEEE Transactions on Power Delivery，1995，10(3):1378-1384.

[33] BARRETT J S,GREEN M A.A Statistical Method for Evaluating Electrical Failures[J].IEEE Transactions on Power Delivery，1994，9(3):1524-1530.

[34] ENGELBRECHT C,GUTMAN I,HARTINGS R.A Practical Implementation of Statistical Principles to Dimension AC Line Insulators with Respect to Contaminated Conditions[J].IEEE Transactions on Power Delivery，2007，22(1):667-673.

[35] ENGELBRECHT C S，HARTINGS R，LUNDQUIST J. Statistical Dimensioning of Insulators with Respect to Pollution Conditions[J]. IEEE Proceedings-Generation Transmission and Distribution，2004，151(3):321-326.

[36] HOULGATE R G，LAMBETH P J，ROBERTS W J. The Performance of Insulators at Extra and Ultra High Voltage in Coastal Environment[C]// Presented at 29th General Session of CIGRE. New York：IEEE，1982:33-40.

[37] RISK F A. A Systematic Approach to High-Voltage Insulator Selection for Polluted Environment[C]//Presented at CIGRE 2th Regional in Arab Countries.New York：IEEE，1997:1019-1023.

[38] YE H,ZHANG J,JI Y M,et al.Contamination Accumulation and Withstand Voltage Characteristics of Various Types of Insulators[C]//Proceedings of the 7th International Conference on Properties and Applications of Properties Materials.Nagoya：ELSEVIER，2003:1019-1023.

[39] TOPALIS F V,GONOS I F,STATHOPULOS I A.Dielectric Behavior of Polluted Porcelain Insulators[J].IEEE Proceedings-Generation Transmission and Distribution，2001，148(4):269-274.