发明名称 一种不对称电力系统暂态稳定的建模和仿真方法
摘要 一种不对称电力系统暂态稳定的建模和仿真方法,采用相分量法和序分量法混合建模的方法进行不对称电力系统的建模,在此基础上,基于边界节点的三相功率平衡条件构造了隐函数形式的边界协调方程,提出基于分解协调的不对称电力系统暂态稳定仿真算法。采用本发明能够解决传统方法无法兼顾高效和灵活性的困难,可方便地处理复杂故障和多重故障计算问题,克服了传统的基于序网连接的故障处理方式的不足。同时,还能很好地解决参数不对称情况下系统的暂稳仿真问题。IEEE标准系统的测试表明,本发明具有较高的准确性和收敛性。
申请公布号 CN103761351B 申请公布日期 2017.01.04
申请号 CN201310469025.1 申请日期 2013.10.09
申请人 广西电网公司;清华大学 发明人 黄秀琼;黄少伟;陈颖;沈沉;陈卫东;白先红;覃惠玲
分类号 G06F17/50(2006.01)I;G06F9/455(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 广西南宁公平专利事务所有限责任公司 45104 代理人 黄永校
主权项 一种不对称电力系统暂态稳定的建模和仿真方法,其特征在于,包括以下步骤:步骤(1)给定一组t<sub>n</sub>时刻边界母线集的状态量<img file="FDA0001047769460000011.GIF" wi="60" he="70" />和<img file="FDA0001047769460000012.GIF" wi="83" he="79" />计算t<sub>n</sub>时刻边界母线集的功率偏差<maths num="0001"><math><![CDATA[<mrow><msub><mi>&Delta;P</mi><msub><mi>t</mi><mi>n</mi></msub></msub><mo>=</mo><msubsup><mi>P</mi><msub><mi>t</mi><mi>n</mi></msub><mi>B</mi></msubsup><mo>+</mo><msubsup><mi>P</mi><msub><mi>t</mi><mi>n</mi></msub><mover><mi>B</mi><mo>~</mo></mover></msubsup><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>P</mi><msub><mi>t</mi><mi>n</mi></msub><mrow><mi>B</mi><mn>1</mn></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>P</mi><msub><mi>t</mi><mi>n</mi></msub><mrow><mi>B</mi><mn>2</mn></mrow></msubsup></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>P</mi><msub><mi>t</mi><mi>n</mi></msub><mrow><mover><mi>B</mi><mo>~</mo></mover><mn>1</mn></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>P</mi><msub><mi>t</mi><mi>n</mi></msub><mrow><mover><mi>B</mi><mo>~</mo></mover><mn>2</mn></mrow></msubsup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001047769460000013.GIF" wi="1142" he="183" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>&Delta;Q</mi><msub><mi>t</mi><mi>n</mi></msub></msub><mo>=</mo><msubsup><mi>Q</mi><msub><mi>t</mi><mi>n</mi></msub><mi>B</mi></msubsup><mo>+</mo><msubsup><mi>Q</mi><msub><mi>t</mi><mi>n</mi></msub><mover><mi>B</mi><mo>~</mo></mover></msubsup><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>Q</mi><msub><mi>t</mi><mi>n</mi></msub><mrow><mi>B</mi><mn>1</mn></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>Q</mi><msub><mi>t</mi><mi>n</mi></msub><mrow><mover><mi>B</mi><mo>~</mo></mover><mn>2</mn></mrow></msubsup></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>Q</mi><msub><mi>t</mi><mi>n</mi></msub><mrow><mover><mi>B</mi><mo>~</mo></mover><mn>1</mn></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>Q</mi><msub><mi>t</mi><mi>n</mi></msub><mrow><mover><mi>B</mi><mo>~</mo></mover><mn>2</mn></mrow></msubsup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001047769460000014.GIF" wi="1166" he="181" /></maths>将<img file="FDA0001047769460000015.GIF" wi="85" he="63" />和<img file="FDA0001047769460000016.GIF" wi="98" he="71" />作为协调目标,得到t<sub>n</sub>时刻的边界协调方程:<maths num="0003"><math><![CDATA[<mrow><msubsup><mi>&Phi;</mi><msub><mi>t</mi><mi>n</mi></msub><mi>B</mi></msubsup><mrow><mo>(</mo><msubsup><mi>y</mi><msub><mi>t</mi><mi>n</mi></msub><mi>B</mi></msubsup><mo>)</mo></mrow><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mi>&Delta;P</mi><msub><mi>t</mi><mi>n</mi></msub></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&Delta;Q</mi><msub><mi>t</mi><mi>n</mi></msub></msub></mrow></mtd></mtr></mtable></mfenced><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001047769460000017.GIF" wi="1045" he="159" /></maths>步骤(2)基于式(3),构造基于迭代求解的不对称系统暂态仿真分解协调求解算法:t<sub>n</sub>=t<sub>0</sub>,t<sub>n+1</sub>=t<sub>0</sub>+△t,协调侧初始化计算环境并设置t<sub>n+1</sub>时刻的边界母线状态<img file="FDA0001047769460000018.GIF" wi="75" he="70" />和<img file="FDA0001047769460000019.GIF" wi="110" he="79" />令<img file="FDA00010477694600000110.GIF" wi="227" he="78" />步骤(3)根据协调侧给定的t<sub>n+1</sub>时刻的边界母线状态求解各自子系统的动态方程,得到边界母线集功率偏差<img file="FDA00010477694600000111.GIF" wi="98" he="63" />和<img file="FDA00010477694600000112.GIF" wi="115" he="63" />并返回协调侧;步骤(4)若<img file="FDA00010477694600000113.GIF" wi="211" he="78" />且<img file="FDA00010477694600000114.GIF" wi="219" he="79" />则转至步骤(5),否则由<img file="FDA00010477694600000115.GIF" wi="101" he="63" />和<img file="FDA00010477694600000116.GIF" wi="115" he="63" />求解<img file="FDA00010477694600000117.GIF" wi="115" he="71" />更新边界母线状态,令<maths num="0004"><math><![CDATA[<mrow><msubsup><mi>y</mi><msub><mi>t</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mi>B</mi></msubsup><mo>=</mo><msubsup><mi>y</mi><msub><mi>t</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mover><mi>B</mi><mo>~</mo></mover></msubsup><mo>=</mo><msubsup><mi>y</mi><msub><mi>t</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mi>B</mi></msubsup><mo>+</mo><msubsup><mi>&Delta;y</mi><msub><mi>t</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mi>B</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00010477694600000118.GIF" wi="1066" he="87" /></maths>然后返回步骤(3);步骤(5)t<sub>n+1</sub>时刻全网一致收敛,t<sub>n+1</sub>=t<sub>n+1</sub>+△t,若t<sub>n+1</sub>≤t<sub>end</sub>,则返回步骤(3)开始下一时步的仿真,否则仿真结束。
地址 530023 广西壮族自治区南宁市兴宁区民主路6号