发明名称 一种基于广域量测的电力系统动态稳定性评估方法
摘要 本发明公开了一种基于广域量测的电力系统动态稳定性评估方法,包括以下步骤:建立电压稳定指标;监测主站,收集拓扑结构变化信息,将数据通过广域通信网发送给监测子站;各节点的监测子站接收电源管理单元测量的该节点的母线电压和源节点电压信息;计算系统电压稳定指标;本发明通过预测负荷节点最大传输功率和监测系统运行点的运行区域,提出了一种监测电压稳定的局部指标方法,并利用广域测量系统,将它应用于在线电压稳定监测,克服了基于参数辨识的指标不能对联络节点进行电压稳定评估及线路指标在线应用时等值电源的电压幅值相对于被监控母线的电压幅值有强耦合的缺陷。
申请公布号 CN103595043B 申请公布日期 2016.06.08
申请号 CN201310612194.6 申请日期 2013.11.26
申请人 国家电网公司;国网山东省电力公司电力科学研究院;国网山东节能服务有限公司 发明人 辛卫东;汪东军
分类号 H02J3/00(2006.01)I;H02J13/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 济南圣达知识产权代理有限公司 37221 代理人 张勇
主权项 一种基于广域量测的电力系统动态稳定性评估方法,其特征是:包括以下步骤:(1)建立电压稳定指标;(2)监测主站,收集拓扑结构变化信息,将数据通过广域通信网发送给监测子站;(3)各节点的监测子站接收电源管理单元测量的该节点的母线电压和源节点电压信息;(4)计算系统电压稳定指标:各监测子站将计算的电压稳定指标通过广域通信网传递给监测主站,监测主站计算系统电压稳定指标;所述步骤(3)具体为:节点j的监测子站接收到电源管理单元测量到的节点j母线电压及源节点电压信息后,按单电源功率传输等值系统模型计算函数sgn x<sub>j</sub>、功率极限Pmax<sub>j</sub>、Qmax<sub>j</sub>和Smax<sub>j</sub>、功率裕度Pmargin<sub>j</sub>、Qmargin<sub>j</sub>和Smargin<sub>j</sub>及电压稳定功率裕度指标;所述步骤(1)建立电压稳定指标的具体为:建立电力系统传输功率模型:将每条电力系统的电力传输线路简化成简单电力系统传输功率模型,其中,U<sub>s</sub>∠θ为电源端电压,线路阻抗为Z<sub>L</sub>=R+jX,负荷侧电压为U<sub>r</sub>∠0,负荷功率为<img file="FDA0000845183690000017.GIF" wi="252" he="74" />则,符合有功功率和无功功率可以分别表示为:<maths num="0001"><math><![CDATA[<mrow><mi>P</mi><mo>=</mo><msub><mi>U</mi><mi>r</mi></msub><mo>&lsqb;</mo><mrow><mo>(</mo><msub><mi>U</mi><mi>S</mi></msub><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi><mo>-</mo><msub><mi>U</mi><mi>r</mi></msub><mo>)</mo></mrow><mfrac><mi>R</mi><mrow><msup><mi>R</mi><mn>2</mn></msup><mo>+</mo><msup><mi>X</mi><mn>2</mn></msup></mrow></mfrac><mo>+</mo><msub><mi>U</mi><mi>S</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi><mfrac><mi>X</mi><mrow><msup><mi>R</mi><mn>2</mn></msup><mo>+</mo><msup><mi>X</mi><mn>2</mn></msup></mrow></mfrac><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000845183690000011.GIF" wi="1678" he="102" /></maths><maths num="0002"><math><![CDATA[<mrow><mi>Q</mi><mo>=</mo><msub><mi>U</mi><mi>r</mi></msub><mo>&lsqb;</mo><mrow><mo>(</mo><msub><mi>U</mi><mi>S</mi></msub><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi><mo>-</mo><msub><mi>U</mi><mi>r</mi></msub><mo>)</mo></mrow><mfrac><mi>X</mi><mrow><msup><mi>R</mi><mn>2</mn></msup><mo>+</mo><msup><mi>X</mi><mn>2</mn></msup></mrow></mfrac><mo>+</mo><msub><mi>U</mi><mi>S</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi><mfrac><mi>R</mi><mrow><msup><mi>R</mi><mn>2</mn></msup><mo>+</mo><msup><mi>X</mi><mn>2</mn></msup></mrow></mfrac><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000845183690000012.GIF" wi="1686" he="103" /></maths>联合两式,消除cosθ和sinθ,可以解得:<maths num="0003"><math><![CDATA[<mrow><msup><msub><mi>U</mi><mi>r</mi></msub><mn>2</mn></msup><mo>=</mo><mrow><mo>(</mo><mi>&beta;</mi><mo>&PlusMinus;</mo><msqrt><mrow><msup><mi>&beta;</mi><mn>2</mn></msup><mo>-</mo><mn>1</mn></mrow></msqrt><mo>)</mo></mrow><msqrt><mrow><mo>(</mo><msup><mi>P</mi><mn>2</mn></msup><mo>+</mo><msup><mi>Q</mi><mn>2</mn></msup><mo>)</mo><mo>(</mo><msup><mi>R</mi><mn>2</mn></msup><mo>+</mo><msup><mi>X</mi><mn>2</mn></msup><mo>)</mo></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000845183690000013.GIF" wi="1429" he="87" /></maths>式中,<maths num="0004"><math><![CDATA[<mrow><mi>&beta;</mi><mo>=</mo><mfrac><mrow><msup><msub><mi>U</mi><mi>S</mi></msub><mn>2</mn></msup><mo>-</mo><mn>2</mn><mi>P</mi><mi>R</mi><mo>-</mo><mn>2</mn><mi>Q</mi><mi>X</mi></mrow><mrow><mn>2</mn><msqrt><mrow><mo>(</mo><msup><mi>P</mi><mn>2</mn></msup><mo>+</mo><msup><mi>Q</mi><mn>2</mn></msup><mo>)</mo><mo>(</mo><msup><mi>R</mi><mn>2</mn></msup><mo>+</mo><msup><mi>X</mi><mn>2</mn></msup><mo>)</mo></mrow></msqrt></mrow></mfrac><mo>&GreaterEqual;</mo><mn>1</mn><mo>,</mo></mrow>]]></math><img file="FDA0000845183690000014.GIF" wi="694" he="165" /></maths>由于<img file="FDA0000845183690000015.GIF" wi="191" he="109" />则当β&gt;1时,式(3)中的高值解对应PV曲线的上半支解,低值解对应PV曲线的下半支解;而当β&lt;1时,系统没有平衡点,当β=1时系统位于PV曲线的拐点;定义符号函数sgn x为:<maths num="0005"><math><![CDATA[<mrow><mi>sgn</mi><mi> </mi><mi>x</mi><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>x</mi><mo>&gt;</mo><mn>1</mn></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>x</mi><mo>&lt;</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000845183690000016.GIF" wi="1655" he="253" /></maths>式中,<maths num="0006"><math><![CDATA[<mrow><mi>x</mi><mo>=</mo><mfrac><mrow><msup><msub><mi>U</mi><mi>r</mi></msub><mn>2</mn></msup></mrow><mrow><msqrt><mrow><msup><mi>P</mi><mn>2</mn></msup><mo>+</mo><msup><mi>Q</mi><mn>2</mn></msup></mrow></msqrt><mo>&CenterDot;</mo><msqrt><mrow><msup><mi>R</mi><mn>2</mn></msup><mo>+</mo><msup><mi>X</mi><mn>2</mn></msup></mrow></msqrt></mrow></mfrac></mrow>]]></math><img file="FDA0000845183690000021.GIF" wi="536" he="164" /></maths>函数sgnx给出了运行点是否工作在PV曲线上半支或下半支的信息;取负荷裕度Pmargin、Qmargin和Smargin分别为Pmargin=Pmax‑P   (5)Qmargin=Qmax‑Q   (6)Smargin=Smax‑S   (7)其中:<img file="FDA0000845183690000022.GIF" wi="350" he="87" />Pmax为无功功率Q固定,增加有功功率P时,PV曲线拐点的功率;Qmax为有功功率P固定,增加负荷无功功率Q时,QV曲线拐点的功率值;Smax为恒定功率因数条件下,增加负荷S时,PV曲线拐点的负荷总量;由系统达到功率传输极限的条件(β=1)及式(5)~(7)可知:在系统功率传输极限点,负荷裕度Pmargin、Qmargin和Smargin同时等于0;因此,度量运行点偏离系统功率传输极限点的距离用其中最小值给出,定义电压稳定指标voltage instability index,VSI为:<maths num="0007"><math><![CDATA[<mrow><msub><mi>I</mi><mrow><mi>V</mi><mi>S</mi><mi>I</mi></mrow></msub><mo>=</mo><mi>sgn</mi><mi> </mi><mi>x</mi><mo>&CenterDot;</mo><mi>min</mi><mrow><mo>(</mo><mrow><mo>|</mo><mfrac><msub><mi>P</mi><mrow><mi>m</mi><mi>arg</mi><mi>i</mi><mi>n</mi></mrow></msub><msub><mi>P</mi><mi>max</mi></msub></mfrac><mo>|</mo><mo>,</mo><mo>|</mo><mfrac><msub><mi>Q</mi><mrow><mi>m</mi><mi>arg</mi><mi>i</mi><mi>n</mi></mrow></msub><msub><mi>Q</mi><mi>max</mi></msub></mfrac><mo>|</mo><mo>,</mo><mo>|</mo><mfrac><msub><mi>S</mi><mrow><mi>m</mi><mi>arg</mi><mi>i</mi><mi>n</mi></mrow></msub><msub><mi>S</mi><mi>max</mi></msub></mfrac><mo>|</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000845183690000023.GIF" wi="1596" he="182" /></maths>
地址 250002 山东省济南市市中区二环南路1号