发明名称 用于变压器局部放电定位的超声波接收阵列的定向方法
摘要 本发明公开了一种用于变压器局部放电定位的超声波接收阵列的定向方法,该方法针对沿着平面阵阵列边沿或者中心线布置的、呈L形或形的超声波接收阵列。本发明将宽带聚焦和阵列扩展相结合,首先对宽带的超声阵列接收信号进行聚焦操作,将局部放电产生的宽带超声信号按时间分段,利用无噪声数据来计算聚焦矩阵,将信号聚焦到同一频率段;再以四阶累积量算法对阵列进行虚拟扩展,虚拟扩展后的阵列的数据协方差矩阵用未扩展的原阵列按顺序形成的四阶累积量矩阵来代替;最后,对虚拟扩展后的阵列利用多重信号分类算法估计局部放电的来波方向。
申请公布号 CN103995221A 申请公布日期 2014.08.20
申请号 CN201410235567.7 申请日期 2014.05.29
申请人 西安交通大学 发明人 罗勇芬;辛晓虎;李彦明
分类号 G01R31/12(2006.01)I 主分类号 G01R31/12(2006.01)I
代理机构 西安通大专利代理有限责任公司 61200 代理人 陆万寿
主权项 用于变压器局部放电定位的超声波接收阵列的定向方法,其特征在于,包括以下步骤:1)利用超声波阵列传感器接收变压器局部放电产生的超声波阵列信号,将接收到的超声波阵列信号分为J个时间段X<sub>0</sub>,X<sub>1</sub>,...,X<sub>J‑1</sub>,计算每个时间段无噪声数据的协方差矩阵,其计算公式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>P</mi><mi>j</mi></msub><mo>=</mo><msub><mi>R</mi><mi>j</mi></msub><mo>-</mo><msubsup><mi>&sigma;</mi><mi>j</mi><mn>2</mn></msubsup><mi>I</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000513440160000011.GIF" wi="1203" he="88" /></maths>式中:j=0,1,...J‑1;R<sub>j</sub>为X<sub>j</sub>的协方差矩阵;<img file="FDA0000513440160000012.GIF" wi="56" he="78" />为R<sub>j</sub>小特征值的和;I为单位矩阵;2)选定分段阵列信号的聚焦矩阵,其计算公式如下:<img file="FDA0000513440160000013.GIF" wi="204" he="80" />其中,j=1,...J‑1  (2)式中:<img file="FDA0000513440160000014.GIF" wi="54" he="67" />和<img file="FDA0000513440160000015.GIF" wi="54" he="71" />分别是P<sub>0</sub>和P<sub>j</sub>的特征矢量;3)用步骤2)得到的聚焦矩阵T<sub>j</sub>左乘X<sub>j</sub>可得聚焦后的每段窄带阵列信号Y<sub>j</sub>(f<sub>0</sub>),其计算公式如下:Y<sub>j</sub>(f<sub>0</sub>)=T<sub>j</sub>X<sub>j</sub> j=1,...J‑1  (3)式中:f<sub>0</sub>为阵列信号的参考频率;4)计算聚焦后每段窄带阵列信号Y<sub>j</sub>(f<sub>0</sub>)的四阶累积量矩阵R<sub>4,j</sub>,并计算其平均值Rcum,其计算公式如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>R</mi><mrow><mn>4</mn><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mi>E</mi><mo>{</mo><msup><mrow><mrow><mo>(</mo><msub><mi>Y</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>&CircleTimes;</mo><msub><mi>Y</mi><mi>j</mi></msub><msup><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>*</mo></msup><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>Y</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>&CircleTimes;</mo><msub><mi>Y</mi><mi>j</mi></msub><msup><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>*</mo></msup><mo>)</mo></mrow></mrow><mi>H</mi></msup></mtd></mtr><mtr><mtd><mo>-</mo><mi>E</mi><mo>{</mo><mrow><mo>(</mo><msub><mi>Y</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>&CircleTimes;</mo><msub><mi>Y</mi><mi>j</mi></msub><msup><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>*</mo></msup><mo>)</mo></mrow><mo>}</mo><mo>&CenterDot;</mo><mi>E</mi><mo>{</mo><msup><mrow><mo>(</mo><msub><mi>Y</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>&CircleTimes;</mo><msub><mi>Y</mi><mi>j</mi></msub><msup><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>*</mo></msup><mo>)</mo></mrow><mi>H</mi></msup><mo>}</mo></mtd></mtr><mtr><mtd><mo>-</mo><mi>E</mi><mo>{</mo><mrow><mo>(</mo><msub><mi>Y</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>)</mo></mrow><msub><mi>Y</mi><mi>j</mi></msub><msup><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>H</mi></msup><mo>)</mo></mrow><mo>}</mo><mo>&CircleTimes;</mo><mi>E</mi><mo>{</mo><msup><mrow><mo>(</mo><msub><mi>Y</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>)</mo></mrow><msub><mi>Y</mi><mi>j</mi></msub><msup><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>H</mi></msup><mo>)</mo></mrow><mo>*</mo></msup><mo>}</mo></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000513440160000028.GIF" wi="1454" he="240" /></maths>式中:E为求矩阵互相关;<img file="FDA0000513440160000029.GIF" wi="65" he="62" />为求矩阵的克罗内克积;*为求矩阵转置;H为求矩阵共轭转置;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>Rcum</mi><mo>=</mo><mfrac><mn>1</mn><mi>J</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>J</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>R</mi><mrow><mn>4</mn><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000513440160000021.GIF" wi="1277" he="148" /></maths>5)对Rcum进行特征分解,得到噪声子空间En,并按式(6)的谱估计公式进行谱峰搜索,得到极大点所对应的来波方向<img file="FDA0000513440160000022.GIF" wi="130" he="73" />即为局部放电信号来波方向,其计算公式如下:<img file="FDA0000513440160000023.GIF" wi="1260" he="141" />式中:<img file="FDA0000513440160000024.GIF" wi="173" he="70" />代表噪声信号在<img file="FDA0000513440160000025.GIF" wi="114" he="64" />方向的空间谱能量的倒数;θ和<img file="FDA0000513440160000026.GIF" wi="52" he="67" />分别代表来波方向的俯仰角和方位角;<img file="FDA0000513440160000027.GIF" wi="167" he="69" />为超声阵列的流行矢量。
地址 710049 陕西省西安市咸宁西路28号