发明名称 变电站-调度中心两级分布式电网的非线性状态估计方法
摘要 本发明涉及一种变电站-调度中心两级分布式电网的非线性状态估计方法,属于电力系统运行和控制技术领域。首先在变电站中利用来自PMU、RTU的多源量测,采用变电站三相无阻抗非线性多源状态估计方法进行变电站状态估计,同步辨识并剔除变电站内开关开合状态坏数据和模拟量坏数据,获得各变电站的高可靠性的三相拓扑结构、高精度的三相状态估计结果及三相不平衡度并实时上传到调度中心,同时没有参与变电站状态估计的变电站选择部分原始量测量上传到调度中心,调度中心利用上传的正序或单相量对电网进行拓扑分析及非线性单相状态估计,有效监视全网三相不平衡程度和非全相运行的情况,并实现了传统状态估计到两级分布式状态估计的过渡。
申请公布号 CN101924364B 申请公布日期 2013.05.01
申请号 CN201010237485.8 申请日期 2010.07.23
申请人 清华大学 发明人 孙宏斌;张伯明;吴文传;郭庆来;李青芯
分类号 H02J3/00(2006.01)I;H02J13/00(2006.01)I;G01R21/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 罗文群
主权项 1.一种变电站-调度中心两级分布式电网的非线性状态估计方法,其特征在于该方法包括以下步骤: (1)设T1为电网的变电站的状态估计周期,T2为电网的调度中心的状态估计周期,T1<T2;采集各变电站中各电压等级下来自电网的远方终端单元和相量量测单元的多源量测数据,来自远方终端单元的量测数据为:变电站内的各三相电压幅值、三相电流幅值、三相有功功率、三相无功功率以及变电站内各开关各相的开关开合状态,来自相量量测单元的量测数据为:变电站内的各三相复电压、三相复电流、三相有功功率和三相无功功率; (2)设电网变电站的各电压等级下各开关三相为闭合,形成由连通的无阻抗开关支路组成的开关岛,在各开关岛内分别进行无阻抗电压状态估计和无阻抗功率状态估计,并根据无阻抗功率状态估计结果分别进行模拟量坏数据的辨识和开关开合状态坏数据辨识,剔除坏数据,根据新的开关开合状态,形成变电站内各电压等级下的新开关岛,设新开关岛有M个,则<img file="FSA00000205924400011.GIF" wi="244" he="120" />其中,K为该变电站内电压等级个数,N<sub>k</sub>为第k个电压等级下开关岛个数,对各新开关岛进行状态估计,得到变电站内各电压等级下各开关的三相开合状态<img file="FSA00000205924400012.GIF" wi="89" he="53" />各节点的三相复电压<img file="FSA00000205924400013.GIF" wi="94" he="59" />各开关上三相有功功率<img file="FSA00000205924400014.GIF" wi="65" he="52" />和无功功率<img file="FSA00000205924400015.GIF" wi="71" he="52" />以及各节点三相注入有功功率<img file="FSA00000205924400016.GIF" wi="248" he="144" />和无功功率<img file="FSA00000205924400017.GIF" wi="284" he="144" />其中<img file="FSA00000205924400018.GIF" wi="209" he="53" />分别为与变电站内的变压器、线路及容抗器相连的出线节点的三相注入有功功率和无功功率,<img file="FSA00000205924400019.GIF" wi="166" he="51" />分别为除出线节点以外的物理节点的三相注入有功功率和无功功率,上标<img file="FSA000002059244000110.GIF" wi="191" he="48" />表示三相;(3)以步骤(2)的变电站内各电压等级下各新开关岛分别为一条母线,则共有M条母线,母线的复电压为<img file="FSA000002059244000111.GIF" wi="95" he="60" />通过求解<img file="FSA000002059244000112.GIF" wi="302" he="173" />得到变电站内第m个新开关岛的母线复电压<img file="FSA000002059244000113.GIF" wi="117" he="61" />其中,<img file="FSA000002059244000114.GIF" wi="86" he="55" />为上述各节点三相复电压<img file="FSA000002059244000115.GIF" wi="68" he="59" />的第i个值,n为第m个新开关岛内的节点数,则母线的a、b线电压相量为<img file="FSA000002059244000116.GIF" wi="94" he="60" /><img file="FSA000002059244000117.GIF" wi="344" he="61" /><img file="FSA000002059244000118.GIF" wi="71" he="60" />为母线a相电压相量,<img file="FSA000002059244000119.GIF" wi="72" he="61" />为母线b相电压相量;(4)根据步骤(2)的变电站三相状态估计结果,通过求解<img file="FSA00000205924400021.GIF" wi="408" he="142" />得到第m个新开关岛内各开关上三相复电流<img file="FSA00000205924400022.GIF" wi="89" he="61" />通过求解<img file="FSA00000205924400023.GIF" wi="425" he="143" />得到第m个新开关岛内各节点的三相注入复电流<img file="FSA00000205924400024.GIF" wi="266" he="143" />其中<img file="FSA00000205924400025.GIF" wi="63" he="60" />为与变电站内的变压器、线路及容抗器相连的出线节点的三相注入复电流,<img file="FSA00000205924400026.GIF" wi="56" he="60" />为除出线节点以外的物理节点的三相注入复电流,出线节点a相的电流相量为<img file="FSA00000205924400027.GIF" wi="87" he="60" />(5)根据步骤(3)中变电站内各母线三相复电压<img file="FSA00000205924400028.GIF" wi="95" he="60" />得到变电站内各母线三相电压的正序分量<img file="FSA00000205924400029.GIF" wi="45" he="60" />和负序分量<img file="FSA000002059244000210.GIF" wi="72" he="60" /><img file="FSA000002059244000211.GIF" wi="543" he="254" />其中,<img file="FSA000002059244000212.GIF" wi="452" he="124" />计算得到变电站内各母线的三相复电压正序分量幅值|U<sub>1</sub>|和负序分量幅值|U<sub>2</sub>|:<img file="FSA000002059244000213.GIF" wi="445" he="193" />其中,U<sub>1real</sub>和U<sub>1imag</sub>为变电站内各母线三相复电压正序分量<img file="FSA000002059244000214.GIF" wi="44" he="61" />的实部和虚部,U<sub>2real</sub>和U<sub>2imag</sub>为负序分量<img file="FSA000002059244000215.GIF" wi="50" he="60" />的实部和虚部,通过求解<img file="FSA000002059244000216.GIF" wi="330" he="138" />得到变电站内各母线的三相电压不平衡度;(6)由步骤(4)变电站内各出线节点的三相注入复电流<img file="FSA000002059244000217.GIF" wi="87" he="60" />计算得到各出线节点三相电流的正序分量<img file="FSA000002059244000218.GIF" wi="32" he="60" />和负序分量<img file="FSA000002059244000219.GIF" wi="57" he="60" /><img file="FSA000002059244000220.GIF" wi="552" he="255" />其中,<img file="DEST_PATH_FSB00000751738000011.GIF" wi="447" he="123" />计算得到各出线节点三相电流的正序分量幅值<img file="DEST_PATH_FSB00000751738000012.GIF" wi="63" he="78" />和负序分量幅值<img file="DEST_PATH_FSB00000751738000013.GIF" wi="89" he="89" /><img file="DEST_PATH_FSB00000751738000014.GIF" wi="487" he="256" />其中,I<sub>1real</sub>和I<sub>1imag</sub>分别为正序分量<img file="DEST_PATH_FSB00000751738000015.GIF" wi="33" he="62" />的实部和虚部,I<sub>2real</sub>和I<sub>2imag</sub>分别为负序分量<img file="DEST_PATH_FSB00000751738000016.GIF" wi="38" he="68" />的实部和虚部,通过求解<img file="DEST_PATH_FSB00000751738000017.GIF" wi="329" he="165" />得到变电站内各出线节点三相电流不平衡度;(7)通过如下公式:<img file="DEST_PATH_FSB00000751738000018.GIF" wi="641" he="407" />计算变电站内各开关三相总有功功率P<sub>CB</sub>和三相总无功功率Q<sub>CB</sub>,以及变电站内各节点的三相注入总有功功率<img file="DEST_PATH_FSB00000751738000019.GIF" wi="243" he="139" />和三相总无功功率<img file="DEST_PATH_FSB000007517380000110.GIF" wi="283" he="167" />其中P<sub>LN</sub>、Q<sub>LN</sub>分别为变电站内各出线节点的三相总有功功率和无功功率,P<sub>nd</sub>、Q<sub>nd</sub>分别为除出线节点以外的物理节点的三相总有功功率和三相总无功功率;(8)计算变电站内各开关的三相总开合状态z<sub>CB</sub>,其中第i个开关的三相总开合状态计算公式如下:<img file="DEST_PATH_FSB000007517380000111.GIF" wi="760" he="267" />其中,0表示开关的三相都闭合,1表示开关的三相都断开,2表示开关的三相开合状态不一致,处于报警状态,并向调度中心发送非全相运行状态的报警信号;(9)当到达电网的变电站的状态估计周期T1时,重复步骤(1)-(9),进行周期性变电站状态估计,当到达电网的调度中心的状态估计周期T2时,进行以下步骤;(10)电网的调度中心从各变电站采集量测量熟数据z<sub>se</sub>,其中,当调度中心需要正序熟数据对电网进行状态估计时,<img file="FSA00000205924400041.GIF" wi="587" he="83" />其中,<img file="FSA00000205924400042.GIF" wi="47" he="59" />表示正序熟数据,当调度中心需要单相熟数据对电网进行状态估计时,<img file="FSA00000205924400043.GIF" wi="649" he="81" /><img file="FSA00000205924400044.GIF" wi="55" he="59" />表示单相熟数据,同时调度中心从各变电站采集以下变电站状态估计结果:变电站的各母线三相电压不平衡度ε<sub>U</sub>、各出线节点三相电流不平衡度ε<sub>I</sub>以及各开关的三相总开合状态z<sub>CB</sub>;(11)调度中心采集电网中不参与变电站状态估计的各变电站原始量测量z<sub>meas</sub>,其中,当调度中心需要正序量测量对电网进行状态估计时,<img file="FSA00000205924400045.GIF" wi="1028" he="85" />其中,<img file="FSA00000205924400046.GIF" wi="83" he="58" />表示正序量测量,|U′<sub>1</sub>|为远方终端单元采集到的变电站内各电气母线电压幅值,|I′<sub>1</sub>|为远方终端单元采集到的各出线节点的电流幅值,P′、Q′分别为远方终端单元采集到的各出线节点的三相总的有功功率和无功功率,<img file="FSA00000205924400047.GIF" wi="48" he="60" />为相量量测单元采集到的各电气母线正序复电压,<img file="FSA00000205924400048.GIF" wi="37" he="60" />分别为相量量测单元采集到的各出线节点的正序复电流,P′<sub>PMU</sub>、Q′<sub>PMU</sub>分别为相量量测单元采集到的各出线节点的三相总有功功率和无功功率,当调度中心需要单相熟数据对电网进行状态估计时,<img file="FSA00000205924400049.GIF" wi="1078" he="84" />其中,<img file="FSA000002059244000410.GIF" wi="83" he="58" />表示单相量测量,|U′<sub>ab</sub>|为远方终端单元采集到的变电站内各电气母线a、b线电压幅值,|I′<sub>a</sub>|为远方终端单元采集到的变电站内各出线节点的a相电流幅值,<img file="FSA000002059244000411.GIF" wi="64" he="60" />为相量量测单元采集到的变电站内各电气母线a、b线电压相量,<img file="FSA000002059244000412.GIF" wi="40" he="60" />为相量量测单元采集到的变电站内各出线节点的a相电流相量;(12)当调度中心的第i个开关的三相总开合状态z<sub>CB,i</sub>等于2时,即为非全相运行状态时,处于报警状态;设置一个三相电流不平衡度阈值,当变电站内各出线节点的三相复电流不平衡度ε<sub>I</sub>大于此阈值时,处于报警状态;设置一个三相电压不平衡度阈值,当变电站内各母线的三相电压不平衡度ε<sub>U</sub>大于此阈值时,处于报警状态;利用各变电站状态估计结果,对电网进行传统的电力系统拓扑分析和状态估计,传统状态估计目标函数如下: J(x)=[z-h(x)]<sup>T</sup>R<sup>-1</sup>]T[z-h(x)] 其中,状态量x=(v<sub>1</sub>,...v<sub>m</sub>,θ<sub>1</sub>,...,θ<sub>m</sub>)<sup>T</sup>表示电网各母线的电压幅值和相角,量测矢量<img file="FSA000002059244000413.GIF" wi="239" he="136" />包括变电站状态估计得到的量测量熟数据z<sub>se</sub>及没有进行状态估计的变电站的原始量测量z<sub>meas</sub>;(13)重复步骤(1)-(13),进行变电站-调度中心两级分布式电网的非线性状态估计。 
地址 100084 北京市海淀区清华园1号