发明名称 局部放电在线检测中的特征因子提取方法
摘要 本发明涉及一种局部放电在线检测中的特征因子提取方法,其特点是,包括如下步骤:通过对放电脉冲时频信息的分析,提出能够表征脉冲时频信息的特征量,然后对脉冲进行分类,分类完成后再计算出每一子类内表征脉冲信号特征的特征因子,以形成用于放电模式识别的特征因子。本发明提出了一种描述放电脉冲时频特征的新方法,通过此方法可以提取表征放电脉冲波形特征的高维特征量,以满足分类算法对系统特征量的需求。
申请公布号 CN103675610B 申请公布日期 2016.06.15
申请号 CN201310456277.0 申请日期 2013.09.29
申请人 国家电网公司;国网宁夏电力公司银川供电公司 发明人 鲍永胜;郝婷
分类号 G01R31/12(2006.01)I 主分类号 G01R31/12(2006.01)I
代理机构 宁夏专利服务中心 64100 代理人 赵明辉
主权项 一种局部放电在线检测中的特征因子提取方法,其特征在于,包括如下步骤:通过对放电脉冲时频信息的分析,提出能够表征脉冲时频信息的特征量,然后对脉冲进行分类,分类完成后再计算出每一子类内表征脉冲信号特征的特征因子,以形成用于放电模式识别的特征因子;其中提出能够表征脉冲时频信息的特征量具体是:21)为便于脉冲群的分类,对于采集到的放电脉冲信号做如下处理<maths num="0001"><math><![CDATA[<mrow><msub><mi>p</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>a</mi><mn>0</mn></msub><mo>,</mo><msub><mi>a</mi><mn>1</mn></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>a</mi><mi>i</mi></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>a</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo><mi>&Delta;</mi><mi>t</mi><mo>,</mo><mo>...</mo><mo>,</mo><mi>&Delta;</mi><mi>t</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo><mi>&Delta;</mi><mi>t</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000902760590000011.GIF" wi="1149" he="190" /></maths>式中,j表示第j个脉冲,n表示脉冲由n个点组成;a<sub>i</sub>为第i个点的时域波形值;Δt为采样时间间隔;Δt(i‑1)为第i个点对应的时间;22)对脉冲波形做傅里叶变换可得<maths num="0002"><math><![CDATA[<mrow><msub><mi>P</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>A</mi><mn>0</mn></msub><mo>,</mo><msub><mi>A</mi><mn>1</mn></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>A</mi><mi>i</mi></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>A</mi><mrow><mi>n</mi><mo>/</mo><mn>2</mn><mo>-</mo><mn>1</mn></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo><mi>&Delta;</mi><mi>f</mi><mo>,</mo><mo>...</mo><mo>,</mo><mi>&Delta;</mi><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo><mi>&Delta;</mi><mi>f</mi><mrow><mo>(</mo><mi>n</mi><mo>/</mo><mn>2</mn><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000902760590000012.GIF" wi="1285" he="215" /></maths>式中,A<sub>i</sub>为第i个点的频谱幅值;Δf(i‑1)为第i个点的频率值;23)通过下式求取信号的时间重心<img file="FDA0000902760590000014.GIF" wi="61" he="69" />和频率重心<img file="FDA0000902760590000015.GIF" wi="100" he="71" /><maths num="0003"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msubsup><mi>T</mi><mn>0</mn><mi>j</mi></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>t</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msubsup><mi>p</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>p</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>F</mi><mn>0</mn><mi>j</mi></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>/</mo><mn>2</mn><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>f</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msubsup><mi>p</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>/</mo><mn>2</mn><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>p</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000902760590000013.GIF" wi="902" he="333" /></maths>24)将信号脉冲的波形作为时间变量的概率密度分布,其傅里叶变换所得各频谱分量的幅值作为频率变量的概率密度分布,则可以引入变量的信息量作为表征脉冲波形特征的特征量,通过下式便可计算出表征放电脉冲波形时频信息的新的高维特征量<maths num="0004"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msubsup><mi>T</mi><mi>j</mi><mi>k</mi></msubsup><mo>=</mo><msup><mrow><mo>{</mo><msub><mi>log</mi><mn>2</mn></msub><mo>&lsqb;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>T</mi><mn>0</mn><mi>j</mi></msubsup><mo>)</mo></mrow><mrow><mn>2</mn><mi>k</mi></mrow></msup><mo>&CenterDot;</mo><msubsup><mi>p</mi><mi>j</mi><mrow><mn>2</mn><mi>k</mi></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><msubsup><mi>p</mi><mi>j</mi><mrow><mn>2</mn><mi>k</mi></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo><mo>}</mo></mrow><mi>k</mi></msup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>F</mi><mi>j</mi><mi>k</mi></msubsup><mo>=</mo><msup><mrow><mo>{</mo><msub><mi>log</mi><mn>2</mn></msub><mo>&lsqb;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>/</mo><mn>2</mn><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>F</mi><mn>0</mn><mi>j</mi></msubsup><mo>)</mo></mrow><mrow><mn>2</mn><mi>k</mi></mrow></msup><mo>&CenterDot;</mo><msubsup><mi>p</mi><mi>j</mi><mrow><mn>2</mn><mi>k</mi></mrow></msubsup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>/</mo><mn>2</mn><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>p</mi><mi>j</mi><mrow><mn>2</mn><mi>k</mi></mrow></msubsup><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo><mo>}</mo></mrow><mi>k</mi></msup></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000902760590000021.GIF" wi="1110" he="319" /></maths>式中,k=2,...,m,对同一类中的所有放电脉冲进行特征提取后,即得脉冲群的特征向量序列<img file="FDA0000902760590000022.GIF" wi="438" he="102" />j=1,2,...,N,N为脉冲个数,通过选取不同的m值就可以得到不同维数的特征向量序列,以满足分类算法对脉冲信号特征量的需求。
地址 100031 北京市西城区西长安街86号