发明名称 一种基于全微分方程的戴维南等值参数跟踪的计算方法
摘要 本发明提供一种求解电力系统戴维南等值参数的方法,将潮流方程的有功功率和无功功率分别对戴维南等值参数求导形成全微分方程组,通过求解以戴维南等值参数作为变量的全微分方程组,得到戴维南等值参数。本发明克服了现有计算方法对两个时步之间戴维南等值参数不变进行假设所带来的严重误差,具有适应性强,使用简单,计算量小的特点,可以应用于需要形成戴维南等值系统的电力系统分析计算和现场装置。
申请公布号 CN101350524B 申请公布日期 2011.08.31
申请号 CN200810222191.0 申请日期 2008.09.11
申请人 中国电力科学研究院 发明人 汤涌;孙华东;易俊
分类号 H02J3/00(2006.01)I;H02J3/16(2006.01)I;G06F17/13(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 北京安博达知识产权代理有限公司 11271 代理人 徐国文
主权项 1.一种基于全微分方程的戴维南等值电阻、电抗及电势的跟踪计算方法,其特征在于:将潮流方程的有功功率和无功功率分别对戴维南等值阻抗及电势求导形成全微分方程组,通过求解全微分方程组得到戴维南等值参数,包括下列步骤:步骤A:监测戴维南等值负荷母线电压幅值、负荷有功功率和无功功率,并将其作为每一个计算时步的已知量;步骤B:将初始时刻的戴维南等值阻抗R<sub>0</sub>+jX<sub>0</sub>作为求取后续一系列时刻戴维南等值阻抗及电势的初始值;步骤C:利用初始时刻的戴维南等值阻抗R<sub>0</sub>+jX<sub>0</sub>以及初始时刻的电压幅值、有功功率和无功功率值,计算得到第一个时间序列的戴维南等值阻抗及电势,利用得到的第一个时间序列的戴维南等值阻抗值,计算求取第二个时间序列的戴维南等值阻抗及电势,以此类推计算得到随着时间变化的一系列戴维南等值电势E<sub>rk</sub>+jE<sub>ik</sub>和戴维南等值阻抗R<sub>k</sub>+jX<sub>k</sub>;在所述步骤C中,提出的全微分方程法的推导过程如下:各项电气参数有如下关系:<img file="FSB00000303641500011.GIF" wi="1183" he="108" />式(1)中,<img file="FSB00000303641500012.GIF" wi="34" he="46" />为戴维南等值电势,<img file="FSB00000303641500013.GIF" wi="32" he="47" />为负荷母线端电压,P+jQ为负荷有功功率和无功功率,R+jX为戴维南等值阻抗, 为<img file="FSB00000303641500015.GIF" wi="32" he="47" />的共轭,不失一般性,令<img file="FSB00000303641500016.GIF" wi="215" he="48" />将上式实部和虚部分开得:V<sup>2</sup>-V×E<sub>r</sub>+P×R+Q×X=0     (2)V×E<sub>i</sub>-P×X+Q×R=0        (3)式(2)、(3)中,E<sub>r</sub>、E<sub>i</sub>分别为<img file="FSB00000303641500017.GIF" wi="35" he="47" />的实部和虚部;由(2)、(3)两式可以解得:<img file="FSB00000303641500018.GIF" wi="1238" he="116" /><img file="FSB00000303641500021.GIF" wi="1293" he="240" />用(4)、(5)两式求P、Q的全微分,得:<img file="FSB00000303641500022.GIF" wi="968" he="139" /><img file="FSB00000303641500023.GIF" wi="1022" he="104" />(6)<img file="FSB00000303641500025.GIF" wi="1019" he="121" /><img file="FSB00000303641500026.GIF" wi="901" he="123" /><img file="FSB00000303641500027.GIF" wi="1021" he="131" /><img file="FSB00000303641500028.GIF" wi="1024" he="110" />(7)<img file="FSB000003036415000210.GIF" wi="942" he="122" /><img file="FSB000003036415000211.GIF" wi="1030" he="121" />式(6)和(7)可以写为:<img file="FSB000003036415000212.GIF" wi="1005" he="117" /><img file="FSB000003036415000213.GIF" wi="1041" he="104" />(8)<img file="FSB000003036415000215.GIF" wi="1025" he="121" /><img file="FSB000003036415000216.GIF" wi="895" he="122" /><img file="FSB000003036415000217.GIF" wi="1017" he="117" /><img file="FSB000003036415000218.GIF" wi="1041" he="103" />(9)<img file="FSB000003036415000220.GIF" wi="949" he="123" /><img file="FSB000003036415000221.GIF" wi="1163" he="121" />即:<img file="FSB00000303641500031.GIF" wi="1739" he="170" /><img file="FSB00000303641500032.GIF" wi="1337" he="129" />(10)<img file="FSB00000303641500034.GIF" wi="1230" he="128" /><img file="FSB00000303641500035.GIF" wi="955" he="75" /><img file="FSB00000303641500036.GIF" wi="1765" he="138" /><img file="FSB00000303641500037.GIF" wi="1276" he="131" />(11)<img file="FSB00000303641500039.GIF" wi="1373" he="128" /><img file="FSB000003036415000310.GIF" wi="952" he="66" />于是,由公式(2)、(3)、(10)、(11)可以得到如下包含六个方程的方程组:<img file="FSB000003036415000311.GIF" wi="1665" he="546" />式(12)中,有<img file="FSB000003036415000312.GIF" wi="676" he="59" />八个未知数,因此必须首先给定其中任意两个未知数的值,或者再给出两个方程,才能得到所有未知数的解,可利用全网信息计算求得初始时刻的戴维南等值阻抗R<sub>0</sub>+jX<sub>0</sub>,然后就可以利用式(12)求得初始时刻以后各个时刻的系统戴维南等值参数。
地址 100192 北京市海淀区清河小营东路15号中国电力科学研究院科技部
您可能感兴趣的专利