发明名称 一种基于马拉算法的电力系统故障录波数据分析方法
摘要 一种基于马拉算法的电力系统故障录波数据分析方法:1.录波数据的采样:每个周期的采样点数取N=128,对应的采样频率fs=Nff=128×50=6400Hz;2.暂态波形的提取;3.尺度分析;4.马拉算法的分解与重构:对信号进行5层多分辨分析,再利用小波变换检测信号的突变点,检测模极大值;5.故障持续时间和故障波形分析。本发明可以大大降低小波变换的计算量,有利于对含有大量信息的电力系统故障信号进行实时处理,很好地分析电力系统故障过程并提取故障特征基因,从而进行故障的诊断以及定位。
申请公布号 CN103018629A 申请公布日期 2013.04.03
申请号 CN201210483763.7 申请日期 2012.11.23
申请人 广东电网公司电力科学研究院;科大智能科技股份有限公司;广东粤恒电力科技有限公司 发明人 高新华;陶维青;余南华;李林;陈炯聪;何刚;李传健;郭晋楠;李瑞;单开;周克林;柳慧超;黄向明
分类号 G01R31/08(2006.01)I 主分类号 G01R31/08(2006.01)I
代理机构 广州知友专利商标代理有限公司 44104 代理人 周克佑
主权项 一种基于马拉算法的电力系统故障录波数据分析方法,其特征是:包括以下步骤:S1录波数据采样每个周期的采样点数取N=128,对应的采样频率fs=Nff=128*50=6400Hz;S2暂态波形提取采用马拉算法使用多贝西滤波器对故障波形进行分解和重构,在第一尺度下将故障录波信号的频谱分离成高频带和低频带;S3尺度分析在第一尺度下,以多贝西4为小波基,得出下式: <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>w</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>w</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>w</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>w</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mn>64</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mi>H</mi> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>x</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>x</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>x</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mi>x</mi> <mrow> <mo>(</mo> <mn>128</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> <mrow> <mi>H</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>w1是第一尺度下的小波系数,h是小波滤波器的系数;同样,第一尺度下,以多贝西4为小波基的马拉算法的近似系数由下式给出: <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>a</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>a</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>a</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>a</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mn>64</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mi>G</mi> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>x</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>x</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>x</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mi>x</mi> <mrow> <mo>(</mo> <mn>128</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> <mrow> <mi>G</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mtd> <mtd> <mi>g</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>a1是第一尺度下的近似系数,g是尺度滤波器的系数;不断递归计算出不同尺度下的小波系数和近似系数;S4马拉算法的分解与重构对信号进行5层多分辨分析,再利用小波变换检测信号的突变点,检测模极大值;模极大值的检测步骤如下:S3‑1对数据窗内的采样序列进行小波分解,得到第一层高频系数和第二层高频系数;S3‑2求第一层高频系数和第二层高频系数的模平均值和模最大值;S3‑3比较模最大值与模平均值,若大于设定的阈值,则认为已检测到突变点,否则认为没有检测到突变点;S3‑4若检测到突变点,则记下突变发生的时间,并进一步进行小波分解和重构,若没有检测到突变点,则视为稳态信号;S5故障和故障持续时间分析如果在第一层和第二层高频系数中检测到了模极大值点,说明存在故障或电能质量扰动,否则未出现故障或电能质量扰动;模极大值点还可以得出故障或扰动的持续时间;由扰动持续时间对其所在的频带进行重构,提取出扰动波形,通过重构的扰动波形对扰动进一步分析:如果是短期或者长期变化扰动,则调用相应的短期或者长期变化扰动识别子程序;如果是暂态扰动,则将扰动波形作为人工神经网络和模糊专家系统的输入信号;如果是负荷扰动,则进一步判断是负荷变化引起的扰动还是负荷本身引起的扰动。
地址 510080 广东省广州市东风东路水均岗8号