发明名称 基于动态无功设备的电力系统中中枢母线电压的控制方法
摘要 本发明涉及一种基于动态无功设备的电力系统中中枢母线电压的控制方法,属于电力系统评估和控制技术领域。首先,监测系统当中中枢节点电压值,根据该电压值计算出参与电压控制的发电机与动态无功补偿设备的控制指令。将指令下发给相应的发电机与动态无功补偿设备。判断是否需要进行动态无功优化,如果不需要,则在时刻T<sub>0</sub>+△T,返回第一步骤,若需要,则根据发电机与动态无功补偿装置的无功功率及上限和下限值,计算动态无功优化的控制指令并下发。发电机与动态无功补偿装置根据指令调整各自的无功功率。本发明方法能够有效的进行暂态电压稳定的预防控制,提高电力系统中动态无功储备水平和扰动后的系统电压的恢复效果。
申请公布号 CN103825281B 申请公布日期 2016.02.24
申请号 CN201410030530.0 申请日期 2014.01.22
申请人 清华大学 发明人 孙宏斌;郭庆来;王彬;张伯明;王旭冉;葛怀畅;吴文传
分类号 H02J3/16(2006.01)I 主分类号 H02J3/16(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 罗文群
主权项 一种基于动态无功补偿设备的电力系统中中枢母线电压的控制方法,其特征在于该控制方法包括以下步骤:(1)从电力系统的调度中心获取中枢母线电压的设定值<img file="FDA0000889310470000011.GIF" wi="111" he="77" />当前采集时刻T<sub>0</sub>的电力系统中中枢母线的电压值V<sub>p</sub>,通过下式计算得到中枢母线电压控制阶段的参与控制的发电机的电压调整量ΔV<sub>g</sub>和参与控制的动态无功补偿设备的电压调整量ΔV<sub>s</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><msub><mi>&Delta;V</mi><mi>g</mi></msub><mo>,</mo><msub><mi>&Delta;V</mi><mi>s</mi></msub></mrow></munder><mo>{</mo><msub><mi>W</mi><mi>p</mi></msub><mo>|</mo><mo>|</mo><mi>a</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>V</mi><mi>p</mi></msub><mo>-</mo><msubsup><mi>V</mi><mi>p</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>C</mi><mrow><mi>p</mi><mi>g</mi></mrow></msub><msub><mi>&Delta;V</mi><mi>g</mi></msub><mo>+</mo><msub><mi>C</mi><mrow><mi>p</mi><mi>s</mi></mrow></msub><msub><mi>&Delta;V</mi><mi>s</mi></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><mo>+</mo><msub><mi>W</mi><mi>q</mi></msub><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>&Theta;</mi><mi>g</mi></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><mo>+</mo><mo>|</mo><mo>|</mo><msub><mi>&Theta;</mi><mi>s</mi></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><mo>)</mo></mrow><mo>}</mo></mrow>]]></math><img file="FDA0000889310470000012.GIF" wi="1486" he="125" /></maths>受约束于:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>V</mi><mi>p</mi><mi>min</mi></msubsup><mo>&le;</mo><msub><mi>V</mi><mi>p</mi></msub><mo>+</mo><msub><mi>C</mi><mrow><mi>p</mi><mi>g</mi></mrow></msub><msub><mi>&Delta;V</mi><mi>g</mi></msub><mo>+</mo><msub><mi>C</mi><mrow><mi>p</mi><mi>s</mi></mrow></msub><msub><mi>&Delta;V</mi><mi>s</mi></msub><mo>&le;</mo><msubsup><mi>V</mi><mi>p</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000889310470000013.GIF" wi="736" he="77" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>Q</mi><mi>g</mi><mi>min</mi></msubsup><mo>&le;</mo><msub><mi>Q</mi><mi>g</mi></msub><mo>+</mo><msub><mi>C</mi><mi>g</mi></msub><msub><mi>&Delta;V</mi><mi>g</mi></msub><mo>+</mo><msub><mi>C</mi><mrow><mi>g</mi><mi>s</mi></mrow></msub><msub><mi>&Delta;V</mi><mi>s</mi></msub><mo>&le;</mo><msubsup><mi>Q</mi><mi>g</mi><mi>max</mi></msubsup></mrow>]]></math><img file="FDA0000889310470000014.GIF" wi="719" he="84" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>Q</mi><mi>s</mi><mi>min</mi></msubsup><mo>&le;</mo><msub><mi>Q</mi><mi>s</mi></msub><mo>+</mo><msub><mi>C</mi><mrow><mi>s</mi><mi>g</mi></mrow></msub><msub><mi>&Delta;V</mi><mi>g</mi></msub><mo>+</mo><msub><mi>C</mi><mi>s</mi></msub><msub><mi>&Delta;V</mi><mi>s</mi></msub><mo>&le;</mo><msubsup><mi>Q</mi><mi>s</mi><mi>max</mi></msubsup></mrow>]]></math><img file="FDA0000889310470000015.GIF" wi="710" he="80" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mo>|</mo><msub><mi>&Delta;V</mi><mi>g</mi></msub><mo>|</mo><mo>&le;</mo><msubsup><mi>&Delta;V</mi><mi>g</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000889310470000016.GIF" wi="288" he="94" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mo>|</mo><msub><mi>&Delta;V</mi><mi>s</mi></msub><mo>|</mo><mo>&le;</mo><msubsup><mi>&Delta;V</mi><mi>s</mi><mi>max</mi></msubsup></mrow>]]></math><img file="FDA0000889310470000017.GIF" wi="289" he="85" /></maths>其中,Θ<sub>g</sub>为发电机无功功率的裕度矢量,该裕度矢量Θ<sub>g</sub>中的第i个分量为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>&Theta;</mi><msub><mi>g</mi><mi>i</mi></msub></msub><mo>=</mo><mfrac><mrow><msub><mi>Q</mi><msub><mi>g</mi><mi>i</mi></msub></msub><mo>+</mo><munder><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><mi>N</mi><mi>g</mi></mrow></munder><mrow><msub><mi>C</mi><msub><mi>g</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></msub><msub><mi>&Delta;V</mi><msub><mi>g</mi><mi>j</mi></msub></msub></mrow><mo>+</mo><munder><mo>&Sigma;</mo><mrow><mi>k</mi><mo>&Element;</mo><mi>N</mi><mi>s</mi></mrow></munder><mrow><msub><mi>C</mi><mrow><msub><mi>gs</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub></mrow></msub><msub><mi>&Delta;V</mi><msub><mi>s</mi><mi>k</mi></msub></msub></mrow><mo>-</mo><msubsup><mi>Q</mi><msub><mi>g</mi><mi>i</mi></msub><mi>min</mi></msubsup></mrow><mrow><msubsup><mi>Q</mi><msub><mi>g</mi><mi>i</mi></msub><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msubsup><mo>-</mo><msubsup><mi>Q</mi><msub><mi>g</mi><mi>i</mi></msub><mi>min</mi></msubsup></mrow></mfrac></mrow>]]></math><img file="FDA0000889310470000018.GIF" wi="918" he="195" /></maths>Θ<sub>s</sub>为动态无功补偿设备无功功率的裕度矢量,裕度矢量Θ<sub>s</sub>的第i个分量为:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>&Theta;</mi><msub><mi>s</mi><mi>i</mi></msub></msub><mo>=</mo><mfrac><mrow><msub><mi>Q</mi><msub><mi>s</mi><mi>i</mi></msub></msub><mo>+</mo><munder><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><mi>N</mi><mi>s</mi></mrow></munder><mrow><msub><mi>C</mi><msub><mi>s</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></msub><msub><mi>&Delta;V</mi><msub><mi>s</mi><mi>j</mi></msub></msub></mrow><mo>+</mo><munder><mo>&Sigma;</mo><mrow><mi>k</mi><mo>&Element;</mo><mi>N</mi><mi>g</mi></mrow></munder><mrow><msub><mi>C</mi><mrow><msub><mi>sg</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub></mrow></msub><msub><mi>&Delta;V</mi><msub><mi>g</mi><mi>k</mi></msub></msub></mrow><mo>-</mo><msubsup><mi>Q</mi><msub><mi>s</mi><mi>i</mi></msub><mi>min</mi></msubsup></mrow><mrow><msubsup><mi>Q</mi><msub><mi>s</mi><mi>i</mi></msub><mi>max</mi></msubsup><mo>-</mo><msubsup><mi>Q</mi><msub><mi>s</mi><mi>i</mi></msub><mi>min</mi></msubsup></mrow></mfrac></mrow>]]></math><img file="FDA0000889310470000019.GIF" wi="1029" he="222" /></maths>其中下标p表示中枢母线,下标g表示发电机,下标s表示动态无功补偿设备;C<sub>pg</sub>是从调度中心获取的中枢母线电压偏差对参与控制的发电机的电压调整量ΔV<sub>g</sub>的灵敏度矩阵;C<sub>ps</sub>是从调度中心获取的中枢母线电压偏差对参与控制的动态无功补偿设备的电压调整量ΔV<sub>s</sub>的灵敏度矩阵;W<sub>p</sub>为中枢母线电压偏差的权重系数,W<sub>p</sub>的取值范围为0.1‑1,W<sub>q</sub>为发电机和动态无功补偿设备的无功功率储备裕度的权重系数,W<sub>q</sub>的取值范围为10‑100,a是中枢母线电压的增益系数,a的取值范围为1‑10,Q<sub>g</sub>是参与控制的发电机当前发出的无功功率,<img file="FDA0000889310470000021.GIF" wi="244" he="79" />分别为参与控制的发电机的无功功率下限与无功功率上限,Q<sub>s</sub>是参与控制的动态无功补偿设备当前发出的无功功率,<img file="FDA0000889310470000022.GIF" wi="252" he="77" />分别为参与控制的动态无功补偿设备发出的无功功率下限与无功功率的上限,C<sub>g</sub>为从调度中心获取的参与控制的发电机的无功功率调整量对参与控制的发电机的电压调整量ΔV<sub>g</sub>的灵敏度矩阵,C<sub>s</sub>为从调度中心获取的参与控制的动态无功补偿设备的无功功率调整量对参与控制的动态无功补偿设备的电压调整量ΔV<sub>s</sub>的灵敏度矩阵,C<sub>gs</sub>为从调度中心获取的参与控制的发电机的无功功率调整量对参与控制的动态无功补偿设备的电压调整量ΔV<sub>s</sub>的灵敏度矩阵,C<sub>sg</sub>为从调度中心获取的参与控制的动态无功补偿设备的无功功率调整量对参与控制的发电机的电压调整量ΔV<sub>g</sub>的灵敏度矩阵,<img file="FDA0000889310470000023.GIF" wi="98" he="79" />为设定的中枢母线电压的下限,<img file="FDA0000889310470000024.GIF" wi="102" he="78" />为设定的中枢母线电压上限,<img file="FDA0000889310470000025.GIF" wi="141" he="79" />为从调度中心获取的参与控制的发电机允许的单步最大调整量,<img file="FDA0000889310470000026.GIF" wi="143" he="71" />为从调度中心获取的参与控制的动态无功补偿设备允许的单步最大调整量,Ns为参与控制的动态无功补偿设备的编号的集合,Ng为参与控制的发电机的编号的集合;将上述计算得到的参与控制的发电机的电压调整量ΔV<sub>g</sub>和参与控制的动态无功补偿设备的电压调整量ΔV<sub>s</sub>作为控制指令下达给相应的发电机和动态无功补偿设备,各发电机和动态无功补偿设备根据控制指令,完成中枢母线电压控制阶段;(2)设定动态无功补偿设备无功功率的裕度阈值为<img file="FDA0000889310470000027.GIF" wi="359" he="71" />将当前动态无功补偿设备的裕度值与设定的裕度阈值比较,若当前动态无功补偿设备的裕度值大于或等于设定的裕度阈值,则经过中枢母线电压两次控制之间的时间间隔ΔT后重复步骤(1),ΔT的取值范围为5至10分钟,若当前动态无功补偿设备的裕度值小于设定的裕度阈值,则进行步骤(3);(3)根据下式计算,得到动态无功补偿设备动态无功功率裕度优化阶段的参与控制的发电机和动态无功补偿设备的电压调整量ΔV<sub>g</sub>和ΔV<sub>s</sub>:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><msub><mi>&Delta;V</mi><mi>g</mi></msub><mo>,</mo><msub><mi>&Delta;V</mi><mi>s</mi></msub></mrow></munder><mo>{</mo><msub><mi>W</mi><mi>s</mi></msub><mo>|</mo><mo>|</mo><msub><mi>&mu;</mi><mi>s</mi></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><mo>+</mo><msub><mi>W</mi><mi>g</mi></msub><mo>|</mo><mo>|</mo><msub><mi>&mu;</mi><mi>g</mi></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><mo>}</mo></mrow>]]></math><img file="FDA0000889310470000028.GIF" wi="558" he="127" /></maths>受约束于:V<sub>p</sub>‑V<sub>band</sub>≤V<sub>p</sub>+C<sub>pg</sub>ΔV<sub>g</sub>+C<sub>ps</sub>ΔV<sub>s</sub>≤V<sub>p</sub>+V<sub>band</sub><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msubsup><mi>Q</mi><mi>g</mi><mi>min</mi></msubsup><mo>&le;</mo><msub><mi>Q</mi><mi>g</mi></msub><mo>+</mo><msub><mi>C</mi><mi>g</mi></msub><msub><mi>&Delta;V</mi><mi>g</mi></msub><mo>+</mo><msub><mi>C</mi><mrow><mi>g</mi><mi>s</mi></mrow></msub><msub><mi>&Delta;V</mi><mi>s</mi></msub><mo>&le;</mo><msubsup><mi>Q</mi><mi>g</mi><mi>max</mi></msubsup></mrow>]]></math><img file="FDA0000889310470000031.GIF" wi="693" he="79" /></maths><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msubsup><mi>Q</mi><mi>s</mi><mi>min</mi></msubsup><mo>&le;</mo><msub><mi>Q</mi><mi>s</mi></msub><mo>+</mo><msub><mi>C</mi><mrow><mi>s</mi><mi>g</mi></mrow></msub><msub><mi>&Delta;V</mi><mi>g</mi></msub><mo>+</mo><msub><mi>C</mi><mi>s</mi></msub><msub><mi>&Delta;V</mi><mi>s</mi></msub><mo>&le;</mo><msubsup><mi>Q</mi><mi>s</mi><mi>max</mi></msubsup></mrow>]]></math><img file="FDA0000889310470000032.GIF" wi="678" he="76" /></maths><maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><mo>|</mo><msub><mi>&Delta;V</mi><mi>g</mi></msub><mo>|</mo><mo>&le;</mo><msubsup><mi>&Delta;V</mi><mi>g</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000889310470000033.GIF" wi="268" he="87" /></maths><maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><mo>|</mo><msub><mi>&Delta;V</mi><mi>s</mi></msub><mo>|</mo><mo>&le;</mo><msubsup><mi>&Delta;V</mi><mi>s</mi><mi>max</mi></msubsup></mrow>]]></math><img file="FDA0000889310470000034.GIF" wi="245" he="76" /></maths>上述目标函数中,μ<sub>s</sub>为动态无功补偿设备的无功功率矢量,表示动态无功补偿设备经控制后无功功率偏离设定值的大小,μ<sub>g</sub>为发电机的无功功率矢量,表示发电机经控制后的无功功率偏离所有发电机无功功率平均值的大小;<maths num="0014" id="cmaths0014"><math><![CDATA[<mrow><msub><mi>&mu;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>Q</mi><mrow><mi>s</mi><mi>i</mi></mrow></msub><mo>+</mo><munder><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><mi>N</mi><mi>s</mi></mrow></munder><mrow><msub><mi>C</mi><msub><mi>s</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></msub><msub><mi>&Delta;V</mi><msub><mi>s</mi><mi>j</mi></msub></msub></mrow><mo>+</mo><munder><mo>&Sigma;</mo><mrow><mi>k</mi><mo>&Element;</mo><mi>N</mi><mi>g</mi></mrow></munder><mrow><msub><mi>C</mi><mrow><msub><mi>sg</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub></mrow></msub><msub><mi>&Delta;V</mi><msub><mi>g</mi><mi>k</mi></msub></msub></mrow><mo>-</mo><msubsup><mi>Q</mi><msub><mi>s</mi><mi>i</mi></msub><mrow><mi>r</mi><mi>e</mi><mi>f</mi><mi>e</mi><mi>r</mi><mi>e</mi><mi>n</mi><mi>c</mi><mi>e</mi></mrow></msubsup></mrow><mrow><msubsup><mi>Q</mi><mrow><mi>s</mi><mi>i</mi></mrow><mi>max</mi></msubsup><mo>-</mo><msubsup><mi>Q</mi><mrow><mi>s</mi><mi>i</mi></mrow><mi>min</mi></msubsup></mrow></mfrac></mrow>]]></math><img file="FDA0000889310470000035.GIF" wi="811" he="165" /></maths><maths num="0015" id="cmaths0015"><math><![CDATA[<mrow><msub><mi>&mu;</mi><mrow><mi>g</mi><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>Q</mi><mrow><mi>g</mi><mi>i</mi></mrow></msub><mo>+</mo><munder><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><mi>N</mi><mi>g</mi></mrow></munder><mrow><msub><mi>C</mi><msub><mi>g</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></msub><msub><mi>&Delta;V</mi><msub><mi>g</mi><mi>j</mi></msub></msub></mrow><mo>+</mo><munder><mo>&Sigma;</mo><mrow><mi>k</mi><mo>&Element;</mo><mi>N</mi><mi>s</mi></mrow></munder><mrow><msub><mi>C</mi><mrow><msub><mi>gs</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub></mrow></msub><msub><mi>&Delta;V</mi><msub><mi>s</mi><mi>k</mi></msub></msub></mrow><mo>-</mo><msubsup><mi>Q</mi><mi>g</mi><mrow><mi>a</mi><mi>v</mi><mi>e</mi><mi>r</mi></mrow></msubsup></mrow><msubsup><mi>Q</mi><mi>g</mi><mrow><mi>s</mi><mi>u</mi><mi>m</mi></mrow></msubsup></mfrac></mrow>]]></math><img file="FDA0000889310470000036.GIF" wi="763" he="167" /></maths><img file="FDA0000889310470000037.GIF" wi="165" he="84" />为动态无功补偿设备的裕度设定值,<img file="FDA0000889310470000038.GIF" wi="106" he="79" />为当前时刻所有发电机的无功功率平均值,<img file="FDA0000889310470000039.GIF" wi="100" he="73" />为当前时刻发电机无功功率总和,W<sub>s</sub>为动态无功补偿设备的无功功率矢量的权重系数,W<sub>s</sub>的取值范围为0.1‑1,W<sub>g</sub>为发电机的无功功率矢量的权重系数,W<sub>g</sub>的取值范围为10‑100;将计算得到的动态无功功率裕度优化阶段的电压调整量ΔV<sub>g</sub>和ΔV<sub>s</sub>作为控制指令下达给参与控制的发电机与动态无功补偿设备,完成动态无功功率裕度优化阶段控制;(4)在时刻T<sub>0</sub>+△T,返回步骤(1)。
地址 100084 北京市海淀区清华园1号