发明名称 一种容性绝缘设备介质损耗角的计算方法
摘要 本发明公开了监测技术领域的一种容性绝缘设备介质损耗角的计算方法。其技术方案是,首先构建Blackman自卷积窗,并采集待测设备的电流和电压信号,并对电流和电压信号的进行离散采样;其次采用构建的Blackman自卷积窗对离散信号序列进行加窗截断,得到加窗截断后离散信号的离散频谱函数;接着对加窗截断后离散信号的离散频谱函数进行分析,得到主瓣宽度与旁瓣特性;最后采用四谱线插值修正理论求解出电压、电流基波信号相位角进而求出介损角。本发明利用Blackman自卷积窗极好的抑制频谱泄露效应,利用四谱线插值进行频谱校正,解决非同步采样下基于FFT理论的介损角测量误差,测量精度高,硬件要求小。
申请公布号 CN103576002B 申请公布日期 2016.01.20
申请号 CN201310557636.1 申请日期 2013.11.11
申请人 华北电力大学(保定) 发明人 王永强;谢军;律方成
分类号 G01R27/26(2006.01)I 主分类号 G01R27/26(2006.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 陈波
主权项 一种容性绝缘设备介质损耗角的计算方法,其特征在于,所述方法具体包括以下步骤:步骤1:构建Blackman自卷积窗;步骤2:采集待测设备的电流和电压信号,并对电流和电压信号进行离散采样,得到电流和电压信号的离散数字信号序列;步骤3:采用步骤1中构建的Blackman自卷积窗对离散信号序列进行加窗截断,并进行快速傅里叶变换FFT,从而得到加窗截断后离散信号的离散频谱函数;步骤4:在步骤3的基础上,对加窗截断后离散信号的离散频谱函数进行分析,得到主瓣宽度与旁瓣特性;步骤5:通过分析加窗截断信号经FFT后离散频谱,比较并找出基波频点准确位置附近里所有离散谱线幅值最大的四根谱线,采用四谱线插值修正理论求解出电压、电流基波信号相位角进而求出介损角;所述步骤3中的具体过程为:步骤301:对电流、电压离散序列加窗截断后,加窗截断信号为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>h</mi><mo>=</mo><mn>0</mn></mrow><mi>H</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>A</mi><mi>h</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mn>2</mn><mi>&pi;</mi><mfrac><mrow><mi>h</mi><mo>&CenterDot;</mo><msub><mi>f</mi><mn>0</mn></msub></mrow><msub><mi>f</mi><mi>s</mi></msub></mfrac><mi>n</mi><mo>+</mo><msub><mi>&phi;</mi><mi>h</mi></msub><mo>)</mo></mrow><msub><mi>W</mi><mrow><mi>B</mi><mi>P</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000799974740000011.GIF" wi="852" he="149" /></maths>其中,A<sub>h</sub>为h次谐波分量幅值;f<sub>s</sub>为采样频率;φ<sub>h</sub>为h次谐波分量初相角;f<sub>0</sub>为基波频率;H为所能采集的最高次谐波次数;N为离散Blackman自卷积窗的长度;W<sub>BP</sub>(n)为P阶Blackman自卷积窗的频率特性函数;步骤302:根据快速傅里叶变换FFT的基本性质,加窗截断信号经FFT后,离散频谱为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>X</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>h</mi><mo>=</mo><mn>1</mn></mrow><mi>H</mi></munderover><mfrac><mn>1</mn><mrow><mn>2</mn><mi>j</mi></mrow></mfrac><mrow><mo>(</mo><msub><mi>A</mi><mi>h</mi></msub><msup><mi>e</mi><mrow><msub><mi>j&phi;</mi><mi>m</mi></msub></mrow></msup><msub><mi>W</mi><mrow><mi>B</mi><mi>P</mi></mrow></msub><mo>(</mo><mrow><mi>k</mi><mo>-</mo><msub><mi>hk</mi><mn>0</mn></msub></mrow><mo>)</mo><mo>+</mo><msub><mi>A</mi><mi>h</mi></msub><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>j&phi;</mi><mi>h</mi></msub></mrow></msup><msub><mi>W</mi><mrow><mi>B</mi><mi>P</mi></mrow></msub><mo>(</mo><mrow><mi>k</mi><mo>+</mo><msub><mi>hk</mi><mn>0</mn></msub></mrow><mo>)</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000799974740000021.GIF" wi="1140" he="143" /></maths>其中,k<sub>0</sub>=f<sub>0</sub>N/f<sub>s</sub>,为基波频点准确位置;步骤303:根据离散傅里叶变换的基本性质,长度为N的P阶Blackman自卷积窗其离散频谱函数为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>W</mi><mrow><mi>B</mi><mi>P</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>{</mo><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mfrac><mi>k</mi><mi>P</mi></mfrac><mi>&pi;</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mfrac><mrow><mi>k</mi><mi>&pi;</mi></mrow><mi>P</mi></mfrac></mrow></msup><munderover><mo>&Sigma;</mo><mrow><mi>h</mi><mo>=</mo><mn>0</mn></mrow><mn>2</mn></munderover><mfrac><mrow><msup><mrow><mo>(</mo><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>h</mi></msup><msub><mi>a</mi><mi>h</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi><mi>k</mi></mrow><mi>N</mi></mfrac></mrow><mrow><mn>2</mn><mi>s</mi><mi>i</mi><mi>n</mi><mo>&lsqb;</mo><mi>&pi;</mi><mrow><mo>(</mo><mfrac><mi>k</mi><mi>N</mi></mfrac><mo>-</mo><mfrac><mi>h</mi><mi>M</mi></mfrac><mo>)</mo></mrow><mo>&rsqb;</mo><mi>s</mi><mi>i</mi><mi>n</mi><mo>&lsqb;</mo><mi>&pi;</mi><mrow><mo>(</mo><mfrac><mi>k</mi><mi>N</mi></mfrac><mo>+</mo><mfrac><mi>h</mi><mi>M</mi></mfrac><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mfrac><mo>}</mo></mrow><mi>P</mi></msup></mrow>]]></math><img file="FDA0000799974740000022.GIF" wi="1324" he="295" /></maths>其中,M为原Blackman窗的长度。
地址 071003 河北省保定市永华北大街619号