发明名称 一种含逆变型分布式电源的配电网故障计算通用方法
摘要 一种含逆变型分布式电源的配电网故障计算通用方法,它首先针对逆变型分布式电源的短路电流特性进行建立次暂态及稳态的短路等效模型,然后进行故障计算:步骤1,输入系统参数,形成系统导纳矩阵;步骤2:将滤波电抗并入分布式电源导纳矩阵;步骤3:形成次暂态电流的三相分量迭代格式;步骤4:计算稳态模型的内电势正序分量;步骤5:形成稳态电流的三相分量迭代格式;步骤6:统一求解方程得到次暂态故障电压三相分量和稳态故障电压三相分量;步骤7:判断方程<img file="DDA0001122361830000011.GIF" wi="322" he="71" />是否收敛,如果收敛则迭代结束,否则返回步骤3继续迭代。本发明提出了一种相序混合通用计算方法,可适用于含逆变型分布式电源的配电网故障计算。
申请公布号 CN106291255A 申请公布日期 2017.01.04
申请号 CN201610859391.1 申请日期 2016.09.28
申请人 国网山东省电力公司电力科学研究院;国家电网公司 发明人 苏建军;张林利;李立生;邵志敏;孙勇;张世栋;李建修;刘合金;樊迪;王昕;孙树敏;李玉敦;李广磊
分类号 G01R31/08(2006.01)I 主分类号 G01R31/08(2006.01)I
代理机构 济南诚智商标专利事务所有限公司 37105 代理人 李修杰
主权项 一种含逆变型分布式电源的配电网故障计算通用方法,其特征是,包括以下过程:一、对逆变型分布式电源进行统一化建模(1)建立短路电流的次暂态模型定义逆变器的虚拟内电势E,因此短路电流次暂态值<img file="FDA0001122361800000015.GIF" wi="43" he="55" />的模型用公式(1)描述:<maths num="0001"><math><![CDATA[<mrow><msup><mover><mi>I</mi><mo>&CenterDot;</mo></mover><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup><mo>=</mo><mfrac><mrow><msub><mover><mi>E</mi><mo>&CenterDot;</mo></mover><mn>0</mn></msub><mo>-</mo><msub><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mi>f</mi></msub></mrow><mrow><mi>R</mi><mo>+</mo><mi>j</mi><mi>X</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122361800000011.GIF" wi="1165" he="148" /></maths>式中,<img file="FDA0001122361800000012.GIF" wi="53" he="71" />表示短路前的逆变器内电势正序分量,<img file="FDA0001122361800000013.GIF" wi="69" he="71" />为逆变器并网短路正序电压,R表示逆变器出口电阻,X为逆变器滤波电抗,j为虚数单位;(2)建立短路电流的稳态模型逆变型分布式电源采用恒功率控制,按照给定的功率指令输出,在有低电压穿越控制要求时,逆变器按照并网技术要求输出设定的无功功率,在逆变器并网电压小于90%时,低电压穿越控制下的无功功率指令为:Q<sub>ref1</sub>=K(U<sub>f</sub>‑U<sub>0</sub>)U<sub>f</sub> (K≥2)   (2)式中,U<sub>0</sub>为故障前逆变器并网正序电压,U<sub>f</sub>为并网短路电压正序分量,K为无功支撑的比例系数;该阶段分布式电源的短路计算模型采用恒功率模型表示为:<maths num="0002"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>P</mi><mo>=</mo><msub><mi>P</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>,</mo><mi>Q</mi><mo>=</mo><msub><mi>Q</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mrow><mo>(</mo><msub><mi>U</mi><mi>f</mi></msub><mo>&gt;</mo><mn>0.9</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>P</mi><mo>=</mo><msub><mi>P</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mo>,</mo><mi>Q</mi><mo>=</mo><msub><mi>Q</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mn>0.2</mn><mo>&lt;</mo><msub><mi>U</mi><mi>f</mi></msub><mo>&lt;</mo><mn>0.9</mn><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122361800000014.GIF" wi="1341" he="158" /></maths>式中,P表示逆变器输出有功功率,Q表示逆变器输出无功功率,P<sub>ref</sub>为有功功率参考值,Q<sub>ref</sub>与Q<sub>ref1</sub>为不同电压情况下的无功功率参考值;短路电流稳态值<img file="FDA0001122361800000016.GIF" wi="27" he="58" />的模型用公式(4)描述:<maths num="0003"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>E</mi><mo>=</mo><mfrac><mrow><mi>Q</mi><mi>X</mi></mrow><msub><mi>U</mi><mi>f</mi></msub></mfrac><mo>+</mo><mfrac><mrow><mi>P</mi><mi>R</mi></mrow><msub><mi>U</mi><mi>f</mi></msub></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mi>&delta;</mi><mo>=</mo><mfrac><mrow><mi>P</mi><mi>X</mi></mrow><msub><mi>U</mi><mi>f</mi></msub></mfrac><mo>-</mo><mfrac><mrow><mi>Q</mi><mi>R</mi></mrow><msub><mi>U</mi><mi>f</mi></msub></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mover><mi>I</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mfrac><mrow><mover><mi>E</mi><mo>&CenterDot;</mo></mover><mo>-</mo><msub><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mi>f</mi></msub></mrow><mrow><mi>R</mi><mo>+</mo><mi>j</mi><mi>X</mi></mrow></mfrac></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122361800000021.GIF" wi="1517" he="476" /></maths>式中,<img file="FDA0001122361800000022.GIF" wi="45" he="54" />表示短路后的内电势,δ为<img file="FDA0001122361800000023.GIF" wi="38" he="54" />的相角,R表示逆变器出口电阻,X为逆变器滤波电抗;二、故障计算过程步骤1,输入系统参数,形成系统导纳矩阵Y<sub>Sys</sub>;步骤2:将滤波电抗并入分布式电源导纳矩阵:<maths num="0004"><math><![CDATA[<mrow><msub><mi>Y</mi><mrow><mi>a</mi><mi>b</mi><mi>c</mi></mrow></msub><mo>=</mo><msub><mi>Y</mi><mrow><mi>S</mi><mi>y</mi><mi>s</mi></mrow></msub><mo>+</mo><mfrac><mn>1</mn><mrow><mi>R</mi><mo>+</mo><mi>j</mi><mi>X</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122361800000024.GIF" wi="1294" he="126" /></maths>其中R与X分别为逆变器出口电阻和电抗;步骤3:形成次暂态电流的三相分量迭代格式:<maths num="0005"><math><![CDATA[<mrow><msub><mi>Y</mi><mrow><mi>a</mi><mi>b</mi><mi>c</mi></mrow></msub><msubsup><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mrow><mi>f</mi><mi>a</mi><mi>b</mi><mi>c</mi><mn>0</mn></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mfrac><msubsup><mover><mi>E</mi><mo>&CenterDot;</mo></mover><mrow><mi>a</mi><mi>b</mi><mi>c</mi></mrow><mn>0</mn></msubsup><mrow><mi>R</mi><mo>+</mo><mi>j</mi><mi>X</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122361800000025.GIF" wi="1293" he="143" /></maths>其中,<img file="FDA0001122361800000026.GIF" wi="118" he="77" />表示第k+1次迭代的次暂态故障电压三相分量,<img file="FDA0001122361800000027.GIF" wi="84" he="71" />表示内电势的三相分量;步骤4:计算稳态模型的内电势正序分量:<maths num="0006"><math><![CDATA[<mrow><msup><mover><mi>E</mi><mo>&CenterDot;</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>=</mo><mfrac><mrow><mi>Q</mi><mi>X</mi></mrow><msubsup><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mi>f</mi><mi>k</mi></msubsup></mfrac><mo>+</mo><mi>j</mi><mfrac><mrow><msub><mi>P</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub><mi>R</mi></mrow><msubsup><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mi>f</mi><mi>k</mi></msubsup></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122361800000028.GIF" wi="1324" he="143" /></maths>其中<img file="FDA0001122361800000029.GIF" wi="89" he="55" />表示第k+1次迭代时的稳态内电势正序分量,<img file="FDA00011223618000000210.GIF" wi="63" he="71" />为第k次迭代时的稳态故障电压正序分量,Q与P<sub>ref</sub>由式(4)计算,R与X分别为逆变器出口电阻和电抗;步骤5:形成稳态电流的三相分量迭代格式:<maths num="0007"><math><![CDATA[<mrow><msub><mi>Y</mi><mrow><mi>a</mi><mi>b</mi><mi>c</mi></mrow></msub><msubsup><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mrow><mi>f</mi><mi>a</mi><mi>b</mi><mi>c</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mfrac><msubsup><mover><mi>E</mi><mo>&CenterDot;</mo></mover><mrow><mi>a</mi><mi>b</mi><mi>c</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mrow><mi>R</mi><mo>+</mo><mi>j</mi><mi>X</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122361800000031.GIF" wi="1278" he="135" /></maths>其中,<img file="FDA0001122361800000032.GIF" wi="87" he="70" />表示第k+1次迭代时的稳态内电势三相分量,<img file="FDA0001122361800000033.GIF" wi="101" he="78" />表示第k+1次迭代的稳态故障电压三相分量;步骤6:统一求解方程:<maths num="0008"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>Y</mi><mrow><mi>a</mi><mi>b</mi><mi>c</mi></mrow></msub><msubsup><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mrow><mi>f</mi><mi>a</mi><mi>b</mi><mi>c</mi><mn>0</mn></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mfrac><msubsup><mover><mi>E</mi><mo>&CenterDot;</mo></mover><mrow><mi>a</mi><mi>b</mi><mi>c</mi></mrow><mn>0</mn></msubsup><mrow><mi>R</mi><mo>+</mo><mi>j</mi><mi>X</mi></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>Y</mi><mrow><mi>a</mi><mi>b</mi><mi>c</mi></mrow></msub><msubsup><mover><mi>U</mi><mo>&CenterDot;</mo></mover><mrow><mi>f</mi><mi>a</mi><mi>b</mi><mi>c</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mfrac><msubsup><mover><mi>E</mi><mo>&CenterDot;</mo></mover><mrow><mi>a</mi><mi>b</mi><mi>c</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mrow><mi>R</mi><mo>+</mo><mi>j</mi><mi>X</mi></mrow></mfrac></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122361800000034.GIF" wi="1318" he="295" /></maths>得到次暂态故障电压三相分量和稳态故障电压三相分量;步骤7:判断方程<img file="FDA0001122361800000035.GIF" wi="323" he="71" />是否成立,即判断方程<img file="FDA0001122361800000036.GIF" wi="322" he="71" />是否收敛,如果成立则迭代结束,否则返回步骤3继续迭代,其中,1≤i≤n,ε为给定的迭代计算误差,i表示各节点编号,n表示总节点数。
地址 250002 山东省济南市市中区望岳路2000号