发明名称 一种基于单频电流传输特性的分布式谐波源辨识方法
摘要 本发明涉及一种基于单频电流传输特性的分布式谐波源辨识方法,其包括以下步骤:1)建立电流谐波源辨识模型;2)导出约束条件;3)构建电流谐波观测矩阵并进行预处理;4)设置用电网络内谐波源最大辨识次数T以及非高斯性阈值F,令k=0;5)令k=k+1,采用FastICA方法求解混矩阵W<sub>h</sub>;6)计算W<sub>h</sub>的逆矩阵W<sub>h</sub><sup>-1</sup>,依据约束条件调整W<sub>h</sub><sup>-1</sup>的相位及行顺序,记为P<sub>h</sub>;若P<sub>h</sub>满足约束条件,则执行步骤7;否则,返回执行步骤5;7)计算P<sub>h</sub>的逆矩阵,更新解混矩阵W<sub>h</sub>=P<sub>h</sub><sup>-1</sup>,并计算辨识结果的非高斯性之和<img file="DDA0000595940820000011.GIF" wi="39" he="42" />;若<img file="DDA0000595940820000012.GIF" wi="35" he="40" />不小于阈值F,输出解混矩阵W<sub>h</sub>和计算谐波源辨识结果,并分析主要谐波源位置,结束迭代;否则执行步骤8;8)当k≤T时,返回步骤5;当k>T时,输出解混矩阵W<sub>h</sub>和非高斯性阈值,并进一步分析网络内主要谐波源的位置。
申请公布号 CN104280612A 申请公布日期 2015.01.14
申请号 CN201410588860.1 申请日期 2014.10.28
申请人 清华大学 发明人 王雪;吴江伟;刘佑达;孙欣尧
分类号 G01R23/16(2006.01)I 主分类号 G01R23/16(2006.01)I
代理机构 北京纪凯知识产权代理有限公司 11245 代理人 徐宁;刘美丽
主权项 一种基于单频电流传输特性的分布式谐波源辨识方法,其包括以下步骤:1)建立电流谐波源辨识模型:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>I</mi><mi>h</mi><mi>S</mi></msubsup><mo>=</mo><msubsup><mi>P</mi><mi>h</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msubsup><mi>I</mi><mi>h</mi><mi>M</mi></msubsup><mo>=</mo><msub><mi>W</mi><mi>h</mi></msub><msubsup><mi>I</mi><mi>h</mi><mi>M</mi></msubsup></mrow>]]></math><img file="FDA0000595940790000011.GIF" wi="400" he="86" /></maths>式中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>I</mi><mi>h</mi><mi>M</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>I</mi><mn>1</mn><mi>m</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>I</mi><mn>2</mn><mi>m</mi></msubsup></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msubsup><mi>I</mi><mi>k</mi><mi>m</mi></msubsup></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msubsup><mi>I</mi><mi>N</mi><mi>m</mi></msubsup></mtd></mtr></mtable></mfenced><mo>,</mo><msubsup><mi>I</mi><mi>h</mi><mi>S</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>I</mi><mrow><mn>1</mn><mi>s</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>I</mi><mrow><mn>2</mn><mi>s</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>I</mi><mi>ks</mi></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>I</mi><mi>Ns</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo><msub><mi>P</mi><mi>h</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>p</mi><mn>11</mn></msub></mtd><mtd><msub><mi>p</mi><mn>12</mn></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>p</mi><mrow><mn>1</mn><mi>N</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>p</mi><mn>21</mn></msub></mtd><mtd><msub><mi>p</mi><mn>22</mn></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>p</mi><mrow><mn>2</mn><mi>N</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>p</mi><mrow><mi>N</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>p</mi><mrow><mi>N</mi><mn>2</mn></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>p</mi><mi>NN</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000595940790000012.GIF" wi="1001" he="402" /></maths><img file="FDA0000595940790000013.GIF" wi="65" he="83" />表示第k个支路上电能信息监测节点测得的第h次电流谐波序列,I<sub>ks</sub>表示第k个监测支路向用电网络注入的第h次电流谐波序列,p<sub>ij</sub>表示从第j个支路对第i个支路的谐波影响系数,W<sub>h</sub>为谐波源辨识的解混矩阵;2)导出电流谐波源辨识模型在低压用电网络内电流谐波源辨识的约束条件;3)依据用电网络内各支路测量到的单频电流谐波信息构建电流谐波观测信号矩阵,对观测信号矩阵进行预处理去除其中的缓变分量,并记为X;4)设置用电网络内谐波源最大辨识次数T以及非高斯性阈值F,令初始辨识迭代次数k=0;5)令k=k+1,采用FastICA求解电流谐波源辨识模型中的解混矩阵W<sub>h</sub>;6)计算解混矩阵W<sub>h</sub>的逆矩阵W<sub>h</sub><sup>‑1</sup>,依据约束条件调整W<sub>h</sub><sup>‑1</sup>的相位及行顺序,记为P<sub>h</sub>,若P<sub>h</sub>满足约束条件,则进入步骤7);否则,返回执行步骤5);7)计算P<sub>h</sub>的逆矩阵,更新解混矩阵W<sub>h</sub>=P<sub>h</sub><sup>‑1</sup>,依据电流谐波源辨识模型分析各监测支路监测对象向用电网络注入的统计独立分量,并计算辨识结果的非高斯性之和<img file="FDA0000595940790000014.GIF" wi="78" he="59" />若<img file="FDA0000595940790000015.GIF" wi="50" he="62" />不小于非高斯性阈值F,输出解混矩阵W<sub>h</sub>,由观测信号矩阵X取代<img file="FDA0000595940790000016.GIF" wi="70" he="83" />估计各监测支路向用电网络的谐波注入量<img file="FDA0000595940790000017.GIF" wi="86" he="83" />分析各支路电流谐波的方差,结合各监测支路电流谐波最大值进而确定主要谐波扰动源所处的支路位置,结束辨识迭代过程;否则进入步骤8);8)当k≤T时,返回步骤5);当k>T时,输出解混矩阵W<sub>h</sub>和当前的非高斯性阈值,并进一步分析用电网络内主要谐波源的所处的支路,辨识出各监测支路电流谐波信号后,分析各电流谐波的方差,结合各监测支路电流谐波的最大值进而确定主要谐波源的位置。
地址 100084 北京市海淀区100084信箱82分箱清华大学专利办公室