[发明专利]一种利用Jeffery差异量的变压器运行状态振声检测方法和系统

权利要求书

1.一种利用Jeffery差异量的变压器运行状态振声检测方法，其特征在于，包括：

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

步骤102生成信号一阶差分序列，具体为：第n个信号一阶差分序列记为其第1个元素为0，第i个元素为其中，n为序列序号，取值范围为n＝1,2,···,N；N为所述信号序列S的长度；s_i为所述信号序列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变换系数为F₁(m₁)，其求取公式为第n个所述信号二阶差分序列的DFT变换系数为F₂(m₂)；其求取公式为其中，m₁为所述信号一阶差分序列的DFT变换系数序号，其取值范围为m₁＝1,2,···,N；m₂为所述信号二阶差分序列的DFT变换系数序号，其取值范围为m₂＝1,2,···,N；

步骤105求取归一化DFT系数，具体为：所述信号一阶差分序列的归一化DFT系数为P₁(m₁)，其计算公式为所述信号二阶差分序列的归一化DFT系数为P₂(m₂)，其计算公式为

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

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

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

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

2.一种利用Jeffery差异量的变压器运行状态振声检测系统，其特征在于，包括：

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

模块202生成信号一阶差分序列，具体为：第n个信号一阶差分序列记为其第1个元素为0，第i个元素为其中，n为序列序号，取值范围为n＝1,2,···,N；N为所述信号序列S的长度；s_i为所述信号序列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变换系数为F₁(m₁)，其求取公式为第n个所述信号二阶差分序列的DFT变换系数为F₂(m₂)；其求取公式为其中，m₁为所述信号一阶差分序列的DFT变换系数序号，其取值范围为m₁＝1,2,···,N；m₂为所述信号二阶差分序列的DFT变换系数序号，其取值范围为m₂＝1,2,···,N；

模块205求取归一化DFT系数，具体为：所述信号一阶差分序列的归一化DFT系数为P₁(m₁)，其计算公式为所述信号二阶差分序列的归一化DFT系数为P₂(m₂)，其计算公式为

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

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

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

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