基于异常值分析的变压器在线监测方法.docVIP

基于异常值分析的变压器在线监测方法 郭基伟，谢敬东，唐国庆 (东南大学电气工程系，南京210096) ??? 摘? 要：从变压器在线监测数据中异常值产生的机理出发，提出了一种利用在线观测值中异常点的发生时间和类型进行报警的变压器在线监测方法。实例表明，这种方法可用于运行条件存在较大差异，难以用统一的标准值来判别状态的设备，存在较高的实用价值。??? 关键词：变压器；在线监测；故障诊断；时间序列；异常值分析 1 变压器在线监测数据中异常值的产生　　变压器是电力系统中最重要的设备之一，一旦发生故障，将造成非常严重的损失。因此变压器的状态监测和故障诊断一直受到人们的高度重视。随着对变压器可靠运行要求的不断提高，传统的离线检测技术已经是难满足要求，因此近年来陆续开发了多种变压器在线监测装置。在线监测手段可以在设备不退出运行的情况下对设备进行连续的监测，及时发现设备的潜伏性故障，因而得到了越来越广泛的应用。　　目前的在线监测装置大部分都是采用一定的标准值作为评判标准，当监测的数据超过评判值时发出报警信号。但是，由于变压器运行的环境各不相同，各种装置的监测原理和数据采集精度也有所差别，因而在某台变压器上适用的标准值在另一台变压器上却往往会发生误报现象，在这种情况下，根据变压器自身监测资料的变化来判断变压器的状态应该是一种更好的选择。　　实际上变压器可以被看作为一个动态系统，在变压器正常运行时，其传递函数一定，观测到的资料可被视为一时间序列［1］；当变压器内部在故障时，传递函数发生变化，因此在观测到的数据中会产生异常点，通过检测和识别这些异常点的发生时间和类型，可以发现变压器内部是否存在故障，作为在线监测装置报警的依据。2 异常值的数学模型　　在变压器的在线监测资料中产生的异常值按照它们各自的性质和产生的机理，可以分为两种类型［2］。??? 1)附加异常值(additive outlier AO)　　这是一种孤立的异常值，通常这种异常值是由测量设备受到外界的干扰或是自身的扰动产生的误差造成的，它们的出现并不波及到邻近的观测值。就在线观测值的时间序列而言，这类异常值是非本质的，它们并不涉及到时间序列内在的相关结构。??? 2)新息异常值(innovation outlier IO)　　这是一种成片出现的异常值，它通过时间序列的相关结构，波及到与之邻近的一系列资料点。这类异常值是由于动态系统的结构发生了变换，从而使观测值的时间序列相关结构发生畸变造成的。　　在变压器在线监测中，AO异常值是由于仪器或人为误差造成的伪资料，分析时应当删除；IO异常值是由于变压器内部结构发生变化造成的，是可能的故障发生点，应当予以检测并作为报警依据。　　设不含异常值的在线监测资料时间序列ARIMA模型为??? 　　其中：｛Zt｝为不含异常值的在线监测资料时间序列；｛｝为零均值，方差e2的白噪声序列；??? ??? 式中：δT(t)表示时刻的冲击性脉冲，δT(T)＝0若t≠T，δT(t)＝1若t＝T；ω表示异常值的大小。??? 在时刻的新息异常值IO可以表示为［3］?? 　　由式(2)和式(3)可知，一个附加异常值AO只影响被观测序列的水平，而一个新息异常值通过由所描述的系统动态特性影响后面所有的观测值。考虑多个异常　的情况，其一般形式为　　　3 异常值检测和识别的迭代方法［］　　为了检测变压器在线监测资料中的异常值并作为报警的依据，必须识别出异常值发生的时间，并识别出异常值的类型。　　假设异常值发生在时刻，测量的在线监测资料序列为yt，t＝1，2，…，n。?? 　　根据最小二乘原理，IO模型中异常值影响ω的最小二乘估计就是时刻T的残差??? ??? ????　　在给定时刻是否存在AO型或IO型异常值的显著性检验，可以转化为检验式(7)或式(8)中的ω＝0还是ω≠0。设H0表示式(7)和式(8)中的ω＝0的假设，H1表示式(7)中ω≠0的假设，H2表示式(8)中ω≠0的假设，各个假设的似然比检验准则为：???? 　　在零假设H0下，λ1,T和λ2,T具有标准正态分布，λ3,T的分布比较复杂，Craig证明了它可以用收敛的贝赛尔函数来表示。验证一个时间序列是否存在IO或AO异常点可以通过检验??? 由于λ3，T的计算比较复杂，为了简化起见，IO和AO异常点的区分可以通过比较λ1,T和λ2,T的绝对值大小得到。如果｜λ1,T｜＞｜λ2,T｜，则将此异常点归IO异常点；如果｜λ1,T｜＞｜λ2,T｜，则将此异常点归AO异常点。　　由以上讨论可以得出在异常值出现的时间和模型参数未知情况下，检验异常值的迭代方法：　　1)假设序列不存在异常值，用常规方法对观测到的序列估计ARIMA模型；　　2)由估计