蘇 寧，李亞錦，張婉瑩，張國新，于大洋，劉英男

（1.海南電網(wǎng)有限責(zé)任公司瓊海供電局，海南 瓊海 571400；2.山東大學(xué)電氣工程學(xué)院，山東 濟(jì)南 250061）

0 引言

人工智能的紅外圖像識別技術(shù)［1］能夠識別紅外測溫圖像中的發(fā)熱點(diǎn)并計(jì)算溫差指標(biāo)，在機(jī)器人巡檢中能夠代替人工發(fā)現(xiàn)缺陷［2］。但針對變電設(shè)備發(fā)熱的診斷是在缺陷發(fā)生后，單純依靠溫差閾值［3］僅能判斷設(shè)備的當(dāng)前狀態(tài)，無法對其發(fā)熱趨勢進(jìn)行跟蹤預(yù)測。環(huán)境溫度和負(fù)載升高，有可能導(dǎo)致設(shè)備絕對溫度的進(jìn)一步升高，缺陷由一般發(fā)展為嚴(yán)重缺陷，但在迎峰度夏、度冬封網(wǎng)期間，設(shè)備不易停電處理。因此需要新的手段對變電設(shè)備熱點(diǎn)溫度、相對溫差進(jìn)行預(yù)測，在迎峰度夏之前進(jìn)行預(yù)見性檢修。

設(shè)備狀態(tài)并不能單純依據(jù)紅外測溫指標(biāo)來評判，還要結(jié)合設(shè)備的電壓等級、廠家型號、運(yùn)行年限、環(huán)境天氣等運(yùn)行因素做出綜合評價(jià)和判斷［4-8］。這種分析需要依賴運(yùn)維檢修人員的經(jīng)驗(yàn)和專業(yè)水平，需要長期積累，并且因人而異，難以實(shí)現(xiàn)標(biāo)準(zhǔn)化，難以保證延續(xù)性。而且在運(yùn)維實(shí)際中，這些多維度、多來源的數(shù)據(jù)往往不完整，存在大量錯(cuò)誤信息，給缺陷的綜合評判帶來困難。

當(dāng)前人工智能的關(guān)鍵技術(shù)已應(yīng)用到設(shè)備運(yùn)維檢修中［9-10］，并取得了一定的診斷效果。針對設(shè)備狀態(tài)預(yù)測方法包括基于時(shí)間序列的回歸分析模型［11］、支持向量機(jī)［12］、多維度評價(jià)［13-14］和深度學(xué)習(xí)［15-17］等算法。但基于時(shí)間序列的回歸模型是以線性模型為基礎(chǔ)，難以滿足變電設(shè)備運(yùn)行狀態(tài)多維非線性特性。支持向量機(jī)雖然具有一定的泛化能力，但隨著變電設(shè)備運(yùn)行狀態(tài)特征因素的增多，將影響模型預(yù)測精度。而基于深度學(xué)習(xí)的預(yù)測方法使模型具有非線性表達(dá)能力，其中長短期記憶網(wǎng)絡(luò)（Long Short Term Memory，LSTM）［18］可以更好地實(shí)現(xiàn)長時(shí)間序列的分析與趨勢預(yù)測。文獻(xiàn)［19］利用LSTM 的方法對變壓器的運(yùn)行狀態(tài)做出了預(yù)測，文獻(xiàn)［20］驗(yàn)證了LSTM 預(yù)測模型在故障時(shí)間序列中具有很強(qiáng)的適用性。上述研究為大數(shù)據(jù)背景下的變電設(shè)備運(yùn)行狀態(tài)評估提供了思路。

在上述背景下，以變電設(shè)備歷史運(yùn)行功率、熱點(diǎn)溫度、環(huán)境溫度等檢測數(shù)據(jù)為基礎(chǔ)，利用LSTM 算法預(yù)測多源因素影響下設(shè)備熱點(diǎn)溫度，實(shí)現(xiàn)設(shè)備熱點(diǎn)溫度發(fā)展趨勢和缺陷嚴(yán)重程度的動(dòng)態(tài)預(yù)測，對變電設(shè)備運(yùn)維檢修方案的制訂提供參考。

1 LSTM算法

采用LSTM 模型建立輸入映射到輸出的內(nèi)在關(guān)聯(lián)。LSTM 是一種遞歸型神經(jīng)網(wǎng)絡(luò)，精確控制LSTM記憶單元上一個(gè)時(shí)間步長到下一個(gè)時(shí)間步長的步距，解決了梯度消失或者爆炸的問題。LSTM單元結(jié)構(gòu)如圖1所示。

圖1 LSTM的拓?fù)浣Y(jié)構(gòu)

圖1中，Yt和Qt是輸入向量和隱藏層向量，Bt是記憶單元的狀態(tài)向量，Bt通過Yt和Qt判斷是否從上一時(shí)間段更新到下一時(shí)間段，jt表示輸入門，gt表示遺忘門，ut表示輸出門。

遺忘門gt的值可以表示為

式中：Wg和cg分別為遺忘門gt的系數(shù)和偏移；α是激活函數(shù)。式（1）計(jì)算出遺忘門gt能夠沿著圖中的信息流逐點(diǎn)將gt和Bt相乘來控制記憶單元中包含的信息量。

輸入門jt的值可以表示為

式中：Wj和cj分別為輸入門jt的系數(shù)和偏移，由神經(jīng)網(wǎng)絡(luò)訓(xùn)練得到。式（2）計(jì)算出輸入門jt能夠沿著圖中的信息流逐點(diǎn)將jt和相乘來控制記憶單元中暫時(shí)包含的信息量。的計(jì)算公式為

式中：Wb和cb分別為～Bt的系數(shù)和偏移，由神經(jīng)網(wǎng)絡(luò)訓(xùn)練得到。

最終記憶單元在此時(shí)刻的狀態(tài)由式（4）決定。

更新后的新值通過輸出門決定當(dāng)前神經(jīng)網(wǎng)絡(luò)的隱藏層，同樣，輸出門ut的表達(dá)式為

式中：Wu和cu分別為輸出門ut的系數(shù)和偏移，由神經(jīng)網(wǎng)絡(luò)訓(xùn)練得到。

隱藏層Qt由輸出門ut和Bt逐點(diǎn)相乘得到，即為

綜上可知，LSTM 網(wǎng)絡(luò)可以動(dòng)態(tài)更新網(wǎng)絡(luò)蘊(yùn)含值信息，并從當(dāng)前時(shí)刻傳遞到下一時(shí)刻。

2 基于LSTM的熱點(diǎn)溫度預(yù)測流程

算法模型具體應(yīng)用流程為：

1）模型輸入輸出變量的選取。根據(jù)熱點(diǎn)溫度影響因素，選取正常相溫度、環(huán)境溫度、季節(jié)、運(yùn)行電流/功率，作為預(yù)測模型的輸入變量，輸出變量為熱點(diǎn)溫度?；跉v史數(shù)據(jù)構(gòu)建數(shù)據(jù)樣本集。將數(shù)據(jù)集數(shù)據(jù)進(jìn)行歸一化，并按照一定比例劃分成訓(xùn)練集和測試集。

2）將訓(xùn)練數(shù)據(jù)輸入LSTM 預(yù)測模型中，通過LSTM 對訓(xùn)練集中的數(shù)據(jù)進(jìn)行建模，確定LSTM 網(wǎng)絡(luò)輸入的節(jié)點(diǎn)數(shù)和隱含層神經(jīng)元個(gè)數(shù)以及輸出結(jié)果，由tanh 和sigmoid 函數(shù)的輸出相乘得到網(wǎng)絡(luò)的輸出熱點(diǎn)溫度。tanh 和sigmoid 激活函數(shù)用于將數(shù)值控制在-1至1之間，可表示為：

在每個(gè)層級，LSTM 模型訓(xùn)練采用隨時(shí)間反向傳播（Back Propagation Through Time，BPTT）算法提升其對時(shí)間序列的學(xué)習(xí)能力，選取對數(shù)似然函數(shù)作為損失函數(shù)。輸入層的神經(jīng)元個(gè)數(shù)選取為4 個(gè)，輸入的時(shí)序變量為｛x1，x2，…，xT｝，其中在采樣時(shí)刻t神經(jīng)元的輸入時(shí)間序列為xt［It（1），It（2），It（3），It（4）］。輸出變量為熱點(diǎn)溫度。神經(jīng)網(wǎng)絡(luò)選擇2 層LSTM，最后連接一個(gè)全連接層后通過一個(gè)激活函數(shù)作為最終的輸出。

訓(xùn)練過程中為評估預(yù)測值和實(shí)際值之間的偏差大小，選取平均絕對誤差作為損失函數(shù)評價(jià)訓(xùn)練效果。以損失函數(shù)值最小為目標(biāo)不斷訓(xùn)練優(yōu)化神經(jīng)網(wǎng)絡(luò)參數(shù)，直到達(dá)到較好效果并保存訓(xùn)練模型。平均絕對誤差計(jì)算方法為

式中：m為樣本數(shù)；yt和?t為t時(shí)刻熱點(diǎn)溫度的實(shí)際值和預(yù)測值。

3）將t 時(shí)刻的正常相溫度、環(huán)境溫度、季節(jié)、運(yùn)行電流狀態(tài)參數(shù)輸入步驟2）中的訓(xùn)練模型，得到t+1 時(shí)刻熱點(diǎn)溫度的預(yù)測值。

4）根據(jù)熱點(diǎn)溫度和相對溫差判斷缺陷性質(zhì)。相對溫差計(jì)算公式為

式中：T1為預(yù)測的熱點(diǎn)溫度；T2正常相溫度；T0為環(huán)境溫度。T2和T0選取歷史同期時(shí)的檢測數(shù)據(jù)。

搜集2015年至2020年的多組設(shè)備同一部位發(fā)熱數(shù)據(jù)進(jìn)行分析，數(shù)據(jù)間隔為1 個(gè)月。根據(jù)檢測時(shí)間將春夏秋冬四個(gè)季節(jié)因素分別轉(zhuǎn)換成1、2、3、4 量化的變量。輸入的樣本數(shù)據(jù)如表1所示。

表1 數(shù)據(jù)樣本

基于LSTM 算法的預(yù)測模型建立在Python3.6 環(huán)境以及Keras、Scikit?learn、Tensorflow 等函數(shù)庫。將訓(xùn)練樣本集中的狀態(tài)監(jiān)測信息作為網(wǎng)絡(luò)的輸入，對網(wǎng)絡(luò)模型進(jìn)行離線訓(xùn)練，首先對參數(shù)進(jìn)行初始化，設(shè)置學(xué)習(xí)步長為0.01，輸入維度為4，輸出層節(jié)點(diǎn)個(gè)數(shù)為1，權(quán)重參數(shù)矩陣所有元素初始化為（-1，1），將誤差限值設(shè)置為1×10-2，設(shè)置迭代次數(shù)為500次，對樣本進(jìn)行離線訓(xùn)練，從而得到熱點(diǎn)溫度預(yù)測模型。熱點(diǎn)溫度預(yù)測模型流程如圖2所示。

圖2 熱點(diǎn)溫度預(yù)測模型流程

以某800 kV 換流站A、B、C 三相歷史數(shù)據(jù)為例，采用監(jiān)測和運(yùn)維數(shù)據(jù)共計(jì)1 745條數(shù)據(jù)，前1 545項(xiàng)數(shù)據(jù)作為訓(xùn)練集，后200 項(xiàng)數(shù)據(jù)作為測試集輸入預(yù)測模型。選取繞組正常相溫度、環(huán)境溫度、整流側(cè)有功功率為狀態(tài)變量，預(yù)測每一相的繞組溫度。

以預(yù)測A相繞組溫度變化為例，訓(xùn)練數(shù)據(jù)集的損失值收斂速度很快，如圖3 所示，迭代20 次之后損失值便接近預(yù)期水平，迭代40 次之后損失值穩(wěn)定在最低水平0.008 4，滿足最低誤差要求，所有訓(xùn)練樣本誤差小于閾值。

圖3 A相繞組溫度預(yù)測損失值

如圖4 所示為A 相繞組溫度真實(shí)值和預(yù)測值的情況，可以看出真實(shí)值和預(yù)測值的波形和走勢基本一致，重合程度高，預(yù)測值的振幅略小于真實(shí)值，尤其當(dāng)真實(shí)值的振動(dòng)幅度較大時(shí)，預(yù)測值難以達(dá)到真實(shí)值的變化速度，當(dāng)真實(shí)值波動(dòng)較為穩(wěn)定時(shí)，曲線貼合程度較高?？傮w來看，預(yù)測具有較高的準(zhǔn)確性。

圖4 A相繞組溫度真實(shí)值和預(yù)測值

3 實(shí)例分析

將訓(xùn)練模型應(yīng)用到變電站設(shè)備熱點(diǎn)溫度趨勢預(yù)測中。2019 年8 月對某220 kV 變電站1 號變壓器進(jìn)行紅外測溫，圖5為紅外檢測圖譜。發(fā)現(xiàn)110 kV側(cè)套管B 相柱頭存在發(fā)熱，最高溫度為71 ℃，相對溫差為69.5%，判斷為一般缺陷。10 月進(jìn)行復(fù)測，最高溫度為73.6 ℃，相對溫差為71.5%，溫度雖然有所上升，但仍未達(dá)到嚴(yán)重缺陷。隨后按照檢測周期繼續(xù)跟蹤檢測，不同時(shí)間的檢測數(shù)據(jù)如表2所示。

圖5 檢測圖譜

分析表2 的檢測數(shù)據(jù)可知，環(huán)境溫度和負(fù)載的降低，導(dǎo)致設(shè)備的絕對溫度降低，但發(fā)熱缺陷仍未消失。且隨著環(huán)境溫度和負(fù)載的升高，相對溫差已達(dá)到79.59%，接近嚴(yán)重缺陷指標(biāo)80%。利用上述訓(xùn)練好的LSTM熱點(diǎn)溫度預(yù)測模型，根據(jù)歷史同期正常相溫度、環(huán)境溫度、季節(jié)、運(yùn)行電流數(shù)據(jù)輸入模型，兩組輸入數(shù)據(jù)如表3 所示。根據(jù)預(yù)測模型最終得到相應(yīng)的熱點(diǎn)溫度，分別為82.3 ℃和94.3 ℃。計(jì)算相對溫差分別為81.5%、58.6%。根據(jù)缺陷嚴(yán)重程度判斷標(biāo)準(zhǔn)，90 ℃≤熱點(diǎn)溫度≤130 ℃或δ≥80%但熱點(diǎn)溫度未達(dá)緊急缺陷溫度值，判斷1 號主變壓器110 kV 側(cè)套管B 相柱頭發(fā)熱屬于嚴(yán)重缺陷。

表2 1號主變壓器110 kV側(cè)套管B相柱頭數(shù)據(jù)樣本

表3 歷史同期檢測數(shù)據(jù)

根據(jù)LSTM 熱點(diǎn)溫度預(yù)測可知，隨著環(huán)境溫度和負(fù)載的升高，缺陷由一般發(fā)展為嚴(yán)重缺陷。而在迎峰度夏期間，不易停電處理。通過對變電設(shè)備熱點(diǎn)溫度、相對溫差進(jìn)行預(yù)測，在迎峰度夏之前可進(jìn)行預(yù)見性檢修。根據(jù)運(yùn)維知識庫可知，該部位發(fā)熱的原因可能是柱頭老化、松動(dòng)，按照檢修工藝對柱頭連接處進(jìn)行打磨，涂抹導(dǎo)電膏，更換并緊固螺栓，消除缺陷。

將上述方法應(yīng)用到其他設(shè)備發(fā)熱預(yù)警上，檢測和預(yù)測數(shù)據(jù)如表4 所示。110 kV 變電站1 號主變壓器35 kV 側(cè)B 相套管接線柱隨著季節(jié)下降，熱點(diǎn)溫度未下降，該缺陷已為嚴(yán)重缺陷，需在秋檢期間盡快消除，否則在迎峰度冬期間產(chǎn)生安全隱患。

表4 110 kV變電站1號主變壓器35 kV側(cè)B相套管接線柱數(shù)據(jù)

圖6 110 kV變電站1號主變壓器35 kV側(cè)B相套管接線柱數(shù)據(jù)

4 結(jié)語

在充分考慮變電設(shè)備發(fā)熱受運(yùn)行工況影響的基礎(chǔ)上，利用LSTM 構(gòu)建變電設(shè)備熱點(diǎn)溫度預(yù)測模型，解決無法跟蹤預(yù)測發(fā)熱趨勢的運(yùn)維現(xiàn)狀。通過熱點(diǎn)溫度預(yù)測模型，可獲取在環(huán)境溫度和負(fù)載升高的工況下熱點(diǎn)溫度狀態(tài)值，根據(jù)缺陷判斷導(dǎo)則判定缺陷性質(zhì)，從而提前消缺。仿真結(jié)果表明，建立的溫度預(yù)測模型可有效的預(yù)測熱點(diǎn)溫度發(fā)展趨勢，輔助現(xiàn)場人員在迎峰度夏之前進(jìn)行預(yù)見性檢修。