6-3-7基于MKPCA的间歇过程故障检测终

PAGEPAGE1基于多向核主元分析的啤酒发酵过程故障诊断模型摘要：针对主元分析故障诊断模型在非线性时变过程中应用的局限性，基于间歇过程的周期性特点，将核变换理论引入非线性空间的数据特征提取中，提出了一种改进的多向核主元分析故障诊断模型，有效地解决了过程数据的非线性问题，保证数据信息抽取的完整性。通过与其他方法的对比实验，结果表明所提出的方法对缓慢时变的间歇过程具有良好的实时性和准确性。关键词：间歇过程故障检测多向核主元分析1引言间歇过程是批次生产的重复过程，广泛应用于生物制药、化工原料、食品等行业，其具有生产过程重复性高、动态特性变化快、建模困难等特点，这导致传统的故障诊断方法难以得到较好的应用效果[1]。主元分析(principalcomponentanalysis,PCA)是多元统计过程监测(multivariatestatisticalprocessmonitoring,MSPM)的重要方法之一，但是PCA在过程建模时假定过程是线性的，这导致在具有强非线性生产过程的在线监测中存在误报率过高的现象[2]。近年来针对间歇过程提出的多向主元分析(MPCA)方法得到了较多的研究，然而MPCA实质上仍是一种线性化建模方法，对复杂的非线性过程在线监控的可靠性和实时性往往也难以保证[3]。针对非线性过程监测的建模问题，Scholkopf等人将核函数理论引入到统计过程监控中，将主元分析(PCA)方法推广到代表非线性领域的高维特征空间，据此发展的KPCA模型可以从数据样本中提取出非线性特征，与PCA算法相比，该方法表现出更优的监测性能[4]。本文针对间歇过程特点，将核函数理论应用于多向主元分析中，提出一种改进的多向核主元分析（MKPCA）过程故障监测算法，并通过啤酒发酵过程的故障检测实验对算法性能进行了验证。2核主元分析（KPCA）核主元分析通过非线性映射将输入集合映射到一个高维特征空间，使数据具有更好的可分性，再对高维空间的映射数据进行PCA处理，得到非线性主元。KPCA不直接计算特征向量，而是将其转化为求核矩阵的特征值和特征向量，避免了在特征空间求特征向量，而数据在特征向量上的投影转换为求核函数的线性组合，大大简化了计算量[5]。首先通过非线性映射函数：，将输入空间，k=1,2,…,M映射到特征空间F：，k=1,2,...,M中，然后在该特征空间中对式(1)的协方差矩阵进行线性主元分析。(1)在特征空间中计算主元，可通过求解式(2)中的特征值和特征向量得到：(2)将每个样本与式(2)作内积，可得式(3)。(3)因为式(2)的所有解均在张成的子空间内，所以存在系数使得式(4)成立。(4)对式(2)、(3)和(4)进行合并，得式(5)。(5)取作为核函数，，可得到式(6)。(6)式中，,其特征向量所对应的特征值为，为了提取主元特征，将投影到上可得到式(7)。，(7)式(7)称为KPCA的第k个主元。3多向核主元分析故障诊断模型对于间歇过程其数据集比连续过程数据集多一维“批量”元素，每批数据都可以看作一个二维数据阵，多批数据则构成了三维数据阵，其中I为批次数目，J为变量数目，K为采样点数。将数据按批次方向展开为，X的每一行均表示一个批次数据,如图1示。图1MKPCA建模三维数据矩阵展开后，数据处理和分析过程等同于KPCA方法[6]。建模步骤如下：(1)对于数据集按批次方向展开成二维数据阵，并对其按式(8)进行标准化。(8)式中：—x(j)的样本均值，S(j)—x(j)的样本标准差。(2)计算核矩阵K，记其元素为，其中：(9)(3)在特征空间中，根据式(10)和(11)对核矩阵进行标定得到。(10)(11)其中：。(4)对核矩阵进行特征值分解,并且使得满足式(12)。(12)(5)对于每一个正常批次的数据x，根据式(7)提取其非线性主元。(6)按式(13)和(14)构建监控统计量和SPE。(13)(14)(7)按式(15)和(16)确定统计量的置信限。(15)其中：n为样本个数，m为主元个数，是检验水平为、自由度为m，n-1时的F分布临界值。(16)其中：为建模所用数据的协方差矩阵的特征值，是当检验水平为时的正态分布临界值，M是全部主元个数，m为主元模型中的主元个数。Inthis:isusedinmodelingofthedatacovariancematrixeigenvalue,iswhenthetestlevelisnormaldistributioncriticalvalues,Misthetotalnumberofprincipalcomponents,masthenumberofprincipalcomponentsinthePCAmodel.运用多向核主元法对间歇过程进行故障检测的步骤如下：Usingmultiwaykernelprincipalcomponentmethodforfaultdetectionofbatchprocess，itsstepsareasfollows:当对批次进行在线监测时，仅可知自批次开始时刻到监测时刻的采样数据。然而，监测过程的测试数据应为完整的批次数据。因此，需要对自监测时刻至批次结束时刻的数据进行估计。针对此问题已经提出了多种方法，本文采用各变量的均值来代替其估计值。Whentheon-linemonitoringofbatch,Onlyknown,thesamplingdatasincebatchmonitoringtimetoMonitoringtime.However,testdataofMonitorprocessshallbethecompletebatchdata.Therefore,needtobeestimatedatasincemonitoredthemomenttotheendofbatchmoment.Thedatasincemonitoredthemomenttotheendofbatchmomentneedtobeestimated.Tosolvethisproblem,avarietyofmethodshavebeenproposed,inthispaper,usingthemeanofeachvariabletoreplacetheestimates.(1)在第k个采样时刻，新的反应批次数据为，展开处理采集到的数据，得到展开后的数据矩阵，对此矩阵依据式(8)进行标准化。(1)Inthefirstksamplingtime,Thenewreactionbatchdatais,processingsampleddatagettheunfoldeddatamatrix,tostandardizethematrixbasedonthistype(8).(2)估计新批次未反应完时刻的数据，补足第一步标准化后的数据矩阵，得到，作为完整的新批次数据。(2)Estimationofthenewbatchdidnotreacttimedata,supplyingthefirststepofthestandardizeddatamatrix,gettingasanewintegritybatchdata.(3)根据式(9)计算测试数据相应的核向量。(3)Accordingtoequation(9)calculationthetestdatacorrespondingkernelvector(4)根据式(17)对核向量作标准化处理得到。(4)Tostandardizekernelvectoraccordingtothetype(17)getting(17)其中：K和在训练时得到，。Amongthem:Kandobtainedduringtraining,(5)根据式(18)提取非线性主元。(5)Accordingtoequation(18)extractnonlinearprincipalcomponent.(18)(6)按式(13)和(14)分别计算测试数据的和SPE统计量，并判断是否超出了各自的置信限。如果出现超出其置信限的情况，则说明过程中出现了故障。(6)Accordingtoformula(13)and(14)respectivelytocalculatethetestdataoftheandSPEstatistics,anddeterminewhetheritbeyondtherespectiveconfidencelimits.Ifthereisaconditionthatbeyonditsconfidencelimits,thenitappearedfailureintheprocess.4实验研究实验采用微型啤酒生产装置，测试数据来自发酵过程监控数据。根据生产运行中各变量的活跃程度和对生产状态的影响，选择温度、压力、液位、糖度、PH值和酒精度6个过程监测变量，这些变量反应了酵母菌菌体生长和发酵产物的合成状况。过程周期15天，每1小时采样1次，每批次采样360次。实验选取12个正常批次的数据建模。由于每一批次数据(为采样次数)的反应时间不同，因此，在将转换成之后，对多于2160列的直接截取到2160列，对不足2160列的批次补零，然后将矩阵排列成形式，进行标准化处理，核函数采用径向基核函数，按93%的累计贡献率提取主成分。其中，MPCA算法的主元数目为2；而MKPCA算法的主元数目为4。可以看出，MKPCA算法所选的主元数目高于MPCA算法所选的主元数目，这是由于前者从高维特征空间中提取主元，而后者从输入空间中提取主元。4.Theexperimentalstudy Thisexperimentuseddeviceforminiaturebeerproduction,testingdatafromthefermentationprocesscontroldata.Accordingtotheactivedegreeofeachvariableintheproductionfunctionandtheinfluenceontheproductionstatus,choosingthetemperature,pressure,liquidlevel,sugardegree,PHvalueandalcoholdegree,sixprocessmonitoringvariables,thesevariableshasbeensynthesizedbythereactionofyeastcellgrowthandthefermentationproducts.15daysasaprocesscycle,sampling1timesevery1hours,eachbatchessamples360times.Theexperimentselected12normalbatchesofdatamodeling.Becauseeachbatchofdata(isthenumberofsampling)ofdifferentreactiontime,Therefore,afterconvertingto,directinginterceptionofmorethan2160to2160,tolessthan2160batchesofzeropadding,thenthematrixisarrangedintheformof,standardtreatment,kernelfunctionusingrbfkernelfunction,Accordingto93%ofthecontributionratetoextractprincipalcomponent.Amongthem,principalcomponentnumberoftheMPCAalgorithmis2;TheprincipalcomponentnumberofMKPCAalgorithmis4.ItisshownthatprincipalcomponentnumberselectedinMKPCAalgorithmishigherthanthatselectedinMPCAalgorithm.Thisisduetotheformerfromhighdimensionalfeaturespacetoextracttheprincipalcomponent,andthelatterfromtheinputspacetoextracttheprincipalcomponent.Figure4PCAstatisticsmonitoringchartFigure4PCASPEstatisticsmonitoringchartFigure4MPCAstatisticsmonitoringchartFigure5MPCASPEstatisticsmonitoringchartFigure4MKPCAstatisticsmonitoringchartFigure4MKPCASPEstatisticsmonitoringchart对啤酒发酵过程进行在线监测，在317-360采样时刻引入压力传感器故障，对测试数据分别采用PCA算法、MPCA算法和MKPCA算法进行在线监测。PCA的和SPE监测结果如图2，3所示。MPCA的和SPE监测结果如图4，5所示。MKPCA的和SPE监测结果如图6，7所示。Theon-linemonitoringofthebeerfermentationprocess,Thepressuresensorfaultwasintroducedin317-360samplingtime,PCAalgorithm,MPCAalgorithmandMKPCAalgorithmwereusedfortheonlinemonitoringofbeerfermentationprocess.ThemonitoringresultsofstatisticsandSPEstatisticsaboutPCAwereshowninFigure2,3.ThemonitoringresultsofstatisticsandSPEstatisticsaboutMPCAwereshowninFigure4,5.ThemonitoringresultsofstatisticsandSPEstatisticsaboutMKPCAwereshowninFigure6,7.实验结果分析：图1中PCA的统计量在故障时刻不能检测出压力传感器故障的存在，并且在第12和34采样时刻还存在着故障误报现象，图2中PCA的SPE统计量在317-360采样时刻能够及时的检测出故障。由于统计量没有检测出过程故障而SPE统计量检测出了过程故障，所以PCA算法不能实现对啤酒发酵过程的监测；从图3、4中可以看出，当采用MPCA算法在线监测时，图3的统计量在317-351采样时刻并没有检测出过程故障，而在352-360采样时刻检测出了过程故障，所以MPCA算法的统计量应用在啤酒发酵过程时存在检测滞后的现象，即不能及时检测出故障。图4的SPE统计量在317-360采样时刻能够及时的检测出了过程故障。同理，MPCA算法也不能及时准确的实现对啤酒发酵过程的在线监测；从图5、6中可以看出，通过引入核函数并结合MPCA算法复合而成的MKPCA算法的统计量和SPE统计量都能及时准确的检测出过程故障，而且不存在误报现象。因此采用MKPCA算法用于啤酒发酵过程的在线监测较PCA算法和MPCA算法可靠。Analysisofexperimentalresults:ThepressuresensorfaultcannotbedetectedinthefaulttimefromthestatisticsofPCAinfigure1，andtherearefaultmisreportingphenomenoninthe12and34samplingtime.Thepressuresensorfaultcanbedetectedinthe317-360samplingtimefromtheSPEstatisticsofPCAinfigure2intime.Becausethepressuresensorfaultcan’tbedetectedfromthestatisticsandthepressuresensorfaultcanbedetectedfromtheSPEstatistics.SoPCAalgorithmcan'tbeusedfortheonlinemonitoringofbeerfermentationprocess.Ascanbeseenfromthefigure3,4,WhenusingMPCAalgorithmon-linemonitoring,Figure3,thestatisticsin317-351thesamplingtimedidn'tdetectprocessfaults,however,theprocessfaultsweremonitoredin352-360samplingtimes,therefore,theapplicationofthestatisticsofMPCAalgorithmfortheonlinemonitoringofbeerfermentationprocessexistthephenomenonofhysteresis.Thatcan'tdetectthefaultintime.Figure4,SPEstatisticsin317-360samplingtimecandetectedtheprocessfaulttimely.Inthesameway,MPCAalgorithmcan'ttimelyandaccuratelyrealizetheon-linemonitoringofthebeerfermentationprocess;Inthefigure5and6,byintroducingkernelfunctionandcombiningtheMPCAalgorithmofcompositeMKPCATstatisticandSPEstatisticsofthealgorithmcanaccuratelyandtimelydetectprocessfaults,andthereisnofalsepositives.Aboveall,UsingMKPCAalgorithmisbetterthanPCAalgorithmandMPCAalgorithm通过实验结果可知，引入非线性核函数能够充分提取过程中存在的非线性信息，有效计算出高维特征空间中的主元。与PCA和MPCA算法相比，MKPCA算法表现出更好的监测性能，更适于对非线性间歇过程进行在线监测。Throughtheexperimentalresultswecanknowthatbyintroducingthenonlinearkernelfunctioncanfullyextractthenonlinearinformationwhichexistedintheprocess,principalcomponentinthehighdimensionalfeaturespacecanbecalculatedeffectively.ComparedwithPCAandMPCAalgorithm,MKPCAalgorithmshowsbettermonitoringperformance,moresuitableforon-linemonitoringofnonlinearbatchprocess.5结论本文针对间歇发酵过程缓慢时变和非线性等特点，利用核理论方法对MPCA算法进行了改进，提出了适用的多向核主元分析故障诊断算法。通过引入非线性核函数，能够充分提取过程中存在的非线性信息，有效计算出高维特征空间中的主元，并将研究结果应用于啤酒发酵过程监测。通过与PCA算法、MPCA算法进行对比实验表明所提出的模型可以有效处理间歇过程批次间存在的非线性属性，获取过程变量间的非线性关系，提高了故障诊断的及时性和准确性。5Conclusion Thisarticlebasedontheintermittentfermentationprocessslowtime-varying,nonlinearandothercharacteristics,usingofkerneltheorymethodimprovedtheMPCAalgorithm,itputsforwardthesuitablemultiwaykernelprincipalcomponentanalysisalgorithmforfaultdiagnosis.Byintroducingnonlinearkernelfunction,tofullyextractthenonlinearinformationwhichexistintheprocess.Effectivelycalculatetheprincipalcomponentinthehighdimensionalfeaturespace,andtheresearchresultscanbeappliedtobeerfermentationprocessmonitoring.ThroughwiththePCAalgorithmandtheMPCAalgorithmcomparativeexperiments,itshowthattheexistenceofthenonlinearpropertywhereamongbatchprocessbatchcanbeeffectivelytreatedbytheprogramsmodelwhathavebeenproposed,obtainingthenonlinearrelationshipamongtheprocessvariables,Improvingthetimelinessandaccuracyoffaultd