发明名称 一种多谐波源辨识方法
摘要 本发明涉及一种新的多谐波源辨识方法,它可为多谐波源系统中扰动源的定位、谐波的治理以及谐波的奖惩实施等提供理论依据,具有广阔的应用前景和良好的社会经济效益。它的步骤为:(1)定义谐波影响指标;(2)求取谐波影响指标;(3)数据的采集和预处理;(4)进行数据的选择和分析。
申请公布号 CN102323494B 申请公布日期 2013.05.15
申请号 CN201110135481.3 申请日期 2011.05.24
申请人 山东电力研究院;国家电网公司 发明人 臧宏志;孙媛媛;杨永军;俞晓冬;李广磊;尹志明
分类号 G01R31/00(2006.01)I;G01R23/16(2006.01)I 主分类号 G01R31/00(2006.01)I
代理机构 济南圣达知识产权代理有限公司 37221 代理人 张勇
主权项 1.一种多谐波源辨识方法,其特征是,它的过程为:(1)定义谐波责任指标假设系统节点X处的谐波问题是由n个“可疑谐波源负荷”共同作用产生,对于h次谐波,测得n个“可疑谐波源负荷”节点的谐波电流I<sub>h1</sub>,I<sub>h2</sub>,I<sub>h3</sub>……I<sub>hn</sub>;节点X处的h次谐波电压<img file="FDA00002952308000011.GIF" wi="74" he="82" />根据系统的谐波阻抗[Z<sub>h</sub>]计算,将节点X处的h次谐波电压分为两部分:一部分由谐波源负荷i产生<img file="FDA00002952308000012.GIF" wi="102" he="82" />另一部分由其他谐波源负荷产生的背景谐波<img file="FDA00002952308000013.GIF" wi="96" he="79" /><img file="FDA00002952308000014.GIF" wi="1874" he="179" />其中,I<sub>hi</sub>为谐波源负荷节点i的h次谐波电流;<img file="FDA00002952308000015.GIF" wi="81" he="78" />为节点i与节点X之间的谐波阻抗;<img file="FDA00002952308000016.GIF" wi="623" he="82" />为除了谐波源负荷节点i之外的其他谐波源负荷节点与节点X之间的谐波阻抗;谐波源负荷i在X处产生的h次谐波电压<img file="FDA00002952308000017.GIF" wi="79" he="82" />对节点X处的h次谐波电压<img file="FDA00002952308000018.GIF" wi="77" he="82" />的贡献通过<img file="FDA00002952308000019.GIF" wi="73" he="80" />在<img file="FDA000029523080000110.GIF" wi="74" he="80" />上的投影表示,谐波源负荷i对X处的h次谐波电压的谐波责任指标<img file="FDA000029523080000111.GIF" wi="124" he="77" />定义为:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>HF</mi><mi>i</mi><mi>X</mi></msubsup><mo>=</mo><mfrac><mrow><mo>|</mo><msubsup><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>hi</mi><mi>X</mi></msubsup><mo>|</mo></mrow><mrow><mo>|</mo><msubsup><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>h</mi><mi>X</mi></msubsup><mo>|</mo></mrow></mfrac><mi>cos</mi><mi>&beta;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中β为<img file="FDA000029523080000113.GIF" wi="216" he="79" />的夹角;(2)求取谐波责任指标对电力系统中关键节点的实际波形进行采样,获取谐波数据,进而计算出谐波责任指标<img file="FDA000029523080000114.GIF" wi="139" he="79" />建立如下方程:Y=Sθ+ε    (3)其中,Y为由<img file="FDA000029523080000115.GIF" wi="177" he="98" />构造的矩阵参数,<maths num="0002"><![CDATA[<math><mrow><mi>Y</mi><mo>=</mo><mo>[</mo><msup><mrow><mo>|</mo><msubsup><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>h</mi><mi>X</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><msup><mrow><mo>|</mo><msubsup><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>h</mi><mi>X</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msup><mrow><mo>|</mo><msubsup><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>h</mi><mi>X</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mi>m</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>,</mo></mrow></math>]]></maths><img file="FDA000029523080000117.GIF" wi="178" he="97" />(k=1,2,…m)为各个采样时刻t<sub>1</sub>,t<sub>2</sub>,……,t<sub>m</sub>的节点X处的h次谐波电压<img file="FDA000029523080000118.GIF" wi="78" he="82" />的幅值;S为由|I<sub>hi</sub>(t<sub>k</sub>)|构造的矩阵参数,<img file="FDA000029523080000212.GIF" wi="808" he="390" />(k=1,2,…m)为各采样时刻节点i的h次谐波电流的幅值;<maths num="0003"><![CDATA[<math><mrow><mi>&theta;</mi><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mrow><mo>|</mo><msub><mover><mi>E</mi><mo>&RightArrow;</mo></mover><mi>h</mi></msub><mo>|</mo></mrow><mn>2</mn></msup></mtd><mtd><mo>-</mo><mn>2</mn><mo>|</mo><msubsup><mi>Z</mi><mi>hi</mi><mi>X</mi></msubsup><mo>|</mo><mo>|</mo><msub><mover><mi>E</mi><mo>&RightArrow;</mo></mover><mi>h</mi></msub><mo>|</mo><mi>COS</mi><mrow><mo>(</mo><msub><mi>&alpha;</mi><mi>eq</mi></msub><mo>)</mo></mrow></mtd><mtd><msup><mrow><mo>|</mo><msubsup><mi>Z</mi><mi>hi</mi><mi>X</mi></msubsup><mo>|</mo></mrow><mn>2</mn></msup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths>α<sub>eq</sub>为<img file="FDA00002952308000023.GIF" wi="215" he="79" />的夹角α的平均值,<img file="FDA00002952308000024.GIF" wi="87" he="97" />为<img file="FDA00002952308000025.GIF" wi="64" he="84" />的幅值,ε是用α<sub>eq</sub>来代替α所产生的误差;<img file="FDA00002952308000026.GIF" wi="80" he="81" />为节点i和X之间的谐波阻抗;在上述方程中,当进行相关性分析时,假设背景谐波<img file="FDA00002952308000027.GIF" wi="61" he="84" />保持不变,使用最小二乘法来估计出未知参数θ,代入(2)得,<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>HF</mi><mi>i</mi><mi>X</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mrow><mo>|</mo><msubsup><mi>Z</mi><mi>hi</mi><mi>X</mi></msubsup><mo>|</mo></mrow><mn>2</mn></msup><msup><mrow><mo>|</mo><msub><mi>I</mi><mi>hi</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><mo>|</mo><msubsup><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>h</mi><mi>X</mi></msubsup><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>|</mo><mo>-</mo><msup><mrow><mo>|</mo><msub><mover><mi>E</mi><mo>&RightArrow;</mo></mover><mi>h</mi></msub><mo>|</mo></mrow><mn>2</mn></msup></mrow><mrow><mn>2</mn><msup><mrow><mo>|</mo><msubsup><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>h</mi><mi>X</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>&times;</mo><mn>100</mn><mo>%</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths><img file="FDA00002952308000029.GIF" wi="202" he="79" />表示的是t<sub>k</sub>时刻谐波源负荷i对节点X的h次谐波电压的谐波责任指标,而一段特定时间内的平均谐波责任指标为:<maths num="0005"><![CDATA[<math><mrow><msubsup><mi>HF</mi><mi>i</mi><mi>X</mi></msubsup><mo>=</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>+</mo><mfrac><mn>1</mn><mi>m</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mfrac><mrow><msup><mrow><mo>-</mo><mo>|</mo><msub><mover><mi>E</mi><mo>&RightArrow;</mo></mover><mi>h</mi></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>|</mo><msubsup><mi>Z</mi><mi>hi</mi><mi>X</mi></msubsup><mo>|</mo></mrow><mn>2</mn></msup><msup><mrow><mo>|</mo><msub><mi>I</mi><mi>hi</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow><mrow><mn>2</mn><msup><mrow><mo>|</mo><msubsup><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>h</mi><mi>X</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mo>&times;</mo><mn>100</mn><mo>%</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>m为测量的总次数,即只要测得可疑谐波源负荷i的h次谐波电流I<sub>hi</sub>和节点X处的h次谐波电压<img file="FDA000029523080000211.GIF" wi="78" he="80" />,根据上述方程计算出负荷i的谐波责任;(3)数据的采集和预处理在电力系统的主要谐波源负荷点安装测量装置,同时采集关键节点的单相电压和电流信号,并去除含有明显噪声和暂态干扰的信号;(4)进行数据的选择和分析找出只有一个“可疑谐波源负荷”的谐波电流有较大变化的时间段,取出满足此要求的时间段的节点电压和电流的数据进行相关性分析,在这个时间段内负荷i的h次谐波电流变化较大而其余负荷的h次谐波电流变化较小,对所选择出的时间段的实测电压和电流数据进行傅里叶分解,得出各变量的基频分量和各次谐波分量,根据公式(3)估计出未知参数θ,再根据公式(4)和(5)计算以百分比形式表示的可疑谐波源负荷的谐波责任指标<img file="FDA00002952308000031.GIF" wi="150" he="77" />
地址 250002 山东省济南市市中区二环南路500号