一种利用Jeffery差异量的变压器运行状态振声检测方法和系统与流程

本发明涉及电力领域，特别是涉及一种变压器运行状态振声检测方法及系统。

背景技术：

随着智能电网的高速发展，电力设备安全稳定运行显得尤其重要。目前，对超高压及以上电压等级的电力设备开展运行状态检测，尤其是对异常状态的检测显得愈加重要和迫切。电力变压器作为电力系统的重要组成部分，是变电站中最重要的电气设备之一，其可靠运行关系到电网的安全。

变压器运行状态检测的基本原理是提取变压器运行中的各个特征量，分析、辨识并跟踪特征量以此监测变压器的异常运行状态。当前变压器运行状态的常用检测方法中，包括检测局部放电的脉冲电流法和超声波检测法、检测绕组变形的频率响应法以及检测机械及电气故障的振动检测法等。这些检测方法主要检测变压器绝缘状况及机械结构状况，其中以变压器振动信号(振声)的检测最为全面，对于大部分变压器故障及异常状态均能有所反应。

虽然变压器振声检测方法在变压器运行状态监测中有着广泛的应用，且技术相对成熟，但是由于振声检测方法利用了变压器发出的振动信号，很容易受到环境噪声的影响，所以此方法在实际工作环境中应用时常常得不到令人满意的结果。

技术实现要素：

如前所述，变压器振声检测方法在变压器运行状态监测中有着广泛的应用，且技术相对成熟，但是由于振声检测方法利用了变压器发出的振动信号，很容易受到环境噪声的影响，所以此方法在实际工作环境中应用时常常得不到令人满意的结果。

本发明的目的是提供一种利用jeffery差异量的变压器运行状态振声检测方法和系统，所提出的方法利用了不同运行状态下变压器振声信号差值与背景噪声差值在jeffery差异量方面的差异，提高了状态监测的性能。所提出的方法具有较好的鲁棒性，计算也较为简单。

为实现上述目的，本发明提供了如下方案：

一种利用jeffery差异量的变压器运行状态振声检测方法，包括：

步骤101获取按时间顺序采集的信号序列s；

步骤102生成信号一阶差分序列，具体为：第n个信号一阶差分序列记为其第1个元素为0，第i个元素为其中，n为序列序号，取值范围为n＝1,2,···,n；n为所述信号序列s的长度；si为所述信号序列s的第i个元素，为所述信号序列s的第|i-1|n个元素，i为元素序号，取值范围为i＝1,2,···,n；||n表示以n为模取余数运算；

步骤103求取信号二阶差分序列，具体为：第n个信号二阶差分序列记为其第1个元素为0；第i个元素为其中，为所述信号序列s的第|i-2|n个元素，

步骤104求取dft变换系数，具体为：第n个所述信号一阶差分序列的dft变换系数为f1(m1)，其求取公式为第n个所述信号二阶差分序列的dft变换系数为f2(m2)；其求取公式为其中，m1为所述信号一阶差分序列的dft变换系数序号，其取值范围为m1＝1,2,···,n；m2为所述信号二阶差分序列的dft变换系数序号，其取值范围为m2＝1,2,···,n；

步骤105求取归一化dft系数，具体为：所述信号一阶差分序列的归一化dft系数为p1(m1)，其计算公式为所述信号二阶差分序列的归一化dft系数为p2(m2)，其计算公式为

步骤106求取中心频率，具体为：所述信号一阶差分序列的中心频率为其计算公式为所述信号二阶差分序列的中心频率为其计算公式为δf为所述信号序列s的采样频率；

步骤107求取jeffery差异量，具体为：第n个jeffery差异量记为hn，其计算公式为其中，dj为jeffery差异因子，其计算公式为m为求和参数；

步骤108求取状态判断阈值，具体为：状态判断阈值记为ε0，其计算公式为其中，k为求和参数；hk为第k个所述jeffery差异量；

步骤109判断变压器运行状态，具体为：判断第n个所述jeffery差异量hn是否大于或者等于所述状态判断阈值ε0；若第n个所述jeffery差异量hn大于或者等于所述状态判断阈值ε0，则在所述信号序列s的第n点处，变压器处于非正常运行状态；若第n个所述jeffery差异量hn小于所述状态判断阈值ε0，则在所述信号序列s的第n点处，变压器处于正常运行状态。

一种利用jeffery差异量的变压器运行状态振声检测系统，包括：

模块201获取按时间顺序采集的信号序列s；

模块202生成信号一阶差分序列，具体为：第n个信号一阶差分序列记为其第1个元素为0，第i个元素为其中，n为序列序号，取值范围为n＝1,2,···,n；n为所述信号序列s的长度；si为所述信号序列s的第i个元素，为所述信号序列s的第|i-1|n个元素，i为元素序号，取值范围为i＝1,2,···,n；||n表示以n为模取余数运算；

模块203求取信号二阶差分序列，具体为：第n个信号二阶差分序列记为其第1个元素为0；第i个元素为其中，为所述信号序列s的第|i-2|n个元素，

模块204求取dft变换系数，具体为：第n个所述信号一阶差分序列的dft变换系数为f1(m1)，其求取公式为第n个所述信号二阶差分序列的dft变换系数为f2(m2)；其求取公式为其中，m1为所述信号一阶差分序列的dft变换系数序号，其取值范围为m1＝1,2,···,n；m2为所述信号二阶差分序列的dft变换系数序号，其取值范围为m2＝1,2,···,n；

模块205求取归一化dft系数，具体为：所述信号一阶差分序列的归一化dft系数为p1(m1)，其计算公式为所述信号二阶差分序列的归一化dft系数为p2(m2)，其计算公式为

模块206求取中心频率，具体为：所述信号一阶差分序列的中心频率为其计算公式为所述信号二阶差分序列的中心频率为其计算公式为δf为所述信号序列s的采样频率；

模块207求取jeffery差异量，具体为：第n个jeffery差异量记为hn，其计算公式为其中，dj为jeffery差异因子，其计算公式为m为求和参数；

模块208求取状态判断阈值，具体为：状态判断阈值记为ε0，其计算公式为其中，k为求和参数；hk为第k个所述jeffery差异量；

模块209判断变压器运行状态，具体为：判断第n个所述jeffery差异量hn是否大于或者等于所述状态判断阈值ε0；若第n个所述jeffery差异量hn大于或者等于所述状态判断阈值ε0，则在所述信号序列s的第n点处，变压器处于非正常运行状态；若第n个所述jeffery差异量hn小于所述状态判断阈值ε0，则在所述信号序列s的第n点处，变压器处于正常运行状态。

根据本发明提供的具体实施例，本发明公开了以下技术效果：

如前所述，变压器振声检测方法在变压器运行状态监测中有着广泛的应用，且技术相对成熟，但是由于振声检测方法利用了变压器发出的振动信号，很容易受到环境噪声的影响，所以此方法在实际工作环境中应用时常常得不到令人满意的结果。

本发明的目的是提供一种利用jeffery差异量的变压器运行状态振声检测方法和系统，所提出的方法利用了不同运行状态下变压器振声信号差值与背景噪声差值在jeffery差异量方面的差异，提高了状态监测的性能。所提出的方法具有较好的鲁棒性，计算也较为简单。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍。显而易见，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明的方法流程示意图；

图2为本发明的系统流程示意图；

图3为本发明的具体实施案例流程示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

图1一种利用jeffery差异量的变压器运行状态振声检测方法的流程示意图

图1为本发明一种利用jeffery差异量的变压器运行状态振声检测方法的流程示意图。如图1所示，所述的一种利用jeffery差异量的变压器运行状态振声检测方法具体包括以下步骤：

步骤101获取按时间顺序采集的信号序列s；

步骤102生成信号一阶差分序列，具体为：第n个信号一阶差分序列记为其第1个元素为0，第i个元素为其中，n为序列序号，取值范围为n＝1,2,···,n；n为所述信号序列s的长度；si为所述信号序列s的第i个元素，为所述信号序列s的第|i-1|n个元素，i为元素序号，取值范围为i＝1,2,···,n；||n表示以n为模取余数运算；

步骤103求取信号二阶差分序列，具体为：第n个信号二阶差分序列记为其第1个元素为0；第i个元素为其中，为所述信号序列s的第|i-2|n个元素，

步骤104求取dft变换系数，具体为：第n个所述信号一阶差分序列的dft变换系数为f1(m1)，其求取公式为第n个所述信号二阶差分序列的dft变换系数为f2(m2)；其求取公式为其中，m1为所述信号一阶差分序列的dft变换系数序号，其取值范围为m1＝1,2,···,n；m2为所述信号二阶差分序列的dft变换系数序号，其取值范围为m2＝1,2,···,n；

步骤105求取归一化dft系数，具体为：所述信号一阶差分序列的归一化dft系数为p1(m1)，其计算公式为所述信号二阶差分序列的归一化dft系数为p2(m2)，其计算公式为

步骤106求取中心频率，具体为：所述信号一阶差分序列的中心频率为其计算公式为所述信号二阶差分序列的中心频率为其计算公式为δf为所述信号序列s的采样频率；

步骤107求取jeffery差异量，具体为：第n个jeffery差异量记为hn，其计算公式为其中，dj为jeffery差异因子，其计算公式为m为求和参数；

步骤108求取状态判断阈值，具体为：状态判断阈值记为ε0，其计算公式为其中，k为求和参数；hk为第k个所述jeffery差异量；

步骤109判断变压器运行状态，具体为：判断第n个所述jeffery差异量hn是否大于或者等于所述状态判断阈值ε0；若第n个所述jeffery差异量hn大于或者等于所述状态判断阈值ε0，则在所述信号序列s的第n点处，变压器处于非正常运行状态；若第n个所述jeffery差异量hn小于所述状态判断阈值ε0，则在所述信号序列s的第n点处，变压器处于正常运行状态。

图2一种利用jeffery差异量的变压器运行状态振声检测系统的结构意图

图2为本发明一种利用jeffery差异量的变压器运行状态振声检测系统的结构示意图。如图2所示，所述一种利用jeffery差异量的变压器运行状态振声检测系统包括以下结构：

模块201获取按时间顺序采集的信号序列s；

模块202生成信号一阶差分序列，具体为：第n个信号一阶差分序列记为其第1个元素为0，第i个元素为其中，n为序列序号，取值范围为n＝1,2,···,n；n为所述信号序列s的长度；si为所述信号序列s的第i个元素，为所述信号序列s的第|i-1|n个元素，i为元素序号，取值范围为i＝1,2,···,n；||n表示以n为模取余数运算；

模块203求取信号二阶差分序列，具体为：第n个信号二阶差分序列记为其第1个元素为0；第i个元素为其中，为所述信号序列s的第|i-2|n个元素，

模块204求取dft变换系数，具体为：第n个所述信号一阶差分序列的dft变换系数为f1(m1)，其求取公式为第n个所述信号二阶差分序列的dft变换系数为f2(m2)；其求取公式为其中，m1为所述信号一阶差分序列的dft变换系数序号，其取值范围为m1＝1,2,···,n；m2为所述信号二阶差分序列的dft变换系数序号，其取值范围为m2＝1,2,···,n；

模块205求取归一化dft系数，具体为：所述信号一阶差分序列的归一化dft系数为p1(m1)，其计算公式为所述信号二阶差分序列的归一化dft系数为p2(m2)，其计算公式为

模块206求取中心频率，具体为：所述信号一阶差分序列的中心频率为其计算公式为所述信号二阶差分序列的中心频率为其计算公式为δf为所述信号序列s的采样频率；

模块207求取jeffery差异量，具体为：第n个jeffery差异量记为hn，其计算公式为其中，dj为jeffery差异因子，其计算公式为m为求和参数；

模块208求取状态判断阈值，具体为：状态判断阈值记为ε0，其计算公式为其中，k为求和参数；hk为第k个所述jeffery差异量；

模块209判断变压器运行状态，具体为：判断第n个所述jeffery差异量hn是否大于或者等于所述状态判断阈值ε0；若第n个所述jeffery差异量hn大于或者等于所述状态判断阈值ε0，则在所述信号序列s的第n点处，变压器处于非正常运行状态；若第n个所述jeffery差异量hn小于所述状态判断阈值ε0，则在所述信号序列s的第n点处，变压器处于正常运行状态。

下面提供一个具体实施案例，进一步说明本发明的方案

图3为本发明具体实施案例的流程示意图。如图3所示，具体包括以下步骤：

步骤301获取按时间顺序采集的信号序列s；

步骤302生成信号一阶差分序列，具体为：第n个信号一阶差分序列记为其第1个元素为0，第i个元素为其中，n为序列序号，取值范围为n＝1,2,···,n；n为所述信号序列s的长度；si为所述信号序列s的第i个元素，为所述信号序列s的第|i-1|n个元素，i为元素序号，取值范围为i＝1,2,···,n；||n表示以n为模取余数运算；

步骤303求取信号二阶差分序列，具体为：第n个信号二阶差分序列记为其第1个元素为0；第i个元素为其中，为所述信号序列s的第|i-2|n个元素，

步骤304求取dft变换系数，具体为：第n个所述信号一阶差分序列的dft变换系数为f1(m1)，其求取公式为第n个所述信号二阶差分序列的dft变换系数为f2(m2)；其求取公式为其中，m1为所述信号一阶差分序列的dft变换系数序号，其取值范围为m1＝1,2,···,n；m2为所述信号二阶差分序列的dft变换系数序号，其取值范围为m2＝1,2,···,n；

步骤305求取归一化dft系数，具体为：所述信号一阶差分序列的归一化dft系数为p1(m1)，其计算公式为所述信号二阶差分序列的归一化dft系数为p2(m2)，其计算公式为

步骤306求取中心频率，具体为：所述信号一阶差分序列的中心频率为其计算公式为所述信号二阶差分序列的中心频率为其计算公式为δf为所述信号序列s的采样频率；

步骤307求取jeffery差异量，具体为：第n个jeffery差异量记为hn，其计算公式为其中，dj为jeffery差异因子，其计算公式为m为求和参数；

步骤308求取状态判断阈值，具体为：状态判断阈值记为ε0，其计算公式为其中，k为求和参数；hk为第k个所述jeffery差异量；

步骤309判断变压器运行状态，具体为：判断第n个所述jeffery差异量hn是否大于或者等于所述状态判断阈值ε0；若第n个所述jeffery差异量hn大于或者等于所述状态判断阈值ε0，则在所述信号序列s的第n点处，变压器处于非正常运行状态；若第n个所述jeffery差异量hn小于所述状态判断阈值ε0，则在所述信号序列s的第n点处，变压器处于正常运行状态。

本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述较为简单，相关之处参见方法部分说明即可。

本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。