大型火電廠主、輔機數(shù)量較多,結(jié)構(gòu)復雜,其作用是完成從熱能到機械能再到電能的轉(zhuǎn)換過程。設(shè)備與設(shè)備之間的耦合性、系統(tǒng)的復雜性等因素決定了火電廠是一個高故障率和故障危害性很大的生產(chǎn)場所,這些故障都將造成重大的經(jīng)濟損失和社會后果,因此應用設(shè)備故障預警及狀態(tài)監(jiān)測技術(shù),掌握關(guān)鍵設(shè)備的正常運行狀態(tài),及時發(fā)現(xiàn)易損部件的工況惡化趨勢,及時進行維護和檢修,對提高設(shè)備的可靠性和安全性十分必要^[1-5]。

目前,發(fā)電設(shè)備預警和預知維修的關(guān)鍵技術(shù)主要有狀態(tài)監(jiān)測技術(shù)、故障診斷與預測技術(shù)和維修分析與決策,而故障模式識別是這些技術(shù)的基礎(chǔ)。本文采用K-means算法與基于EM（期望最大）的高斯混合模型（Gaussian Mixture Model,GMM）算法相結(jié)合的聚類分析算法進行故障模式識別,實踐結(jié)果證明該算法能夠很好地用于發(fā)電設(shè)備智能診斷預測。

 1 算法介紹

1.1 K-means算法

K-means算法,也被稱為K-平均或K-均值算法,是一種得到廣泛應用的聚類分析方法^[6]。該算法的主要思想是：將各個聚類子集內(nèi)所有數(shù)據(jù)樣本的均值作為該聚類的代表點,通過迭代的方法逐次更新各聚類中心的值,把數(shù)據(jù)集劃分為不同的類別,最終目標是使評價聚類性能準則的函數(shù)達到最優(yōu),使生成的每個聚類（簇）內(nèi)緊湊,類間獨立^[7-10]。

假設(shè)要把樣本集X={x₁, x₂, x₃,…, x_n}分為K個類別,算法的過程描述如下：

算法：K-means算法,劃分并計算基于簇中對象的平均值;

輸入：簇的數(shù)目K和包含n個對象的數(shù)據(jù)庫;

輸出：K個簇,使平方誤差總和最小。

算法步驟：

1）為每個聚類確定一個初始聚類中心,形成K個初始聚類中心;

2）將樣本集中的樣本按照最小距離原則分配到最鄰近的聚類;

3）利用每個聚類中的樣本均值作為新的聚類中心;

4）重復步驟2和步驟3,直到聚類不再變化;

5）迭代結(jié)束,得到K個聚類。

K-means算法的優(yōu)點是算法快速、簡單,對大數(shù)據(jù)集有較高的效率并且是可伸縮的。缺點是需要根據(jù)初始聚類中心來確定一個初始劃分,然后對初始劃分進行優(yōu)化。在K-means 算法中,多維空間相似性度量基于歐氏距離進行計算,并不能準確反映多維空間點中的相似情況。

1.2 高斯混合模型（GMM）

1.2.1 多維空間單高斯模型（SGM）

假設(shè)有一組在多維空間（假設(shè)維度為D）的點X_i,i=1,2,…,n,若這些點的分布近似橢球狀,則可以用高斯密度函數(shù)g(X,μ,Σ)來描述產(chǎn)生這些點的密度函數(shù)：

 （2）

由于x_i為己知點,可通過p(X;μ,Σ)來求解μ和Σ。這類問題通常以極大似然估計的方法來進行求解,即在p(X;μ,Σ)為最大值條件下（樣本點x_i已經(jīng)發(fā)生,故可認為p(X;μ,Σ)樣本x_i發(fā)生的概率最大）,求得μ和Σ。根據(jù)極值理論,可以通過對μ和Σ分別求導的方式求解。

1.2.2 高斯混合模型（GMM）

SGM樣本點的分布為橢球狀。而實際中,大部分樣本點分布不可能是或近似橢球狀,因此引用了高斯混合模型（GMM）。GMM是一種常用的描述混合密度函數(shù)分布的模型,可視為多個高斯分布以一定比例混合而成^[11-14],GMM模型可描述為：

Fig.1 Solution program