发明名称 含风电场的注入空间电压稳定域局部边界求解方法
摘要 本发明属于电力系统技术领域,提供一种含风电场的注入空间电压稳定域局部(IVSR)边界求解方法,用于风速难以准确预测情况下电网电压稳定性的分析与在线安全监控。该方法首先,利用潮流追踪与双层调度模型获得与风电场紧密关联的调控机组,用以平衡风速变化所引起的风电场输出功率的波动;进而,利用模态分析获得关键发电机节点,用于IVSR的有效降维;最后,对可能的风速和风电场输出功率区间进行分段并行计算,并通过定向微扰来获得对应风速下的IVSR的局部边界。本发明可快速获得不同风速下IVSR的局部边界,可用于包含风电场的电力系统在线电压稳定监控,具有很好的工程应用前景。
申请公布号 CN101895130A 申请公布日期 2010.11.24
申请号 CN201010247868.3 申请日期 2010.08.08
申请人 天津大学 发明人 贾宏杰;穆云飞;余晓丹
分类号 H02J3/38(2006.01)I 主分类号 H02J3/38(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 程毓英
主权项 1.一种含风电场的注入空间电压稳定域局部边界求解方法,包括下面的步骤:第一步:确定风电场分布,并确定风电场的功率变动补偿方案(1)定义如下集合:系统风电场集合Ω<sub>WS</sub>={G<sub>w,1</sub>,G<sub>w,2</sub>,...G<sub>w,s</sub>}风电场直接紧密关联负荷集合Ω<sub>WL</sub>={L<sub>w,1</sub>,L<sub>w,2</sub>,...L<sub>w,t</sub>}与Ω<sub>WL</sub>紧密关联电源集合Ω<sub>CS</sub>={G<sub>c,1</sub>,G<sub>c,2</sub>,...G<sub>c,u</sub>}风电场非直接关联紧密负荷集合Ω<sub>OL</sub>={L<sub>o,1</sub>,L<sub>o,2</sub>,...,L<sub>o,v</sub>}风电场非直接紧密关联电源集合Ω<sub>OS</sub>={G<sub>o,1</sub>,G<sub>o,2</sub>,...G<sub>o,z</sub>},设发电机i有功输出功率P<sub>G,i</sub>派送给负荷k的派送因子为:k<sub>Gi-Lk</sub>=P<sub>Gi-Lk</sub>/P<sub>G,i</sub>,负荷k有功功率P<sub>L,k</sub>从发电机i汲取功率的汲取因子为:k<sub>Lk-Gi</sub>=P<sub>Gi-Lk</sub>/P<sub>L,k</sub>,其中:P<sub>Gi-Lk</sub>为发电机i和负荷k之间传输的有功功率的总和;(2)确定风速变化引起出力调整的补偿方案:先确定由风电场Ω<sub>WS</sub>出力调整引起的Ω<sub>CS</sub>和Ω<sub>OS</sub>的输出功率的变化量:假设负荷Ω<sub>WL</sub>、Ω<sub>OL</sub>在一定时期内保持不变,而Ω<sub>WS</sub>中第w座风电场受风速影响,其出力变化为ΔP<sub>Gw</sub>,则Ω<sub>CS</sub>中第i台发电机送往Ω<sub>WL</sub>负荷的输出功率调整量为:<img file="FDA0000024218420000011.GIF" wi="660" he="86" />i∈Ω<sub>CS</sub>,设Ω<sub>CS</sub>中第i台发电机供给Ω<sub>WL</sub>负荷比重为γ<sub>i</sub>,则供给Ω<sub>OL</sub>中负荷的比重为1-γ<sub>i</sub>,为此Ω<sub>CS</sub>中第i台发电机为弥补Ω<sub>WS</sub>中风电场功率的波动而做出的输出功率调整量为:ΔP<sub>Gi</sub>=ΔP′<sub>Gi</sub>/γ<sub>i</sub>,i∈Ω<sub>CS</sub>,其中,γ<sub>i</sub>由如下关系求得:<img file="FDA0000024218420000012.GIF" wi="399" he="86" />再求出由于Ω<sub>CS</sub>中发电机出力调整导致Ω<sub>OL</sub>中第k个负荷供电不平衡而需要Ω<sub>OS</sub>中第j台发电机输出功率调整量为:<img file="FDA0000024218420000013.GIF" wi="823" he="189" />(3)确定风电场出力随机微小波动补偿方案设定一个固定阈值,在Ω<sub>CS</sub>中选择γ<sub>i</sub>大于该固定阈值的机组来跟踪风电场随机功率波动,称其为自动发电控制AGC机组集合Ω<sub>CSA</sub>,则Ω<sub>CSA</sub>中机组i的功率调整量为:<img file="FDA0000024218420000014.GIF" wi="775" he="142" />Ω<sub>OL</sub>中的不平衡功率将由Ω<sub>OS</sub>中的平衡机来补偿;第二步:利用模态分析法选择对注入空间静态电压稳定域IVSR的边界有决定性作用的关键发电机节点,设关键发电机的有功注入向量为:P<sub>k</sub>=[P<sub>k,1</sub>,P<sub>k,2</sub>,...,P<sub>k,m</sub>];第三步:确定每个风电场的风速,进而确定各风电场的出力,形成风电场输送到电网中的有功功率P<sub>Gw</sub>及无功功率Q<sub>Gw</sub>向量;第四步:对有功功率P<sub>Gw</sub>及无功功率Q<sub>Gw</sub>向量P<sub>Gw</sub>进行双向微扰当微扰量为ΔP时,按如下方式对P<sub>Gw,i</sub>(1≤i≤s)实施微扰:<img file="FDA0000024218420000015.GIF" wi="305" he="63" /><img file="FDA0000024218420000016.GIF" wi="305" he="59" />式中:上标+和-分别表示正向微扰和负向微扰,微扰所产生的微小功率不平衡量通过Ω<sub>CSA</sub>中的机组按照第一步确定的风电场出力随机微小波动补偿方案进行平衡,Q<sub>Gw</sub>的扰动量由Q<sub>Gw</sub>和P<sub>Gw</sub>和风速之间的映射关系决定,将对P<sub>Gw</sub>和Q<sub>Gw</sub>实施微扰后系统的运行点记为<img file="FDA0000024218420000017.GIF" wi="98" he="71" />以它为初始点,确定系统在P<sub>k</sub>注入空间中的电压稳定临界点,计为<img file="FDA0000024218420000018.GIF" wi="99" he="53" />按上述方式对P<sub>Gw</sub>和Q<sub>Gw</sub>中所有的风电场依次实施正向微扰,可得正向微扰对应的临界点集:<img file="FDA0000024218420000021.GIF" wi="507" he="53" />同样,对P<sub>Gw</sub>和Q<sub>Gw</sub>所有风电场依次实施负向微扰,可得负向微扰对应的临界点集:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>X</mi><msub><mi>P</mi><mi>k</mi></msub><mo>-</mo></msubsup><mo>=</mo><mo>{</mo><msubsup><mi>x</mi><mrow><msub><mi>P</mi><mi>k</mi></msub><mo>,</mo><mn>1</mn></mrow><mrow><mo>&prime;</mo><mo>-</mo></mrow></msubsup><mo>,</mo><msubsup><mi>x</mi><mrow><msub><mi>P</mi><mi>k</mi></msub><mo>,</mo><mn>2</mn></mrow><mrow><mo>&prime;</mo><mo>-</mo></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>x</mi><mrow><msub><mi>P</mi><mi>k</mi></msub><mo>,</mo><mi>m</mi></mrow><mrow><mo>&prime;</mo><mo>-</mo></mrow></msubsup><mo>}</mo><mo>;</mo></mrow></math>]]></maths>第五步:边界超平面的求解域修正;(1)利用<img file="FDA0000024218420000023.GIF" wi="86" he="53" />通过求解m个线性方程得到IVSR边界的近似超平面:<img file="FDA0000024218420000024.GIF" wi="397" he="97" />(2)利用<img file="FDA0000024218420000025.GIF" wi="59" he="49" />可得对应的边界超平面:<img file="FDA0000024218420000026.GIF" wi="350" he="97" />(3)采用如下公式加以修正:<img file="FDA0000024218420000027.GIF" wi="154" he="99" /><img file="FDA0000024218420000028.GIF" wi="232" he="95" />,式中:α<sub>i</sub>代表超平面系数,上标”+”表示利用<img file="FDA0000024218420000029.GIF" wi="86" he="53" />通过求解m个线性方程组得到IVSR边界的近似超平面的边界系数,上标“-”表示利用<img file="FDA00000242184200000210.GIF" wi="59" he="49" />可得对应的边界超平面的系数,<img file="FDA00000242184200000211.GIF" wi="345" he="101" /><img file="FDA00000242184200000212.GIF" wi="453" he="58" />表示初始点<img file="FDA00000242184200000213.GIF" wi="63" he="45" />对应的临界点,并得到修正后的边界超平面H<sub>P</sub>为:<img file="FDA00000242184200000214.GIF" wi="258" he="101" />第六步:判断是否遇到终止条件,若是则终止计算,否则继续,终止条件是监控系统调度员人为设定的监控计算终止条件,或者当所有稳定域局部边界计算完成后的终止条件;第七步:当风电场风速变化时,首先根据第一步确定的风速变化引起出力调整的补偿方案进行补偿,然后转第三步,继续。
地址 300072 天津市南开区卫津路92号天津大学