发明名称 基于功率振荡次数评估系统阻尼比的方法
摘要 本发明公开了一种基于功率振荡次数评估系统阻尼比的方法。该方法基于广域测量系统和电力系统分析综合程序软件提供的数据生成互联电网的小干扰线性化方程,解方程得到各支路功率振荡的时域表达式;结合各支路功率振荡的时域表达式和有效扰动的判定方法求得在区域间振荡频率范围内两区域间各联络线的有效扰动次数和互联电网总有效扰动次数的理论比值与区域间阻尼比的关系式;由广域测量系统获取区域间各联络线每日有效扰动次数和互联电网总每日有效扰动次数,计算得到两者之间的实测比值;绘制关系式对应的曲线图,获取区域间阻尼比的评估值。本发明针对WAMS历史数据能够评估系统阻尼比的大小及变化情况,为运行调度人员提供参考依据。
申请公布号 CN104459409B 申请公布日期 2017.04.12
申请号 CN201410770911.2 申请日期 2014.12.12
申请人 国家电网公司;国网湖北省电力公司经济技术研究院;华中科技大学 发明人 许汉平;陈伟彪;杨东俊;刘子全;魏聪;文劲宇
分类号 G01R31/00(2006.01)I 主分类号 G01R31/00(2006.01)I
代理机构 武汉荆楚联合知识产权代理有限公司 42215 代理人 王健
主权项 一种基于功率振荡次数评估系统阻尼比的方法,其特征在于,该方法包括以下步骤:(a)在有H个区域、L条区域间联络线、N条支路、n个节点、其中发电机节点m个的互联电网中装设广域测量系统,该电网在所有支路上均装设同步相量测量单元;(b)在广域测量系统获取所有节点的电压幅值和相角以及有效扰动次数的信息,通过电力系统分析综合程序软件得到互联电网中各发电机的参数和所有线路的参数,将上述信息和参数生成互联电网的小干扰线性化方程;(c)求解小干扰线性化方程,得到各支路功率振荡的时域表达式<maths num="0001"><math><![CDATA[<mrow><msub><mi>&Delta;P</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>V</mi><mrow><mi>j</mi><mn>0</mn></mrow></msub><mi>sin</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>-</mo><msub><mi>&theta;</mi><mrow><mi>j</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow></mrow><msub><mi>X</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mfrac><msub><mi>&Delta;V</mi><mi>i</mi></msub><mo>+</mo><mfrac><mrow><msub><mi>V</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mi>sin</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>-</mo><msub><mi>&theta;</mi><mrow><mi>j</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow></mrow><msub><mi>X</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mfrac><msub><mi>&Delta;V</mi><mi>j</mi></msub><mo>+</mo><mfrac><mrow><msub><mi>V</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><msub><mi>V</mi><mrow><mi>j</mi><mn>0</mn></mrow></msub><mi>cos</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>-</mo><msub><mi>&theta;</mi><mrow><mi>j</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow></mrow><msub><mi>X</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mfrac><mrow><mo>(</mo><msub><mi>&Delta;&theta;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&Delta;&theta;</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001175553450000011.GIF" wi="1619" he="142" /></maths>其中:Δθ<sub>i</sub>和Δθ<sub>j</sub>分别表示第i节点和j个节点的电压相角变化量,ΔV<sub>i</sub>和ΔV<sub>j</sub>分别表示第i节点和j个节点的电压幅值的变化量,θ<sub>i0</sub>和θ<sub>j0</sub>分别表示第i节点和j个节点的稳态电压相角,V<sub>i0</sub>和V<sub>j0</sub>分别表示第i节点和j个节点的稳态电压幅值,X<sub>ij</sub>表示第i节点和j个节点之间的支路电抗值,上述i=1,2,…,n;j=1,2,…,n;(d)在区域A<sub>s</sub>和区域A<sub>k</sub>间有g条联络线,这两个区域间的振荡频率在[f<sub>1</sub>,f<sub>2</sub>]之间,s、k=1,2,…,H,s≠k,结合步骤(c)中各支路功率振荡的时域表达式和有效扰动的判定方法求得在[f<sub>1</sub>,f<sub>2</sub>]的频率范围内联络线l<sub>p</sub>的有效扰动次数q<sub>p</sub>和互联电网总有效扰动次数Q的理论比值Q<sub>p</sub>与区域A<sub>s</sub>和区域A<sub>k</sub>间阻尼比ξ的关系式<maths num="0002"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>p</mi></msub><mo>=</mo><msub><mi>q</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>)</mo></mrow><mo>/</mo><mi>Q</mi><mrow><mo>(</mo><mi>&xi;</mi><mo>)</mo></mrow><mo>=</mo><mo>(</mo><mrow><mn>1</mn><mo>-</mo><mfrac><mn>1</mn><msqrt><mrow><mn>2</mn><mi>&pi;</mi></mrow></msqrt></mfrac><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mi>&infin;</mi></mrow><mrow><mo>(</mo><mfrac><mrow><msub><mi>f</mi><mrow><mi>m</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow></mrow></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&mu;</mi><mn>1</mn></msub></mrow><msub><mi>&sigma;</mi><mn>1</mn></msub></mfrac><mo>)</mo></mrow></msubsup><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msup><mi>t</mi><mn>2</mn></msup><mn>2</mn></mfrac></mrow></msup><mi>d</mi><mi>t</mi></mrow><mo>)</mo><mo>/</mo><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mn>1</mn><msqrt><mrow><mn>2</mn><mi>&pi;</mi></mrow></msqrt></mfrac><mrow><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mi>&infin;</mi></mrow><mrow><mo>(</mo><mfrac><mrow><msub><mi>F</mi><mrow><mi>m</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&mu;</mi><mn>1</mn></msub></mrow><msub><mi>&sigma;</mi><mn>1</mn></msub></mfrac><mo>)</mo></mrow></msubsup><mrow><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msup><mi>t</mi><mn>2</mn></msup><mn>2</mn></mfrac></mrow></msup><mi>d</mi><mi>t</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow>]]></math><img file="FDA0001175553450000012.GIF" wi="1578" he="143" /></maths>其中:<img file="FDA0001175553450000013.GIF" wi="830" he="110" />M为有效扰动阈值,K<sub>l(p)</sub>(ξ)为第p条联络线功率振荡幅值在t=0时刻的初始值,μ<sub>1</sub>、σ<sub>1</sub>分别为冲击扰动面积F服从的正态分布的均值和方差;(e)通过广域测量系统获取历史数据中区域A<sub>s</sub>和区域A<sub>k</sub>间在[f<sub>1</sub>,f<sub>2</sub>]的频率范围内各联络线每日有效扰动次数qa<sub>w</sub>和互联电网总每日有效扰动次数Q',计算得到区域A<sub>s</sub>和区域A<sub>k</sub>间各联络线每日有效扰动次数和互联电网总每日有效扰动次数的实测比值QA<sub>w</sub>=qa<sub>w</sub>/Q',w=1,2,…,g;(f)将步骤(d)的关系式绘制成区域A<sub>s</sub>和区域A<sub>k</sub>间各联络线对应的曲线图,共g个,并根据步骤(e)求得的g个联络线每日有效扰动次数qa<sub>w</sub>和互联电网总每日有效扰动次数Q'的实测比值QA<sub>w</sub>查图得到g个比值对应的阻尼比ξ<sub>p</sub>,p=1,2,…,g;取<img file="FDA0001175553450000014.GIF" wi="256" he="125" />作为区域A<sub>s</sub>和区域A<sub>k</sub>间阻尼比的评估值。
地址 100031 北京市西城区西长安街86号