发明名称 基于三阶收敛牛顿法VSC-HVDC潮流计算方法
摘要 本发明公布了一种基于三阶收敛牛顿法的VSC-HVDC潮流计算方法。首先推导了VSC-HVDC的稳态潮流模型,并在此基础上推导了具有三阶收敛速度的牛顿迭代法的多变量矩阵求解格式,将它应用VSC-HVDC潮流计算;对IEEE5-57节点测试系统等4个算例进行了仿真测试,结果表明,该算法具有良好的收敛特性,并且在达到同样精度要求的情况下,它较之经典牛顿法需要较少的迭代次数。由于在每步迭代中充分利用了Jacobian矩阵三角分解的因子表,提高了潮流计算的速度。
申请公布号 CN102545252B 申请公布日期 2014.05.07
申请号 CN201110433307.7 申请日期 2011.12.22
申请人 河海大学 发明人 孙国强;卫志农;季聪;韦延方;杨雄;袁阳;陆子刚;张伟;刘玉娟;陈凡;潘春兰
分类号 H02J3/36(2006.01)I;G06F17/16(2006.01)I 主分类号 H02J3/36(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 许方
主权项 1.一种基于三阶收敛牛顿法的VSC-HVDC潮流计算方法,其特征在于包括以下步骤:(1)获得电力系统的网络参数,包括:母线编号、名称、负荷有功、负荷无功、补偿电容,输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗;(2)程序初始化,包括:对状态量设置初值,节点次序优化,形成节点导纳矩阵,设置误差精度ε;(3)恢复迭代计数器:k=1;(4)计算功率不平衡量f(x<sup>(k)</sup>);(5)计算系统的Jacobian矩阵;(6)解下述方程组,求得<img file="FDA0000434523990000011.GIF" wi="142" he="60" />并得到修正后的<img file="FDA0000434523990000012.GIF" wi="108" he="61" /><maths num="0001"><![CDATA[<math><mrow><mi>J</mi><mrow><mo>(</mo><msup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>)</mo></mrow><mi>&Delta;</mi><msup><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>=</mo><mo>-</mo><mi>f</mi><mrow><mo>(</mo><msup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msup><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>=</mo><msup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>+</mo><mi>&Delta;</mi><msup><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup></mrow></math>]]></maths>(7)判断<img file="FDA0000434523990000015.GIF" wi="187" he="93" />&lt;ε是否成立,若是,停止计算,输出结果,若否,置k=k+0.5,0.5表示完成三阶牛顿法中的半次迭代,返回步骤(4);(8)计算功率不平衡量<img file="FDA0000434523990000016.GIF" wi="355" he="73" />(9)解下述方程组,求得Δx<sup>(k)</sup>,并得到修正后的x<sup>(k)</sup>:<maths num="0003"><![CDATA[<math><mrow><mi>J</mi><mrow><mo>(</mo><msup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>)</mo></mrow><mi>&Delta;</mi><msup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>=</mo><mo>-</mo><mi>f</mi><mrow><mo>(</mo><msup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>+</mo><mi>&Delta;</mi><msup><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>)</mo></mrow></mrow></math>]]></maths>x<sup>(k+1)</sup>=x<sup>(k)</sup>+Δx<sup>(k)</sup>(10)判断<img file="FDA0000434523990000018.GIF" wi="188" he="93" />&lt;ε是否成立,若是,停止计算,输出结果,若否,置k=k+0.5,0.5表示完成三阶牛顿法中的半次迭代,返回步骤(4)。
地址 210098 江苏省南京市西康路1号