李　莉，熊　煒，何　杰 ，袁旭峰，鄒曉松

(1. 貴州電力試驗研究院，貴州 貴陽550001；2. 貴州大學(xué) 電氣工程學(xué)院，貴州 貴陽550025)

基于改進(jìn)灰色馬爾可夫模型的設(shè)備故障率預(yù)測

李莉1，熊煒2，何杰2，袁旭峰2，鄒曉松2

(1. 貴州電力試驗研究院，貴州 貴陽550001；2. 貴州大學(xué) 電氣工程學(xué)院，貴州 貴陽550025)

摘要：電力設(shè)備故障率具有時變性、隨機性、回退等特點，預(yù)測難度大,因此在灰色GM(1,1)模型的基礎(chǔ)上，采用模糊C均值聚類方法對GM(1,1)模型擬合值的誤差序列進(jìn)行狀態(tài)劃分；通過計算誤差序列的狀態(tài)轉(zhuǎn)移概率矩陣，建立了電力設(shè)備故障率的灰色馬爾可夫預(yù)測模型。該模型既考慮了GM(1,1)模型較強的處理單調(diào)數(shù)列的特性，又計及了通過狀態(tài)轉(zhuǎn)移概率矩陣的變換提取數(shù)據(jù)隨機波動響應(yīng)的特點，避免了最大概率狀態(tài)不為實際狀態(tài)而出現(xiàn)最差的預(yù)測結(jié)果現(xiàn)象。通過實例證明，基于模糊C均值聚類的灰色馬爾可夫模型預(yù)測結(jié)果優(yōu)于傳統(tǒng)的GM(1,1)模型和基于K均值聚類的灰色馬爾可夫預(yù)測模型，具有較高的預(yù)測精度。

關(guān)鍵詞：電力設(shè)備故障率；灰色理論；馬爾可夫鏈；模糊C均值聚類

中圖分類號：TM774

文獻(xiàn)標(biāo)識碼:A

DOI:10.3969/j.issn.1672-0792.2015.08.004

收稿日期：2015-05-21。

基金項目：貴州省科技廳基金(黔科合J字[2011]2060號)。

作者簡介：李莉( 1968-)，女，高級工程師，主要從事電力系統(tǒng)可靠性研究,E-mail:420034562@qq.com。

Abstract：Based on the characteristics of time variation and randomness of power equipment failure rate, this paper proposes a state division method for fitted values error series of GM (1,1) model by fuzzy C means clustering method based on grey GM (1,1) model. By computing the state transition probability matrix of error series, a grey Markov prediction model is established. This model analyzes high ability of GM (1,1) model processing monotone series then explores random fluctuation responses when extracting data through transformation of state transition probability matrix. So the worst prediction results when the maximum probability state is not the actual state can be avoided. The actual case study validates that the grey Markov forecasting model based on fuzzy C means clustering method has better performance than conventional GM (1,1) model and grey Markov forecasting model based on K means clustering, and has higher forecast precision.

Keywords：failure rate of power equipment; grey theory; Markov chain; fuzzy C means clustering

0引言

電力設(shè)備故障率預(yù)測是電力系統(tǒng)可靠性評估及狀態(tài)檢修決策優(yōu)化的基礎(chǔ)，而電力設(shè)備故障率與設(shè)備所處的地理位置、氣候條件、運行條件等眾多因素相關(guān)，呈現(xiàn)出時變性和隨機性的特點，以及檢修后故障率回退現(xiàn)象，導(dǎo)致電力設(shè)備故障率的預(yù)測難度增大。

目前電力設(shè)備故障率預(yù)測方法主要包括：基于歷史數(shù)據(jù)的方法、基于設(shè)備役齡和健康狀態(tài)的方法和基于人工智能的方法。其中基于歷史數(shù)據(jù)的預(yù)測方法采用統(tǒng)計方法進(jìn)行建模，要求具備大量的歷史數(shù)據(jù)，且數(shù)據(jù)不能有周期性變化或突變[1，2]?；谠O(shè)備役齡和健康狀態(tài)的故障率預(yù)測模型，需要對設(shè)備的健康狀態(tài)進(jìn)行連續(xù)跟蹤[3～5]?；谌斯ぶ悄艿碾娏υO(shè)備故障預(yù)測方法主要是利用BP,RBF等人工神經(jīng)網(wǎng)絡(luò)方法研究設(shè)備故障率變化規(guī)律，需要大量的統(tǒng)計樣本，預(yù)測結(jié)果具有很大的主觀性和隨機性[6，7]。因此，如何提高電力設(shè)備故障率的預(yù)測精度，是電力系統(tǒng)可靠性評估及開展?fàn)顟B(tài)檢修的研究熱點。

灰色模型對于序列較短并具有明確上升趨勢的數(shù)據(jù)來說，具有較好的預(yù)測效果，但對于隨機波動性較大的數(shù)據(jù)序列預(yù)測精度較低，而馬爾可夫理論通過狀態(tài)轉(zhuǎn)移概率描述各隨機因素的影響以及各狀態(tài)之間轉(zhuǎn)移的內(nèi)在規(guī)律，可以彌補灰色模型的不足[8～11]。為此，本文提出了基于模糊C均值聚類的灰色馬爾可夫模型對電力設(shè)備故障率進(jìn)行預(yù)測，以便能充分利用歷史數(shù)據(jù)，提高隨機波動性較大的數(shù)列預(yù)測精度，拓寬灰色理論的應(yīng)用范圍。

1灰色GM(1,1)模型

(1)設(shè)電力設(shè)備故障率原始數(shù)據(jù)序列為：

(1)

其中：

(2)對原始數(shù)列進(jìn)行一次累加：

(2)

其中

(3)生成緊鄰均值序列：

(3)

其中：

(4)建立GM(1,1)模型白化方程：

(4)

其中：a和b是待確定參數(shù)。

(5)GM(1,1)模型灰色微分方程為：

(5)

(6)對應(yīng)的時間響應(yīng)方程為：

(6)

(7)參數(shù)a和b的最小二乘估計為：

(7)

式中：

(8)對預(yù)測結(jié)果進(jìn)行累減還原得到預(yù)測數(shù)列：

(8)

其中：k=1,2,…,n-1

(9)計算預(yù)測數(shù)列與原始數(shù)列的殘差：

(9)

式中：

由上可見，GM(1,1)模型是將原始數(shù)列經(jīng)一次累加后進(jìn)行預(yù)測的，由于一次累加數(shù)列具有單調(diào)性，因此GM(1,1)模型適用于預(yù)測指數(shù)變化規(guī)律的數(shù)據(jù)，對隨機波動型數(shù)據(jù)則預(yù)測誤差較大。

2灰色馬爾可夫鏈預(yù)測模型

本文采用的灰色馬爾可夫預(yù)測基本思路，是通過GM(1,1)模型的殘差序列計算狀態(tài)轉(zhuǎn)移矩陣，根據(jù)狀態(tài)轉(zhuǎn)移矩陣對GM(1,1)模型結(jié)果進(jìn)行修正，以提高預(yù)測結(jié)果的準(zhǔn)確性。具體步驟為：

(1)將殘差序列E劃分為n個狀態(tài)，任一狀態(tài)可表示為?i(i=1,2,…,n)。

為提高灰色馬爾可夫鏈的預(yù)測精度，使殘差序列狀態(tài)劃分更科學(xué)、更有效，本文采用模糊C均值聚類算法確定殘差序列E的分布狀態(tài)，其算法表述為：

(10)

(11)

(12)

且滿足約束條件

(13)

式中：m>1；c是聚類數(shù)量；n是聚類空間的樣本數(shù)；U={μij}表示1個c×n維的矩陣，μij為樣本點ej隸屬于第i類的程度；J=(U,V)是類內(nèi)誤差的加權(quán)平方和目標(biāo)函數(shù)。dij=‖νi-ej‖為第i個聚類中心Vi與第j個數(shù)據(jù)點ej間的歐氏距離，F(xiàn)CM算法通過迭代式使目標(biāo)函數(shù)最小化，從而達(dá)到最佳聚類。

(2)根據(jù)殘差序列E的狀態(tài)分布，計算n步轉(zhuǎn)移概率矩陣：

(14)

(4)根據(jù)未來不同時刻殘差所處的狀態(tài)，對灰色預(yù)測結(jié)果進(jìn)行修正得到預(yù)測區(qū)間。依據(jù)文獻(xiàn)[11]中灰色馬爾可夫預(yù)測方法的討論，本文采用式(15)計算殘差序列的預(yù)測值，以避免馬爾可夫過程中最大概率狀態(tài)不為實際狀態(tài)而出現(xiàn)最差的預(yù)測結(jié)果現(xiàn)象。

(15)

(5)更新原始數(shù)據(jù)，重復(fù)上述步驟，直至完成預(yù)測目標(biāo)。

3算例分析

以某電網(wǎng)2011年6月～2013年12月共31個月220 kV輸電線路故障率為原始數(shù)據(jù)序列，如圖1所示。從圖中可看出，220 kV輸電線路故障率具有時變性、隨機性和回退等特點，不能采用時間序列法進(jìn)行準(zhǔn)確預(yù)測。

圖1　220 kV輸電線路故障率數(shù)據(jù)序列

本文以前26個月的數(shù)據(jù)為基礎(chǔ)數(shù)據(jù)，建立灰色馬爾可夫預(yù)測模型進(jìn)行預(yù)測和誤差分析。

(1)GM(1,1)預(yù)測

采用GM(1,1)預(yù)測模型對220 kV輸電線路故障率進(jìn)行預(yù)測，時間響應(yīng)方程為：

x(1)(k+1)=[x(0)(1)-132.716]e0.0026t+132.716

擬合曲線如圖2所示。

圖2　220 kV輸電線路故障率GM(1,1)預(yù)測結(jié)果

從圖2可看出， GM(1,1)模型殘差也具有時變、隨機分布等特點，能體現(xiàn)設(shè)備故障率的變化規(guī)律。

(2)對殘差序列進(jìn)行狀態(tài)劃分

為進(jìn)行對比分析，本文分別采用K均值聚類和模糊C均值聚類對灰色預(yù)測數(shù)據(jù)殘差序列進(jìn)行狀態(tài)劃分，當(dāng)狀態(tài)數(shù)量為4時，殘差序列聚類中心如表1所示，狀態(tài)分布情況如圖3所示。

表1　兩種聚類方法獲得的殘差序列狀態(tài)中心

圖3　GM(1,1)預(yù)測數(shù)據(jù)殘差序列狀態(tài)劃分結(jié)果

(3)灰色馬爾可夫預(yù)測模型

分別對上述2種狀態(tài)劃分結(jié)果，采用馬爾可夫過程對殘差序列進(jìn)行擬合，擬合誤差如表2所示，擬合曲線如圖4所示。

表2　26個原始數(shù)據(jù)的灰色馬爾可夫擬合誤差

圖4　26個原始數(shù)據(jù)的擬合結(jié)果與實際值的比較

從表2和圖4可看出，2種聚類方法的灰色馬爾可夫擬合結(jié)果趨勢與實際符合，但基于模糊C均值聚類的灰色馬爾可夫擬合結(jié)果較基于K均值聚類的誤差更小，擬合結(jié)果更精確。

(4)灰色馬爾可夫預(yù)測結(jié)果

利用基于模糊C均值聚類的灰色馬爾可夫預(yù)測模型及ARMA模型[12]分別對后5年設(shè)備故障率進(jìn)行預(yù)測，預(yù)測結(jié)果如圖5所示，誤差如表3所示。

圖5　預(yù)測結(jié)果與實際值的比較

表3　灰色馬爾可夫模型及ARMA模型預(yù)測誤差

從表3和圖5可看出，當(dāng)原始數(shù)據(jù)隨機波動性增大時，基于模糊C均值聚類的灰色馬爾可夫預(yù)測結(jié)果趨勢與實際基本相符，均方差較基于K均值聚類的灰色馬爾可夫預(yù)測、ARMA模型預(yù)測小，可見，具有較高的預(yù)測精度。

4結(jié)論

由于電力設(shè)備故障率具有時變性、隨機性、回退等特點，本文提出基于模糊C均值聚類的灰色馬爾可夫預(yù)測方法對電力設(shè)備故障率進(jìn)行預(yù)測，通過預(yù)測結(jié)果的比較，該模型預(yù)測優(yōu)于傳統(tǒng)的GM(1,1)模型、基于K均值聚類的灰色馬爾可夫預(yù)測模型和ARMA預(yù)測模型，具有較高的預(yù)測精度，為電力系統(tǒng)可靠性的評估及電力設(shè)備狀態(tài)檢修的開展提供了可靠的科學(xué)依據(jù)。但在該模型中，如何引入氣候、檢修后故障率回退等因素的影響，提高設(shè)備故障率非等時間間隔序列的預(yù)測精度是下一步研究的重點。

參考文獻(xiàn):

[1]茹斌，張?zhí)靷ィ跤钚?基于小波去噪及ARMA模型的故障率預(yù)測方法研究[J].測控技術(shù),2014,33(10):43-46,50.

[2]顧潔，杜漸,秦杰，等.基于全壽命周期成本的變電站設(shè)備檢修模糊規(guī)劃模型[J].電力系統(tǒng)自動化,2014,38(1):44-50.

[3]孫鵬，陳紹輝，張彩慶. 基于Marquardt法參數(shù)估計的變電設(shè)備壽命周期故障率評估[J]. 電力系統(tǒng)保護(hù)與控制，2012，40(1)：85-90.

[4]潘樂真，張焰，俞國勤，等.狀態(tài)檢修決策中的電氣設(shè)備故障率推算[J].電力自動化設(shè)備，2010，30(2)：91-94.

[5]劉若愚.佛山電網(wǎng)輸變電設(shè)備的狀態(tài)檢修[D].廣州：華南理工大學(xué)，2014.

[6]胡澤江，劉宗兵，束洪春.基于BP神經(jīng)網(wǎng)絡(luò)的設(shè)備故障率獲取[J].云南水力發(fā)電，2008，24(1)：85-88.

[7]楊森，孟晨，王成.基于改進(jìn)灰色神經(jīng)網(wǎng)絡(luò)的故障預(yù)測方法研究[J].計算機應(yīng)用研究，2013，30(12):3625-3628.

[8]劉思峰，謝乃明，黨耀國，等. 灰色系統(tǒng)理論及其應(yīng)用[M].北京：科學(xué)出版社，2010.

[9]章偉，鄧院昌. 基于灰色-馬爾可夫鏈的短期風(fēng)速及風(fēng)電功率預(yù)測[J]. 中國電力，2013，46(2)：98-102.

[10]吳文可，文福拴，薛禹勝，等. 基于馬爾可夫鏈的電力系統(tǒng)連鎖故障預(yù)測[J]. 電力系統(tǒng)自動化，2013，37(5)：29-37.

[11]李小芳，基于改進(jìn)灰色馬爾可夫模型的城市污水量和COD總量預(yù)測[D].天津：天津大學(xué)，2012.

[12]Moradkhani A, Haghifam M R, Mohammadzadeh M. Failure rate modelling of electric distribution overhead lines considering preventive maintenance [J], Generation Transmission & Distribution IET, 2014, 8(6):1028-1038.

Prediction of the Equipment Failure Rate Based on Improved Grey Markov Model

Li Li1，Xiong Wei2，He Jie2，Yuan Xufeng2，Zou Xiaosong2

(1. Guizhou Institute of Power Equipment Test and Research, Guiyang 550001, China; 2. Electrical Engineering School, Guizhou University, Guiyang 550025, China)

歡迎投稿,歡迎訂閱

投稿網(wǎng)址：http://www.dlkxygc.com