发明名称 一种采用SVC作为无功补偿设备的无功规划方法
摘要 本发明公开一种采用SVC进行考虑电网大扰动后系统暂态过程稳定性约束的电力系统无功规划方法。通过不同的电压稳定原则生成无功补偿设备选点指标矩阵,采用模糊聚类进行无功补偿装置选点布置,得到适应SVC运行性能的无功补偿布点方案;通过建立基于不同运行方式下无功补偿优化模型,基于复杂高阶系统设备动态模型,利用时域仿真及优化算法得到相应候选无功补偿点SVC容量配置最优结果,从而实现SVC无功规划。本发明的无功规划方法根据模糊聚类理论、动态电力系统理论、时域仿真技术及内点法优化综合算法,解决了考虑系统发生故障后的暂态过程稳定性约束的无功规划问题,能充分发挥SVC动态响应特性,有效提高了电力系统无功补偿运行效率。
申请公布号 CN104701864A 申请公布日期 2015.06.10
申请号 CN201510097259.7 申请日期 2015.03.04
申请人 东南大学 发明人 王玉荣;陈昊;李庆昌
分类号 H02J3/18(2006.01)I 主分类号 H02J3/18(2006.01)I
代理机构 江苏永衡昭辉律师事务所 32250 代理人 王斌
主权项 一种采用SVC作为无功补偿设备的无功规划方法,其特征在于:具体步骤如下:步骤1:基于系统故障排序法,在高峰负荷、平均负荷水平下根据公式(1)、(2)、(3)分别计算最严重故障后的暂态过程电压稳定性指标:系统暂态电压跌落:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>V</mi><mi>dip</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>V</mi><mn>0</mn></msub><mo>-</mo><msub><mi>V</mi><mi>s</mi></msub></mrow><msub><mi>V</mi><mn>0</mn></msub></mfrac><mo>=</mo><mfrac><mi>&Delta;V</mi><msub><mi>V</mi><mn>0</mn></msub></mfrac><mo>&times;</mo><mn>100</mn><mo>%</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000677577180000011.GIF" wi="1167" he="142" /></maths>其中,V<sub>s</sub>为短期稳定过程中节点电压最低幅值,V<sub>0</sub>为节点电压初始幅值;暂态过程中电压偏移指标:<img file="FDA0000677577180000012.GIF" wi="1392" he="316" />暂态过程中低电压持续时间指标:<img file="FDA0000677577180000013.GIF" wi="1438" he="176" />其中,其中,t<sub>cl</sub>为故障清除时间,t<sub>f</sub>为仿真终止时间;步骤2:根据步骤1得到的高峰负荷、平均负荷水平下的暂态过程中的各个电压稳定性指标,采用模糊聚类方法,分别得到SVC在系统高峰负荷、平均负荷水平下的的候选无功补偿的节点i位置集合N<sub>C</sub>;步骤3:建立SVC容量规划最优化模型:选择三阶或以上发电机模型、一阶励磁系统模型及混合负荷模型:三阶发电机模型:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mrow><mi>d&delta;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac><mo>=</mo><mrow><mo>(</mo><mi>&omega;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>&omega;</mi><mi>B</mi></msub></mtd></mtr><mtr><mtd><msub><mi>T</mi><mi>j</mi></msub><mfrac><mrow><mi>d&omega;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac><mo>=</mo><msub><mi>T</mi><mi>m</mi></msub><mo>-</mo><msub><mi>T</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>D</mi><mrow><mo>(</mo><mi>&omega;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>T</mi><mrow><mi>d</mi><mn>0</mn></mrow><mo>&prime;</mo></msubsup><mfrac><mrow><msubsup><mi>dE</mi><mi>q</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac><mo>=</mo><msub><mi>E</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>E</mi><mi>q</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>X</mi><mi>d</mi></msub><mo>-</mo><msubsup><mi>X</mi><mi>d</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><msub><mi>I</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000677577180000014.GIF" wi="762" he="356" /></maths>其中,δ(t)为q轴相对于同步旋转坐标系的角位移;ω<sub>B</sub>与ω(t)分别为转子额定转速及转子t时刻角速度;T<sub>j</sub>与T′<sub>d0</sub>为时间常数,T<sub>m</sub>与T<sub>e</sub>(t)分别为机械力矩与机电力矩;E<sub>f</sub>,E′<sub>q</sub>分别为定子励磁电动势及电机q轴瞬变电动势;X<sub>d</sub>与X′<sub>d</sub>分别为d轴同步电抗与d轴瞬变电抗;I<sub>d</sub>为d轴电流;其机端电压方程为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>V</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>x</mi><mi>q</mi></msub><msub><mi>I</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>r</mi><mi>a</mi></msub><msub><mi>I</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>V</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>E</mi><mi>q</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>x</mi><mi>d</mi><mo>&prime;</mo></msubsup><msub><mi>I</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>r</mi><mi>a</mi></msub><msub><mi>I</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000677577180000021.GIF" wi="559" he="143" /></maths>其中,V<sub>d</sub>、V<sub>q</sub>(I<sub>d</sub>、I<sub>q</sub>)分别为d轴电压及q轴电压(电流);x<sub>q</sub>、r<sub>a</sub>分别为q轴同步电抗及绕组电阻;一阶励磁系统模型:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>T</mi><mi>A</mi></msub><mfrac><mrow><msub><mi>dE</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac><mo>=</mo><msub><mi>K</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>V</mi><mi>ref</mi></msub><mo>-</mo><msub><mi>V</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><msub><mi>E</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>E</mi><mrow><mi>f</mi><mn>0</mn></mrow></msub></mrow>]]></math><img file="FDA0000677577180000022.GIF" wi="720" he="123" /></maths>其中,T<sub>A</sub>为时间常数;K<sub>A</sub>为控制系数;V<sub>ref</sub>及V<sub>m</sub>分别为节点电压参考值及幅值;E<sub>f0</sub>为励磁电动势初始值;以SVC无功投入成本最小为条件作出目标函数:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>min</mi><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>b</mi></msub></munderover><mo>[</mo><msub><mi>C</mi><mrow><mi>fix</mi><mo>_</mo><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>C</mi><mrow><mi>i</mi><mo>_</mo><mi>v</mi><mn>1</mn></mrow></msub><mi>max</mi><mrow><mo>(</mo><msubsup><mi>Q</mi><mi>SVCi</mi><mo>+</mo></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>C</mi><mrow><mi>i</mi><mo>_</mo><mi>v</mi><mn>2</mn></mrow></msub><mi>max</mi><mrow><mo>(</mo><mo>|</mo><msubsup><mi>Q</mi><mi>SVCi</mi><mo>-</mo></msubsup><mo>|</mo><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000677577180000023.GIF" wi="1029" he="140" /></maths>其中,C<sub>fix_i</sub>为无功投资固定成本;C<sub>i_v1</sub>为容性无功可变成本;C<sub>i_v2</sub>为感性无功可变成本;<img file="FDA0000677577180000024.GIF" wi="98" he="64" />为节点i容性SVC补偿容量;<img file="FDA0000677577180000025.GIF" wi="97" he="62" />为感性的SVC补偿容量;静态安全约束条件为<maths num="0006" id="cmaths0006"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>P</mi><mi>Gi</mi></msub><mo>-</mo><msub><mi>P</mi><mi>Li</mi></msub><mo>-</mo><msub><mi>V</mi><mi>i</mi></msub><munder><mi>&Sigma;</mi><mi>j&omega;i</mi></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><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mi>Gi</mi></msub><mo>+</mo><msub><mi>Q</mi><mi>ci</mi></msub><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>Q</mi><mi>Li</mi></msub><mo>-</mo><msub><mi>V</mi><mi>i</mi></msub><munder><mi>&Sigma;</mi><mi>j&omega;i</mi></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><mn>0</mn></mtd></mtr><mtr><mtd><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>V</mi><mrow><mi>i</mi><mi>min</mi></mrow></msub><mo>&le;</mo><msub><mi>V</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>V</mi><mrow><mi>i</mi><mi>max</mi></mrow></msub></mtd><mtd><mi>i</mi><mo>&Element;</mo><msub><mi>N</mi><mi>b</mi></msub></mtd></mtr></mtable></mfenced></mtd></mtr><mtr><mtd><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>P</mi><mrow><mi>Gi</mi><mi>min</mi></mrow></msub><mo>&le;</mo><msub><mi>P</mi><mi>Gi</mi></msub><mo>&le;</mo><msub><mi>P</mi><mrow><mi>Gi</mi><mi>max</mi></mrow></msub></mtd><mtd><mi>i</mi><mo>&Element;</mo><msub><mi>N</mi><mi>G</mi></msub></mtd></mtr></mtable></mfenced></mtd></mtr><mtr><mtd><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>Q</mi><mrow><mi>Gi</mi><mi>min</mi></mrow></msub><mo>&le;</mo><msub><mi>Q</mi><mi>Gi</mi></msub><mo>&le;</mo><msub><mi>Q</mi><mrow><mi>Gi</mi><mi>max</mi></mrow></msub></mtd><mtd><mi>i</mi><mo>&Element;</mo><msub><mi>N</mi><mi>G</mi></msub></mtd></mtr></mtable></mfenced></mtd></mtr><mtr><mtd><mo>-</mo><msub><mi>S</mi><mrow><mi>ij</mi><mi>max</mi></mrow></msub><mo>&le;</mo><msub><mi>S</mi><mi>ij</mi></msub><mo>&le;</mo><msub><mi>S</mi><mrow><mi>ij</mi><mi>max</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mrow><mi>SVCi</mi><mi>min</mi></mrow></msub><mo>&le;</mo><msub><mi>Q</mi><mi>SVCi</mi></msub><mo>&le;</mo><msub><mi>Q</mi><mrow><mi>SVCi</mi><mi>max</mi></mrow></msub></mtd></mtr><mtr><mtd><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>y</mi><mi>i</mi></msub><mo>=</mo><msub><mi>N</mi><mi>c</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000677577180000026.GIF" wi="1084" he="772" /></maths>约束条件中,P<sub>Gi</sub>,P<sub>Li</sub>分别节点i处的发电机有功出力和负荷有功需求;Q<sub>Gi</sub>,Q<sub>Li</sub>分别节点i处的发电机无功出力和负荷无功需求;Q<sub>ci</sub>为节点i的无功补偿容量;N<sub>b</sub>为节点集合,N<sub>G</sub>为系统发电机集合,N<sub>C</sub>为SVC补偿节点集合;G<sub>ij</sub>与B<sub>ij</sub>分别为节点导纳矩阵的电导与电纳;θ<sub>ij</sub>为节点i与节点j的相角差;V<sub>i</sub>,V<sub>imax</sub>,V<sub>imin</sub>分别为节点i电压幅值、电压幅值最大值及最小值;P<sub>imax</sub>,P<sub>imin</sub>分别为发电机节点i有功输出最大值及最小值;Q<sub>imax</sub>,Q<sub>imin</sub>分别为发电机节点i无功输出最大值及最小值;Q<sub>SVCi</sub>、Q<sub>SVCimax</sub>及Q<sub>SVCimin</sub>分别为无功补偿候选节点i无功补偿容量、无功补偿最大值及最小值;S<sub>ij</sub>,S<sub>ijmax</sub>分别为线路ij的视在功率及其最大值;y<sub>i</sub>为布尔变量,y<sub>i</sub>=1表示节点i为候选无功补偿点;暂态过程中电压稳定约束包括:短期电压跌落约束:<img file="FDA0000677577180000031.GIF" wi="790" he="162" />短期电压稳定过程中低电压持续时间约束:负荷节点电压在0.75V<sub>0</sub>≤V<sub>s</sub>(t)≤0.8V<sub>0</sub>范围的持续时间不能超过20周期,即0.4s;达到新的稳态过程后电压幅值约束:V<sub>i</sub>≥0.95V<sub>0i</sub> i=1,...,N<sub>b</sub>;步骤4:采用隐式梯形积分法,通过微分方程差分化方法,将优化模型中的微分方程转变为代数差分化方程;根据时域仿真法及内点优化法,根据目标函数以及约束条件,分别计算得到在系统高峰负荷、平均负荷运行模式下SVC的各个候选无功规划位置的无功补偿容量;步骤5:根据步骤2、4在不同负荷水平下得到的SVC选点因素,结合时域仿真分析,校核SVC最终需要补偿的位置及容量。
地址 210096 江苏省南京市四牌楼2号