发明名称 一种变电站三相无阻抗非线性多源状态估计方法
摘要 本发明涉及一种变电站三相无阻抗非线性多源状态估计方法,属于电力系统运行和控制技术领域。首先进行量测采集,然后将变电站各电压等级各相各无阻抗开关支路形成连通的开关岛,在各开关岛内分别进行无阻抗电压状态估计和无阻抗功率状态估计,再根据功率估计值分别进行坏数据的辨识和开关状态错误辨识,用删除坏数据后的量测量重新进行状态估计。本发明的优点是利用了变电站内高度冗余的多源三相量测信息,对无阻抗开关支路建立量测方程进行三相状态估计,开关状态辨识和遥测坏数据辨识实现解耦,同时将坏数据和开关状态错误解决在变电站层,从而提高了全网状态估计的可靠性。
申请公布号 CN101958543B 申请公布日期 2012.10.10
申请号 CN201010140811.3 申请日期 2010.04.02
申请人 清华大学 发明人 孙宏斌;张伯明;吴文传;郭庆来;李青芯
分类号 H02J3/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 罗文群
主权项 1.一种变电站三相无阻抗非线性多源状态估计方法,其特征在于该方法包括以下步骤:(1)以T为采集周期,对变电站内各电压等级进行扫描,采集当前时间t<sub>(0)</sub>时变电站中各电压等级内来自数据采集与监控系统和相量测量系统的量测数据,来自数据采集与监控系统的量测数据有:变电站内各开关上流过的三相电流幅值<img file="FSB00000831357100011.GIF" wi="79" he="51" />三相有功功率量测<img file="FSB00000831357100012.GIF" wi="124" he="56" />和三相无功功率量测<img file="FSB00000831357100013.GIF" wi="147" he="56" />变电站内各节点的三相电压幅值<img file="FSB00000831357100014.GIF" wi="122" he="50" />三相注入电流幅值<img file="FSB00000831357100015.GIF" wi="84" he="59" />三相注入有功功率量测<img file="FSB00000831357100016.GIF" wi="122" he="62" />和三相注入无功功率量测<img file="FSB00000831357100017.GIF" wi="154" he="59" />以及变电站内各开关各相的开关状态量测<img file="FSB00000831357100018.GIF" wi="94" he="51" />来自相量测量系统的量测数据有:变电站内各节点的三相复电压<img file="FSB00000831357100019.GIF" wi="111" he="51" />三相注入复电流<img file="FSB000008313571000110.GIF" wi="70" he="59" />三相注入有功功率量测<img file="FSB000008313571000111.GIF" wi="125" he="59" />和三相注入无功功率量测<img file="FSB000008313571000112.GIF" wi="162" he="57" />其中,上标<img file="FSB000008313571000113.GIF" wi="206" he="46" />分别表示第<img file="FSB000008313571000114.GIF" wi="29" he="35" />相量测数据;(2)初始化相位,令<img file="2010101408113100001DEST_PATH_IMAGE002.GIF" wi="43" he="15" />(3)初始化变电站内第k个电压等级第<img file="FSB000008313571000116.GIF" wi="28" he="35" />相各开关和刀闸的状态为闭合,则变电站内第k个电压等级第<img file="FSB000008313571000117.GIF" wi="29" he="35" />相各闭合的开关支路为无阻抗开关支路;(4)定义变电站内第k个电压等级第<img file="FSB000008313571000118.GIF" wi="29" he="35" />相无阻抗开关支路形成N个连通的开关岛;(5)对变电站内第k个电压等级内第<img file="FSB000008313571000119.GIF" wi="28" he="36" />相第n个开关岛进行无阻抗电压状态估计,得出第<img file="FSB000008313571000120.GIF" wi="30" he="35" />相第n个开关岛复电压估计值<img file="FSB000008313571000121.GIF" wi="228" he="149" />上式中幅值<img file="FSB000008313571000122.GIF" wi="1000" he="301" />相角<img file="FSB000008313571000123.GIF" wi="513" he="318" />其中,m<sub>1</sub>是第<img file="FSB000008313571000124.GIF" wi="29" he="37" />相第n个开关岛内节点复电压量测<img file="FSB000008313571000125.GIF" wi="81" he="52" />的个数,m<sub>2</sub>是节点电压幅值量测<img file="FSB000008313571000126.GIF" wi="97" he="50" />的个数,<img file="FSB000008313571000127.GIF" wi="111" he="76" />和<img file="FSB000008313571000128.GIF" wi="86" he="51" />分别是上述节点复电压量测<img file="FSB000008313571000129.GIF" wi="81" he="51" />中第i个量测的幅值和相角,<img file="FSB000008313571000130.GIF" wi="102" he="58" />是上述节点电压幅值量测<img file="FSB000008313571000131.GIF" wi="92" he="59" />中第i个量测,w<sub>magli</sub>、w<sub>angi</sub>和w<sub>mag2i</sub>分别是上述<img file="FSB000008313571000132.GIF" wi="135" he="80" /><img file="FSB00000831357100021.GIF" wi="89" he="52" />和<img file="FSB00000831357100022.GIF" wi="99" he="52" />的权重;(6)设定一无阻抗电压状态估计误差阈值,定义第k个电压等级内第<img file="FSB00000831357100023.GIF" wi="28" he="37" />相第n个开关岛内电压量测向量为:<img file="FSB00000831357100024.GIF" wi="266" he="150" />其中<img file="FSB00000831357100025.GIF" wi="82" he="51" />为所述的节点复电压量测,<img file="FSB00000831357100026.GIF" wi="89" he="51" />为所述的节点电压幅值量测,若第i个电压量测<img file="FSB00000831357100027.GIF" wi="55" he="57" />与上述开关岛复电压估计值<img file="FSB00000831357100028.GIF" wi="58" he="51" />之间正则化残差值最大且该差值大于上述设定的无阻抗电压状态估计误差阈值,则该量测<img file="FSB00000831357100029.GIF" wi="55" he="58" />为坏数据,删除该坏数据,重复步骤(5)-(6),直到各电压量测向量中没有坏数据,且第n个开关岛电压估计幅值为上述<img file="FSB000008313571000210.GIF" wi="84" he="51" />相角为上述<img file="FSB000008313571000211.GIF" wi="78" he="50" />(7)计算上述第k个电压等级的第<img file="FSB000008313571000212.GIF" wi="30" he="39" />相第n个开关岛内各节点注入视在功率平方伪量测<img file="FSB000008313571000213.GIF" wi="77" he="66" />及各开关支路视在功率平方伪量测<img file="FSB000008313571000214.GIF" wi="96" he="59" />视在功率平方伪量测的求解如下:<img file="FSB000008313571000215.GIF" wi="339" he="145" />其中,<img file="FSB000008313571000216.GIF" wi="54" he="59" />为第<img file="FSB000008313571000217.GIF" wi="28" he="36" />相第n个开关岛内各节点注入电流幅值量测,<img file="FSB000008313571000218.GIF" wi="50" he="52" />为第<img file="FSB000008313571000219.GIF" wi="28" he="37" />相第n个开关岛内各开关支路电流幅值量测,<img file="FSB000008313571000220.GIF" wi="77" he="51" />为上述第n个开关岛电压幅值估计值;(8)以第<img file="FSB000008313571000221.GIF" wi="30" he="36" />相第n个开关岛内无阻抗开关支路上有功功率<img file="FSB000008313571000222.GIF" wi="51" he="52" />和无功功率<img file="FSB000008313571000223.GIF" wi="51" he="58" />为状态量<img file="FSB000008313571000224.GIF" wi="235" he="146" />设上述第k个电压等级内第<img file="FSB000008313571000225.GIF" wi="29" he="37" />相第n个开关岛内状态量初值<img file="FSB000008313571000226.GIF" wi="83" he="45" />为相应的各无阻抗开关支路上有功功率和无功功率量测,即<img file="FSB000008313571000227.GIF" wi="537" he="147" />定义第<img file="FSB000008313571000228.GIF" wi="29" he="36" />相第n个开关岛内线性量测向量为:<img file="FSB000008313571000229.GIF" wi="307" he="450" />非线性量测向量为:<img file="FSB000008313571000230.GIF" wi="254" he="149" />则上述第k个电压等级第<img file="FSB000008313571000231.GIF" wi="30" he="36" />相第n个开关岛内量测向量<img file="FSB000008313571000232.GIF" wi="216" he="146" />(9)设迭代次数的初始值count=0,将上述<img file="FSB000008313571000233.GIF" wi="53" he="48" />作为已知量,<img file="FSB000008313571000234.GIF" wi="87" he="51" />作为状态量初值,对各无阻抗开关支路建立量测方程进行无阻抗功率状态估计,得到第k个电压等级内第<img file="FSB000008313571000235.GIF" wi="29" he="35" />相第n个开关岛内各无阻抗开关支路上功率估计值<img file="FSB00000831357100031.GIF" wi="103" he="47" />无阻抗功率状态估计迭代方程如下:<img file="FSB00000831357100032.GIF" wi="822" he="238" />其中<img file="FSB00000831357100033.GIF" wi="78" he="39" />是无阻抗功率状态估计迭代修正量,W是变电站自动化系统中量测的权重矩阵,<img file="FSB00000831357100034.GIF" wi="187" he="57" />为非线性部分量测函数,雅可比矩阵<img file="FSB00000831357100035.GIF" wi="213" he="155" />通过量测方程得到,其中线性部分量测方程为:<img file="FSB00000831357100036.GIF" wi="1132" he="463" />上式中A<sub>KCL</sub>是第k个电压等级内第<img file="FSB00000831357100037.GIF" wi="29" he="35" />相第n个开关岛无阻抗电网中的节点-支路关联矩阵,I为单位矩阵,r<sub>l</sub>为线性部分量测误差矢量,非线性部分量测方程为:<img file="FSB00000831357100038.GIF" wi="335" he="109" />上式中,r<sub>nl</sub>是非线性部分量测的误差矢量,<img file="FSB00000831357100039.GIF" wi="138" he="51" />为非线性部分量测函数:<img file="FSB000008313571000310.GIF" wi="894" he="177" />上式中,<img file="FSB000008313571000311.GIF" wi="127" he="57" />分别为状态量<img file="FSB000008313571000312.GIF" wi="61" he="39" />的有功功率部分和无功功率部分,P<sub>inj</sub>,Q<sub>inj</sub>分别为由开关上功率状态量<img file="FSB000008313571000313.GIF" wi="48" he="40" />求得的各节点注入有功功率和无功功率:<img file="FSB000008313571000314.GIF" wi="353" he="165" />计算得到非线性部分量测方程的雅可比矩阵H<sub>nl</sub>:<img file="FSB000008313571000315.GIF" wi="617" he="111" />(10)使迭代次数count=count+1,设定一个无阻抗功率状态估计迭代精度ε,当上述状态估计迭代修正量<img file="FSB000008313571000316.GIF" wi="156" he="46" />时,重复步骤(9)-(10),当<img file="FSB000008313571000317.GIF" wi="153" he="41" />时,迭代结束,设 定一迭代次数门槛值countmax,若迭代次数count<countmax,状态估计收敛,得到第<img file="FSB00000831357100041.GIF" wi="29" he="37" />相第n个开关岛内各无阻抗开关支路上功率状态量最终估计结果<img file="FSB00000831357100042.GIF" wi="339" he="145" />并进入步骤(11),若count≥countmax,状态估计发散,重复步骤(1)-(10);(11)根据上述计算结果,求解式<img file="FSB00000831357100043.GIF" wi="357" he="145" />以得到估计后的第<img file="FSB00000831357100044.GIF" wi="29" he="36" />相量测量的估计值<img file="FSB00000831357100045.GIF" wi="71" he="51" />设定一个无阻抗功率状态估计误差阈值,将上述第k个电压等级第<img file="FSB00000831357100046.GIF" wi="30" he="43" />相第n个开关岛内量测向量<img file="FSB00000831357100047.GIF" wi="42" he="42" />中第i个量测值<img file="FSB00000831357100048.GIF" wi="44" he="52" />与上述估计后的第<img file="FSB00000831357100049.GIF" wi="29" he="36" />相第i个开关上的量测量估计值<img file="FSB000008313571000410.GIF" wi="56" he="52" />进行比较,若两者正则化残差值最大且该差值大于设定的无阻抗功率状态估计误差阈值,则该量测<img file="FSB000008313571000411.GIF" wi="44" he="50" />为坏数据,删除该坏数据,置该量测的权重w<sub>i</sub>=0,重复步骤(9)-(11),直到各量测量中没有坏数据;(12)设定一个估计阈值,将上述第k个电压等级内第<img file="FSB000008313571000412.GIF" wi="31" he="37" />相第n个开关岛内各无阻抗开关支路上的视在功率估计值<img file="FSB000008313571000413.GIF" wi="458" he="81" />与估计阈值进行比较,若第i个开关视在功率估计值x<sub>cbi</sub>大于估计阈值,则该开关的估计状态为闭合状态,若第i个开关视在功率估计值x<sub>cbi</sub>小于估计阈值,则该开关的估计状态为断开状态;将比较得到的第i个开关的估计状态与上述第i个开关开合状态量测<img file="FSB000008313571000414.GIF" wi="80" he="52" />进行比较,若开关的估计状态与上述开合状态量测值不一致,则该开关的开合状态量测<img file="FSB000008313571000415.GIF" wi="81" he="52" />为坏数据,将开合状态量测<img file="FSB000008313571000416.GIF" wi="80" he="52" />替代为开关的估计状态,若开关的估计状态与上述开合状态量测一致,则保持开关状态量测不变,重复步骤(5)-(12),直到所有开关岛状态估计完成则转入步骤(13);(13)根据上述开关岛状态估计结果,重新确定第k个电压等级第<img file="FSB000008313571000417.GIF" wi="29" he="36" />相内各开关状态,开关闭合的支路为无阻抗开关支路,重复步聚(4)-(13),重新形成开关岛并进行第二次状态估计,完成第二次状态估计后,进入步骤(14);(14)上述步骤完成对<img file="FSB000008313571000418.GIF" wi="34" he="34" />相状态估计后,令<img file="FSB000008313571000419.GIF" wi="168" he="48" />重复步骤(3)-(14),对b,c两相分别进行状态估计,完成三相状态估计后,转入步骤(15);(15)当此电压等级状态估计完成之后,进入下一个电压等级进行状态估计,重复步骤(2)-(15),直到完成所有电压等级状态估计;(16)将当前采集时间t<sub>(1)</sub>与t<sub>(0)</sub>进行比较,当t<sub>(1)</sub>-t<sub>(0)</sub>大于上述采集周期T时,重复步骤(1)-(16)。 
地址 100084 北京市海淀区清华园1号