張雷，王光華，李金鑠，耿宏賢，戴志輝

(1.國(guó)網(wǎng)河北省電力有限公司保定供電分公司，河北 保定 071000；2.華北電力大學(xué)電氣與電子工程學(xué)院，河北 保定 071003)

0 引言

繼電保護(hù)是電力系統(tǒng)的第一道防線，有效利用繼電保護(hù)運(yùn)行數(shù)據(jù)對(duì)繼保設(shè)備狀態(tài)進(jìn)行評(píng)估，是保持電網(wǎng)安全穩(wěn)定運(yùn)行的有效途徑[1—2]。

針對(duì)繼保設(shè)備的狀態(tài)評(píng)估，國(guó)內(nèi)外做了大量研究。從繼保設(shè)備的原理、硬件、配置等方面進(jìn)行分析、建模、仿真試驗(yàn)，然后對(duì)其可靠性進(jìn)行分析，為設(shè)備運(yùn)行狀態(tài)評(píng)估提供理論依據(jù)，也為運(yùn)維人員提供了提高繼保設(shè)備可靠性的方法和策略參考[3—4]。例如文獻(xiàn)[5]詳細(xì)介紹了繼保設(shè)備可靠性研究現(xiàn)狀和發(fā)展趨勢(shì)，并指出設(shè)備可靠運(yùn)行是保證電網(wǎng)穩(wěn)定的重要因素。文獻(xiàn)[6]指出設(shè)備本身、信息系統(tǒng)和人為因素是保護(hù)系統(tǒng)風(fēng)險(xiǎn)的主要來(lái)源，并進(jìn)一步分析了繼保設(shè)備故障風(fēng)險(xiǎn)對(duì)一次系統(tǒng)的影響。文獻(xiàn)[7]提出了一種基于變權(quán)重模糊綜合評(píng)價(jià)的繼保設(shè)備狀態(tài)評(píng)價(jià)方法，其指標(biāo)更加完善，但需要足夠多的失效樣本，而繼保設(shè)備本身可靠性較高，失效樣本數(shù)量較少，因此難以滿足精度要求。對(duì)此，文獻(xiàn)[8]在三參數(shù)威布爾分布模型的基礎(chǔ)上結(jié)合灰色模型，在小樣本條件下對(duì)繼保設(shè)備進(jìn)行壽命評(píng)估，取得了精度較高的計(jì)算結(jié)果。文獻(xiàn)[9]從擴(kuò)充樣本的角度出發(fā)，使用蒙特卡洛法進(jìn)行抽樣，提高了小樣本情況下評(píng)估結(jié)果的穩(wěn)定性。文獻(xiàn)[10]考慮了繼保設(shè)備缺陷數(shù)據(jù)的隨機(jī)截尾特征，基于兩參數(shù)威布爾分布，采用極大似然估計(jì)(maximum likelihood estimation，MLE)實(shí)現(xiàn)參數(shù)估計(jì)，但未考慮區(qū)間刪失數(shù)據(jù)的影響。整體上，目前繼保設(shè)備運(yùn)行狀態(tài)評(píng)估方法很少計(jì)及刪失數(shù)據(jù)的特點(diǎn)，且對(duì)包含刪失數(shù)據(jù)的保護(hù)運(yùn)行“大數(shù)據(jù)”未充分挖掘，造成了信息損失、可靠性評(píng)估結(jié)果偏頗和分析模型假設(shè)過(guò)多的問(wèn)題。圍繞電力大數(shù)據(jù)的分析與利用，國(guó)內(nèi)外學(xué)者已進(jìn)行了大量研究[11—12]。但在大數(shù)據(jù)背景下，對(duì)繼保設(shè)備的失效數(shù)據(jù)進(jìn)行挖掘、分析，從而為判別繼保設(shè)備狀態(tài)、評(píng)估其運(yùn)行水平提供可行途徑的相關(guān)研究，還較少見(jiàn)到。研究基于大數(shù)據(jù)的繼保設(shè)備可靠性分析與運(yùn)行狀態(tài)評(píng)估方法具有重要意義。

文中結(jié)合繼保設(shè)備運(yùn)行數(shù)據(jù)特點(diǎn)，首先基于期望-最大化(expectation-maximization，EM)算法對(duì)失效模型的參數(shù)進(jìn)行估計(jì)，并計(jì)算得到設(shè)備的可靠性指標(biāo)；其次，對(duì)比各失效期內(nèi)不同估計(jì)方法得到的模型參數(shù)精度，驗(yàn)證了方法處理刪失數(shù)據(jù)的有效性；最后，通過(guò)算例驗(yàn)證了利用文中方法規(guī)劃設(shè)備檢修周期的可行性。

1 大數(shù)據(jù)背景下失效數(shù)據(jù)特點(diǎn)分析

長(zhǎng)期實(shí)踐經(jīng)驗(yàn)表明，設(shè)備失效率與累計(jì)運(yùn)行時(shí)間之間的函數(shù)關(guān)系可用浴盆曲線表示。設(shè)備的全壽命周期由早期失效期、偶然失效期與老化失效期組成[13]。對(duì)于繼保設(shè)備，出廠前會(huì)經(jīng)過(guò)充分的測(cè)試以排除隱患，現(xiàn)場(chǎng)安裝后，運(yùn)維人員也會(huì)對(duì)設(shè)備進(jìn)行嚴(yán)密的調(diào)試。因此，可認(rèn)為繼保設(shè)備投運(yùn)后其運(yùn)行狀態(tài)便處于偶然失效期。繼保設(shè)備的失效數(shù)據(jù)一般是在其投入運(yùn)行之后收集的現(xiàn)場(chǎng)數(shù)據(jù)。

受觀測(cè)條件與監(jiān)測(cè)手段的限制，現(xiàn)場(chǎng)收集到的失效數(shù)據(jù)可能包含4種類(lèi)型。

(1)精確數(shù)據(jù)。繼保設(shè)備的自檢功能日益完善，自檢周期非常短，因此自檢發(fā)現(xiàn)缺陷的時(shí)間即可認(rèn)為是缺陷發(fā)生的準(zhǔn)確時(shí)間。此外，在周期較短的巡視過(guò)程中的缺陷發(fā)現(xiàn)時(shí)間也可認(rèn)定是缺陷發(fā)生時(shí)間。這類(lèi)失效數(shù)據(jù)稱為精確數(shù)據(jù)，如圖1中編號(hào)為1、2的失效數(shù)據(jù)，線段的左端點(diǎn)表示設(shè)備的投運(yùn)時(shí)間，右端點(diǎn)表示缺陷發(fā)生時(shí)間。

圖1 原始失效數(shù)據(jù)Fig.1 Raw failure data

(2)區(qū)間刪失數(shù)據(jù)。部分隱蔽性較強(qiáng)的缺陷無(wú)法自檢查出，無(wú)法確定具體發(fā)生時(shí)間，如專(zhuān)業(yè)巡視時(shí)發(fā)現(xiàn)的“復(fù)歸按鈕脫落”“備用線圈斷線”等缺陷，只可推斷出發(fā)生在本次巡視前、最近一次巡視之后。文中稱此類(lèi)無(wú)法確定缺陷發(fā)生的確切時(shí)間，只能推斷時(shí)間區(qū)間的失效數(shù)據(jù)為區(qū)間刪失數(shù)據(jù)。如圖1中編號(hào)為3、4的失效數(shù)據(jù)，線段的左端點(diǎn)表示設(shè)備的投運(yùn)時(shí)間，尖括號(hào)之間的線段表示缺陷發(fā)生時(shí)間所在的區(qū)間。

(3)右刪失數(shù)據(jù)。在實(shí)際工作情況下收集繼保設(shè)備的失效數(shù)據(jù)，往往伴隨著各種偶然事件。在觀察結(jié)束之前，有可能會(huì)出現(xiàn)樣本設(shè)備中途退出運(yùn)行的情形。設(shè)備在退出運(yùn)行之前始終保持正常工作，只能得到設(shè)備保持正常工作的時(shí)長(zhǎng)，此類(lèi)數(shù)據(jù)稱為右刪失數(shù)據(jù)。如圖1中編號(hào)為5、6的失效數(shù)據(jù)，線段的左端點(diǎn)表示設(shè)備的投運(yùn)時(shí)間，右端點(diǎn)表示設(shè)備退出運(yùn)行的時(shí)間。

(4)隨機(jī)截尾數(shù)據(jù)。繼保設(shè)備的可靠性水平較高，因此直到觀察結(jié)束，大部分設(shè)備仍能保持正常工作。此類(lèi)直至觀察結(jié)束仍未出現(xiàn)失效的樣本數(shù)據(jù)稱為截尾數(shù)據(jù)[10]，如圖1中編號(hào)為7、8、9的失效數(shù)據(jù)，線段的左端點(diǎn)表示設(shè)備的投運(yùn)時(shí)間，右端點(diǎn)表示觀察結(jié)束的時(shí)間，即截尾時(shí)間ttr。

因各繼保設(shè)備的投運(yùn)時(shí)間不同，所以當(dāng)?shù)竭_(dá)ttr時(shí)，各繼保設(shè)備的運(yùn)行時(shí)間也不盡相同。如果把各樣本設(shè)備的投運(yùn)時(shí)刻移動(dòng)到同一計(jì)時(shí)起點(diǎn)，這些設(shè)備的失效數(shù)據(jù)將表現(xiàn)出隨機(jī)截尾的特點(diǎn)，見(jiàn)圖2。

圖2 計(jì)時(shí)起點(diǎn)對(duì)齊后的失效數(shù)據(jù)Fig.2 Failure data after starting time is aligned

針對(duì)繼保設(shè)備失效數(shù)據(jù)的上述特點(diǎn)，常用的最小二乘法等算法已無(wú)法滿足設(shè)備失效模型參數(shù)估計(jì)的需求。文中結(jié)合EM算法，利用帶有刪失特性的失效數(shù)據(jù)，估計(jì)繼保設(shè)備失效模型的參數(shù)。

2 繼保設(shè)備失效分布模型

失效分布模型能反映設(shè)備所處運(yùn)行階段及設(shè)備的運(yùn)行狀態(tài)。采用指數(shù)分布模型和威布爾分布模型，通過(guò)計(jì)算分布模型的參數(shù)，估計(jì)設(shè)備的可靠度R、故障率λ、平均無(wú)故障時(shí)間tMTBF等指標(biāo)。其中，指數(shù)分布反映偶然失效期的失效特征，常用于擬合具有恒定失效率的設(shè)備壽命分布，運(yùn)行的繼保設(shè)備多處于該失效期。威布爾分布則可通過(guò)不同的形狀和尺度參數(shù)反映早期失效期(失效率遞減)、偶然失效期和損耗失效期(失效率呈加速趨勢(shì)，如繼保設(shè)備的老化期)特征，具備較強(qiáng)的靈活性。

2.1 指數(shù)分布

指數(shù)分布模型[14]的可靠性指標(biāo)如下。

可靠度函數(shù)：

R(t)=e-λt

(1)

式中：t為設(shè)備失效時(shí)間。

故障概率密度函數(shù)：

(2)

失效率：

λ(t)=λ

(3)

平均無(wú)故障時(shí)間：

tMTBF=1/λ

(4)

2.2 威布爾分布

威布爾分布模型[15—17]可分成三參數(shù)和兩參數(shù)2種形式。兩參數(shù)模型的函數(shù)表達(dá)式如下。

可靠度函數(shù)：

R(t)=e-(t/η)k

(5)

式中：η>0為尺度參數(shù)；k>0為形狀參數(shù)。

故障概率密度函數(shù)：

(6)

失效率：

λ(t)=(k/η)(t/η)k-1

(7)

平均無(wú)故障時(shí)間：

tMTBF=ηΓ(1+1/k)

(8)

式中：Γ(·)為伽馬函數(shù)。

3 繼保設(shè)備運(yùn)行狀態(tài)評(píng)估

3.1 EM算法

EM算法的基本思想是：當(dāng)失效數(shù)據(jù)樣本中存在刪失數(shù)據(jù)時(shí)，若已知分布模型的參數(shù)，則可以結(jié)合已知樣本數(shù)據(jù)估計(jì)刪失數(shù)據(jù)的數(shù)學(xué)期望。反之，若已知?jiǎng)h失數(shù)據(jù)的值，則可利用MLE估計(jì)分布模型的參數(shù)[18]。

EM算法流程簡(jiǎn)述如下。(1)輸入包括刪失數(shù)據(jù)的樣本數(shù)據(jù)；(2)對(duì)分布模型的參數(shù)賦初值，并設(shè)定可接受的最大參數(shù)估計(jì)誤差；(3)執(zhí)行E步：將上次迭代得到的參數(shù)估計(jì)結(jié)果代入分布模型，估計(jì)刪失數(shù)據(jù)的概率分布，并計(jì)算刪失數(shù)據(jù)服從此概率分布下對(duì)數(shù)似然函數(shù)的數(shù)學(xué)期望；(4)執(zhí)行M步：對(duì)參數(shù)尋優(yōu)，使上述期望達(dá)到最大，將尋優(yōu)結(jié)果作為本次參數(shù)估計(jì)的結(jié)果并計(jì)算估計(jì)誤差；(5)若參數(shù)估計(jì)誤差小于可接受的最大誤差，計(jì)算結(jié)束，否則返回步驟(3)。

3.2 基于EM算法的指數(shù)分布(EM-exp)模型參數(shù)估計(jì)

指數(shù)分布的故障概率密度函數(shù)如式(2)所示，其對(duì)數(shù)似然函數(shù)的數(shù)學(xué)期望為：

(9)

式中：n為樣本總數(shù)；ti為第i個(gè)樣本的真實(shí)失效時(shí)間；z為發(fā)生數(shù)據(jù)刪失的樣本集；θ，θ′分別為不同分布模型中本次和上次迭代得到的參數(shù)；對(duì)指數(shù)分布，λ′，λ分別為上次和本次估計(jì)的故障率參數(shù)；Q(·)為似然函數(shù)的數(shù)學(xué)期望；E(·)為ti在服從上次估計(jì)的故障率參數(shù)下對(duì)數(shù)似然函數(shù)的數(shù)學(xué)期望。

當(dāng)Q(θ|θ′)取得極大值時(shí)，有：

(10)

根據(jù)EM算法的收斂性質(zhì)，當(dāng)?shù)螖?shù)ω→+∞時(shí)，λ收斂于定值，即：

λ=λ′ω→+∞

(11)

迭代次數(shù)足夠多時(shí)，將式(10)中的各數(shù)學(xué)期望項(xiàng)展開(kāi)，并將式(11)代入式(10)，得：

(12)

采用數(shù)值方法求解上述方程即可得到指數(shù)分布模型的參數(shù)估計(jì)值。

3.3 基于EM算法的威布爾分布(EM-wbl)模型參數(shù)估計(jì)

威布爾分布的故障概率密度函數(shù)如式(6)所示。為簡(jiǎn)化計(jì)算，令μ=1/ηk，則式(6)可改寫(xiě)為：

(13)

威布爾分布模型對(duì)數(shù)似然函數(shù)的數(shù)學(xué)期望為：

(14)

式中：k′,μ′分別為上次迭代所得威布爾分布的形狀參數(shù)和尺度參數(shù)估計(jì)結(jié)果；k,μ為本次迭代估計(jì)的參數(shù)。

當(dāng)Q(θ|θ′)取得極大值時(shí)，有：

(15)

(16)

當(dāng)?shù)螖?shù)ω→+∞時(shí)，參數(shù)k,μ將收斂于定值，即：

(17)

迭代次數(shù)足夠多時(shí)，將式(15)、式(16)中的各數(shù)學(xué)期望項(xiàng)展開(kāi)，并將式(17)代入式(15)、式(16)，得：

(18)

(19)

求解上述隱函數(shù)方程組，即可求得威布爾分布參數(shù)k,μ的數(shù)值解。

傳統(tǒng)MLE在處理刪失數(shù)據(jù)時(shí)只計(jì)及了隨機(jī)截尾數(shù)據(jù)，不能處理區(qū)間刪失數(shù)據(jù)，會(huì)影響繼保設(shè)備運(yùn)行狀態(tài)評(píng)估的準(zhǔn)確性。而基于EM算法的失效分布模型，在考慮區(qū)間刪失數(shù)據(jù)的基礎(chǔ)上，可通過(guò)對(duì)刪失數(shù)據(jù)的還原處理，得到理想的參數(shù)估計(jì)值并據(jù)此計(jì)算設(shè)備的可靠度、故障率、平均無(wú)故障時(shí)間等指標(biāo)，從而實(shí)現(xiàn)繼保設(shè)備運(yùn)行狀態(tài)的評(píng)估。

4 算例分析

為檢驗(yàn)?zāi)Ｐ偷恼_性，利用計(jì)算機(jī)生成帶有刪失特性的失效數(shù)據(jù)進(jìn)行分析。首先，使用蒙特卡羅法生成服從指數(shù)分布和威布爾分布的精確失效數(shù)據(jù)[19—21]，并以此為基礎(chǔ)生成帶有刪失特性的樣本集。為模擬區(qū)間刪失數(shù)據(jù)，隨機(jī)選取部分樣本進(jìn)行刪失處理。在失效數(shù)據(jù)的左右兩側(cè)分別擴(kuò)充出刪失區(qū)間，兩側(cè)的刪失區(qū)間長(zhǎng)度服從參數(shù)為ζ2,σ2的正態(tài)分布，且相互獨(dú)立。

為模擬右刪失數(shù)據(jù)，在非區(qū)間刪失樣本中隨機(jī)選取部分樣本數(shù)據(jù)，令其刪失時(shí)間tc服從(t0,ttr)上的均勻分布，若刪失發(fā)生在失效之后，則重新分配小于失效時(shí)間的tc。其中，t0為投運(yùn)時(shí)間。

為模擬失效數(shù)據(jù)的隨機(jī)截尾特征，令計(jì)時(shí)起點(diǎn)為t0=0，各樣本設(shè)備的投運(yùn)時(shí)間服從參數(shù)為ζ1,σ1的正態(tài)分布。若設(shè)備的投運(yùn)時(shí)間在計(jì)時(shí)起點(diǎn)之前，則重新分配大于t0的投運(yùn)時(shí)間。失效時(shí)間大于ttr的樣本數(shù)據(jù)即為截尾數(shù)據(jù)。

4.1 偶然失效區(qū)

繼保設(shè)備工作于偶然失效期時(shí)，其失效率基本保持不變。使用計(jì)算機(jī)生成服從指數(shù)分布的失效數(shù)據(jù)，樣本容量為100。生成失效數(shù)據(jù)時(shí)所用參數(shù)如表1所示。

表1 參數(shù)表(偶然失效期)Table 1 Parameters table (random failure)

分別使用EM-exp模型、EM-wbl模型和傳統(tǒng)基于MLE的指數(shù)分布模型(MLE-exp)進(jìn)行參數(shù)估計(jì)。由于MLE-exp方法不能處理區(qū)間刪失數(shù)據(jù)，所以在使用此類(lèi)數(shù)據(jù)時(shí)，選取缺陷發(fā)現(xiàn)時(shí)間，即刪失區(qū)間的右端點(diǎn)作為設(shè)備的失效時(shí)間。參數(shù)估計(jì)結(jié)果如表2所示，表中p1為使用計(jì)算機(jī)生成數(shù)據(jù)時(shí)，精確數(shù)據(jù)在所有非截尾數(shù)據(jù)中所占比例；p2為區(qū)間刪失數(shù)據(jù)所占比例；p3為右刪失數(shù)據(jù)所占比例；tcal為計(jì)算耗時(shí)。

表2 偶然失效區(qū)參數(shù)擬合結(jié)果Table 2 Parameter titting results of random failures

由表2可得：

(1)刪失數(shù)據(jù)所占比例較小時(shí)，與EM-exp模型相比，MLE-exp模型的參數(shù)估計(jì)耗時(shí)最少，且2種模型得到的λ的估計(jì)值與真值之間的差距較小。因此，對(duì)于失效數(shù)據(jù)為精確數(shù)據(jù)的繼保設(shè)備，宜采用MLE-exp模型進(jìn)行參數(shù)估計(jì)。

(2)隨著刪失數(shù)據(jù)增多，3種模型得到的參數(shù)估計(jì)值與真值之間的差距不斷增大，此時(shí)EM-exp模型準(zhǔn)確性最好。因此，對(duì)于失效數(shù)據(jù)中含刪失數(shù)據(jù)的繼保設(shè)備，宜采用EM-exp模型進(jìn)行參數(shù)估計(jì)。

(3)對(duì)處于偶然失效期的設(shè)備，EM-wbl模型的參數(shù)估計(jì)結(jié)果精度較差。因此，對(duì)于具有恒定失效率或定期維護(hù)的繼保設(shè)備，宜采用EM-exp模型進(jìn)行參數(shù)估計(jì)。

4.2 損耗失效區(qū)

使用計(jì)算機(jī)生成服從威布爾分布的失效數(shù)據(jù)，樣本容量為100。生成數(shù)據(jù)時(shí)所用參數(shù)，如表3所示。

表3 參數(shù)表(損耗失效區(qū))Table 3 Parameters table (wearout failure)

指數(shù)分布模型不適用于失效率隨時(shí)間不斷變化的情形，所以本節(jié)只討論威布爾分布模型的參數(shù)估計(jì)精度。分別使用EM-wbl模型和傳統(tǒng)的基于MLE的威布爾分布模型(MLE-wbl)進(jìn)行參數(shù)估計(jì)，估計(jì)結(jié)果如表4所示。

表4 損耗失效區(qū)參數(shù)擬合結(jié)果Table 4 Parameter fitting results of wearout failures

由表4可知，當(dāng)樣本數(shù)據(jù)全為精確數(shù)據(jù)時(shí)，2種方法計(jì)算得到的參數(shù)k,η以及tMTBF的值均相同。隨著刪失數(shù)據(jù)所占比例不斷上升，2種方法的估計(jì)誤差均有所增加，但EM-wbl模型的估計(jì)結(jié)果更接近準(zhǔn)確值，準(zhǔn)確度更高。

5 實(shí)例分析

5.1 數(shù)據(jù)預(yù)處理

選取某市同一型號(hào)的繼保設(shè)備50臺(tái)，記設(shè)備投運(yùn)時(shí)間為t0，停止觀察的時(shí)間為截尾時(shí)間ttr[22—23]。對(duì)于直到觀察結(jié)束仍未發(fā)生缺陷的設(shè)備，其累計(jì)運(yùn)行時(shí)長(zhǎng)tacc為：

tacc=ttr-t0

(20)

對(duì)于沒(méi)有發(fā)生故障，但在觀察期間退出運(yùn)行的設(shè)備，記其停運(yùn)時(shí)間為tc，其累計(jì)運(yùn)行時(shí)長(zhǎng)為：

tacc=tc-t0

(21)

對(duì)于發(fā)生缺陷的設(shè)備，若缺陷發(fā)現(xiàn)方式為自檢，則發(fā)現(xiàn)缺陷的時(shí)間即為缺陷發(fā)生時(shí)間tf，設(shè)備的失效時(shí)間記為：

t=tf-t0

(22)

若發(fā)現(xiàn)缺陷的方式為巡視過(guò)程，則設(shè)備的失效時(shí)間區(qū)間(t-,t+)的計(jì)算方法為：

(23)

式中：T為設(shè)備歷次巡視時(shí)間，取其中距缺陷發(fā)現(xiàn)時(shí)間最近的一次巡視時(shí)間作為設(shè)備失效區(qū)間的左端點(diǎn)，缺陷發(fā)現(xiàn)時(shí)間為失效區(qū)間的右端點(diǎn)。

整理后的部分失效數(shù)據(jù)如表5所示，觀察期間，50臺(tái)設(shè)備中共有9臺(tái)發(fā)生失效，6臺(tái)設(shè)備的失效時(shí)間為精確數(shù)據(jù)，3臺(tái)設(shè)備的失效時(shí)間為區(qū)間刪失數(shù)據(jù)，無(wú)設(shè)備在正常運(yùn)行情況下中途退出運(yùn)行。

表5 繼保設(shè)備失效數(shù)據(jù)Table 5 Failure data of relay protection equipment

5.2 計(jì)算分析

將經(jīng)預(yù)處理的失效數(shù)據(jù)代入EM-wbl模型中，得到的參數(shù)估計(jì)值為：

(24)

將參數(shù)值代入失效模型中，可得出此型號(hào)繼保設(shè)備的可靠性指標(biāo)如下。

可靠度：

R(t)=e-(t/60 641)1.717 2

(25)

故障概率密度：

f(t)=1.050 6×10-8t0.717 2e-(t/60 641)1.717 2

(26)

失效率：

λ(t)=1.050 6×10-8t0.717 2

(27)

平均無(wú)故障時(shí)間：

tMTBF=54 071

(28)

可靠度、故障概率密度、失效率隨時(shí)間變化的曲線如圖3—圖5所示。

圖3 可靠度曲線Fig.3 Reliability function curve

圖4 故障概率密度曲線Fig.4 Fault probability density curve

圖5 失效率曲線Fig.5 Failure rate curve

如圖3所示，當(dāng)R(t)=0.98時(shí)，得到繼保設(shè)備運(yùn)行時(shí)間為6 250 h，實(shí)際運(yùn)行情況中設(shè)備首次失效時(shí)間為7 283 h，即通過(guò)EM算法得到的參數(shù)估計(jì)值可以很好的擬合設(shè)備可靠度曲線。形狀參數(shù)k>1，由圖4、圖5可知，故障概率密度及失效率隨運(yùn)行時(shí)間的增加逐漸增大，反映了繼保設(shè)備處于損耗期的特征。平均無(wú)故障時(shí)間tMTBF可用于檢修周期制定，通常情況下?tīng)顟B(tài)檢修周期Tm=0.1tMTBF，即建議狀態(tài)檢修周期為5 407 h。

綜上可知，當(dāng)存在刪失數(shù)據(jù)時(shí)，采用EM-wbl得到的參數(shù)估計(jì)值可更精確地?cái)M合設(shè)備實(shí)際的失效情況，再根據(jù)繼保設(shè)備可靠度、失效率、故障概率密度等指標(biāo)隨時(shí)間變化的曲線預(yù)判其壽命分布，達(dá)到評(píng)價(jià)繼保設(shè)備運(yùn)行水平，預(yù)測(cè)設(shè)備運(yùn)行年限的目的。

此外，上述分析的失效數(shù)據(jù)來(lái)自同類(lèi)型設(shè)備，因此得出的結(jié)論適用于該類(lèi)型設(shè)備。但所提模型對(duì)不同型號(hào)設(shè)備具有普適性。

6 結(jié)語(yǔ)

首先分析了繼保設(shè)備失效數(shù)據(jù)的特點(diǎn)，總結(jié)了現(xiàn)場(chǎng)收集到的失效數(shù)據(jù)類(lèi)型：精確數(shù)據(jù)、區(qū)間刪失數(shù)據(jù)、右刪失數(shù)據(jù)和隨機(jī)截尾數(shù)據(jù)。然后在失效數(shù)據(jù)樣本包含刪失數(shù)據(jù)的情況下，對(duì)比分析了偶然失效期和老化失效期內(nèi)不同估計(jì)方法得到的模型參數(shù)精度。證明在偶然失效期內(nèi)，EM-exp模型得到的參數(shù)估計(jì)精度更高；在老化失效期內(nèi)，EM-wbl模型能以更高的精度擬合繼保設(shè)備的失效特性曲線。

最后通過(guò)對(duì)某型號(hào)設(shè)備數(shù)據(jù)處理，結(jié)合失效分布模型，分析了該繼保設(shè)備可靠度、失效率、故障概率密度等可靠性指標(biāo)的時(shí)變特征。利用設(shè)備可靠度曲線，可提前判斷設(shè)備可靠度，降低到閥值所需的時(shí)間，并及時(shí)發(fā)出告警提醒工作人員加強(qiáng)巡視。利用失效率曲線，可得設(shè)備的故障風(fēng)險(xiǎn)隨時(shí)間變化的規(guī)律，用于安排檢修計(jì)劃。利用故障概率密度曲線，可計(jì)算設(shè)備的平均無(wú)故障時(shí)間，為合理安排服役年限提供參考。下一步將針對(duì)不同類(lèi)型刪失數(shù)據(jù)，對(duì)比研究不同算法的參數(shù)估計(jì)精度。