发明名称 一种变压器有载分接开关运行状态的在线监测方法
摘要 一种变压器有载分接开关运行状态的在线监测方法,属信号监测领域。其实时采集变压器箱壁的机械振动信号;使用延迟坐标法对振动信号x(t)进行相空间重构;确定重构相空间的指定区域V;计算空间相点的密度分布;计算有载分接开关动作前后空间密度曲线的相关系数ρ;若相关系数ρ>0.8,则判断有载分接开关运行状态正常;若相关系数ρ<0.8,则判断有载分接开关运行状态异常。其通过对有载分接开关动作过程中的振动信号进行实时监控,能够实现对变压器有载分接开关工作状态的高效、准确地判断。可广泛用于各种变压器有载分接开关的在线监测与故障诊断领域。
申请公布号 CN103048036B 申请公布日期 2016.06.29
申请号 CN201210525518.8 申请日期 2012.12.07
申请人 上海市电力公司;国家电网公司 发明人 傅坚;林嘉杨;王丰华;金之俭;周翔
分类号 G01H1/00(2006.01)I;G01M13/00(2006.01)I 主分类号 G01H1/00(2006.01)I
代理机构 上海三和万国知识产权代理事务所(普通合伙) 31230 代理人 蔡海淳
主权项 一种变压器有载分接开关运行状态的在线监测方法,包括采集变压器有载分接开关表面的振动信号,使用延迟坐标法对振动信号时间序列x(t)进行相空间重构和确定重构相空间的指定区域V,所述的使用延迟坐标法对振动信号时间序列x(t)进行相空间重构,为<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mi>X</mi><mo>(</mo><mn>1</mn><mo>)</mo><mo>=</mo><mo>&lsqb;</mo><mi>x</mi><mo>(</mo><mn>1</mn><mo>)</mo><mo>,</mo><mi>x</mi><mo>(</mo><mn>1</mn><mo>+</mo><mi>&tau;</mi><mo>)</mo><mo>,</mo><mo>...</mo><mo>,</mo><mi>x</mi><mo>(</mo><mn>1</mn><mo>+</mo><mo>(</mo><mi>m</mi><mo>-</mo><mn>1</mn><mo>)</mo><mi>&tau;</mi><mo>)</mo><mo>&rsqb;</mo></mtd></mtr><mtr><mtd><mi>X</mi><mo>(</mo><mn>2</mn><mo>)</mo><mo>=</mo><mo>&lsqb;</mo><mi>x</mi><mo>(</mo><mn>2</mn><mo>)</mo><mo>,</mo><mi>x</mi><mo>(</mo><mn>2</mn><mo>+</mo><mi>&tau;</mi><mo>)</mo><mo>,</mo><mo>...</mo><mo>,</mo><mi>x</mi><mo>(</mo><mn>2</mn><mo>+</mo><mo>(</mo><mi>m</mi><mo>-</mo><mn>1</mn><mo>)</mo><mi>&tau;</mi><mo>)</mo><mo>&rsqb;</mo></mtd></mtr><mtr><mtd><mo>...</mo><mo>...</mo></mtd></mtr><mtr><mtd><mi>X</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>=</mo><mo>&lsqb;</mo><mi>x</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>,</mo><mi>x</mi><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo><mo>,</mo><mo>...</mo><mo>,</mo><mi>x</mi><mo>(</mo><mi>k</mi><mo>+</mo><mo>(</mo><mi>m</mi><mo>-</mo><mn>1</mn><mo>)</mo><mi>&tau;</mi><mo>)</mo><mo>&rsqb;</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000848411510000011.GIF" wi="1006" he="301" /></maths>式中,m为嵌入维数,τ为延迟时间,这k个向量形成了一个重构的相空间;其特征在于所述的在线监测方法包括下列步骤:(1)将采集到的振动信号进行抗混迭数字滤波和高速缓存,然后通过高速总线传输至数据分析模块;(2)在使用延迟坐标法对振动信号时间序列x(t)进行相空间重构和确定重构相空间的指定区域V后,计算空间相点的密度分布;(3)计算有载分接开关动作前后空间密度曲线的相关系数ρ;若相关系数ρ&gt;0.8,则判断有载分接开关运行状态正常;若相关系数ρ&lt;0.8,则判断有载分接开关运行状态异常;其所述的变压器有载分接开关运行状态的在线监测方法,以通过对变压器有载分接开关表面振动信号的实时监控,来直接判断变压器有载分接开关的工作状态,便于操作人员及时发现变压器有载分接开关的异常,从而根据异常情况对变压器有载分接开关进行及时地检修,可大大降低变压器有载分接开关的故障损坏率;其中,其所述的嵌入维数m和延迟时间τ的计算过程为:(2‑1)计算振动信号时间序列x(t)的标准差σ;(2‑2)将振动信号时间序列x(t)分成t个不相交的子序列,每个子序列的长度l为l=N/t,此处有t=6;式中,N为振动信号时间序列的均值;(2‑3)对每个子序列,计算三个统计量<img file="FDA0000848411510000021.GIF" wi="312" he="86" />和S<sub>cor</sub>(t),为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>S</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>16</mn></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>2</mn></mrow><mn>5</mn></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>4</mn></munderover><mi>S</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><msub><mi>r</mi><mi>j</mi></msub><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000848411510000022.GIF" wi="597" he="166" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&Delta;</mi><mover><mi>S</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>2</mn></mrow><mn>5</mn></munderover><mi>&Delta;</mi><mi>S</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000848411510000023.GIF" wi="518" he="159" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>S</mi><mrow><mi>c</mi><mi>o</mi><mi>r</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>&Delta;</mi><mover><mi>S</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mo>|</mo><mover><mi>S</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo></mrow>]]></math><img file="FDA0000848411510000024.GIF" wi="533" he="94" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>S</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>r</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>t</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>t</mi></munderover><mo>&lsqb;</mo><msub><mi>C</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>r</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>C</mi><mi>s</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mn>1</mn><mo>,</mo><mi>r</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000848411510000025.GIF" wi="916" he="155" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>S</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>N</mi><mo>,</mo><mi>r</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>t</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>t</mi></munderover><mo>&lsqb;</mo><msub><mi>C</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>N</mi><mo>/</mo><mi>t</mi><mo>,</mo><mi>r</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>C</mi><mi>s</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mn>1</mn><mo>,</mo><mi>N</mi><mo>/</mo><mi>t</mi><mo>,</mo><mi>r</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>,</mo><mi>N</mi><mo>&RightArrow;</mo><mi>&infin;</mi></mrow>]]></math><img file="FDA0000848411510000026.GIF" wi="1478" he="153" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>C</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>N</mi><mo>,</mo><mi>r</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>2</mn><mrow><mi>k</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><munder><mo>&Sigma;</mo><mrow><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>k</mi></mrow></munder><mi>H</mi><mrow><mo>(</mo><mi>r</mi><mo>-</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo><mi>r</mi><mo>&gt;</mo><mn>0</mn></mrow>]]></math><img file="FDA0000848411510000027.GIF" wi="1046" he="159" /></maths>d<sub>ij</sub>=||X(i)‑X(j)||,其中的d<sub>ij</sub>为欧氏距离;<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>r</mi><mo>&GreaterEqual;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>r</mi><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000848411510000028.GIF" wi="430" he="174" /></maths>为Heaviside阶跃函数式中,r为半径,其取值范围为1~4;(2‑4)寻找<img file="FDA0000848411510000029.GIF" wi="110" he="87" />的第一个零点或<img file="FDA00008484115100000210.GIF" wi="151" he="86" />的第一个极小值点,第一个零点或第一个极小值点对应的时间即为延迟时间τ;(2‑5)寻找S<sub>cor</sub>(t)的最小值点,对应的时间为振动信号时间序列x(t)的第一个整体最大值时间窗口τ<sub>w</sub>,计算嵌入维数m,为m=int(τ<sub>w</sub>/τ+1)式中,int为取整函数;其所述的空间相点密度分布的计算过程为:等间隔划分计算空间[2R<sub>a</sub>,3R<sub>a</sub>],此处等间隔分为5份,其中的2R<sub>a</sub>为圆心,3R<sub>a</sub>为圆域的半径;计算重构相空间各个相点距离原点的欧氏距离,统计相点落在各个空间的点数,得到各个划分空间的相点分布密度;其所述的相关系数ρ的计算公式为:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mi>&rho;</mi><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><mover><mi>X</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>Y</mi><mi>i</mi></msub><mo>-</mo><mover><mi>Y</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow><mrow><msqrt><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><mover><mi>X</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><msqrt><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msub><mi>Y</mi><mi>i</mi></msub><mo>-</mo><mover><mi>Y</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></mfrac><mo>.</mo></mrow>]]></math><img file="FDA0000848411510000031.GIF" wi="837" he="334" /></maths>
地址 200122 上海市浦东新区源深路1122号