发明名称 一种基于同步相角测量装置的电力系统多区域分布式状态估计方法
摘要 本发明公开了一种基于同步相角测量装置的电力系统多区域分布式状态估计方法,利用来自SCADA系统和PMU的量测数据实时快速的估计出电网各节点的电压幅值和相角值并将各个区域的状态信息进行实时交互。本发明采用基于地域特点的互不重叠解耦策略,将一个大电力系统解耦为几个子系统,并通过一个协调系统来整合各个子系统的状态;多区域分布式状态估计方法的采用,将大系统状态估计变为一系列局部小区域状态估计,大大提高了计算速度。为系统提供实时准确的电压、相角等量测信息,保证系统获得更高的量测冗余度进而提高状态估计的精度。本发明的区域交互管理更利于新能源安全有效的接入,以满足未来智能电网发展要求。
申请公布号 CN103248043B 申请公布日期 2015.01.28
申请号 CN201310157443.7 申请日期 2013.04.28
申请人 西南交通大学 发明人 张葛祥;赵俊博
分类号 H02J3/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 成都信博专利代理有限责任公司 51200 代理人 张澎
主权项 一种基于同步相角测量装置的电力系统多区域分布式状态估计方法,利用来自SCADA系统和PMU(Phasor Measurement Unit,PMU)的量测数据实时快速的估计出电网各节点的电压幅值和相角值并将各个区域的状态信息进行实时交互,包含如下的处理步骤:步骤1:基于地域特点的互不重叠多区域电力系统解耦在一个含有n条接线的电力系统中,按照地区分为r个不重叠的区域S<sub>i</sub>,每个S<sub>i</sub>有n<sub>i</sub>条接线,并且和连接线相连;解耦之后,每一个区域都由自己的本地控制中心控制,该本地控制中心负责估计自己的状态,并通过通信线路连接到一个协调的控制中心;步骤2:读取各个区域的网络数据,包括网络结构和线路阻抗,并形成节点导纳矩阵和支路‑节点关联矩阵;步骤3:系统量测和PMU配置;系统的量测包括电压幅值量测、有功和无功潮流量测、有功和无功注入潮流量测,量测方程如下:●电压幅值:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>V</mi><mi>l</mi><mi>mea</mi></msubsup><mo>=</mo><msub><mi>V</mi><mi>l</mi></msub><mo>+</mo><msub><mi>e</mi><msub><mi>V</mi><mi>l</mi></msub></msub></mrow>]]></math><img file="FDA0000601362170000011.GIF" wi="250" he="82" /></maths>●有功和无功潮流:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>lm</mi><mi>mea</mi></msubsup><mo>=</mo><msub><mi>P</mi><mi>lm</mi></msub><mo>+</mo><msub><mi>e</mi><msub><mi>P</mi><mi>lm</mi></msub></msub></mrow>]]></math><img file="FDA0000601362170000012.GIF" wi="275" he="76" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>Q</mi><mi>lm</mi><mi>mea</mi></msubsup><mo>=</mo><msub><mi>Q</mi><mi>lm</mi></msub><mo>+</mo><msub><mi>e</mi><msub><mi>Q</mi><mi>lm</mi></msub></msub></mrow>]]></math><img file="FDA0000601362170000013.GIF" wi="297" he="81" /></maths>●有功和无功注入潮流:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>l</mi><mi>mea</mi></msubsup><mo>=</mo><mrow><mo>(</mo><msup><msub><mi>V</mi><mi>l</mi></msub><mn>2</mn></msup><msub><mi>g</mi><mi>l</mi></msub><mo>+</mo><munder><mi>&Sigma;</mi><mrow><mi>m</mi><mo>&Element;</mo><mi>a</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><msub><mi>P</mi><mi>lm</mi></msub><mo>)</mo></mrow><mo>+</mo><munder><mi>&Sigma;</mi><mrow><mi>m</mi><mo>&Element;</mo><mi>b</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><msub><mi>P</mi><mi>lm</mi></msub><mo>+</mo><msub><mi>e</mi><msub><mi>P</mi><mi>l</mi></msub></msub></mrow>]]></math><img file="FDA0000601362170000014.GIF" wi="724" he="170" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>Q</mi><mi>l</mi><mi>mea</mi></msubsup><mo>=</mo><mrow><mo>(</mo><msup><msub><mrow><mo>-</mo><mi>V</mi></mrow><mi>l</mi></msub><mn>2</mn></msup><msub><mi>b</mi><mi>l</mi></msub><mo>+</mo><munder><mi>&Sigma;</mi><mrow><mi>m</mi><mo>&Element;</mo><mi>a</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><msub><mi>Q</mi><mi>lm</mi></msub><mo>)</mo></mrow><mo>+</mo><munder><mi>&Sigma;</mi><mrow><mi>m</mi><mo>&Element;</mo><mi>b</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow></munder><msub><mi>Q</mi><mi>lm</mi></msub><mo>+</mo><msub><mi>e</mi><msub><mi>Q</mi><mi>l</mi></msub></msub></mrow>]]></math><img file="FDA0000601362170000015.GIF" wi="757" he="170" /></maths>其中,步骤1中r个不重叠的区域S<sub>i</sub>内的注入有功和无功潮流量测函数方程分别为:P<sub>lm</sub>=V<sub>l</sub><sup>2</sup>(g<sub>lm</sub>+g<sub>slm</sub>)‑V<sub>l</sub>V<sub>m</sub>(g<sub>lm</sub>cos(θ<sub>l</sub>‑θ<sub>m</sub>)+b<sub>lm</sub>sin(θ<sub>l</sub>‑θ<sub>m</sub>))Q<sub>lm</sub>=‑V<sub>l</sub><sup>2</sup>(b<sub>lm</sub>+b<sub>slm</sub>)‑V<sub>l</sub>V<sub>m</sub>(g<sub>lm</sub>sin(θ<sub>l</sub>‑θ<sub>m</sub>)‑b<sub>lm</sub>cos(θ<sub>l</sub>‑θ<sub>m</sub>));l是属于区域S<sub>i</sub>的一条接线;<img file="FDA0000601362170000021.GIF" wi="590" he="71" />和<img file="FDA0000601362170000022.GIF" wi="96" he="70" />分别是接线l和m处电压幅值、有功和无功潮流、有功和无功注入潮流的量测值;V<sub>l</sub>,V<sub>m</sub>分别是接线l和m的电压幅值;θ<sub>l</sub>和θ<sub>m</sub>分别是接线l和m的相角值;<img file="FDA0000601362170000023.GIF" wi="367" he="70" />和<img file="FDA0000601362170000024.GIF" wi="58" he="69" />分别是电压幅值、有功和无功潮流、有功和无功注入潮流的量测误差;g<sub>lm</sub>+jb<sub>lm</sub>是支路l‑m串联导纳值;其中g<sub>lm</sub>为支路l‑m串联电导值,b<sub>lm</sub>为支路l‑m串联电纳值;g<sub>slm</sub>+jb<sub>slm</sub>是支路l‑m并联导纳值;其中g<sub>slm</sub>为支路l‑m并联电导值,b<sub>slm</sub>为支路l‑m并联电纳值;g<sub>l</sub>和b<sub>l</sub>分别是连接到线l的并联电导和电纳值;a(l)为连到l上并且属于区域S<sub>i</sub>的所有线的集合;b(l)为连到l上并且属于不同于区域S<sub>i</sub>的区域S<sub>j</sub>的所有线的集合,i≠j;●PMU配置原则为了保证整个系统的可观测性,本发明中对每一个子系统都至少配置一个PMU并且PMU都配置在子系统与子系统之间的连接线上;步骤4:分区域状态估计对各子区域分别采用最小二乘法进行状态估计,分区域状态估计模型如下:z<sub>i</sub>=h<sub>i</sub>(x<sub>i</sub>)+e<sub>i</sub>,i=1,2,...rz<sub>c</sub>=h<sub>c</sub>(x)+e<sub>c</sub>其中,z<sub>i</sub>是区域S<sub>i</sub>的m<sub>i</sub>×1维内部节点或本地量测向量;z<sub>c</sub>是m<sub>c</sub>×1维边界量测向量;<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>x</mi><mi>i</mi></msub><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>&theta;</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000601362170000025.GIF" wi="176" he="157" /></maths>是区域S<sub>i</sub>中2n<sub>i</sub>×1维本地状态向量,包含n<sub>i</sub>个电压值和n<sub>i</sub>个相角值;<img file="FDA0000601362170000026.GIF" wi="348" he="101" />是全系统的状态向量;h<sub>i</sub>(x<sub>i</sub>),h<sub>c</sub>(x)为基于基尔霍夫电压或电流定律的非线性向量量测函数;e<sub>i</sub>,e<sub>c</sub>为高斯随机误差向量;步骤5:PMU协调采用PMU的实时准确量测来对系统的各子区域状态信息进行线性协调,状态估计在协调水平变为线性,并且线性状态估计模型为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>z</mi><mi>p</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>z</mi><msub><mi>V</mi><mi>p</mi></msub></msub></mtd></mtr><mtr><mtd><msub><mi>z</mi><msub><mi>I</mi><mi>p</mi></msub></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>I</mi></mtd></mtr><mtr><mtd><mi>Y</mi></mtd></mtr></mtable></mfenced><mo>&CenterDot;</mo><mi>x</mi></mrow>]]></math><img file="FDA0000601362170000027.GIF" wi="375" he="167" /></maths>其中,<img file="FDA0000601362170000031.GIF" wi="63" he="62" />和<img file="FDA0000601362170000032.GIF" wi="54" he="62" />分别为PMU的电压和电流量测向量;I表示关于电压向量的量测单位矩阵;Y表示关于电流量测向量的线性矩阵,由节点和支路导纳组成;步骤6:分布式不良数据处理采用协调系统不良数据处理和子区域不良数据处理两个阶段处理,即:1.计算协调系统自身的稀疏逆矩阵,完成协调系统不良数据的处理并将结果传输到子区域的状态估计器中;2.在接收到协调系统传过来的协调系统自身稀疏逆矩阵后,接着进行子区域状态估计并计算各个子区域的稀疏逆矩阵,进而完成子区域不良数据处理;步骤7:若所有的状态估计都收敛,不良数据处理结束,则输出各个节点的电压幅值和相角值并且输出不良数据的辨识结果,否则转步骤4。
地址 610031 四川省成都市二环路北一段111号西南交通大学科技处