主权项 |
用于变压器局部放电定位的超声波接收阵列的定向方法,其特征在于,包括以下步骤: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>σ</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>⊗</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>⊗</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>⊗</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>·</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>⊗</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>⊗</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>Σ</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" />为超声阵列的流行矢量。 |