1、书书书第 卷第期复 杂 系 统 与 复 杂 性 科 学 年 月 文章编号：（）；：基于复杂网络理论的电力网络抗毁性分析郭明健，高岩（上海理工大学系统科学系，上海 ）摘要：抗毁性分析是电网安全研究的核心内容之一，因传统分析方法无法有效分析故障产生的过程，对研究抗毁性存在局限性，基于复杂网络理论研究电力网络的抗毁性并以上海市崇明区为例实证分析。对电力网络随机攻击和蓄意攻击，得出攻击后网络效率和最大连通子图数量变化，并提出网络效率变化率作为评估抗毁性的参数。根据仿真结果，提出一种基于实时接近中心性优先攻击策略的分段式保护方案，提高电力网络的抗毁性，以保证电网安全。关键词：电力系统；抗毁性；复杂网络；
2、网络效率；接近中心性中图分类号：文献标识码：收稿日期：；修回日期：基金项目：国家自然科学基金（）第一作者：郭明健（），男，江苏镇江人，硕士，主要研究方向为电力系统，复杂网络。通信作者：高岩（），男，黑龙江五常人，博士，教授，主要研究方向为电力系统需求侧管理，系统分析与优化等。，（，）：，：；引言电能是现代社会中最重要、最方便的能源。发电厂将其它形式的能量转换成电能，电能经过变压器和不同电压等级的输电线路输送并分配给用户，再通过各种用电设备转换成用户需要的能量。生产、输送、分配和消费电能的各种电气设备连接在一起而组成的整体称为电力系统。现代电力系统是一个复杂系统，随着其规模的不断扩大，维持电网稳
3、定难度不断提升。近些年，国内外发生了多起大型停电事故，例如 年美国纽约城区因复 杂 系 统 与 复 杂 性 科 学 年 月某变电站老旧变压器起火，导致整个城区大型停电事故发生；年中国台湾因工作人员疏忽，致使供气中断导致电厂跳机，进而发生全台大断电事故。这些停电事故会给社会经济和人民生活带来长期的不利影响。近年来，人们对于智能电网的需求侧响应及系统能源管理的研究已经逐步完善，然而对于其安全性及稳定性的研究还不够成熟。因此电力系统的抗毁性分析不仅具有理论研究意义，更具有保证电网安全稳定的现实价值。网络抗毁性是指遭受随机性或确定性攻击时，网络保持其连通性的能力。网络抗毁性可分为静态抗毁性和动态抗毁性
4、，其区别在于故障节点或边是否会导致其他节点或边的故障。抗毁性分析是电力系统中识别易受攻击的脆弱组件的关键问题。由于电力系统的复杂性，基于电气介数的度量方法难以系统分析故障的传播过程，复杂网络理论利用网络拓扑结构及属性可有效分析电力系统抗毁性。复杂网络理论源于图论的研究，随着其理论发展，现已广泛应用于社交网络、生物科学、电力系统 等领域。对于电力系统的抗毁性分析，构建拓扑脆弱性指数是一种流行的评估网络脆弱性的方法，文献 基于拓扑度量来分析电力系统网络，例如平均最短路径长度，出入度等参数。然而这类方法不能揭示失效节点的传播机制，文献 基于故障链的级联故障图提出一种评估输电网络抗毁性的方法；等 基于
5、凸优化及复杂网络提出一种模块化和高度可扩展的算法以便进行安全性评估；文献 基于介数中心性采用蓄意攻击得到边的抗毁性参数识别关键设施。本文首先基于复杂网络理论，将变电站、发电厂等站点看作节点，将传输电缆看作边构建电力网络拓扑图；接着模拟针对该网络的随机攻击和蓄意攻击，并提出网络效率变化率作为抗毁性分析的评估参数；最后以上海市崇明区电力网络做实例分析，进而针对性地提出解决方案，以确保电网安全稳定运行。模型构建相较于其他工业系统，电力系统最大的特点是不易存储，即其生产、输送和消费在同一时刻实现。电力系统的基本要求是保证安全优质地向用户供电，但风电、新能源汽车、光伏、特高压等新能源也将给新型电力系统稳
6、定运行带来巨大挑战。新型电力系统结构示意图如图所示。图新型电力系统的结构示意图 新型电力系统中稳定运行受到光伏或风电发电量的影响，光电风电功率不确定性高、抗扰性弱，易造成电力系统常规机组切除，极大约束了电力系统对新能源的消纳水平。其主体结构包含发电厂、变电站、输配电线路。为了便于研究，在构建电力网络模型时进行如下简化：）不考虑功率、频率、电压、电流、磁链、电动势等动态参数的影响；）不考虑输入输出设备差异，将任一发电厂、变电站视为相同类型的节点；）不考虑输电配电线路的差异，将节点间的连接视为相同的边。基于前文将电力系统看作一个网络，将不同类型的变电站、发电厂等站点看作同类节点，将输配电线路看作边
7、构建电力拓扑网络。电力网络拓扑图表示为（，），其中集合包含节点，即变电站、发电厂等站点；表示网络的边集合，代表节点间的传输电缆。此外，其邻接矩阵为 。其中 时表示节点与节点相连接，时两节点不相连，该模型的加权邻接矩阵 且定义为 节点与节点相连其他（）其中，为节点与节点边上的权值，为了简化网络，均以带入计算。参数选取和评估指标选取参数选取研究网络抗毁性时，首先要确定枢纽节点，其中介数中心性、接近中心性反映节点重要程度的指标。介数是使用最广泛的中心性指标，当网络中移除具有较大介数中心性的点或边时，会使得网络处于高中断的风险中。节点的接近中心性反映某一节点与其他节点之间的接近程度，用一个节点到所有其
8、他节点的最短距离的平均值的倒数来表示接近中心性，该节点距离其他节点越近，它的接近性中心性越大，它也可以有效反映基于时间序列第 卷第期郭明健，等：基于复杂网络理论的电力网络抗毁性分析产生的递归网络之间的紧密程度。本文将介数中心性和接近中心性作为判断节点重要程度的参数。节点介数是指一个网络里通过节点的最短路径条数，节点的介数中心性可表示为（）（）其中，为连接和的最短路径的数量，为其中经过节点的最短路径的数量。节点的接近中心性可表示为（）（）其中，为节点到其余各点的平均距离，平均距离的倒数就是接近中心度。评估指标选取评估网络的抗毁性，即评估重要节点遭遇攻击后网络的连通性。网络在受到攻击时，其连通性越
9、高崩溃概率就越小，因此网络连通性是衡量网络崩溃程度的重要指标。网络的连通性主要体现在故障点数量、连通子图数量、最大连通子图的比率 以及网络连接效率上，故选用它们作为评估指标。其中故障点数量，即为网络受到攻击后，失效节点的数量。连通子图数量是指节点遭受攻击失效后划分为不同连通子图的数量。网络遭受攻击致部分节点失效后，完整网络被分割成多个互不联系的子网，其中弧边连接数及节点数最多的子网为最大连通子图，最大连通子图的节点数与原网络节点数的比值定义为最大连通子图的比率。越大，该网络遭受攻击时被破坏的程度越低。（）网络传输效率可以用网络连接效率衡量，即网络效率，节点和节点之间的网络效率为 ，（）其中，为
10、节点和节点之间的最短距离。整个网络的连通效率定义为所有节点的平均网络效率，即（），（），（）网络在遭受攻击致部分节点失效后，网络效率也会变化，本文提出网络效率变化率作为评估指标。假设被攻击前网络效率为，攻击后网络效率为，攻击前后节点数变化为。将攻击前后的网络效率变化与节点数量的变化的比值定义为网络效率变化率。越大表示其节点重要性越高，其表达式为 （）电力网络的攻击方法通常研究抗毁性时会采用两种攻击模式，分别为随机攻击和蓄意攻击。随机攻击是对网络中的节点或边进行一定概率的随机攻击和破坏。通常展现为非人为故意破坏造成的故障，例如硬件故障以及自然灾害造成的传输设备损坏，属于随机攻击。蓄意攻击是指对网
11、络的进攻是按照策略进行的。若网络连接中起关键作用的节点或边被人为破坏，此时攻击方已获知点与边的信息，清楚地知道网络中的重要节点和边，并按照亟定顺序攻击节点或边。根据静态参数的不同选用两种蓄意攻击策略：网络节点度优先攻击策略和实时接近中心性攻击策略。网络节点度优先攻击策略的攻击过程：）计算网络中所有节点的度，按照其节点度将所有节点排序；）找到度最大的节点并从网络中删除，攻击过程中计算网络效率；）按照）中的排序依次攻击，直至所有节点删除。实时接近中心性攻击策略的攻击过程：）计算网络中所有节点的接近中心性参数；）找到接近中心性最大的节点，并对其发起攻击使其失效；）重新计算剩余节点的接近中心性；）重复
12、）和），直至所有节点失效。在现有的基于介数中心性攻击策略 的研究基础上，本文提出实时攻击的方式，即移除最高接近中心性的节点后重新计算剩余网络的接近中心性，循环直至网络中所有节点被删除的攻击方式。采用实时接近中心性攻击方式时，给予网络反应的时间更短，攻击更具有侵略性。相较节点度优先攻击策略，实时接近中心性的策略更具有实时性，且接近中心性也更能反映节点的重要程度。复 杂 系 统 与 复 杂 性 科 学 年 月实证分析实证分析以上海市崇明区为例，因为：）崇明区包含着工厂，商场及农场住等不同的设施，整个系统的负荷多元化；）崇明区具备风能、太阳能、生物质能、潮汐能、地热能等种可再生能源，种类多，且可再生
13、能源装机容量占最高负荷的比重超过。随着技术及机制创新，崇明区将打造成风光储一体、农光互补和能源互联网的集成示范，针对它的研究具有前瞻性；）崇明区电力网络规模适中且具有典型性。以往对于电力系统抗毁性分析的文章缺少现实数据，我们根据崇明区 年电网地理接线图构建电力网络拓扑图如图所示，网络含有 个节点及 条边。图中的变电站，发电厂等站点均以英文缩写形式表示，英文字母的大小代表站点度的大小，且站点布局与现实布局不关联。接下来基于 软件进行了随机攻击和蓄意攻击的仿真实验，研究图所示拓扑图，得到仿真结果。通过复杂网络参数可知图的网路效率为 ，连通子图数为，最大连通子图比率为 。通过随机攻击实验和蓄意攻击实
14、验分析传输网络的抗毁性。由于蓄意攻击时，攻击方很难获取边的准确位置，因此仅采用攻击节点的攻击方式，本文蓄意攻击包括网络节点度优先攻击和实时接近中心性攻击两种策略。随机攻击仿真结果在电力网络遭受随机攻击时，都会使得站点无法正常工作，此类状况一般不会导致多个节点同时失效。采用逐一随机攻击的方式，网络效率以及最大连通子图数量同被攻击节点数的关系分别如图和图所示。图中纵坐标前后分别代表攻击站点和网络效率变化率，通过分析网络效率变化率，进而判断节点在网络中的重要性。图上海市崇明区电力网络拓扑图 图随机攻击策略网络效率图 图随机攻击策略最大连通子图变化 网络被随机攻击后变化程度如图、图及表所示。由图、图和
15、表可知节点 （长兴站），（堡北），（陈家镇），（团结站），（中双港），（长江站）在遭受攻击后对整个网络的联通影响较大，即该类节点为枢纽节点，需要加强保护。表随机攻击下枢纽节点网络效率变化率 站点网络效率变化率（）蓄意攻击仿真结果采用节点度优先攻击策略和实时接近中心性优先攻击策略对网络展开蓄意攻击，分别按照制定的策略选择 个点对网络采取攻击，其中实时接近中心性优先攻击策略是作者根据前人 的研究改进提出的一种蓄意攻击的策略。相对于随机攻击的偶然性，蓄意攻击类似于黑客攻击，具有针对性及破坏性。通过比较网络的变化程度得出网络的重要节点及变化规律，实验仿真结果如图及图所示。由图和图可知崇明区电力网络遭遇
16、数次攻击后，两种策略下前三次攻击后网络连通效率下降较多，由此可见做好前期应对并及时处理在保护电网中作用巨大。对比两图，采用节点度优先攻击策略且网络受到 次攻击以后，网络效率下降到 ，失效点数目为，整个网络几乎处于瘫痪状态；当采用实时接近中心性优先攻击策略时，遭受到 次攻击后，网络效率下降至 ，失效点数目为，连通子图数量第 卷第期郭明健，等：基于复杂网络理论的电力网络抗毁性分析为，此时网络已瘫痪。不难发现网络效率变化速率越大，整个网络的连通性受到的影响就会越大，网络遭到的破坏性也就越大。图节点度优先攻击策略网络效率图 图实时接近中心性优先攻击策略网络效率图 将实验现象同文献 及文献 中的方法进行比较，可以得出以下几点规律：）无论是采用文献 中故障链的级联故障图，还是采用文献 中同时进行节点攻击及边攻击，或是实时接近中心性优先攻击策略时，针对枢纽节点的攻击会导致网络中 以上节点失效，崩溃程度大。）文献 中仅选用节点的出入度作为判断枢纽节点依据，依此选中的枢纽节点在遭到攻击后，网络效率下降速率几乎没有变化，无法有效判断枢纽节点的重要程度。文献 中选用节点的介数中心性同边的介数中心性的乘积作为