电力系统网络拓扑结构识别

学院毕业设计（论文）题目：电力系统网络拓扑结构识别学生姓名：学号：学部（系）：机械与电气工程学部专业年级：电气工程及其自动化指导教师：职称或学位：教授目录摘要3ABSTRACT4一绪论61.1课题背景及意义61.2研究现状61.3本论文研究的主要工作7二电力系统网络拓扑结构72.1电网拓扑模型72.2拓扑模型的表达102.3广义乘法与广义加法112.4拓扑的传递性质12三矩阵方法在电力系统网络拓扑的应用143.1网络拓扑的基本概念14规定14定义15连通域的分离153.2电网元件的等值方法17厂站级两络拓扑17元件级网络拓扑173.3矩阵方法与传统方法的比较17四基于关联矩阵的网络拓扑结构识别方法研究184.1关联矩阵184.1.1算法18定义19算法基础194.2拓扑识别204.3主接线拓扑辨识原理224.4算法的简化与加速274.5流程图274.5.1算法流程图27节点编号的优化28消去中间节点和开关支路28算法的实现294.6分布式拓扑辨识法304.7举例和扩展30五全文总结32参考文献32致34摘要电力系统拓扑分析是电力能量流（生产、传输、使用）流动过程中，对用于转换、保护、控制这一过程的元件（在电力系统分析中认为阻抗近似为0的元件）状态的分析，目的是形成便于电网分析与计算的模型，它界于EMS底层和高层之间。就调度自动化而言，底层信息（如SCADA）是拓扑分析的基础，高层应用（如状态估计、安全调度等［1］）是拓扑分析的目的。可见，电力系统在实时运行中，这些元件的状态变化决定了运行方式的变化。如何依据厂站实时信息，快速、准确地跟踪这些变化，是实现电力系统调度自动化过程中基础而关键的工作［2］。拓扑分析在电力系统调度自动化中如此重要的地位，至少应该作到如下几点。（1） 拓扑分析的正确性：对任何情形下的运行方式，由元件状态的状况，针对各种电气接线关系，如单、双母线接线及旁路母线、3／2接线、角型接线等，均能进行正确的处理，当然这必须在实时信息可靠前提下才能实现。（2） 拓扑分析的直观性：大规模电力系统的拓扑结构是复杂的，由此拓扑分析本身就是对这一复杂网络的简化，因此其结果的直观性就很重要。如元件状态（运行、停运）标识，不同电压等级的区分等。（3） 拓扑分析的实时性：由拓扑分析的目的可知，拓扑分析必须是快速的，必须满足对实时决策与控制的要求。拓扑分析的通用性：运行方式变，电网结构就变，也即拓扑结构变，由此在拓扑数据的存储、模型表达等诸多方面都应该考虑其开放性、可扩展性及可维护性等。综上，电力系统网络拓扑分析的目的是明确的，同时也显现电力系统网络拓扑分析有一定的难度。关键词：电力系统；关联矩阵；拓扑分析；网络ABSTRACTPowersystemsassociatedtopologyistheelectricalenergy,transportstream(production,use)flow,forconversion,protectionandcontroltheprocessponent(thepowersystemsintheanalysisconsideredtheimpedanceapproximationtotheponents)stateanalysis,theaimistofacilitateanalysisandcalculationofthenetworkmodel,anditboundedontheems.Thebottomandtop.Intheautomation,informationastheunderlying(scada)topologyisthebasisofanalyzinghighapplications,suchasstatesecurity(estimate,thewaiting[1]topologyisthepurposeofanalysis.therefore,theelectricitysysteminreal-timeoperation,theseponentsofstatechangestotherunway.howthechangesonplantsstandreal-timeinformation,rapidandpreciseinfollowingthesechangesisthepowersystemsandautomatedschedulersthebasisofthework[2].Topologyanalysisintheelectricitysysteminthedeploymentofautomationissoimportantposition,atleastshouldbeasfollows.topology：toanyanalysisofthevalidityoftheoperation,theelementsofthestateofthesituationintheelectricalwiring,suchasasingle,doublebuseswiringandotherbuses3/2connection,thetypeofoperator,wecorrectlyhandled,ofcoursethismustbeinareal-timebasisofreliableinformationwillberealized.topologyanalysis,visualization：large-scalepowersystemsoftheassociatedtopologystructureareplexandtheassociatedtopologyanalysisistheplexnetwork,theimmediacyoftheresultisveryimportant.ifelementsstate(run,stopped)identity,thevoltagebetweendifferentlevels,etc.topologyanalysistimely：topologyanalysisbytheendoftheassociatedtopologyanalysismustbequick;wemustmeettherealdecisionsandcontrol.topologyanalysisuniversality：runway,andnetworkconstruction,alsotheassociatedtopologystructure,thedatastoredinthetopology,modelofexpressionandsoonshouldconsideritopen,scalabilityandcanmaintainsuch.Thepowersystem,networktopologythatthepurposeofanalysisisclear,andalsothepowersystemofnetworktopologyanalysisofdifficulty.Keywords：powersystem；incidencematrix；topologyidentification；network绪论1.1课所谓电力系统网络拓扑结构指的实际上是电力系统网络的各发电厂，变电所和开关站的布局，以及连接他它们的各级电压电力线路的连接方式。在电网发展初期，电网规模较小，电源布局对电网结构起重要作用。随着系统规模的不断扩，尤其是互联大电网的形成，电厂的作用相对弱化，于是电力系统网络主结构的规划设计变得尤为重要。电网互联，是各国电业工业发展的的客观规律，是世界各国电力发展的必然趋势。如此庞大的电网中，电网拓扑结构无疑直接决定着电力系统是否稳定，是否存在安全隐患，能否在意外发生的第一时间解决故障等等。从电网的发展中，为了谋求更多的经济效益和系统运行的稳定性，大电网取代了小电网。经济上大电网可以在最大的地理环境获得最好的能源利用,发挥大电网互联的错峰调峰、水火互济、跨流域补偿调节、互为备用和调节余缺等联网效益,实现网间功率交换,在更大围优化能源配置方式。同时,在安全上大电网承受扰动的能力比小电网显著加强,大电网因事故导致大停电的概率明显减小。在这种情况下，为了考量系统中设备的随机故障和负荷的不确定性，有了量化的电网风险评估。研究不同的电网系统拓扑结构，对于电网运行减少风险有着重要的意义[1。]研究现状随着电网状态估计技术的发展，电力系统拓扑结构分析方法得到了专家和学者的广泛重视，传统的电力系统拓扑分析方法一般将拓扑结构表述为链表关系，

用图论中的搜索技术，如深度优先搜索法和广度优先搜索法分析节点的连通性。这种方法一般需要建立反映拓扑结构的链表，通过处理链表实现拓扑分析。由于在电网的实际运行过程中，状态频繁发生变化的开关占少数，因此将追踪技术引入拓扑分析中，仅在开关状态发生改变时进行局部拓扑分析，可以减少拓扑分析的计算量。图论搜索虽易于理解，但较繁琐，不少学者在此基础上进行了更深入地研究和改进。在给出厂站、网络拓扑结构等概念后，独立进行厂站拓扑结构分析和网络拓扑结构分析，并引入稀疏、分块处理等技术进一步提高网络拓扑结构分析的效率。本论文研究的主要工作本篇论文的主要工作是：熟悉电力系统网络拓扑结构识别在电网风险评估中的运用。了解基于关联矩阵的电网拓扑辨识具体算法流程如何实现。使用关联矩阵的方法，写出算法流程图，用matlab编写相应的计算程序。电力系统网络拓扑结构电力系统网络拓扑结构2.1电网拓扑模型电力系统中的各类设备除输电线路外都集中于发电厂和变电站，厂站设备和各种输电线路的相互连接构成了电力系统网络拓扑。厂站包含的一次设备主要有：发电机、开关变压器、隔离开关、电抗器等。母线一般有单母线、单母线分段、双母线及双母线带旁母、倍半接线等形式；变压器根据其结构又可分成双绕组变压器和三绕组变压器。在各种接线形式下，断路器两边一般设置有隔离开关，断路器和隔离开关串用来连接母线、进出线路和变压器等。元件：电力系统一次设备集合中的一个元素，称为元件。元件按照其结构可以分为：单端点元件，只有一端和电网连接的设备，如发电机组、用电负荷、并联补偿器、调相机等；双端点元件，有两端和电网连接的设备，如断路器、隔离开关、输电线路、串联补偿器、双绕组变压器等；多端点元件，有多个端点和电网连接的设备，如三绕组变压器。在实际计算中，多端点元件可以根据端点的连接情况等值为多个双端点元件．按照元件的性质，元件又分为：无阻抗元件，一般将用于转换和控制电力系统运行方式的元件，称为(近似)无阻抗元件(如断路器、隔离开关等)；有阻抗元件，用于电能转换与传输的元件称为有阻抗元件(如输电线路、变压器等)。厂站：由若干元件连成的区域中，不包含任何输电线路元件的整体，若有输电线路仅含有输电线路元件的一个端点，这样的区域称为厂站。网络：由厂站拓扑分析后的逻辑节点和有阻抗元件构成的集合称为网络，网络是与厂站对应的。电气节点：元件之间的连接点称为电气节点，包含电气连接点和物理母线，所有设备通过电气节点连接在一起。逻辑节点：由无阻抗元件直接连接在一起形成的电气节点连通片称为逻辑节

点，逻辑节点都集中在厂站．系统节点：一个逻辑节点也称为系统节点，所有逻辑节点的总数是系统的最大节点号．子系统：由有阻抗元件连接在一起的系统节点的连通片称为子系统，子系统由网络拓扑分析确定。基于上述定义，本文的电网拓扑模型由厂站拓扑和网络拓扑两层构成，其中所有的电气节点、除输电线路外的所有元件都集中于厂站，厂站间的节点编号互不关联、彼此独立，厂站由输电线路连接构成电力系统的拓扑模型。(2)开关闭合图2-1网络模型(2)开关闭合图2-1网络模型1,2,3…电气节点号),(2)…阻抗元件号图2-1是基于本文定义的含有三个厂站的电网拓扑模型。2.2拓扑模型的表达电力系统主接线图可以由一个节点和边的集合G来描述:G=(V,E(s))2-1其中矿为节点集合，E为边集合，ss为边的开断状态（边赋权）。G在表示厂站拓扑时，V=｛电气节点｝,E=｛厂站双端元件｝；G在表示网络拓扑时，V=｛逻辑节点｝,E=｛有阻抗双端元件｝。也就是说，G由边、点及边赋V的三元集合构成，是一个无向图。在拓扑分析中，边赋Vs只有0和1两种赋值，边赋Vs=O表示断开，边赋Vs=1表示连通。对于—个节点数为n的拓扑结构用关联矩阵A作数学表达，矩阵的行，列号对应节点集合V,矩阵A的元素aj表示点i与点j间的连通性，对应边及其赋V集合E（s），具体表示为：a=ij1,i=j0,i主j,i电j(i=1,2, ,n;j=1,2,,n)2-2s,i主j,igj显然，关联矩阵是对称矩阵，其元素反映网络节点位置及节点间的直接连接关系。以图2-1所示的网络模型中的变电站2为例，该站有6个电气节点，4个无阻抗元件，—个有阻抗元件（标号为7），对有阻抗元件视为断开，形成初步关联矩阵为：sA=s0sA=s001ss00ss0ss010s010s0100s(s=0,1)2-3节点之间没有电气连接的取值为0,有电气连接的取值根据支路的状态s确定。变电站2中变压器支路在厂站拓扑中处理为断开，开关支路只有2-4支路断开，则当前状态下变电站的拓扑结构可以由下式关联矩阵表述。1011001100002-401A=1001000001

00001同样的方法可以确定其它变电站的关联矩阵及网络拓扑的矩阵。2.3广义乘法与广义加法关联矩阵的元素值(1或者0)表示节点间的连通关系，故对其值的运算属于布尔代数的运算．基于此本文利用广义乘和广义加两个运算规则，如下：x&x=min(x,x)nc12122-5x㊉x=max(x,x)1212其中，x,x表示关联矩阵元素。®和㊉这两种运算，满足以下运算规则：12加法的交换律和结合律x㊉x=x㊉x1221(x㊉x)㊉x=x㊉(x㊉x)123123乘法的结合律(x®x)®x=x®(x®x)123123乘法与加法间的分配律(x㊉x)®x=x®x㊉x®x1231323(4)加法有幺元0，乘法有幺元10㊉x=x㊉0=x1®x=x®1=x(5)加法的等幂律x®x=x若nxn矩阵A和B均为布尔代数矩阵，则其乘法和加法运算分别为TOC\o"1-5"\h\zAB二A®B£(a®b) 2-6jiikLj=1 」A㊉B=「a+b] 2-7jkjk其中，(j沁1,2,n;k=1,2.n;)2.4拓扑的传递性质...拓扑结构中点与点间的连通关系是可以通过边传递的，也就是说连通关系是具有传递性的。电力网络主接线关联矩阵法拓扑分析就是根据这种连通的传递性质确定拓扑结构中任意两点的连通牲。例如点j与点i相连，同时点i又与点k相连,则点j与点k也一定是相连的(无论点j与点k是否直接相连)，该牲质可用广义代数运算表示为：a=a㊉(a®a)=s㊉(1®1)=1;(s=0,1) 2-8jkjk jiik例如式(2-4)中，a=0说明图2-1变电站节点1和节点2无连接关系，但由于12a=s=1且a=s=1则节点1和节点2可以通过节点3间接相连。式(2-6)中a表13 32 jk示点j与点k的直接连接关系，即节点通过1阶支路的连通性；a'表示点j与点k通jk过2阶(及以下)支路的连通性，其中2阶支路的连通是通过节点的间接连通。这里的支路均为“简单路径”，即排除了有局部回路的支路。一个节点数为n的拓扑图中，最多通过n-1阶支路(简单路径)传递即可确定任意两点(点j与点k)之间的连通牲，只要拓扑结构中点j与点k存在可以连接的支路，经传递牲运算后必定满足a'=1。jk

对式(2-4)进行矩阵广义乘法和加法计算，得111101110011110A2=1011000001000对式(2-4)进行矩阵广义乘法和加法计算，得111101110011110A2=1011000001000010000112-9可以看出，矩阵点1-点2，点3-点4对应位置的元素由0变为1，说明变电站这两对点存在着连通的2阶路径，即支路1-3-2和支路4-1-3。继续计算矩阵广义乘法和加法计算，得1A3=11111001100110011001100001100112-10同样可发现，矩阵点2-点4对应位置的元素由0变为1，说明变电站存在这一对节点的3阶连通路径，即路径4-1-3-2。再进行计算，得1111111011110A4=111101111000001000010000112-11可发现A4对比A3不再有元素值发生变化，说明变电站2最高阶的路径是3阶路径。至此，变电站2任意两点之间的连通状态就确定了。矩阵方法在电力系统网络拓扑的应用网络拓扑的基本概念现代电力系统是由多个不同类型的元件所构成的结构复杂的网络系统，随着电网结构的不断发展，许多在电力生产实践中遇到的问题必须运用计算机进行分析和求解，做为分析电力系统潮流分布、稳定分析、人工智能专家系统等问题的根据将已知的电力系统网络变成计算机可以识别的网络结构，就必须运用网络拓扑技术加以解决。根据电力发展，建立怎样的网络拓扑结构，对运算结果的精确性及网络拓扑结构的可扩展性有很大影响．因此选择什么样的系统模型建立网络拓扑结构，是进行电网分析计算、解决复杂问题的关键。传统的网络拓扑技术是采用链路连接实现的方法．在计算机中实现电力系统网络拓扑结构，这种方法以已知节点出发，按照各节点问是否存在相关参数来建立网络拓扑结构．逐次形成整个电网的网络拓扑结构。该方法在电网结构发生一定的变化时，扩展性显得不够灵活，需要修正计算机程序中的许多容，运行维护较复杂。下面介绍采用设立网络矩阵的方法．实现电网结构的网络拓扑模型。规定设网络节点数为N,节点依次编号，起始节点为1,终止节点为N；节点i和节点j之间有直接连线时•元素A等于1,元素A等于1；p ji矩阵元素A等于0;a节点i和节点j之间没有直接连线时，元素A等于0,元素A等于0；q ij根据以上3条规定构成矩阵A,矩阵A有如下特点：是N阶方阵；是主对角线元素等于0的对角矩阵；是元素大部分为0的稀疏矩阵。3.1.2定义节点i与节点j连通：节点i与节点j至少有一条直接连线；连通域：某一节点集合中，从任意一个节点出发，每次经过一条连线到达另一节点，总可以找到该集合中所有的节点；非连通域：某一节点集合S和另一节点集合P,从集合P中的任意节点出发，无论经过多少条连线，都找不到集合S中的节点。连通域的分离N阶方阵A所代表的N个节点•可能构成一个连通域•也可能构成M个连通域(M小于等N)。把在同一个连通域中的所有节点都找出来是问题的关键，举例说明分离连通域的方法。图3-1网络节点连接情况例：某一网络有5个节点(a，b，c，d.e),它们的连接情况如图3-1。从网络图中可知节点(a,b,d)构成一个连通域，节点(c,e)构成一个连通域。下面利用状态矩阵A得到这个结论。根据规定，这5个节点构成的网络形成的矩阵TOC\o"1-5"\h\z'0 1 0 0 0'10 0 10A是0 0 0 0 10 10 0 0、0 0 1 0 0丿先从矩阵的第1行开始查找，第2列有1,记下列号，该行投有I了，转入第2

步，把刚才记下的列号当作行号，即查找第2行，第I列、第4列有I,记下列号，因为第I行已查过，故从第4行查找，第2列有1.因为第2行已查过，不能转到其他行，停止查找，得到区域（1，2，4）；同样从第3行开始查找，重复以上过程得到区域（3，5）。查找过程用计算机程序流程图表示如图3-2。图3-2程序流程图流程图中X（N,P）中存放得到的各区域的节点号，N是区域的个数，P是每x1x1个区域中节点个数。此程序已用VB语言在AcceSs数据库下调试通过。3.2电网元件的等值方法3.2.1厂站级两络拓扑把一个厂站（包括不同电压等级的母线、变压器、发电机等）看成整体当作一个节点，厂站i与厂站j之间有一条联络线且运彳丁，则A=A=1，有两条运行的联络ijji线，则A=A=2。若停一条线路，相应的矩阵元素减1,利用上面原理很容易判定ijji某操作是否将造成电网解列。3.2.2元件级网络拓扑厂站外部同上。厂站部元件的等值方法如下：每1条母线是1个节点；母联开关看作联系母线节点的联线；两卷变压器看作联系不同电压等级母线节点的联线；三卷变压器的等值，需要增加一个虚拟节点，与该变压器相连不同电压等级的母线节点与该虚拟节点之间增加一条联线。整个电网用上述方法构成矩阵A，操作电网中任意元件将改变矩阵A中元素的值，利用上面原理很容易判定某操作是否将造成电网解列。矩阵方法与传统方法的比较矩阵方法与传统的网络拓扑方法比较，当矩阵阶数较大时，直接使用存在参数较多的缺点。因此，首先应进行一定的结构处理，即可采用此方法。

利用矩阵结构的网络拓扑技术，可以实现各种复杂结构的电力网在电子计算机中形成网络拓扑结构。当电网扩建新厂、站时，根据增加节点在电网中的位置，建立相应的结构数据，对于计算机中运行的原程序并不需要修改传统方法则在这方面存在一定的局限性。因此，矩阵方法具有很强的可扩展性。同时，建立网络拓扑矩阵与电网接线结构相互对应，具有运用灵活和修改方便的优点。四基于关联矩阵的网络拓扑结构识别方法研究4.1关联矩阵4.1.1算法网络拓扑的基本形态大致可以分为链状结构和树状结构，如图4-1所示，准确地讲，应该是通过二者的组合和演化而来．将多个链状进行组装就是树状，将树状拆分就变成多个链状。(a)链状结构 (b)树状结构图4-1网络拓扑基本结构在实际的网络辨识过程中，总是选用其中某个节点作为起始点，相应地称之为根，而那些末端节点则称之为树叶或输入节点，其他的就构成树枝或树干，这全是为了方便描述而形象化的一种定义。定义设网络包含的节点集合为N=(n,n,n,,n)，其中m为有限值。选定其节123m点n作为起始点(根)，当n对节点n进行查询时，n的应答信息途经节点集合Lsskk•••中的节点，或者说n的应答信息经由节点集合L后抵达n,则定义L中的节点与nksk相关联•若L中有P个节点，则p<m,并且nek称之为自相关，显然n与所有Nks的节点相关联，而与末端节点只有自相关。定义一个mxm阶的矩阵A,其元素a只有“0”或“1”这2种赋值■当节点niji与节点n相关联时，a=1,否则a=0,这个矩阵A就称之为N的关联矩阵。jijij根据定义,网络关联矩阵有以下特殊的结构和特性：链状网络的关联矩阵就是一个上三角矩阵；树状(分枝状)网络的关联矩阵同样具有上三角矩阵的特征,但只有第一行和对角元素为“1”,其余元素均为“0”；起始点n不同，网络N的关联矩阵A就不同。s算法基础如前所述,网络拓扑辨识过程就是找出网络节点及其与位置的相互关系．网络拓扑辨识算法就是,找出一种能够提取出网络节点及其位置关系信息的数学模型和方法。能够大致反映网络拓扑形状或形态的方法称之为预估算法或评估算法,这样的算法有很多,而能够准确地反映网络拓扑及其结构的算法称之为辨识算法。设网络N有m个节点，其关联矩阵A是一个mxm阶的矩阵，为了方便提取网络节点及其位置信息，可定义一个m维全“I”列向量s,并做一个简单的运算As（若网络N为图4-1（a）所示的链状网络），则TOC\o"1-5"\h\zAs=[m,m-1,2,1]t （1）若网络N为图1b所示树状（分枝状）网络，则•••As=[m,1,1,1]T （2）将图4-1（b）所示网络N上下倒置，变成倒树状的网络，则•••As=[2,2,2,1]T （3）很明显，这3个量不仅数值意义明确，而且完全能够反映出各自网络的特点•••和结构．可见，根据关联矩阵的特点，三角状矩阵可直接反映网络节点及其位置的相互关系．也可以说，只要找出了一个网络的关联矩阵，就已经得出了网络的拓扑结构图。另外，参照其他的拓扑预估方法，如常用的分级和多级结构划分方法，式（1）和式（2）仍能明确地划分出等级，而且使得2种不同的结构用一种数学形式表达出来，更进一步得出的拓扑不再是一簇拓扑树形成的“森林”，而是惟一准确的拓扑树．4.2拓扑识别实际的网络不可能是标准形状的，显然关联矩阵也有所不同，不可能是标准的三角状矩阵。由于起始点的选择有所不同，得到的关联矩阵也有所不同，但是再复杂的拓扑也是由图4-1所示的标准拓扑元素组合而成的，其关联矩阵也能反映其拓扑关系。首先应当解决的问题是，如何规关联矩阵，使之成为三角状矩阵。之所以关联矩阵不再是三角状矩阵，并不是因为关联矩阵本身不再具备三角状特征，而是因为每个节点的编排顺序不可能是按图1所示的网络节点顺序编排，这正说明节点顺序需要重新编排，或者说理顺各节点的顺序，关联矩阵自然就变成了三角状矩阵。假设一个树状网络N的关联矩阵为非三角状的mxm阶矩阵A,定义一个m维全“1”列向量s,并做运算TOC\o"1-5"\h\zAs=[r,r, ,r]t (4)12m其中r,r,,r，不再满足r>r> >r的关系或规律，但按由大至小的顺1 2 m 1 2 m•••序重新排列r,r,,r，后得到r',r',,r'，，使得满足r'>r'> >r'，按1 2 m 1 2 m 1 2 m••••••照同样的顺序对A进行变换，得到关联矩阵A，则•••••• •••A's=[r',r', ,r']T (5)1 2 m此时的A'所反映的拓扑关系丝毫未变，但具有三角状特性并不一定具有标•••准形状•所以，有必要设法找出标准形状的关联矩阵。如前所述,不论是标准的链状结构还是树状结构,每一个末端节点在关联矩阵中除了自相关元素为“1”之外，该行的其余元素均为“0”。另外，据前所述不论什么形状的网络拓扑，均可以拆分成链状和树状。所以，原则上讲只要将关联矩阵A'中的标准矩阵至上而下地提取出来，就足以完全得出各级的拓扑结构，最终获得网络N的拓扑结构。此外，已拆分过的网络拓扑可以进一步拆分、组合成链状结构的拓扑，特别是可以拆分成q(假设有q个末端节点)条以起始节点为头、末端节点为尾的链状结构。换言之，只要在关联矩阵A'中找出这q个标准三角矩阵，就找出了相应的链状结构，也相当于得出了网络N的结构。很显然，只要将式(4)向量中数值为“1”的元素所对应的A'的列向量提取出来，则这个列向量中非“0”元素对应的节点就是该链状结构所包含的节点，这些节点的排列顺序就是它们在链状结构中的顺序。如果需要，可以从A'提取出相应的标准三角矩阵，如此重复q次，就可以得出网络N的拓扑结构。4.3主接线拓扑辨识原理根据网络拓扑理论,对于一个任意的拓扑网络,可以用节点2支路关联矩阵来描述其拓扑结构,而对于一个电力系统主接线图,可以把它抽象成为一个拓扑图来描述。如对图4-1(a)所示的主接线,把主接线的节点作为拓扑图的节点,把开关元件作为拓扑图的支路(当开关闭合时该支路连通,开关断开时该支路断开),对于单个变电站,把母线的每一进出线连接点也作为节点,节点-支路关联矩阵A[a]的每个元素表示意义如下:ija表示节点i与支路j的关联值(联通性),当节点i与支路j相联时a=1,ijij否则a=0。显然，当支路j的开关断开时，支路j与原来与其相联的2节点不再相ij联,因而与这2节点的关联值应为0。图4-1(a)是当所有开关都闭合时的状态,此时节点2支路关联矩阵(称为原始节点2支路关联矩阵)为11100000000000111001000001001000A=00010010000010010000010010000010001000001如果断开其中的某些开关,如图4-1(b)所示,此时对应的节点2支路关联矩阵A(称之为当前节点2支路关联矩阵或简称为节点2支路关联矩阵)为001000000000000100000100000010000000000101000000000100000001定义开关状态矢量S二［s］,与开关j的状态相对应。当开关闭合时，S=1,jjj开关断开时，s=0。这样，S=［0,0,1,1,1,0,1,0,1］。A就是An的每一行与S的j0各个对应元素进行“与”运算后得到的。用同样的方法,也可以定义支路2节点关联矩阵B二［b］，显然，A=BT。ij①②⑥⑦⑧①②⑥⑦⑧a）③④⑤①②⑥⑦⑧（b）a）③④⑤①②⑥⑦⑧（b）□表示开关合位，■表示开关开位。3图4-23开关接线型式2对电力网络主接线的辨识,实际就是要找出网络中节点间的连通关系。显然,这种连通关系是可以传递的,即①如果节点与j支路相联，同时，支路j又与节点k相联，则节点i与节点k是相联的；②如果节点i与节点k相联，而节点k又与节点I相联，则节点与节点也是相联的。定义节点2节点连通矩阵（它与网络图论中的节点2节点关联矩阵有区别）C=[C]。当节点i与节点j连通时，c=1,不连通时c=0,显然C是对称矩阵。ijijij以上连通关系的传递性质可以表示如下:若a=1,b=1,则c=aQb=1；ijjkikijjk若c=1,c=1,则c=cQc=1。ikklilikkl对于具有m个节点，n条支路的网络，定义以下矩阵乘法运算：C=A•B（1）其中c=n（ab） （2）ijikkjQ表示“与”运算；U表示“或”运算。那么当A为节点-支路关联矩阵，B为支路-节点关联矩阵时。c将表示节点iji与节点j通过任一支路的关联情况。当然，只要节点i与节点j有一条支路相联,则c=1。C=[c]表示了节点与节点之间的连通性，称为节点-节点连通矩阵。这ijij时的节点-节点连通矩阵仅仅表示了节点之间的直接连通性质,把它称为1级节点-节点连通矩阵，并记为C⑴。由于连通性的传递性质，可以通过用1级节点-节点连通矩阵C⑴进行上面定义的矩阵乘法运算，得到2级节点-节点连通矩阵C（2）=C（i）-C（i）（3）2级节点-点连通矩阵02）在1级节点-节点连通矩阵C⑴的基础上。运用连通关系的传递性，把节点之间的部分间接连通关系也表示出来了。用C（2）再自乘得

到C(3),,直到C(n)=C(n-1)-C(n-1)。这时,通过传递,所有连通的节点之间的关联值都是1了,而不连通的节点间关联值都为0。如对图4-2(b)C(i)=A-B=A-C(i)=A-B=A-At=1000100001000101001001000001001010001000011001000001001001000001C(2)=C(i)-C(i)=C(2)=C(i)-C(i)=1000100001100101011001000001001010001000011001010001001000000101C⑶=C(2)-C(2)=C⑶=C(2)-C(2)=1000100001100101011001000001001010001000011001010001001000000101再重复进行以上运算，将发现C不再变化。这时已把网络中的所有连通关系都表示出来了。由C可见，节点连接成以下几组:节点1、5为一组；节点2、3、6、8为一组；节点4、7为一组。4.4算法的简化与加速利用B二At化简在式中(1)由于B二At，所以式(2)可以转化为nnc= (ab)= (aa) (4)ijikkjikjk日n k=1n即C>4是A的第i行与第j行按位求“与”后的指，如果其中有一位为1,则c=1;ijij利用C的对称性C(2)=C⑴•C⑴等运算过程，实际上就是节点间联通关系的传递过程。犹豫C的对称性，c=m(cc)=m(cc) (5)ij ikkj ikjk即c®由C的第i行与第j行按位求“与”后再“或”如果其中有一位为1,则c=ijij1；同时，由于在节点连通关系的传递过程中，如果节点i与节点j在传递前是连通的，则在传递运算后仍然是连通的。因此，如果c在传递运算前是1，则c无ijij需再进行传递运算。为了充分利用前一步连通性传递结果，更进一步加快传递速度，在对C的自乘过程中，可以直接对c中为0的元素进行运算，并用运算结果直接修改该元素。4.5流程4.5流程4.5.1算法流程实际上,在求连通关系的过程中,只要进行以下过程即可c二 (a a),i二1ij ikjkkmic (c c),i=1ij ik jk图4-3算法框图本算法的框图如图4-3所示。4.5.2节点编号的优化为了加快连通关系的传递速度,应把连接开关较多的节点编号在前,如母线节点。因为,这种传递运算实际是从低节点号开始的,编号较小的节点参加传递运算的次数较多,所以,把连接支路较多的节点编号在前有利于加快传递速度。4.5.3消去中间节点和开关支路对网络进行拓扑辨识中,不需要“过多关心”的中间节点和开关支路可以消去,从而减少节点和支路数目,如图4-4所示。图4-2的网络实际已经把隔离开关支路和中间节点消去了。K2③K1④K3②K2③K1④K3②■0 \ 7—©_V-0①KeqK=K-K-Keq 1 2 32图4-4消去中间节点和支路K=K-K1212K =K・K1313算法的实现用c语言实现用软件编程实现以上算法时，虽然矩阵A和C的存储量分别是mXn和mxm,但用c语言编程时可以充分利用C语言的位运算功能，把一个字与16位的位结构作为共用体。这16位与16条支路(或节点)对应，一个mxn矩阵可以用mx(n/16+1)个字表示。另外，矩阵的行间“与”运算用字作按位“与”运算，这样运算速度可以大大加快，使用32位或64位计算机效果更好。用汇编语言实现由于本算法很简单，用汇编语言编程实现也很方便，笔者已用MCS28051单片机开发了能辨识64节点x64支路的专用辨识模块［14用硬件实现本算法的最大特点是能用简单的ASIC实现。4.6分布式拓扑辨识法电力网络是由许多发电厂和变电站组成的一个多电压等级网络，在每一个发电厂和变电站安装一辨识模块，再在主站安装系统辨识模块，可以实现分层分布式辨识。4.7举例和扩展如前所述，采用关联矩阵算法进行网络拓扑辨识，可以惟一地计算出（辨识）网络拓扑结构，而且可以分级、分区域进行，就是说对一个大型的网络辨识，可以将其分为很多小的局部或区域网络来处理。不仅如此，在很多现有的网络拓扑辨识（评估）算法中规定网络拓扑不能构成环状（并联支路）［15，16］，但在现实网络，尤其是在大型网络中，环状结构是不可避免的。下面以一个包含环状结构的网络拓扑辨识为例，来说明本文方法的优越性和有效性。图4-5含环状网络拓扑结构图4-5所示的网络关联矩阵为(1111111\1011011001011A=000111000011000001020304⑹050610000111丿7其中，矩阵右侧的自然数表示节点排列顺序号。As二[7,4,3,3,2,1,3]T经标准变形后，则式(6)转换为(1111111\1011001001001A=000101000011000001121314(7)1715(0000001丿6可见其还是一个三角矩阵，但却是中空的。中空部分正是反映了3个支路构成的环路部分，也正好是树状结构和倒树状结构的组合。按前述的方法，节点6对应的AS的元素为“1”仍然可以认为是一条链状网络，所不同的是该链路中出现了并联支路。将节点5、节点6构成的短链与主链断开，就又可以用标准的处理模式来处理了。至此说明，本文方法不仅有效，而且适用于环路或并联支路的网络拓扑辨识，运算量也仅仅与节点数772成正比。算法流程如图4-6所示。图4-6算法流程五全文总结基于关联矩阵的电网拓扑辨识算法有以下特点算法简单、速度快,可以用多种方法实现。在用硬件实现时,很少占用CPU的时间。可以很方便地实现分层分布式的拓扑辨识,特别适用于电力系统的实时辨识。网络关联矩阵拓扑辨识方法准确、快速、运算量小。通过适当的分级和分解可以将大型的网络辨识问题分解为局部的小型网络拓扑辨识问题，也可以将整个网络的拓扑辨识问题，根据需要化解为部分区域的网络辨识，以适应于多种网络的拓扑辨识。参考文献[1]KambaleP，MackauerJ．J．TheDispatcherTrainingSimulatorfor［2］ 于尔铿，广一周京阳，等•能量管理系统(EMS).：科学，1998［3］ 印永华，郭强，申洪，华东．大电网发展的经验和借鉴中国电力科学研究院［4］ 周家启,霞。电力系统风险评估方法和应用