发明名称 IEC框架下的邻近谐波间谐波分离测量方法
摘要 本发明公开了一种IEC框架下的邻近谐波间谐波分离测量方法,包括以下步骤,步骤一,对电网信号进行离散采样,步骤二,对采样值按照IEC标准进行10周波加Hanning窗DFT/FFT频谱变换,得到频谱;步骤三,对步骤二得到的频谱进行乘旋转相位因子处理得到新频谱;步骤四,对新频谱求取相邻谱线的矢量和,以抵消其余分量对其的旁瓣干扰,用谱线相消插值法求解对应的频率、频率偏移值和幅值相位;步骤五,在频域上将基波分量和其余分量在关注谱线处的频谱泄漏值剔除;步骤六,求解谐波频谱值得到h次谐波参数;求解h次谐波邻近的间谐波参数。解决了有限的非同步采样数据时无法精确测量两邻近的谐波、间谐波参数的这一难题。
申请公布号 CN104122443B 申请公布日期 2017.02.15
申请号 CN201410379753.8 申请日期 2014.08.04
申请人 国家电网公司;江苏省电力公司;江苏省电力公司电力科学研究院 发明人 袁晓冬;李群;陈兵;杨洪耕;史明明
分类号 G01R23/16(2006.01)I 主分类号 G01R23/16(2006.01)I
代理机构 南京纵横知识产权代理有限公司 32224 代理人 董建林
主权项 一种IEC框架下的邻近谐波间谐波分离测量方法,其特征在于:包括以下步骤,步骤一,对电网信号进行离散采样,采样值为g(nΔt)<sub>N</sub>为公式(1),<img file="FDA0001123122640000011.GIF" wi="1589" he="149" />其中,Δt为单位采样间隔,M为采样分量个数,f<sub>i</sub>、A<sub>i</sub>、<img file="FDA0001123122640000012.GIF" wi="46" he="54" />为采样分量的频率、幅值和相位,n=0、1、…、N‑1,N为采样窗口长度;步骤二,对步骤一中的采样值g(nΔt)<sub>N</sub>按照IEC标准进行10周波加Hanning窗DFT/FFT频谱变换,得到频谱G(k),<img file="FDA0001123122640000013.GIF" wi="1294" he="143" />其中,β<sub>i</sub>=f<sub>i</sub>/Δf为归一化频率,Δf为对应采样窗口的长度为N的频率分辨率,W(*)为Hanning窗谱函数,G(k)对应谱线k处的频谱值;步骤三,对步骤二得到的频谱G(k)进行乘旋转相位因子处理得到新频谱<img file="FDA0001123122640000014.GIF" wi="142" he="71" />步骤四,对新频谱<img file="FDA0001123122640000015.GIF" wi="120" he="78" />求取相邻谱线的矢量和,以抵消其余分量对其的旁瓣干扰,用谱线相消插值法求解基波分量和其余分量各自对应的频率f<sub>i</sub>、频率偏移值δ<sub>i</sub>和幅值相位<img file="FDA0001123122640000016.GIF" wi="152" he="71" />所述其余分量为不包含邻近谐波间谐波之外的分量;步骤五,在频域上将基波分量和其余分量在关注谱线k<sub>i</sub>'处的频谱泄漏值<img file="FDA0001123122640000017.GIF" wi="150" he="78" />剔除获得频谱<img file="FDA0001123122640000018.GIF" wi="563" he="79" />k<sub>i</sub>'为对应频率分量f<sub>i</sub>'的最大谱线序号;步骤六,根据基波频率偏移值δ<sub>1</sub>计算对应谐波频率偏移值δ<sub>h</sub>,利用三根和 五根谱线相消插值法在频谱<img file="FDA0001123122640000021.GIF" wi="147" he="79" />上建立间谐波频率偏移值δ<sub>i</sub>'相等方程式,求解谐波频谱值<img file="FDA0001123122640000022.GIF" wi="99" he="78" />得到h次谐波参数;频域上分离出谐波参数后,利用步骤四中的谱线相消插值法求解h次谐波邻近的间谐波参数,其中h取值为与频率f<sub>i</sub>'的间谐波邻近的谐波次数;求解谐波频谱值<img file="FDA0001123122640000023.GIF" wi="95" he="79" />得到h次谐波参数的过程为,(a)根据基波频率偏移值δ<sub>1</sub>计算对应谐波频率偏移值δ<sub>h</sub>,δ<sub>h</sub>=β<sub>h</sub>‑k<sub>h</sub>    (14)其中,β<sub>h</sub>=f<sub>h</sub>/Δf=h(T+δ<sub>1</sub>),β<sub>h</sub>为第h次谐波对应的归一化频率,f<sub>h</sub>为谐波频率,T为分析周波个数,k<sub>h</sub>为对应谐波频率为f<sub>h</sub>的最大谱线序号;(b)在频谱<img file="FDA0001123122640000024.GIF" wi="147" he="78" />上取间谐波峰谱线k<sub>i</sub>'及其左右各两根谱线k<sub>i</sub>'±1、k<sub>i</sub>'±2,每根谱线对应的频谱包含谐波频谱和间谐波频谱,间谐波频谱<img file="FDA0001123122640000025.GIF" wi="229" he="78" />为,<img file="FDA0001123122640000026.GIF" wi="1315" he="80" />式中<img file="FDA0001123122640000027.GIF" wi="527" he="79" />为h次谐波峰谱线频谱,r=0、±1、±2,d=k<sub>i</sub>'‑k<sub>h</sub>为间谐波、谐波峰谱线距离,<img file="FDA0001123122640000028.GIF" wi="1154" he="159" />(c)利用三谱线比值<sub>3</sub>r<sub>i</sub>和五谱线比值<sub>5</sub>r<sub>i</sub>中的间谐波频率偏移值δ<sub>i</sub>'相等得,<img file="FDA0001123122640000029.GIF" wi="682" he="151" />乘旋转相位因子后间谐波频谱相位值保持不变,公式(16)简化为:<img file="FDA00011231226400000210.GIF" wi="922" he="159" />式中,<img file="FDA00011231226400000211.GIF" wi="906" he="79" /><img file="FDA00011231226400000212.GIF" wi="906" he="78" /><img file="FDA0001123122640000031.GIF" wi="1193" he="79" /><img file="FDA0001123122640000032.GIF" wi="1190" he="79" />α<sub>11</sub>=5α<sub>d‑1</sub>+6α<sub>d</sub>+α<sub>d+1</sub>,α<sub>12</sub>=α<sub>d‑1</sub>+6α<sub>d</sub>+5α<sub>d+1</sub>,α<sub>21</sub>=α<sub>d‑2</sub>+3α<sub>d‑1</sub>+3α<sub>d</sub>+α<sub>d+1</sub>,α<sub>22</sub>=α<sub>d‑1</sub>+3α<sub>d</sub>+3α<sub>d+1</sub>+α<sub>d+2</sub>,(d)求解公式(17),得h次谐波峰谱线频谱<img file="FDA0001123122640000033.GIF" wi="123" he="78" /><img file="FDA0001123122640000034.GIF" wi="1134" he="151" />(e)分离出的h次谐波的幅值相位<img file="FDA0001123122640000035.GIF" wi="176" he="74" /><img file="FDA0001123122640000036.GIF" wi="917" he="159" />
地址 100761 北京市西城区西长安街86号