范思遐,吳 斌,李友釗
(1.上海電機(jī)學(xué)院 商學(xué)院,上海 201306;2.上海電機(jī)學(xué)院 電氣學(xué)院,上海 201306)
根據(jù)國(guó)家能源局在《風(fēng)電發(fā)展“十三五”規(guī)劃》做出的指示,未來(lái)風(fēng)電行業(yè)不僅有著總裝機(jī)量的要求,其機(jī)組的運(yùn)維優(yōu)化也被提上議程。作為海上風(fēng)電度電成本的主要影響因素之一,降低運(yùn)維成本的重要性不言而喻。盡管風(fēng)電機(jī)組失效故障與停機(jī)時(shí)間獨(dú)立分布,但海上風(fēng)機(jī)50%的停機(jī)維修活動(dòng)是主要由于齒輪箱故障失效導(dǎo)致[1]。因此對(duì)齒輪箱進(jìn)行健康管理是風(fēng)電場(chǎng)運(yùn)維活動(dòng)的一個(gè)重要環(huán)節(jié)。
目前,對(duì)齒輪箱維護(hù)多為無(wú)計(jì)劃或者低效率計(jì)劃的維護(hù),一般是待齒輪箱的狀態(tài)衰退至一定程度后才進(jìn)行維修。這類維護(hù)策略雖然在傳統(tǒng)的維修上作出了適應(yīng)性改進(jìn),便于開展,然而沒有彈性的維護(hù)計(jì)劃易造成資源浪費(fèi),增加度電成本。因此,尋找“過(guò)度維護(hù)”與“維護(hù)不足”之間平衡點(diǎn)顯得尤為重要。國(guó)內(nèi)外也發(fā)展了很多針對(duì)齒輪箱等設(shè)備的預(yù)防性維護(hù)策略研究[2-6];有學(xué)者提出隨役齡增加,維護(hù)期逐漸變短的維護(hù)策略[7-8];由于海上風(fēng)機(jī)維護(hù)會(huì)隨階段不同具有不同的維護(hù)策略,應(yīng)論證一種符合階段特點(diǎn)的維護(hù)策略。
利用搭建的風(fēng)機(jī)模擬平臺(tái)進(jìn)行退化狀態(tài)評(píng)估實(shí)驗(yàn),劃分出四個(gè)退化區(qū)間,以此建立了分階段的機(jī)會(huì)維護(hù)模型,并給出了各部件的具體維護(hù)動(dòng)作及單位時(shí)間維護(hù)成本,驗(yàn)證了分階段機(jī)會(huì)維護(hù)的合理性。
按照維修經(jīng)驗(yàn),風(fēng)機(jī)齒輪箱的維護(hù)策略大致可以分為三種:反應(yīng)性維護(hù),預(yù)防性維護(hù),預(yù)測(cè)性維護(hù)。
在實(shí)際風(fēng)電場(chǎng)運(yùn)維活動(dòng)中,針對(duì)不同部件安排各自檢修周期具有不可操作性,作為多部件構(gòu)成的復(fù)雜系統(tǒng),考慮到部件之間存在著維護(hù)經(jīng)濟(jì)相關(guān)性,故將風(fēng)機(jī)齒輪箱的幾個(gè)關(guān)鍵部件(A,B,C,D…)作為一個(gè)維護(hù)組,假設(shè)其間發(fā)生失效的概率是互不相關(guān)的,則機(jī)會(huì)維護(hù)的表述與假設(shè)應(yīng)如下:
(1)齒輪箱各部件的狀態(tài)是可以實(shí)時(shí)監(jiān)測(cè)的,當(dāng)不可回復(fù)故障發(fā)生時(shí)立即停機(jī)維修(最小修復(fù)),否則在預(yù)定維護(hù)時(shí)間點(diǎn)時(shí)進(jìn)行維護(hù)(i表示第i次維護(hù)活動(dòng)),部件的失效率函數(shù)符合Weibull比例失效模型,且部件間不存在拆卸順序相關(guān)性及隨機(jī)相關(guān)性;
(2)部件A作為主要維護(hù)部件,當(dāng)其可靠度衰退到維護(hù)閾值Rp時(shí),對(duì)齒輪箱開展維護(hù)活動(dòng),其它部件獲得了同期進(jìn)行維護(hù)的機(jī)會(huì):部件可靠度RA(i)處在機(jī)會(huì)維護(hù)上限Rom內(nèi),則對(duì)其進(jìn)行維護(hù);
(3)在維護(hù)時(shí)間點(diǎn),部件A的維護(hù)動(dòng)作包括不完全維護(hù)(flag=1)與完全維護(hù)(flag=2),為決定采用何種維護(hù)動(dòng)作,引入維護(hù)回報(bào)率的概念:當(dāng)前維護(hù)動(dòng)作的投入成本C與機(jī)組可靠度提升ΔR的比值ROIA(i)來(lái)判斷;維護(hù)組中其他部件B,C,D…的維護(hù)動(dòng)作則包括三種,分別為不完全維護(hù)、完全維護(hù)、不維護(hù)(flag=3),判斷條件與部件A類似[9];
(4)限制條件,即維護(hù)次數(shù)上限Nm,以及置換閾值Rrep。
在傳統(tǒng)的策略考慮中,并沒有對(duì)齒輪箱退化狀態(tài)進(jìn)行分段,更多是依靠經(jīng)驗(yàn)常識(shí)對(duì)設(shè)備的健康狀況進(jìn)行劃分,缺乏科學(xué)性。
經(jīng)驗(yàn)數(shù)據(jù)表明,風(fēng)電機(jī)組齒輪箱在初始安裝階段故障率最高,而隨著役齡增加,故障率急速下降,直至度過(guò)試運(yùn)營(yíng)期,其失效率h(t)比較穩(wěn)定的處于一個(gè)較低的水平上,因此可用失效率的積分,即可靠度R(t)作為齒輪箱狀態(tài)指標(biāo);而到了后期長(zhǎng)時(shí)間的運(yùn)行導(dǎo)致齒輪箱各部件發(fā)生磨損、疲勞、老化等耗損現(xiàn)象嚴(yán)重,部件的狀態(tài)開始急速衰退,使得失效率又開始上升[10]。故在齒輪箱試運(yùn)行期以后的時(shí)間段,開展多個(gè)部件分階段機(jī)會(huì)性維護(hù)動(dòng)作,其退化趨勢(shì)在每個(gè)退化階段內(nèi)有所差異。具體的維護(hù)示意圖,如圖1所示。
圖1 風(fēng)電機(jī)組齒輪箱分階段維護(hù)示意圖Fig.1 Phased Maintenance Diagram of Wind Turbine Gearbox
退化矩陣描述了齒輪箱在整個(gè)壽命區(qū)間內(nèi)的退化特征,為有效將這些特征向量劃分到不同的分類器中,采用K-means聚類分析法將其分別劃分出(3~7)個(gè)聚類數(shù)目,然后計(jì)算各數(shù)目下的Scat-Sep(最小類內(nèi)緊密度—最大類間分離度)指標(biāo),對(duì)應(yīng)該指標(biāo)最小的聚類劃分即為最佳退化階段數(shù)目。
假設(shè)分類器用 C[k]={C1,C2,…,Ck}表示,其中 k 為總聚類數(shù)目,則類內(nèi)緊密度指標(biāo)Scat(C[k])用類內(nèi)各數(shù)據(jù)之間歐氏距離的平方和來(lái)描述,類間分離度指標(biāo)Sep(C[k])用類與類之間的距離來(lái)描述,即類別間點(diǎn)對(duì)的平均距離:
式中:k—第k個(gè)分類器;
xm,xn—類Ci內(nèi)的任意兩個(gè)數(shù)據(jù)點(diǎn);
y—類Cj內(nèi)的數(shù)據(jù)點(diǎn);
考慮到兩者的取值范圍并無(wú)差異,故組合參數(shù)α取為1,組合平衡函數(shù)為:
式中:M—兩指標(biāo)平衡系數(shù),其值等于最大聚類(每個(gè)類別只含一個(gè)樣本點(diǎn))時(shí)的類間分離度指標(biāo),或者最小聚類(唯一聚類)時(shí)的類內(nèi)緊密度指標(biāo)。
模型建立的關(guān)鍵是確定各部件的維護(hù)動(dòng)作,根據(jù)動(dòng)作不同將部件分為兩類,一類是主維護(hù)部件A,另一類是機(jī)會(huì)維護(hù)部件B,C,D,...,引入維護(hù)動(dòng)作標(biāo)記符號(hào)如下:
3.2.1 單位時(shí)間維護(hù)成本
以風(fēng)機(jī)齒輪箱階段性預(yù)防維護(hù)策略所跨時(shí)間段內(nèi)單位時(shí)間維護(hù)費(fèi)用Cut最小為決策依據(jù),介紹計(jì)算Cut的方法。
(1)維護(hù)停機(jī)損失Cs:維護(hù)活動(dòng)造成的發(fā)電量損失可根據(jù)同期平均風(fēng)機(jī)功率Pave、電價(jià)p以及停機(jī)時(shí)間ts決定,將等待時(shí)間tw也加入到停機(jī)時(shí)間中[6]。從而有下列式:
式中:γ—風(fēng)電機(jī)組同期發(fā)電容量因子;
τi—第i次預(yù)防維護(hù)動(dòng)作所用時(shí)間。
(2)維護(hù)成本:齒輪箱的維護(hù)成本包括三個(gè)部分,部件置換成本Crep、不完全維護(hù)成本Cinc以及最小修復(fù)成本Cminr。理論上來(lái)說(shuō),預(yù)防維護(hù)的投入資源越多其維護(hù)效果就越好,維護(hù)后部件的可靠度就越有保障,因此不完全維護(hù)成本應(yīng)該與部件可靠度函數(shù)的該變量 ΔRi有關(guān),則:Cinc(i)=K(ΔRi)μ
總維護(hù)成本Cme就可以用下式表示:
(3)維護(hù)耗材成本與運(yùn)輸?shù)跹b成本Ctrn。維護(hù)(包括最小修復(fù)與不完全修復(fù))與置換所用到的運(yùn)輸成本,及運(yùn)維人員工資一并算入運(yùn)輸?shù)跹b成本中。模型的決策依據(jù)為單位維護(hù)時(shí)間內(nèi)成本最小,故有:
3.2.2 具體維護(hù)動(dòng)作的判斷
維護(hù)回報(bào)率可用下式計(jì)算:
式中:CA,B,C,…(i)—部件 A,B,C,…第 i次維護(hù)動(dòng)作的維護(hù)成本;ΔRall(i)—部件第i次維護(hù)動(dòng)作執(zhí)行后風(fēng)電機(jī)組的可靠度改善情況。
計(jì)算flag置不同值時(shí)的維護(hù)成本,從而求得對(duì)應(yīng)的維護(hù)回報(bào)率ROI,根據(jù)下式判斷維護(hù)動(dòng)作:
為探究風(fēng)電機(jī)組齒輪箱退化過(guò)程,搭建實(shí)驗(yàn)?zāi)M平臺(tái),如圖2所示。其中齒輪箱的傳動(dòng)比為1:1.36;0.5kW 的磁粉制動(dòng)器作為負(fù)載;變頻調(diào)速三相異步電機(jī)作為動(dòng)力源;變頻器分辨率達(dá)到0.01Hz;齒輪箱振動(dòng)信號(hào)的采集利用MPS-140401信號(hào)采集卡。三相電機(jī)頻率設(shè)定f=20Hz,即轉(zhuǎn)動(dòng)一周期需0.1s,MPS采樣率設(shè)定為16000Hz,每半小時(shí)采集一次數(shù)據(jù),每次持續(xù)采樣時(shí)間為30s,共采集了900組有效數(shù)據(jù),每組數(shù)據(jù)切取三個(gè)振動(dòng)樣本。對(duì)這2700組信號(hào)進(jìn)行小波包分解后,提取能量譜作為齒輪箱的退化矩陣,然后用Scat-Sep評(píng)價(jià)指標(biāo)對(duì)各聚類數(shù)目下的特征向量進(jìn)行分析計(jì)算,其結(jié)果如表1所示??梢钥闯鳊X輪箱的退化過(guò)程為四個(gè)階段時(shí),平衡函數(shù)取到了最小值 0.5311。
圖2 實(shí)驗(yàn)臺(tái)設(shè)置Fig.2 Experimental Platform
表1 不同聚類數(shù)目下各評(píng)價(jià)指標(biāo)Tab.1 Evaluation Indexes Under Different Cluster Numbers
4.2.1 模型參數(shù)設(shè)置
以某兆瓦級(jí)風(fēng)電機(jī)組為研究對(duì)象,將傳動(dòng)系統(tǒng)的主軸、主軸軸承、齒輪副、從動(dòng)軸軸承四個(gè)部件納為一個(gè)機(jī)會(huì)維護(hù)組,其初始參數(shù),如表2所示。分四個(gè)階段進(jìn)行維護(hù)。完全維護(hù)成本相關(guān)系數(shù)K直接取Crep,部件的運(yùn)輸?shù)跹b費(fèi)用Ctrn取15萬(wàn)元與25萬(wàn)元(維護(hù)動(dòng)作與置換動(dòng)作),多部件吊裝運(yùn)輸費(fèi)用增加系數(shù)k取0.13,預(yù)防維護(hù)平均等待時(shí)間取72h,平均預(yù)防性維護(hù)時(shí)間15h,平均完全維護(hù)時(shí)間40h,停機(jī)損失發(fā)電量電價(jià)0.9元/(kW·h),同期發(fā)電量容容量因子 γ 取 0.7、0.6、0.5、0.4 四個(gè)等級(jí),狀態(tài)回復(fù)因子 λ 取 0.75、0.7、0.6與0.5四個(gè)等級(jí),役齡回退因子a=i/(3i+18),可靠度遞減因子b=(18i+1)/(19i+1),不完全維護(hù)成本相關(guān)因子μ取0.8、0.75、0.7、0.7??紤]到維護(hù)閾值變化使得維護(hù)間隔不確定,故每個(gè)維護(hù)階段留有一定的時(shí)間裕量 ΔT=50d,取各階段初始維護(hù)閾值為 Rp=[0.950.930.90 0.85],時(shí)間跨度為 2179d、1002d、1166d、1128d。
表2 機(jī)會(huì)維護(hù)組各部件初始參數(shù)Tab.2 Initial Parameters for Each Part of The Opportunity Maintenance Group
4.2.2 階段性機(jī)會(huì)維護(hù)時(shí)間間隔及動(dòng)作判斷
部件A的維護(hù)時(shí)間間隔以及所處維護(hù)階段,如表3所示。從表中可以看出,在維護(hù)后期,階段性的動(dòng)態(tài)維護(hù)閾值很好的解決了維護(hù)間隔短這個(gè)問題,該結(jié)果與實(shí)際風(fēng)場(chǎng)每年(1~2)次的維護(hù)計(jì)劃更為符合。
表3 部件A階段性預(yù)防維護(hù)時(shí)間間隔Tab.3 Part A Phased Preventive Maintenance Time Interval
表4 各部件分階段預(yù)防性維護(hù)動(dòng)作Tab.4 Phased Preventive Maintenance of Each Component
階段性機(jī)會(huì)維護(hù)策略下的各部件預(yù)防性維護(hù)動(dòng)作,如表4所示。對(duì)應(yīng)的單位時(shí)間機(jī)會(huì)維護(hù)成本為(1.0504e+03)元/天,較不分階段的機(jī)會(huì)維護(hù)策略的單位時(shí)間維護(hù)成本(1.5110e+03)元/天有了很大的降低。
為了更直觀的比較齒輪箱階段性機(jī)會(huì)維護(hù)策略與傳統(tǒng)機(jī)會(huì)維護(hù)策略兩者之間維護(hù)動(dòng)作及經(jīng)濟(jì)性的差異,選取兩種維護(hù)策略第6次預(yù)防性維護(hù)為例,比較兩者之間的差異以便說(shuō)明機(jī)會(huì)維護(hù)動(dòng)作的決策過(guò)程,判斷結(jié)果,如表5所示。
表5 兩種維護(hù)策略第6次維護(hù)動(dòng)作對(duì)比Tab.5 The Sixth Maintenance Actions Comparison of Two Maintenance Strategies
表中:非階段性機(jī)會(huì)維護(hù)策略的維護(hù)閾值設(shè)為0.93。
在不考慮階段性機(jī)會(huì)維護(hù)策略中,部件D在第6次預(yù)防維護(hù)開始前可靠度為0.9531,故沒有達(dá)到進(jìn)行機(jī)會(huì)維護(hù)的條件,因此不進(jìn)行維護(hù)動(dòng)作;而在階段性的維護(hù)策略中,部件D的可靠度為0.9533,雖然與非階段性維護(hù)策略中的部件可靠度近似相等,但是其可靠度RD=0.9533卻處在進(jìn)行機(jī)會(huì)維護(hù)的范圍內(nèi)(Rom=0.965),因此進(jìn)行機(jī)會(huì)維護(hù)動(dòng)作;同時(shí),兩個(gè)不同的維護(hù)策略中,前者的故障最小修復(fù)成本較前者更高;從各部件的維護(hù)回報(bào)率也可以看出,分階段的預(yù)防性機(jī)會(huì)維護(hù)策略的ROI更小,這意味著投入相同的維護(hù)資源機(jī)組可以更好的改善可靠度,部件的到了充分的維護(hù),特別是到了維護(hù)后期,很好的平衡了傳統(tǒng)維護(hù)策略的“維護(hù)不足”與固定維護(hù)閾值機(jī)會(huì)維護(hù)策略后期的“過(guò)度維護(hù)”問題。
針對(duì)海上風(fēng)電機(jī)組齒輪箱的實(shí)際衰退特性,提出一種基于動(dòng)態(tài)維護(hù)閾值的機(jī)會(huì)維護(hù)策略,表現(xiàn)為整個(gè)維護(hù)策略是按階段進(jìn)行展開的,動(dòng)態(tài)維護(hù)閾值和固定維護(hù)閾值的差異就在于前者更好的適應(yīng)了齒輪箱隨役齡增加可靠度加速下降的特性,改善了部件的維護(hù)回報(bào)率。