发明名称 一种计及远程电压控制模式的节点类型扩展连续潮流算法
摘要 本发明公开了一种计及远程电压控制模式的节点类型扩展连续潮流算法,包括以下步骤:根据节点类型扩展潮流可解性判据设定模拟远程电压控制模式的P节点和PQV节点;建立节点类型扩展连续潮流模型,利用预测-校正算法求解静态电压稳定临界点;给出基于节点类型扩展连续潮流法计算远程电压控制模式下电力系统静态稳定临界点过程中发电机无功越限的处理方式;给出基于节点类型扩展连续潮流的电压稳定临界点分岔类型识别方法。本发明适用于远程电压控制模式下电力系统的静态电压稳定分析。
申请公布号 CN103218531B 申请公布日期 2016.03.09
申请号 CN201310129484.5 申请日期 2013.04.15
申请人 河海大学 发明人 赵晋泉;周超;陈刚
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 夏雪
主权项 一种计及远程电压控制模式的节点类型扩展连续潮流算法,其特征在于,包括以下步骤:步骤A、根据节点类型扩展潮流可解性判据设定模拟远程电压控制模式的P节点和PQV节点;步骤B、建立节点类型扩展连续潮流模型,利用预测‑校正算法求解静态电压稳定临界点,具体过程如下:B1、建立节点类型扩展连续潮流数学模型,其具体数学表达式如下:发电机有功出力以及负荷有功和无功的变化的规律用如下数学表达式来模拟:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>P</mi><mi>Gi</mi></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>Gi</mi><mn>0</mn></mrow></msub><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>&lambda;</mi><msub><mi>K</mi><mi>PGi</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>P</mi><mi>Li</mi></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>Li</mi><mn>0</mn></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>&lambda;</mi><msub><mi>K</mi><mi>PLi</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mi>Li</mi></msub><mo>=</mo><msub><mi>Q</mi><mrow><mi>Li</mi><mn>0</mn></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>&lambda;</mi><msub><mi>K</mi><mi>QLi</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000870828190000011.GIF" wi="488" he="258" /></maths>式中:P<sub>Gi</sub>、P<sub>Li</sub>、Q<sub>Li</sub>分别为节点i的发电机发出有功、负荷有功、负荷无功,P<sub>Gi0</sub>、P<sub>Li0</sub>、Q<sub>Li0</sub>分别为基态潮流时节点i的发电机发出有功、负荷有功、负荷无功,K<sub>PGi</sub>、K<sub>PLi</sub>、K<sub>QLi</sub>分别为节点i的发电有功增长系数和负荷的有功增长系数、无功增长系数,λ为负荷参数;计及发电机和负荷增长参数的节点类型扩展连续潮流方程为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&Delta;</mi><msub><mi>P</mi><mi>i</mi></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>Gi</mi><mn>0</mn></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>&lambda;</mi><msub><mi>K</mi><mi>PGi</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>P</mi><mrow><mi>Li</mi><mn>0</mn></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>&lambda;</mi><msub><mi>K</mi><mi>PLi</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>V</mi><mi>i</mi></msub><munder><mi>&Sigma;</mi><mrow><mi>i</mi><mo>&Element;</mo><mi>j</mi></mrow></munder><msub><mi>V</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>G</mi><mi>ij</mi></msub><mi>sin</mi><msub><mi>&theta;</mi><mi>ij</mi></msub><mo>+</mo><msub><mi>B</mi><mi>ij</mi></msub><mi>cos</mi><msub><mi>&theta;</mi><mi>ij</mi></msub><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>S</mi><mi>p</mi></msub></mrow>]]></math><img file="FDA0000870828190000012.GIF" wi="1424" he="102" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&Delta;</mi><msub><mi>Q</mi><mi>i</mi></msub><mo>=</mo><msub><mi>Q</mi><mi>Gi</mi></msub><mo>-</mo><msub><mi>Q</mi><mrow><mi>Li</mi><mn>0</mn></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>&lambda;</mi><msub><mi>K</mi><mi>QLi</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>V</mi><mi>i</mi></msub><munder><mi>&Sigma;</mi><mrow><mi>i</mi><mo>&Element;</mo><mi>j</mi></mrow></munder><msub><mi>V</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>G</mi><mi>ij</mi></msub><mi>cos</mi><msub><mi>&theta;</mi><mi>ij</mi></msub><mo>-</mo><msub><mi>B</mi><mi>ij</mi></msub><mi>sin</mi><msub><mi>&theta;</mi><mi>ij</mi></msub><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>S</mi><mi>q</mi></msub></mrow>]]></math><img file="FDA0000870828190000013.GIF" wi="1131" he="91" /></maths>式中:Q<sub>Gi</sub>为节点i的发电机输出无功,ΔP<sub>i</sub>、ΔQ<sub>i</sub>分别为节点i的有功偏差、无功偏差,V<sub>i</sub>、V<sub>j</sub>分别为节点i、j的电压幅值,θ<sub>j</sub>为节点i,j间的电压相角差,G<sub>j</sub>、B<sub>j</sub>分别为导纳矩阵中第(i,j)个元素的实部与虚部,S<sub>p</sub>为1p、3p类节点的集合,S<sub>q</sub>为1q、3q类节点的集合;将以上的方程用向量表示,节点类型扩展连续潮流的模型简写如下:f(x,λ)=0式中f为节点类型扩展潮流方程,x包含系统中的母线电压幅值和相角;B2、求解节点类型扩展连续潮流的预测步:首先由下式求解出切向量[dθ<sub>1p</sub>,dV<sub>1q</sub>,,dθ<sub>2p</sub>,dV<sub>2q</sub>,dλ]<sup>T</sup>:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mfenced open='' close=''><mtable><mtr><mtd><msup><mi>J</mi><mo>&prime;</mo></msup></mtd><mtd><mfrac><mrow><mo>&PartialD;</mo><mi>f</mi></mrow><mrow><mo>&PartialD;</mo><mi>&lambda;</mi></mrow></mfrac></mtd></mtr></mtable></mfenced></mtd></mtr><mtr><mtd><msub><mi>e</mi><mi>h</mi></msub></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>d&theta;</mi><mrow><mn>1</mn><mi>p</mi></mrow></msub></mtd></mtr><mtr><mtd><mi>d</mi><msub><mi>V</mi><mrow><mn>1</mn><mi>q</mi></mrow></msub><mo>/</mo><msub><mi>V</mi><mrow><mn>1</mn><mi>q</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>d&theta;</mi><mrow><mn>2</mn><mi>p</mi></mrow></msub></mtd></mtr><mtr><mtd><mi>d</mi><msub><mi>V</mi><mrow><mn>2</mn><mi>q</mi></mrow></msub><mo>/</mo><msub><mi>V</mi><mrow><mn>2</mn><mi>q</mi></mrow></msub></mtd></mtr><mtr><mtd><mi>d&lambda;</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>&PlusMinus;</mo><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000870828190000021.GIF" wi="499" he="347" /></maths>式中:e<sub>h</sub>为行向量,其中第h个元素为1,其余元素为0,其中±1中的正负号取决于切向量中相应分量的符号,dθ<sub>1p</sub>、dθ<sub>2p</sub>分别为1p、2p类节点相角的切向量,dV<sub>1q</sub>、dV<sub>2q</sub>分别为1q、2q类节点电压幅值切向量,dλ为负荷参数求导,J′为f对x求导得出的雅可比矩阵;切向量得到后,按照给定的步长控制参数σ,通过下式即可获得预测值:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>x</mi><mrow><mo>&prime;</mo><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msup></mtd></mtr><mtr><mtd><msup><mi>&lambda;</mi><mrow><mo>&prime;</mo><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msup></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup></mtd></mtr><mtr><mtd><msup><mi>&lambda;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup></mtd></mtr></mtable></mfenced><mo>+</mo><mi>&sigma;</mi><mfenced open='[' close=']'><mtable><mtr><mtd><mi>dx</mi></mtd></mtr><mtr><mtd><mi>d&lambda;</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000870828190000022.GIF" wi="463" he="134" /></maths>式中:x<sup>(k)</sup>为1p、2p类节点的相角的当前解和1q、2q类节点电压幅值的当前解组成的向量,x′<sup>(k+1)</sup>为相应的预测值;<img file="FDA0000870828190000023.GIF" wi="73" he="51" />为负荷参数的当前的解,<img file="FDA0000870828190000024.GIF" wi="108" he="51" />为相应的预测值;B3、采用局部参数化策略,避免节点类型扩展连续潮流在临界点奇异,参数化的方程如下式:x<sub>h</sub>‑x′<sub>h</sub>‑Δs=0式中:x<sub>h</sub>为当前解向量中的第h个分量,x′<sub>h</sub>为预测解向量中第h个分量,Δs为步长;其中h的取值如下式所示,取电压跌落情况最严重的节点,h的取法:h={h||dx<sub>h</sub>|=||dx||<sub>∞</sub>}式中:dx<sub>h</sub>为dx的第h个元素;B4、用下式求解校正步:以预测步中得到的预测值作为解的初始值,采用牛顿法或者拟牛顿法进行求解:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>x</mi><mi>h</mi></msub><mo>-</mo><msubsup><mi>x</mi><mi>h</mi><mo>&prime;</mo></msubsup><mo>-</mo><mi>&Delta;s</mi><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000870828190000025.GIF" wi="318" he="137" /></maths>步骤C、给出基于节点类型扩展连续潮流法计算远程电压控制模式下电力系统静态稳定临界点过程中发电机无功越限的处理方式;步骤D、给出基于节点类型扩展连续潮流的电压稳定临界点分岔类型识别方法。
地址 211100 江苏省南京市江宁开发区佛成西路8号