发明名称 一种谐波污染用户的概率潮流评估方法
摘要 本发明提出一种谐波污染用户的概率评估方法,可以更准确地实现谐波干扰用户注入污染的运行评估或预评估。该项技术最大的优点是采用用户基波功率概率分布和谐波特性概率分布参数进行概率型谐波潮流计算,得到该用户对电网的谐波污染的概率分布特征,而不仅仅是某种特殊运行工况下的谐波污染结果,其结论更全面更合理。
申请公布号 CN103778341A 申请公布日期 2014.05.07
申请号 CN201410039925.7 申请日期 2014.01.27
申请人 福州大学 发明人 邵振国;苏文博;张嫣;李兆祥
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 福州元创专利商标代理有限公司 35100 代理人 蔡学俊
主权项 1.一种谐波污染用户的概率潮流评估方法,其特征在于包括下步骤:步骤S01:建立节点注入功率和谐波频谱特性概率模型,生成样本集<maths num="0001"><![CDATA[<math><mrow><mover><mi>P</mi><mo>~</mo></mover><mo>=</mo><mo>{</mo><msub><mi>P</mi><mi>ij</mi></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>}</mo><mo>,</mo><mover><mi>Q</mi><mo>~</mo></mover><mo>=</mo><mo>{</mo><msub><mi>Q</mi><mi>ij</mi></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>}</mo><mo>,</mo><mover><mi>K</mi><mo>~</mo></mover><mo>{</mo><msubsup><mi>K</mi><mi>i</mi><mi>h</mi></msubsup><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>}</mo><mo>,</mo></mrow></math>]]></maths>其中,<img file="FDA0000462990620000012.GIF" wi="164" he="86" />分别表示节点注入有功功率和无功功率的样本集,P<sub>ij</sub>、Q<sub>ij</sub>分别表示第i个节点的第j个样本;<img file="FDA0000462990620000013.GIF" wi="52" he="70" />表示谐波频谱的样本集,<img file="FDA0000462990620000014.GIF" wi="71" he="81" />表示第i个节点的第h次谐波含有率样本;N为样本数,n为节点数,记平衡节点为1;步骤S02:从所述样本集中取一组计算数据<img file="FDA0000462990620000015.GIF" wi="252" he="115" />步骤S03:采用P-Q解耦进行谐波潮流计算:<maths num="0002"><![CDATA[<math><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Delta;P</mi></mtd></mtr><mtr><mtd><mi>&Delta;Q</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>H</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>L</mi></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Delta;&theta;</mi></mtd></mtr><mtr><mtd><mi>&Delta;U</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,ΔU为节点电压幅值偏差量矩阵,Δθ为节点电压相角偏差量矩阵。解之可得式(2);Δθ=H<sup>-1</sup>ΔP(2)ΔU=L<sup>-1</sup>ΔQ上式中的功率偏差量通过式(3)考虑谐波功率对基波功率的修正量;<maths num="0003"><![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><mi>&Delta;P</mi><mo>=</mo><msubsup><mi>P</mi><mi>i</mi><mi>S</mi></msubsup><mo>-</mo><msubsup><mi>P</mi><mi>i</mi><mn>1</mn></msubsup><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>2</mn></mrow><mi>h</mi></munderover><msup><msub><mi>P</mi><mi>i</mi></msub><mi>k</mi></msup></mtd></mtr><mtr><mtd><mi>&Delta;Q</mi><mo>=</mo><msubsup><mi>Q</mi><mi>i</mi><mi>S</mi></msubsup><mo>-</mo><msubsup><mi>Q</mi><mi>i</mi><mn>1</mn></msubsup><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>2</mn></mrow><mi>h</mi></munderover><msubsup><mi>Q</mi><mi>i</mi><mi>k</mi></msubsup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(3)中,<img file="FDA0000462990620000018.GIF" wi="201" he="86" />表示节点i的注入有功功率和无功功率给定值;<img file="FDA0000462990620000019.GIF" wi="197" he="84" />分别表示节点i的k谐波有功功率和无功功率计算值;h表示计算的最高次谐波。在基波潮流解的基础上,谐波功率根据式(4)计算;<maths num="0004"><![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><msup><msub><mi>P</mi><mi>i</mi></msub><mi>k</mi></msup><mo>=</mo><msubsup><mi>R</mi><mi>ii</mi><mi>k</mi></msubsup><msubsup><mi>K</mi><mi>i</mi><mi>k</mi></msubsup><mo>[</mo><msup><mrow><mo>(</mo><msubsup><mi>a</mi><mi>i</mi><mn>1</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msubsup><mi>b</mi><mi>i</mi><mn>1</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo></mtd></mtr><mtr><mtd><msubsup><mi>Q</mi><mi>i</mi><mi>k</mi></msubsup><mo>=</mo><msubsup><mi>X</mi><mi>ii</mi><mi>k</mi></msubsup><msubsup><mi>K</mi><mi>i</mi><mi>k</mi></msubsup><mo>[</mo><msup><mrow><mo>(</mo><msubsup><mi>a</mi><mi>i</mi><mn>1</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msubsup><mi>b</mi><mi>i</mi><mn>1</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,<img file="FDA0000462990620000021.GIF" wi="50" he="83" />和<img file="FDA0000462990620000022.GIF" wi="50" he="82" />分别为节点i每次迭代过程中的基波电流的实部和虚部,其计算如式(5)所示;<img file="FDA0000462990620000023.GIF" wi="70" he="84" />和<img file="FDA0000462990620000024.GIF" wi="74" he="84" />分别为k次谐波导纳矩阵的逆矩阵的对应元素的实部和虚部;<img file="FDA0000462990620000025.GIF" wi="81" he="84" />为节点i的k次谐波电流含有率;谐波功率只在有谐波注入电流的节点其值存在,其他节点的计算值为0;<maths num="0005"><![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><msubsup><mi>a</mi><mi>i</mi><mn>1</mn></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>G</mi><mi>ij</mi></msub><msub><mi>U</mi><mi>j</mi></msub><mi>cos</mi><msub><mi>&theta;</mi><mi>j</mi></msub><mo>-</mo><msub><mi>B</mi><mi>ij</mi></msub><msub><mi>U</mi><mi>j</mi></msub><mi>sin</mi><msub><mi>&theta;</mi><mi>j</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>b</mi><mi>i</mi><mn>1</mn></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>G</mi><mi>ij</mi></msub><msub><mi>U</mi><mi>j</mi></msub><mi>sin</mi><msub><mi>&theta;</mi><mi>j</mi></msub><mo>+</mo><msub><mi>B</mi><mi>ij</mi></msub><msub><mi>U</mi><mi>j</mi></msub><mi>cos</mi><msub><mi>&theta;</mi><mi>j</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>根据式(1)到式(4)进行迭代求解直到前后两次电压差满足收敛条件即可完成当次迭代过程;max(|ΔU<sub>i</sub>|)≤ε    (6)根据谐波源特性计算得到谐波注入电流,利用谐波节点电压方程(7)以及支路电流方程(8)即可求得谐波节点电压和支路电流;<maths num="0006"><![CDATA[<math><mrow><msubsup><mi>U</mi><mi>n</mi><mi>h</mi></msubsup><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>Y</mi><mi>n</mi><mi>h</mi></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msubsup><mi>I</mi><mi>n</mi><mi>h</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msubsup><mi>I</mi><mi>b</mi><mi>h</mi></msubsup><mo>=</mo><msubsup><mi>Y</mi><mi>b</mi><mi>h</mi></msubsup><msup><mi>A</mi><mi>T</mi></msup><msubsup><mi>U</mi><mi>n</mi><mi>h</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,<img file="FDA0000462990620000029.GIF" wi="310" he="84" />表示h次谐波的节点电压向量、节点注入电流向量和支路电流向量;<img file="FDA00004629906200000210.GIF" wi="195" he="84" />表示h次谐波节点导纳矩阵和支路导纳矩阵;A表示节点支路关联矩阵;步骤S04:将节点功率和谐波特性的样本集重复代入式(1)到式(4)的谐波潮流迭代模型,每一组的节点功率和谐波特性都能得到一组对应的谐波电压和谐波电流计算结果,用n表示节点数,N表示样本数,N1表示支路数,<img file="FDA00004629906200000211.GIF" wi="189" he="81" />表示h次谐波节点电压计算结果和支路电流计算结果,则计算结果可表示如下;<maths num="0008"><![CDATA[<math><mfenced open='' close='' separators=''><mtable><mtr><mtd><msubsup><mover><mi>U</mi><mo>~</mo></mover><mi>n</mi><mi>h</mi></msubsup><mo>=</mo><mo>{</mo><msubsup><mi>U</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>h</mi></msubsup><mo>,</mo><mi>i</mi><mo>=</mo><mn>2,3</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>}</mo></mtd></mtr><mtr><mtd><msubsup><mover><mi>I</mi><mo>~</mo></mover><mi>b</mi><mi>h</mi></msubsup><mo>=</mo><mo>{</mo><msubsup><mi>I</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>h</mi></msubsup><mo>,</mo><mi>i</mi><mo>=</mo><mn>2,3</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mn>1</mn><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>}</mo></mtd></mtr></mtable><mo>.</mo></mfenced></math>]]></maths>
地址 350108 福建省福州市闽侯县上街镇大学城学园路2号福州大学新区