融合壽命數(shù)據(jù)和退化數(shù)據(jù)的防噴閥剩余壽命預測

陳秀榮，李 娟，于加舉

(1.青島農(nóng)業(yè)大學 理學與信息科學學院,山東 青島266109；2.青島農(nóng)業(yè)大學 機電工程學院，山東 青島 266109)

融合壽命數(shù)據(jù)和退化數(shù)據(jù)的防噴閥剩余壽命預測

陳秀榮1，李 娟2，于加舉1

(1.青島農(nóng)業(yè)大學 理學與信息科學學院,山東 青島266109；2.青島農(nóng)業(yè)大學 機電工程學院，山東 青島 266109)

為了解決目前利用Wiener過程對產(chǎn)品進行壽命預測時，由于考慮個體差異至少需要對模型的三個參數(shù)進行估計從而導致計算量較大的問題，提出了利用個體方差和總體方差的相關關系以減少待估計參數(shù)個數(shù)的壽命預測方法。利用歷史壽命數(shù)據(jù)和實時退化數(shù)據(jù)，采用Bayes估計和EM(expectation maximization)算法得到性能退化信息的參數(shù)值，從而得出防噴閥的剩余壽命的概率密度函數(shù)及相關分布。仿真結(jié)果表明所提方法不僅減少了待估計參數(shù)的個數(shù)，使計算過程更簡單有效，而且有效提高了剩余壽命的預測精度。

防噴閥；剩余壽命預測；Wiener 過程；可靠性；EM算法

海洋石油鉆井平臺是勘探開發(fā)海上油氣的主要途徑，其運行狀況直接影響到油氣產(chǎn)量、人員與環(huán)境的安全[1-2]。其中,防噴閥是防止海洋鉆井平臺漏油的最后一道屏障，是海洋鉆井平臺的關鍵子設備之一，對于油氣的安全生產(chǎn)、控制生態(tài)環(huán)境的污染等具有重大意義[3]。因此準確評估防噴閥的健康狀態(tài)和可靠性具有重要的實際意義，壽命預測[4]技術(shù)是確保系統(tǒng)安全可靠運行的重要手段。

現(xiàn)有的壽命預測技術(shù)可分為[5]：基于機理分析的方法、基于數(shù)據(jù)驅(qū)動的方法、機理與數(shù)據(jù)驅(qū)動融合的方法。其中,基于數(shù)據(jù)驅(qū)動的壽命預測方法受到廣泛關注，如神經(jīng)網(wǎng)絡[6]、支持向量機[7]等。在基于數(shù)據(jù)驅(qū)動的壽命預測方法中，Wiener過程是性能退化建模中的常用模型[8]。一些學者在假設元件是線性Wiener退化過程的前提下，提出了一些行之有效的方法，如文獻[9]提出了一種基于退化建模的備件需求預測方法，得到了備件需求量的預測分布； 文獻[10]根據(jù)歷史退化數(shù)據(jù)利用極大似然法估計參數(shù)，得到了陀螺儀的剩余壽命分布。對于非線性的退化過程大多是將其線性化，再用Wiener過程進行退化建模，如文獻[11]提出一種時間尺度轉(zhuǎn)換方法，并研究基于時間尺度轉(zhuǎn)換的Wiener過程在非線性退化建模中的應用。文獻[12-13]在考慮測量誤差的情況下，采用Wiener過程對性能退化過程進行了建模，得到相應產(chǎn)品的壽命分布。

文獻[9-13]等研究均沒有考慮同類設備個體之間的差異。為此，一些學者研究了同類設備個體之間存在差異情況下的壽命預測問題，如文獻[14-15]將維納過程中的漂移系數(shù)看成隨機變量，并得出了相應產(chǎn)品的壽命分布函數(shù)；文獻[16]將維納過程中的漂移系數(shù)和擴散系數(shù)都看成隨機變量，研究了橋橫梁的壽命預測。文獻[14-16]等主要利用元件的退化數(shù)據(jù)進行壽命預測，沒有利用已有的同類設備或元件的壽命數(shù)據(jù)。為此，文獻[17]融合了同類產(chǎn)品的歷史壽命數(shù)據(jù)和該產(chǎn)品的性能退化數(shù)據(jù)得到參數(shù)的驗后分布，并利用Bayes方法估計得到該產(chǎn)品的剩余壽命；文獻[18]提出了一種融合壽命數(shù)據(jù)退化數(shù)據(jù)的非線性剩余壽命估計方法，實現(xiàn)了對陀螺儀的壽命預測。

在已有的壽命預測方法研究中，均將退化過程中漂移系數(shù)的方差和總體方差視為兩個獨立的參數(shù)。但對于某個單體而言，其個體的退化規(guī)律和同類產(chǎn)品的退化規(guī)律方差具有一定的相關性。為此，本研究以海洋平臺防噴閥作為研究對象，研究個體方差和同類產(chǎn)品總體方差之間存在線性相關時的剩余壽命預測方法,采用Bayes方法融合產(chǎn)品的歷史壽命數(shù)據(jù)和性能退化數(shù)據(jù)，從而達到減少待估計參數(shù)的數(shù)量、簡化計算和提高預測精度的目的。

1 隨機的Wiener過程建模

假設防噴閥的性能退化過程是Wiener過程，則其性能退化過程可以用下式描述[17]:

X(t)=X(0)+v t+σB(t)。

(1)

式中：X(t)是元件在時刻t的性能退化量，X(0)為初始的退化量，為便于分析，令X(0)=0，v表示漂移系數(shù)，σ表示已知的擴散系數(shù)，B(t)是標準的布朗運動?；陔S機過程首達時間的概念，防噴閥在t時刻的首達時間T定義為：

(2)

其中：δ是一個常量，表示防噴閥的失效閾值。則失效時間T的壽命分布為逆高斯分布[14]，若v已知，則其概率密度函數(shù)為：

(3)

(4)

(5)

可靠度函數(shù)R(t)為：

(6)

2 參數(shù)估計與剩余壽命預測

當產(chǎn)品的性能退化過程參數(shù)v已知時，給定當前退化量Δx，由Wiener過程的獨立性增量性質(zhì)可知，其剩余壽命分布的密度為：

(7)

由式(7)可知，只要估計出性能退化參數(shù)中的v，就可以預測產(chǎn)品的剩余壽命。但對于防噴閥來說，實際測量的退化數(shù)據(jù)較少，僅利用這些退化數(shù)據(jù)估計v精度不夠。但在其歷史實驗或者工業(yè)實踐中存在一些壽命數(shù)據(jù)，因此采用Bayes方法融合產(chǎn)品的歷史壽命數(shù)據(jù)和性能退化數(shù)據(jù)，以提高參數(shù)v的估計精度。

2.1 先驗分布的確定

(8)

采用EM算法對式(8)的未知參數(shù)進行估計，具體的計算步驟如下：

(9)

(10)

(11)

解上面的方程組可得遞推公式：

(12)

2.2 后驗分布的確定

由v的假設，可知其驗前分布為：

(13)

其中：v0=μ(0)和η0由前面的EM算法估計得到。假設時刻t0=0,t1,…,tk所測得的產(chǎn)品的性能退化量是Y0,Y1,…,Yk。記ΔYj=Yj-Yj-1，Δtj=tj-tj-1,j=1,2,…,k。由Bayes方法可知：

(14)

(15)

(16)

式(15)中的參數(shù)ηk,vk可根據(jù)上述方法遞推得到：

(17)

v的Bayes估計為：

(18)

將v的估計值代入式(7)，則可預測產(chǎn)品的剩余壽命并實現(xiàn)參數(shù)的在線更新。相對于參數(shù)的離線預測，在線估計可以對防噴閥的剩余壽命進行實時估計。

3 仿真研究

從圖1中可以看出：防噴閥的性能退化數(shù)據(jù)的軌跡大致呈線性單調(diào)變化，故可采用Wiener過程進行建模。從圖2可以看出，隨著監(jiān)測點個數(shù)的增多，估計的平均壽命逐漸趨近于實際值，說明本文提出的估計方法是可行的。為了說明本文方法的有效性，采用另外兩種方法進行比較，將基于歷史壽命數(shù)據(jù)進行剩余壽命估計的方法記為M1，基于自身退化數(shù)據(jù)進行剩余壽命估計的方法記為M2，將融合歷史壽命數(shù)據(jù)和自身退化數(shù)據(jù)進行壽命估計的方法記為M3。得到的防噴閥的剩余壽命概率密度函數(shù)如圖4所示，可以看出利用本文所提方法進行預測得到的剩余壽命密度函數(shù)更高更窄，預測效果更好。

圖1 防噴閥的循環(huán)次數(shù)退化軌跡

圖2 防噴閥的剩余壽命概率密度函數(shù)

圖3 防噴閥在三個時間點處的可靠度曲線

圖4 三種方法得到的剩余壽命的概率密度函數(shù)

4 結(jié)論

考慮了同類設備的個體差異性，將隨機的Wiener過程引入到防噴閥的性能退化建模中，由歷史統(tǒng)計數(shù)據(jù)得到總體的擴散系數(shù)，而個體在退化過程中的擴散系數(shù)和總體的擴散系數(shù)呈線性相關，這樣可以減少待估計參數(shù)的個數(shù)，使得計算過程更簡單有效。利用EM算法和Bayes方法估計分布參數(shù)的值，進而得到防噴閥的剩余壽命的概率密度函數(shù)及相關分布，并實現(xiàn)了參數(shù)的在線更新。進一步將研究含有測量誤差情況下的元器件的壽命預測問題。

[1]MOJTAHEDIA,LOTFOLLAHIYAGHINMA,HASSANZADEHY,etal.Arobustdamagedetectionmethoddevelopedforoffshorejacketplatformsusingmodifiedartificialimmunesystemalgorithm[J].ChinaOceanEngineering,2012,26(3):379-395.

[2]ZHANGBL,LIUYJ,MAH,etal.Discretefeedforwardandfeedbackoptimaltrackingcontrolforoffshoresteeljacketplatforms[J].OceanEngineering,2014,91:371-378.

[3]祁德慶,李帥,曹巧麗，等.海洋平臺損傷診斷研究[J].力學季刊,2008,29(2):329-335.QIDeqing,LIshuai,CAOQiaoli,etal.Damagediagnosisresearchonoffshoreplatform[J].ChineseQuarterlyofMechanics,2008,29(2):329-335.

[4]劉恩朋,楊占才,靳小波.國外故障預測與健康管理系統(tǒng)開發(fā)平臺綜述[J].測控技術(shù),2014,33(9):1-4.LIUEnpeng,YANGZhancai,JINXiaobo.Overviewofforeignprognosticsandhealthmanagementsystemdevelopmentplatform[J].Measurement&ControlTechnology,2014,33(9):1-4.

[5]王兆強,胡昌華,司小勝，等.一種基于維納過程和貝葉斯估計的獨立部件剩余壽命實時預報[C]//第25屆中國控制與決策會議論文集，沈陽：東北大學出版社,2013：4808-4902.

[6]ELGHAZELW,BAHJ,GUYEUXC,etal.Dependabilityofwirelesssensornetworksforindustrialprognosticsandhealthmanagement[J].ComputersinIndustry,2015,68(C):1-15

[7]LIUDT,ZHOUJH,PANDW,etal.Lithium-ionbatteryremainingusefullifeestimationwithanoptimizedrelevancevectormachinealgorithmwithincrementallearning[J].Measurement,2015,63:143-151.

[8]JARDINEAKS,LIND,BANJEVICD.Areviewonmachinerydiagnosticsandprognosticsimplementingcondition-basedmaintenance[J].MechanicalSystemsandSystemsandSignalProcessing,2006,20(7):1483-1510.

[9]司小勝,胡昌華,周東華.基于退化建模的備件需求預測[C]//第25屆中國控制和決策會議論文集，沈陽：東北大學出版社.2013：4817-4820.

[10]司小勝,胡昌華,張琪，等.不確定退化測量數(shù)據(jù)下的剩余壽命估計[J].電子學報,2015,43(1):30-35.SIXiaosheng,HUChanghua,ZHANGQi,etal.Estimatingremainingusefullifeunderuncertaindegradationmeasurements[J].ActaElectronicaSinica,2015,43(1):30-35.

[11]王浩偉,徐廷學,賀英政.Wiener過程在非線性退化建模中的應用[J].科技導報,2014,32(2):64-67.WANGHaowei,XUTingxue,HEYingzheng.ApplicationofWienerprocessesonmodelingnonlineardegradation[J].Science&TechnologyReview,2014,32(2):64-67.

[12]司小勝,胡昌華,李娟，等.具有不確定測量的非線性隨機退化系統(tǒng)剩余壽命預測[J].上海交通大學學報,2015,49(6):855-860.SIXiaosheng,HUChanghua,LIJuan,etal.Remainingusefullifepredictionofnonlinearstochasticdegradingsystemssubjecttouncertainmeasurements[J].JournalofShanghaiJiaoTongUniversity,2015,49(6):855-860.

[13]WHITMOREGA,SCHENKELBERGF.ModellingaccelerateddegradationdatausingWienerdiffusionwithatimescaletransformation[J].LifetimeDataAnalysis,1997,3(1):27-45.

[14]唐圣金,郭曉松,司小勝，等.基于維納過程的衛(wèi)星用光纖陀螺剩余壽命預測[J].紅外與激光工程,2013,42(12):3347-3352.TANGShengjin,GUOXiaosong,SIXiaosheng,etal.RemainingusefullifepredictionofFOGsusedforsatellitebasedonWienerprocess[J].InfraredandLaserEngineering,2013,42(12):3347-3352.

[15]SIXS,WANGGW,HUCH,etal.Awiener-process-baseddegradationmodelwitharecursivefilteralgorithmforremainingforremainingusefullifeestimation[J].MechanicalSystemsandSignalProcessing,2013,35(1-2):219-37.

[16]WANGX.Wienerprocesseswithrandomeffectsfordegradationdata[J].JournalofMultivariateAnalysis,2010,101(2):340-351.

[17]彭寶華,周經(jīng)倫,孫權(quán)，等.基于退化與壽命數(shù)據(jù)融合的產(chǎn)品剩余壽命預測[J].系統(tǒng)工程與電子技術(shù),2011,33(5):1073-1078.PENGBaohua,ZHOUJinglun,SUNQuan,etal.Residuallifetimepredicationofproductsbasedonfusionofdegradationdataandlifetimedata[J].SystemsEngineeringandElectronics,2011,33(5):1073-1078.

[18]裴洪，胡昌華，司小勝，等.融合壽命數(shù)據(jù)與退化數(shù)據(jù)的剩余壽命估計方法[J].電光與控制,2016,23(9):90-95.PEIHong,HUChanghua,SIXiaosheng,etal.Remainingusefullifeestimationbycombininglifetimedatawithdegradationdata[J].ElectronicsOptics&Control,2016,23(9):90-95.

(責任編輯：李 磊)

Remaining Lifetime Prediction of Blowout Preventer Valve Based on Fusion of Lifetime Data and Degradation Data

CHEN Xiurong1, LI Juan2, YU Jiaju1

(1. College of Science and Information, Qingdao Agricultural University, Qingdao, Shandong 266109, China;2. College of Mechanical and Electrical Engineering, Qingdao Agricultural University, Qingdao, Shandong 266109, China)

To reduce the number of the parameters which need to be estimated in the Wiener model, a new lifetime predication approach was proposed. In virtue of the relevant relation between the population variance and the individual variance, the number of the estimated parameters could be reduced by means of the proposed approach. By the aid of the history data and real-time degradation data, Bayesian estimation and expectation maximization (EM) algorithm were used to obtain parameters that reflect the performance degradation information. Thus, the probability density function and distribution of the remaining lifetime were obtained for the blowout preventer valve according to the estimated parameters. The simulation results show that the proposed approach can not only reduce the number of the estimated parameters and make the calculation simpler, but it can also improve the prediction accuracy of the remaining lifetime.

blowout preventer valve; remaining lifetime prediction; Wiener process; reliability; EM algorithm

2017-02-15

國家自然科學基金項目(41674037；61374126；61379029；61403223)；山東省自然科學基金項目(ZR2013FM021)

陳秀榮(1978—)，女，山東成武人，副教授，主要從事剩余壽命預測及故障診斷方面的研究. 李 娟(1969—)，女，山東海陽人，教授，博士，主要從事為故障診斷與壽命預測,智能監(jiān)測與控制方面的研究，本文通信作者.E-mail：lijuan291@sina.com

TB114.3

A

1672-3767(2017)05-0023-06

10.16452/j.cnki.sdkjzk.2017.05.004