基于隶属度函数的模糊综合评价法的电网脆弱性评估方法

基于隶属度函数的模糊综合评价法的电网脆弱性评估方法
【专利摘要】本发明涉及一种基于隶属度函数的模糊综合评价法的电网脆弱性评估方法，包括以下步骤：1)定义电网静态安全性分析、动态安全性分析、电网拓扑结构安全性分析以及可靠性分析指标；2)基于隶属度函数进行脆弱性模糊综合评价。本发明面向负荷的连续可靠供电目标，从电网的电压、潮流、结构、可靠性等多方面综合衡量电网的脆弱性水平，脆弱性评分更加全面也更具参考价值，提高了电网的脆弱性状态评估水平。
【专利说明】
基于隶属度函数的模糊综合评价法的电网脆弱性评估方法
技术领域
[0001] 本发明涉及一种电网脆弱性评估方法，尤其是涉及一种基于隶属度函数的模糊综 合评价法的电网脆弱性综合评估方法。
【背景技术】
[0002] 电网脆弱性评估是从整体、系统和宏观的角度来分析系统的脆弱程度和潜在的薄 弱环节以及关键的系统元件，电力部门可以通过加强对脆弱环节及关键设备的监控来提高 系统的鲁棒性、降低系统潜在的风险和脆弱性。如何对电网的脆弱性进行全面而综合的评 估一直是一个备受关注的课题。
[0003] 对于电网脆弱性的评估大致分为以下几个方面。第一类是基于复杂系统理论的电 力系统脆弱性评估。其基本思路是以复杂网络理论为基础，从电网拓扑结构角度分析故障 的传播机理，寻找电网本身固有的脆弱性。它一般是基于复杂网络研究领域的小世界网络 以及现实世界网络研究成果，研究聚类系数以及特征路径长度等指标与电网脆弱性的关系 来评价网络的脆弱性。该类方法假设电网是一个无权无向网络，并且忽略支路阻抗的影响， 与实际不符，应加以改进。
[0004] 第二类是基于概率论的电力系统脆弱性评估。该类方法计算反映电网安全状态的 低电压风险、过负荷风险及综合风险等指标，针对电力系统电压安全和潮流分布的安全问 题，建立风险指标与电网脆弱性之间的关系，通过运行状态实时反映系统的脆弱程度。该类 方法可以在电网重负荷运行状态下辨识出电网的脆弱线路、脆弱节点，通过电网运行方式 的改善降低系统停电的概率，使系统摆脱临界状态，避免大停电事故的发生。但是这些方法 都是针对电力系统某个方面或某个部分的安全问题，并不能完全反映系统的全局安全状 ??τ 〇
[0005] 第三类方法是基于系统分析的电力系统脆弱性评估。该方法用来评估电网的动态 安全性，该方法下的脆弱性用来描述系统在正常运行情况或各种随机因素的作用下，系统 承受扰动或故障的能力和系统不能保持正常工作的发展趋势及其影响，它反映了某个或某 些参数发生变化时对系统安全稳定裕度的影响及其变化趋势。并通过脆弱指标的构建，对 系统当前安全水平及其变化趋势进行分析，由此确认电网中的脆弱环节。造成电网脆弱的 原因是多样的，该方法只能从一个侧面反映电网的脆弱特性，评估结果的参考价值有限。

【发明内容】

[0006] 本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种评价手段全面 有效、指标参考价值高、综合反映电网的静态及动态安全性能，考虑电网结构安全性及电网 可靠性的电网脆弱性评估综合方法。
[0007] 本发明的目的可以通过以下技术方案来实现：
[0008] -种基于隶属度函数的模糊综合评价法的电网脆弱性评估方法，包括以下步骤：
[0009] 1)完成与电网脆弱性相关的静态安全性、动态安全性、拓扑结构稳定性以及可靠 性的脆弱性指标的计算；
[0010] 2)基于隶属度函数计算电网脆弱性的模糊综合评分。
[0011] 所述步骤1)中，电网脆弱性的分析分为电网静态安全性分析、动态安全性分析、电 网拓扑结构安全性分析以及可靠性分析。各脆弱性分析的具体方法为：
[0012] 1.电网静态安全分析
[0013] 通过SCADA(数据采集与监视控制系统)采集的实时量测，由状态估计可获得电网 当前的状态信息，进而可通过静态安全评估指标的计算，得到电网当前状态下的潮流和电 压运行水平。本文定义的静态安全评估指标如下：
[0014] 1)电网支路平均潮流越限量AQB0(average quantity of branch overload)
[0015]
[0016]
[0017] 式中：NBQ为系统中过载的支路集合;QB0x为支路x的潮流越限量，表示支路x潮流越 限的严重程度，其值越大，越限越严重;L x为支路X的潮流(标么值，以下同），Lx,max为支路X允 许的最大传输有功功率。
[0018] AQB0表示当前状态下潮流越限支路平均越限的严重程度，其值越大，表明电网潮 流越限越严重。
[0019] 2)专路命裁7tC平ALRF(averaee level of branch flow)
[0020]
[0021]式中：N_为系统中未过载支路的集合;Lx为支路x的潮流，Lx,max为支路x允许的最 大传输有功功率。
[0022] ALBF表示当前状态下潮流未越限支路的平均负载水平，其值越大，表明当前电网 负载越重。
[0023] 3)电网节点电压平均越限量AQNVV(average quantity of node voltage violation)
[0024]
[0025]
[0026] 式中：Nw为系统中当前状态下电压越限的节点数;QNVVZ为节点z的电压越限量，表 示节点z电压越限的严重程度，其值越大，越限越严重;V z为节点z的电压，Vz,max为节点z电压 的上限值，Vz,_为节点z电压的下限值。
[0027] AQNVV表示当前状态下电压越限节点平均越限的严重程度，其值越大，表明系统电 压越限越严重。
[0028] 4)电网电压运行水平ALNV(average level of node voltage)
[0029]
[0030] 式中:Nnw为当前状态下电压未越限的节点数，VZ为节点z的电压。
[0031] ALNV表示当前状态下电压未越限节点的平均电压值，其值越大，表明电网电压水 平越高。
[0032] 5)低电压风险指标分析
[0033] 根据单一故障安全检测法则(N-1原则），网络安全运行的要求是网络中任意一条 线路发生故障时，系统的各项运行指标仍能满足要求。当采用随机概率的方法评估系统脆 弱性时，容易忽略发生概率极小但引发严重后果的故障，因此定义风险指标为事故发生的 概率与事故产生后果的乘积。
[0034] 从事故统计数据可以看出电力系统发生事故(这里只考虑三相短路故障）的概率 基本符合泊松分布，即
[0035]
[0036] 式中是可能事故Ei的发生率，P(Ei)表示事故Ei的发生概率。
[0037] 低电压风险反映的是系统发生的N-1故障造成系统中母线电压下降的可能性和危 害程度。当母线电压为l.Opu时，低电压严重度为0;随着电压降低幅值的增加，低电压严重 度也增加。
[0038] 设发生故障E,时母线j的电压为VB^(标幺值），定义母线j的低电压损失值为1^， 则有：
[0039]
[0040] 考虑到各母线重要程度的不同，电压偏移量相同的情况下其后果有所区别，因此， 定义母线j的低电压严重度为：
[0041]
[0042] 其中，
[0043] αΒ」= aB-aj X aB-bj X aB-Cj
[0044] 式中，αΒ^为母线重要程度因子，€^"、<^^、€^。」分别对应于该母线地位、所连设备 类型、电压等级权重因子，当不计及某因素影响时，可取对应权重值为1，某影响因素影响程 度越大，对应权重值越大于1，具体取值参见附录。
[0045] Ν-1故障时，对应于事故Ei的低电压严重度为故障后系统母线电压低电压严重度 之和：
[0046]
[0047] N-1故障系统低电压风险指标的计算公式为：
[0048]
[0049] 式中，Xt表示系统t时刻的状态，Rv(Ei|Xt)表示在t时刻发生事故Ei时的低电压风 险。6)过负荷风险指标分析
[0050] 线路过负荷风险反映的是系统发生N-1故障导致线路过载的可能性和危害程度。 当线路有功潮流小于或者等于传输功率极限Pmax的90%时，过负荷严重度为0;随着线路传 输功率的增加，对应线路的过负荷严重度增大，且增加速率相应变快。
[0051] 设发生故障电力系统在线小干扰主导特征模式识别及强相关机组选择方法
[0052] 时线路j的有功功率为Pij(实际有功与传输功率极限的比值），定义线路j的过负荷 损失值为wPij，则有：
[0053]
[0054]考虑到线路重要程度的不同，等量过载的情况下其后果有所区别，因此，定义线路 j的过负荷严重度为：
[0055]
[0056] 其中，
[0057] Pj=Paj X^bj X^cj
[0058] 式中也为线路重要程度因子，知、知、知分别对应于该线路介数、线路电压等级、 线路所供负荷等级权重因子，当不计及某因素影响时，可取对应权重值为1，某影响因素影 响程度越大，对应权重值越大于1。
[0059] N-1故障时对应于事故Ei的过负荷严重度为故障后系统未开断线路过负荷严重度 之和：
[0060]
[0061 ] N-1故障系统过负荷风险指标的计算公式为：
[0062]
? i
[0063] 式中，Xt表示系统t时刻的状态，RP(Ei|Xt)表示在t时刻发生事故Ei时的过负荷风 险。
[0064] 通过对系统以上6个指标的计算，完成电网静态安全分析。
[0065] 2.电网动态安全分析
[0066] 电网的动态安全分析包括两部分内容，分别是电网小扰动稳定性指标分析和电网 故障切除时间裕度指标分析。
[0067] 1)电网小扰动稳定性指标分析
[0068]采用特征分析法进行多机系统的小干扰动态稳定分析，采用发电机三阶模型，考 虑励磁系统与PSS。
[0069]电力系统的动态特性可以由一组非线性微分方程组和一组非线性代数方程组描 述：
[0070]
[0071] 在平衡点xo附近线性化，把各变量表示为其初始值与微增量之和：
[0072] Xi = Xi〇+ Δ xi
[0073] 将所得方程组在平衡点XQ附近展开成泰勒级数，并略去各微增量的二次及高次 项，得：
[0074]
[0075] 定义状态向量和非状态向量如下：
[0076] AX=[xi,x2,---,xm]T
[0077] AY=[xm+i,xm+2,---,xn]T
[0078] 将上式写成矩阵形式：
[0079]
[0080]在上式中消去非状态变量得线性化方程：
[0081]
[0082] 即得到描述线性系统的状态方程，其中A为ηΧη维系数矩阵，称为该系统的状态矩 阵。
[0083] 对于由状态方程描述的线性系统，其小干扰稳定性由状态矩阵的所有特征值决 定。由李雅普诺夫第一稳定性定理，如果线性化方程的所有特征根都具有负实部，则原系统 的平衡状态是渐近稳定的；如果至少具有一个根具有正实部，则原系统的平衡状态不稳定； 如果特征根有实部为零的情况，则原系统平衡状态为临界情况。因此，分析系统的小干扰稳 定性问题，可以归结为求解状态矩阵Α的全部特征值的问题。
[0084] 复数特征值总是以共辄对的形式出现，每一对对应一个振荡模式。对于一对复数 特征值：
[0085] λ ' = σ ' ± j ω '
[0086] 振荡频率(Hz)由下式给出：
[0087]
[0088] 阻尼比由下式给出：
[0089]
[0090] 当ξ<0,则该模式是不稳定的；当ξ = 0,则该模式处于稳定边界；当ξ>0,则该模式 是稳定的，且ξ越大，该模式稳定阻尼越强。
[0091 ]使用PSD-SSAP进行小干扰动态稳定分析时，全特征值(QR)法会计算得到状态矩阵 的全部特征值，与对应于每个特征值的机电回路相关比Pl(由左特征向量和右特征向量 决定hPi用于从所有特征根中选出和某些变量强相关的特征根。实际电力系统中，若某个 特征根λ'i满足
[0092]
[0093] 则λ ' ,为低频振荡模式，又称机电模式。
[0094] -般，η机系统存在(η-1)个机电模式，阻尼最低或者最危险的模式即存在于机电 模式中，只有(η-1)对机电振荡模式的特征值均具有负实部时系统才被认为是稳定的。因 此，选择机电模式中的最小阻尼比作为电网小扰动稳定分析指标
[0095] ξπι=ηι?η(ξι,ξ2, ··· ,^η-ι)
[0096] 2)电网故障切除时间裕度指标分析
[0097]电力系统分析综合程序(PSASP)是由中国电力科学研究院研发的电力系统分析程 序。主要用于电力系统规划设计人员确定经济合理、技术可行的规划设计方案，供运行调度 人员确定系统运行方式、分析系统事故、寻求反事故措施，供科研人员研究新设备、新元件 投入系统等新问题以及高等院校用于教学和研究。这里我们运用其暂态稳定计算功能，来 计算电网的故障临界切除时间CCT。
[0098]利用PSASP的常用系统元件模型搭建目标电网络，给定网络需要的参数，进行CCT 的计算。对于输电支路数为nCCT(双回路按照两条支路计算，不计发电机到母线的支路)的电 网络其具体计算流程如下：
[0099] Stepl:对电网络所有输电支路进行编号，从1到nCCT，并令线路编号i = l，迭代次数 j = l，设定迭代误差ε〇 = 0.001;
[0100] Step2:在支路i的50%处发生三相短路故障，并切除故障，给定故障切除时间为to (to值参考正常的故障临界切除时间给定一个较大值），令<
[0101] Step3:在故障切除时间为?/时，利用PSASP进行仿真计算电网络各节点的电压曲 线；
[0102] Step4:观察各个节点电压曲线的电压有无明显降低，是否发生电压失稳，若某个 节点发生电压失稳，令〖/+1 =〖/./2_，j = j+Ι转到Step3;若电压无明显降低，未发生电压失稳， f+1 = (g +/广）/2 ,判断〇/+1 -//) <￡'u.是否成立，若成立，令j = j+Ι，转到Step5，若不成立，令 j = j+Ι，转到Step3;
[0103] Step5:记录支路i的故障临界切除时间为& =i/ .，若i<nccT，令i = i+l，转到Step2; 若i=ncxT，仿真结束。
[0104] 计算得到各线路的故障临界切除时间^后，在给定各条线路的额定故障切除时间 V i时，计算系统各线路的故障切除时间裕度为：
[0105]
[0106] 选取所有线路故障切除时间裕度的最小值作为电网的故障切除时间裕度指标 ICCT，当其满足动态稳定性要求时，所有线路均能满足要求：
[0107]
[0108] ICCT越大，则故障切除具有的安全裕度越大，电网的稳定性能越好。
[0109] 通过对系统以上2个指标的计算，完成电网动态安全分析。
[0110] 3.电网拓扑结构安全性分析
[0111] 对于电网拓扑结构的分析，我们以凝聚度来表示节点的重要度。以往习惯上以节 点度数作为衡量节点重要度的标准，然而，在许多实际网络中的某些关键节点并不一定具 有较大的度数，因此采用节点收缩后的网络凝聚度来评估网络中的节点重要度具有更高的 合理性。网络中节点i的凝聚度的计算方法为：
[0112]
[0113] 式中:m。为节点i的收缩网络中节点的数目，h为节点i收缩网络中节点之间的平均 最短路径，dmin, gh是用边的数目表示的网络中任意两节点g和h之间的最短距离，V表示节点i 收缩网络中所有节点组成的集合。
[0114] 节点i的凝聚度越高，则i节点越重要。该网络各节点之间的联系越强，则该网络受 扰动后其崩溃的可能性越大，电网越脆弱。计算网络的平均凝聚度作为网络拓扑结构稳定 性的衡量指标，平均凝聚度指标为：：
[0115]
[0116] 式中：Nc为电网络中节点的数目。
[0117] 平均凝聚度越大，网络脆弱性越高。
[0118] 4.电网可靠性分析
[0119] 根据电网的停电情况，可计算得到电网的供电可靠性指标，供电可靠性指标衡量 了电网的持续供电能力。
[0120] 为了衡量每个由系统供电的用户在单位时间内所遭受到的平均停电次数，计算 SAIFI (系统平均停电频率指标)为：
[0121]
么1
[0122] 式中，Ni为节点i的用户数，Asusi为节点i的永久故障率，Nsus_i为受永久故障影响的 用户数，Ντ为电网全部用户数，Ντ= ΣΚ。
[0123] 为了衡量系统供电能力的可靠性的大小，计算RS(供电可靠率指标)为：
[0124]
[0125] 式中，Ni为节点i的用户数，U^d+i为节点i的年平均停电时间(小时/年）。
[0126] 完成4个方面11个指标的计算后，进行所述步骤2)中基于隶属度函数的模糊综合 评价法的电网脆弱性评价。电网系统的一个状态即是一个方案，电网脆弱性评价是对不同 方案的11个指标模糊综合评价计算最后的综合评分Zl，Zl既能充分利用客观数据所蕴含的 信息，又能反映数据在实际应用时的重要程度。
[0127] 1.计算脆弱性评估指标的权重向量
[0128] 假定有m个待评价的方案，每个待评价方案有N个评价指标，在脆弱性评估过程中， 各指标权重的确定是一个非常重要的环节。目前，对于各指标权重的确定方法主要分为主 观赋权法和客观赋权法。考虑到判断矩阵的一致性与求解精度和收敛速度的问题，采用主 观赋权的三标度法求取指标权重。
[0129] (1)将各指标重要度进行两两比较，三标度法建立互补型模糊判断矩阵，各标度含 义如下表所示：
[0130] 表1标度法参照表
[0131]
[0132] 表中6、心表示作比较的两个指标，分别表示两指标的重要程度。根据表 中含义将各指标两两进行比较，形成互补型主观权重模糊判断矩阵F=(f^)NX N:
[0133]
[0134] (2)将互补型模糊判断矩阵转变为模糊一致性判断矩阵R=(rij)NXN:
[0139] ri表示模糊判断矩阵F第i行的和。[0140] (3)根据和行归一化方法求解得到排序向量初值w(()):
[0135]
[0136]
[0137]
[0138]
[0141]
[0142]
[0143]
[0144]
[0145] t表示模糊一致性判断矩阵R第i行的和，r_表示的和，wf1即表示第i个指 标占所有指标和的权重的初值。
[0146] (4)将模糊一致性判断矩阵转变成互反型矩阵E=(eij)NXN:
[0150] (5)选取排序向量w(())作为迭代初始值，确定误差ε;[0151] (6)求取特征向量的初值w(())的无穷范数| |w(())| μ，置k=l:
[0147]
[0148]
[0149]
[0152]
[0153] (7)求特征向量迭代值w(k):
[0154]
[0155]
[0156] (8)进行判断
[0157] 若11 |w(k)| μ-| μΗ" μι |<ε，所得到的向量w(k)就是目标的权重向量，迭代结 束;否则置k = k+l，转步骤7)，继续计算。
[0158] 迭代结束后所得到的权重向量即为最终的权重向量：
[0159]
[0160] 2.确定指标集与评语集
[0161 ] m个待评价方案，N个评价指标。其中第i个方案的指标集为Ui = {uii，Ui2，…，UiN}，1 <i<m，其中为第i个方案的第j个指标数值；
[0162 ]对指标有1种评语，评语集V = {vi，V2，…，vi}，此处1 = 5，V = {优，良，中，合格，差};
[0163] 3.评价矩阵标准化
[0164] 由评价指标集可得m个评价方案的评价矩阵U:
[0165]
[0166] 其中，1彡i彡m，l彡j彡N;
[0167] 评价方案的评价指标分为效益型和成本型两类:效益型指标是指该指标相对其他 方案越大说明该方案越优秀的指标，而成本型指标是指该指标相对其他方案越小说明方案 越优秀的指标。脆弱性评价的各项指标中，系统静态安全性指标、动态安全风险指标、电网 可靠性指标为脆弱性的成本型指标;电网小扰动稳定性指标、电网故障切除时间裕度指标、 电网平均凝聚度指标为脆弱性的效益型指标。将U标准化后得到标准评价矩阵C:
[0168]
[0169] 叫为Ulj标准化后元素，标准化方法为：
[0170]
[0171] 其中，效益型指标按上述公式第1式处理，成本型指标按上述公式第2式处理；
[0172] 4.构造隶属度函数模糊评价矩阵
[0173] 对于m个待评价方案，第i个方案标准化后的评价指标集为，…，clN}，l < i ，其第j个指标在评价集V上的模糊子集可通过其对V的隶属函数来计算。
[0174]等腰三角形隶属度函数：
[0175]
[0176] 式中：rij(vk)为第i个方案的第j个指标相对于评语vk的隶属度;pk、qk、sk为对应于 Vk的常数，qk的值对应于第k个评语而取定。
[0177] 等腰梯形隶属度函数：
[0178]
[0179] 式中:rij(vk)为第i个方案的第j个指标相对于评语vk的隶属度;pk、qk、sk、Uk为对应 于Vk的常数，qk、Vk的值对应于第k个评语而取定。
[0180] 高斯隶属度函数：
[0181] 根据
[0182]
[0183]
[0184]
[0185] 式中：参数u用于确定曲线的对称轴，参数〇用于确定曲线的跨度。rij(Vk)为第i个 方案的第j个指标相对于评语Vk的隶属度;u k为对应于vk的常数。
[0186] 由此，基于三种隶属度函数，可得出第i个方案的模糊评价矩阵分别为 (3)
Ο
[0187]
[0188] 5.求取综合评价模糊子集
[0189]第i个方案的综合评价集msv上的模糊子集，通过下式计算：
[0190] Bi=wc〇Ri= {bii,bi2,bi3,bi4,bi5}
[0191] 式中，其中we为综合权重集为第i个方案的模糊评价矩阵，算子。采用M( ·，+ )模 型，则有：
[0192]
[0193] 对Bi进行归一化处理：
[0194] Bi
[0195] 可得到第i个方案的模糊综合评价结果为：
[0196]
[0197] 其中，= 为第i个方案相对于评语vk的隶属度，表示方案i在多大程度 上可以被评语Vk描述。
[0198] 6.利用模糊综合评价结果对方案进行评分
[0199] 给评价结果的每一个评语赋予一个分值而将其量化，可得评分集_ = …4}。 对于评语集V= {:优，良，中，合格，差} (1 = 5)，设定将其量化后的评分集f = !90,80,70,60,50丨.
[0200] 方案i的评价结果向量
[0201] i
[0202] Zlk=iikVk
[0203] 式中，zlk表示方案i相对于评语vk的得分。那么，方案i的综合得分为电力系统在线 小干扰主导特征模式识别及强相关机组选择方法。
[0204]
[0205] m个方案的综合评价矩阵
[0206] 等腰三角形隶属度函数、等腰梯形隶属度函数和高斯型隶属度函数的综合评价矩 阵分别）
[0207] 7.组合系数加权求取三种函数的综合评价结果向量
[0208] 三种隶属度函数在评价方案的过程中各有优势，在评价结果方面有一致性也有差 异性。采用向量相似度原理，分别求得不同隶属度函数下评价结果向量之间的向量相似度， 相似度大的两个评价向量在计算综合评价向量时应予以更多考虑。采用组合系数加权的方 法进行综合评价向量的计算。
[0209] 对于任意两个_隹向量父=(1142，~411)，￥=(71，72，~，711)，其内积表示为：
[0210] [Χ,Υ] =Xiyi+X2y2+---+Xnyn [0211]向量的范数表示为：
[0213] " " l[0214] 向量的夹角表示为：
[0212]
[0215]
[0216] 则两个向量的长度相似度.α .为：
[0217]
[0218]两个向量的方向相似度β为：
[0219]
[0220]两个向量的总相似度γ为：
[0221] γ (χ,Υ)=α(Χ,Υ) · β(Χ,Υ)
[0222] 对于方案i，等腰三角形隶属度函数得到的综合评价结果向量
等腰梯形隶属度函数得到的综合评价结果向量 ^
j，高斯隶属度函数得 到的综合评价结果向量
]，设基于等腰三角形-等腰梯形-高斯隶属度 函数的综合评价结果向量Zf计算公式如下：
[0223]
[0224]式中：m、n2、n3为组合系数，求取方法如下：
[0225] 计算、Z!2)、中任一评价结果向量与其他向量的总相似度，得到：
[0226] γ =( γ 1, γ2, γ3)
[0227] 其中
示向量Zf1 1 2 3)与其他两个向量相似度总和。
[0228] 组合系数向量：
[0229]
[0230] 其中
，表示向量与其余向量的相对相似程 度，组合系数大的向量，在计算三种函数的综合评价结果向量时会给予更多的考量。
[0231] 与现有技术相比，本发明具有以下优点：
[0232] (1)本发明建立了基于隶属度函数的模糊综合评价电网脆弱性的评估方法，该方 法综合考虑了电网的静态稳定、动态稳定、拓扑结构稳定以及可靠性四个方面的指标，使得 对电网脆弱性的评价更加全面，得到的脆弱性评价结果可信度更高。
[0233] (2)本发明通过隶属度函数模糊综合计算得到脆弱性评分，该方法不仅能够保持 一定的计算精度，同时能够实现单方案的评估，具有良好的工程应用前景。
【附图说明】
[0234] 图1为本发明的流程示意图；
[0235] 图2为本发明实施例中三机九节点的实验系统的示意图。
【具体实施方式】
[0236] 下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案 为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于 下述的实施例。
[0237] 如图1所示，本实施例提供一种基于隶属度函数的模糊综合评价法的电网脆弱性 综合评估方法，包括步骤：
[0238] S1、完成与电网脆弱性相关的静态安全性、动态安全性、拓扑结构稳定性以及可靠 性的脆弱性指标的计算。
[0239] S101、电网静态安全分析
[0240] 1)计算电网支路平均潮流越限量AQB0:
[0241]
[0242]
[0243]式中：Νβο为系统中过载的支路集合;QB0x为支路X的潮流越限量，表示支路X潮流越 限的严重程度，其值越大，越限越严重;Lx为支路X的潮流(标么值，以下同），Lx,max为支路X允 许的最大传输有功功率。
[0244] 2)计算支路负载水平ALBF:
[0245]
[0246] 式中：N_为系统中未过载支路的集合;Lx为支路X的潮流，Lx,max为支路X允许的最 大传输有功功率。
[0247] 3)计算电网节点电压平均越限量AQNVV:
[0248]
[0249]
[0250] 式中：Nw为系统中当前状态下电压越限的节点数;QNVVZ为节点z的电压越限量，表 示节点z电压越限的严重程度，其值越大，越限越严重;V z为节点z的电压，Vz,max为节点z电压 的上限值，Vz,_为节点z电压的下限值。
[0251] 4)计算电网电压运行水平ALNV:
[0252]
[0253] 式中:Nnw为当前状态下电压未越限的节点数，VZ为节点z的电压。
[0254] 5)计算低电压风险指标Rv(E|Xt):
[0255] 从事故统计数据可以看出电力系统发生事故（只考虑三相短路故障）的概率基本 符合泊松分布，BP
[0256]
[0257] 式中是可能事故Ei的发生率，P(Ei)表示事故Ei的发生概率。
[0258]发生故障E1时母线j的电压为VB^(标幺值），母线j的低电压损失值为w VB_1J:
[0259]
[0260] 母线j的重要程度因子αΒ」：
[0261 ] Ctg j -〇B_aj X CtB bj X CtB cj
[0262] 其中，(^^、(^^^分别对应于该母线地位^斤连设备类型冰压等级权重因子， 当不计及某因素影响时，可取对应权重值为1，某影响因素影响程度越大，对应权重值越大 于1。
[0263] 母线j的低电压严重违
[0264]
[0265] 对应于事故低电压严重度为故障后系统母线电压低电压严重度之和SsKEO:
[0266] }-
[0267] N-1故障系统低电压风险指标的计算公式为：
[0268]
[0269]式中，Xt表示系统t时刻的状态，Rv(Ei|Xt)表示在t时刻发生事故Ei时的低电压风 险。
[0270] 6)计算过负荷风险指标RP(E|Xt):
[0271] 发生故障Ei时线路j的有功功率为Pij(实际有功与传输功率极限的比值），线路j的 过负荷损失值为WPij :
[0272]
[0273]线路j的重要程度因子氏：
[0274] Pj = Paj X^bj X^cj
[0275] 其中，知、知、知分别对应于该线路介数、线路电压等级、线路所供负荷等级权重 因子，当不计及某因素影响时，可取对应权重值为1，某影响因素影响程度越大，对应权重值 越大于1。
[0276] 线路j的过负荷严重度SP(Wij):
[0277]
[0278] 对应于事故Ei的过负荷严重度为故障后系统未开断线路过负荷严重度之和SsP (Ei)：
[0279] J
[0280] N-1故障系统过负荷风险指标的计算公式为：
[0281]
[0282] 式中，Xt表示系统t时刻的状态，RP(Ei|Xt)表示在t时刻发生事故Ei时的过负荷风 险。
[0283] S102、电网动态安全分析
[0284] 1)计算电网小扰动阻尼比指标
[0285] 在电网的稳定运行点处进行电网潮流计算，将描述系统动态特性的非线性方程在 平衡点（电网稳定运行点）x〇附近线性化，其中发电机采用三阶模型，考虑励磁系统与PSS的 作用，最终得到描述系统动态特性的状态矩阵A nXn(n机系统）。
[0286] Αηχη的一对复数特征根对应着一个振荡模式：
[0287] λ'=σ'土 jco'
[0288] 振荡频率(Hz):
[0289]
[0290] 阻尼比：
[0291]
[0292] 根据每个振荡模式Μ的振荡频率fzl与对应的机电回路相关比Pl，按下述条件选出 (n-1)个机电模式。
[0293]
[0294] 求解出(n-1)个机电模式的阻尼比…，，选取最小阻尼比作为电网小扰动 稳定分析指标
[0295] ξπι=ηι?η(ξι,ξ2, ··· ,^η-ι)
[0296] 2)计算电网的故障切除时间裕度指标
[0297] 利用PSASP的常用系统元件模型搭建目标电网络，给定网络需要的参数，进行CCT 的计算。对于输电支路数为nCCT(双回路按照两条支路计算，不计发电机到母线的支路)的电 网络其具体计算流程如下：
[0298] Stepl:对电网络所有输电支路进行编号，从1到nCCT，并令线路编号i = l，迭代次数 j = l，设定迭代误差ε〇 = 0· 0001;
[0299] Step2:在支路i的50%处发生三相短路故障，并切除故障，给定故障切除时间为to (to值参考正常的故障临界切除时间给定一个较大值），令=A;
[0300] Step3:在故障切除时间为?/_时，利用PSASP进行仿真计算电网络各节点的电压曲 线；
[0301] Step4:观察各个节点电压曲线的电压有无明显降低，是否发生电压失稳，若某个 节点发生电压失稳，令(产1=(: //2，」=」+1转到5丨叩3;若电压无明显降低，未发生电压失稳， C+1 ? +d/.2，判断(?/+1 - () < 是否成立，若成立，令j = j+ι，转到Step5，若不成立，令 j = j+Ι，转到Step3;
[0302] Step5:记录支路i的故障临界切除时间为（若i<nccT，令i = i+l，转到Step2; 若i=ncxT，仿真结束。
[0303]计算得到各线路的故障临界切除时间^后，在给定各条线路的额定故障切除时间 1^ = 0.1 s时，计算系统各线路的故障切除时间裕度，选取最小值作为电网故障切除时间裕 度指标>τ·
[0304]
一 … J
[0305] S103、电网拓扑结构安全性分析
[0306]计算电网络中节点i的凝聚度为：
[0307] L g，辦
[0308] 式中，m。为节点i的收缩网络中节点的数目，1 i为节点i收缩网络中节点之间的平均 最短路径，dmin,gh是用边的数目表示的网络中任意两节点g和h之间的最短距离，V表示节点i 收缩网络中所有节点组成的集合。
[0309]计算网络的平均凝聚度作为网络拓扑结构稳定性的衡量指标，平均凝聚度指标 为：
[0310]
[0311]式中，N为电网中所有节点数目。
[0312] S104、电网可靠性分析
[0313]分别计算SAIFI(系统平均停电频率指标）、.RS.(供电可靠率指标)如下：
[0314]
[0315]
[0310]式中，Ni为节点i的用户数，Asusi为节点i的永久故障率，Nsus_i为受永久故障影响的 用户数，Ντ为电网全部用户数，Ντ= ΣNi，ULciad_i为节点i的年平均停电时间（小时/年）。
[0317] S2、基于熵权模糊算法计算电网的脆弱性评分
[0318] S201、计算脆弱性评估指标的权重向量
[0319] 假定有m个待评价的方案，每个待评价方案有11 (N= 11)个评价指标。
[0320] 1)三标度法建立互补型模糊判断矩阵F=(fij)NXN:
[0321]
[0322] 其中：
[0323]
[0324] 式中，#、//分别表示第i、j个指标的重要程度。
[0325] 2)将互补型模糊判断矩阵转变为模糊一致性判断矩阵R=(r^)NXN:
[0326]
[0327]
[0328]
[0329]
[0330] 式中，ri表示模糊判断矩阵F第i行的和。
[0331 ] 3)根据和行归一化方法求解得到排序向量初值w(()):
[0332]
[0333]
[0334]
[0335]
[0336] 式中，t表示模糊一致性判断矩阵R第i行的和，r_表示Nfk的和，wf即表示第 i个指标占所有指标和的权重的初值。
[0337] 4)将模糊一致性判断矩阵转变成互反型矩阵E=(e^)NXN:
[0338]
[0339]
[0340]
[0341] 5)选取排序向量w(())作为迭代初始值，确定误差ε;
[0344] 7)求特征向量迭代值w(k):
[0342] 6)求取特征向量的初值w(())的无穷范数I |w(())| μ，置k = l:[0343]

3D L W …」
[0345]
[0346]
[0347] 8)进行判断
[0348] 若11 |w(k)| μ-1 μΜ? μι<ε，所得到的向量w(k)就是目标的权重向量，迭代结束； 否则置k = k+l，转7)，继续计算。
[0349 ]迭代结束后所得到的权重向量即为最终的权重向量：
[0350]
[0351] S202、确定指标集与评语集
[0352] m个待评价方案，11个评价指标。其中第i个方案的指标集为Ui={uii，Ui2,…，um}， 1彡i彡m，其中m j (1 < j < 11)为第i个方案的第j个指标数值。
[0353] 对指标有5种评语，评语集V= {:优，良，中，合格，差}，评分集V = {90s80,70,6〇,50;。
[0354] S203、评价矩阵标准化
[0355] 由评价指标集可得m个待评价方案的评价矩阵U(N=11):
[0356]
[0357] 评价方案的评价指标分为效益型和成本型两类:效益型指标是指该指标相对其他 方案越大说明该方案越优秀的指标，而成本型指标是指该指标相对其他方案越小说明方案 越优秀的指标。脆弱性评价的各项指标中，系统静态安全性指标、动态安全风险指标、系统 平均停电频率指标、电网平均凝聚度指标为脆弱性的成本型指标;电网小扰动稳定性指标、 电网故障切除时间裕度指标、电网平均凝聚度指标、供电可靠率指标为脆弱性的效益型指 标。将U标准化后得到标准评价矩阵C:
[0358]
[0359] 叫为叫标准化后元素，标准化方法为：
[0360]
[0361]其中，效益型指标按上述公式第1式处理，成本型指标按上述公式第2式处理。
[0362] S204、构造隶属度函数模糊评价矩阵
[0363]对于m个待评价方案，第i个方案标准化后的评价指标集为，…，clN}，l < i ，其第j个指标在评价集V上的模糊子集可通过其对V的隶属函数来计算。
[0364]等腰三角形隶属度函数：
[0365]
[0366] 式中：rij(vk)为第i个方案的第j个指标相对于评语vk的隶属度;pk、qk、sk为对应于 Vk的常数，qk的值对应于第k个评语而取定。
[0367] 等腰梯形隶属度函数：
[0368]
[0369] 式中:rij(vk)为第i个方案的第j个指标相对于评语vk的隶属度;pk、qk、sk、Uk为对应 于Vk的常数，qk、Vk的值对应于第k个评语而取定。
[0370] 高斯隶属度函数：
[0371]
[0372] 式中：参数u用于确定曲线的对称轴，参数〇用于确定曲线的跨度。rij(Vk)为第i个 方案的第j个指标相对于评语vk的隶属度;u k为对应于vk的常数。
[0373] 由此，基于三种隶属度函数，可得出第i个方案的模糊评价矩阵分别为 (3) Ο
[0374]
[0375] S205、求取综合评价模糊子集
[0376] 第i个方案的综合评价集msv上的模糊子集，通过下式计算：
[0377] Bi=wc〇Ri= {bii,bi2,bi3,bi4,bi5}
[0378] 式中，其中/为综合权重集为第i个方案的模糊评价矩阵，算子〇采用M( ·，+ )模 型， 则有：
[0379]
[0380] 对Bi进行归一化处理：
[0381]
k，l
[0382] 可得到第i个方案的模糊综合评价结果为：
[0383]
[0384] S206、利用模糊综合评价结果对方案进行评分
[0385] 对于评语集V= {:优，良，中，合格，差}，设定将其量化后的评分集兮=丨()(),80,7(),60,50;,
[0386] 方案i的评价结果向量
[0387] Zi= [zii,Zi2, ··· ,Ζ?δ]
[0388]
[0389] 式中，zlk表示方案i相对于评语vk的得分。那么，方案i的综合得分为
[0390] m个方案的综合评价矩阵ZKZiZvsZj。
[0391]等腰三角形隶属度函数、等腰梯形隶属度函数和高斯型隶属度函数的综合评价矩 阵分别爻
[0392] S207、组合系数加权求取三种函数的综合评价结果向量，得到最终评分
[0393]利用向量相似度原理，对综合评价结果向量计算中的各向量组合系数进行求取。 以方案i的等腰三角形隶属度函数综合评价矩K
与等腰梯形隶属度 函数综合评价矩阵
3例，计算其相似度，并计算Zf5的组合系数。[0394]长度相似度α为：
[0396][0397] 方向相似度β为：
[0395]
[0398]
[0399] 其中，
[0400]
[0402] 两个向量的总相似度Y为：
[0401]
[0403]
[0404] 同样的方法计算得Zf)与Zf)的向量相似度。
[0405] Z?)与其他向量相似度总矛丨同理可得Ζ{2)与Zf)对
应的γ 2与γ 3。
[0406] ζΡ的组合系数
，同理可得ζΡ与zf)的组合系数η2与η 3。
[0407] 因此基于等腰三角形-等腰梯形-高斯隶属度函数的方案i的综合评价结果向量
[0408]为了验证基于隶属度函数模糊综合评价的脆弱性评估方法的有效性，以美国西部 电网三机九节点系统做为试验系统。设置节点1为平衡节点，节点2、节点3为PV节点，其余节 点为PQ节点，计算系统潮流，并根据该三机九节点系统在120%负荷与80%负荷运行状态下 给定系统运行的脆弱性指标最大与最小状态阀值，同时计算待评估系统各项脆弱性指标 值，其中，各线路额定故障切除时间给定1/ i = 0.1 s。计算结果如下：
[0409]表2电网脆弱性评估指标值
[0410]
[0411] 给定主观权重模糊判断矩阵为：
[0412]
[0413] 按照S2的基于隶属度函数的模糊综合评价的步骤进行当前系统脆弱性评分计算， 得到：
[0414] 基于等腰三角形隶属度函数、等腰梯形隶属度函数、高斯隶属度函数的评分计算 结果分别为67.9、67.95和66.03,组合系数向量为? = (〇·371,0.338,0.291):，则对于该电网系 统进行脆弱性评估最终的得分为67.37。该结果表明系统的脆弱性程度较低，系统的坚强程 度在及格线以上。该结果体现了本文方法的有效性。
[0415] 该评分综合体现了电网系统静态、动态、网络结构、可靠性多个方面的指标，全面 衡量了电网的脆弱性状态，并且衡量了电网的脆弱性程度，具有很高的参考价值。以上所述 的具体实施例仅为说明本发明的实现效果，并不用以限制本发明。凡在本发明所提出的方 法的基本思路和框架之内所作的任何非实质性的修改、转换和改进，均应包含在本发明的 保护范围之内。
【主权项】
1. 一种基于隶属度函数的模糊综合评价法的电网脆弱性评估方法，其特征在于，所述 方法包括W下步骤： 1) 定义并计算与电网脆弱性相关的静态安全性指标、动态安全性指标、拓扑结构稳定 性指标W及可靠性的脆弱性指标； 2) 基于隶属度函数的模糊综合评价法进行电网的脆弱性评分。2. 根据权利要求1所述的基于隶属度函数的模糊综合评价法的电网脆弱性评估方法， 其特征在于，所述静态安全性指标包括如下： 21) 电网支路平均潮流越限量AQBO:式中：Nbo为系统中过载的支路集合，QBOx为支路X的潮流越限量，Lx为支路X的潮流标么 值，W下同，Lx,max为支路X允许的最大传输有功功率； 22) 支路负载水平ALBF:式中：NnBO为系统中未过载支路的集合，Lx为支路X的潮流，Lx,max为支路X允许的最大传 输有功功率； 23) 电网节点电压平均越限量AQNVV:式中：Ννν为系统中当前状态下电压越限的节点数，QNVVz为节点Z的电压越限量，Vz为节 点Z的电压，Vz,"ax为节点Z电压的上限值，Vz,"in为节点Z电压的下限值； 24) 电网电压运行水平ALNV:式中:Nnvv为当前状态下电压未越限的节点数，Vz为节点Z的电压； 25) 低电压风险指柄；Rv巧I Xt):式中：Xt是t时刻系统的运行状态;P化i)为不确定事故Ei发生的概率;Ssv化i)是发生事 故El时系统母线电压低电压严重度;Rv巧i|Xt)是发生事故El时系统低电压风险值； 26) 过负荷风险指标Rp巧|Xt):式中：Xt是t时刻系统的运行状态;P化i)为不确定事故Ei发生的概率;SsP化i)是发生事 故Ei时系统线路过负荷严重度;Rp化i|Xt)是发生事故Ei时系统过负荷风险值。3. 根据权利要求1所述的基于隶属度函数的模糊综合评价法的电网脆弱性评估方法， 其特征在于，所述动态安全性评估指标如下： 31) 小扰动最小阻尼比指标Cm: ξιη=η?η(ξι,ξ2,... ,ξη-1) 式中：ξι，ξ2，···，ξη-汾别为η-1个低频振荡模式的阻尼比，η为网络中发电机数； 32) 故障切除最小时间裕度指标IccT: 各线路的故障临界切除时间为ti，在给定各条线路的额定故障切除时间ti'时，计算拥 有nccT条输电线路的电网系统的故障切除最小时间裕度指标为：4. 根据权利要求3所述的基于隶属度函数的模糊综合评价法的电网脆弱性评估方法， 其特征在于，所述的电网各线路的故障临界切除时间为ti，计算方法为： 对于输电支路数为nccT双回路按照两条支路计算，不计发电机到母线的支路的电网络 其指标具体计算流程如下： Stepl:对电网络所有输电支路进行编号，从巧IjnccT，并令线路编号i = l，迭代次数j = 1，设定迭代误差ε〇 = 0.0001; Step2:在支路i的50%处发生Ξ相短路故障，并切除故障，给定故障切除时间为to(to值 参考正常的故障临界切除时间给定一个较大值），令皆'=? ; Step3:在故障切除时间为?/时，利用PSASP进行仿真计算电网络各节点的电压曲线； Step4:观察各个节点电压曲线的电压有无明显降低，是否发生电压失稳，若某个节点 发生电压失稳，令r，=i//2, j = j + l转到Step3;若电压无明显降低，未发生电压失稳， //"=化'+^尸)/2,判断的"-钓 <卸是否成立，若成立，令j = j+l，转到steps,若不成立，令 j = j+1，转到 St 巧3; Step5:记录支路i的故障临界切除时间为?,.=岭.，.若i <ncCT，令i = i+1，转到St巧2;若i = nccT，计算过程结束。5. 根据权利要求1所述的基于隶属度函数的模糊综合评价法的电网脆弱性评估方法， 其特征在于，所述拓扑结构稳定性指标如下： 51)网络的平均凝聚度指标 电网络的各节点凝聚度及平均凝聚度指标的计算公式为：式中：化为网络中节点的数目，m。为节点i的收缩网络中节点的数目，li为节点i收缩网 络中节点之间的平均最短路径，cUin,gh是用边的数目表示的网络中任意两节点g和h之间的 最短距离，V表示节点i收缩网络中所有节点组成的集合，αι，α2，···，αΝ为电网络中各节点的 疑聚度。6. 根据权利要求1所述的基于隶属度函数的模糊综合评价法的电网脆弱性评估方法， 其特征在于，所述可靠性的脆弱性指标如下： 61) 系统平均停电频率指标SAIFI:式中，Ν功节点i的用户数，Asus_功节点i的永久故障率，Nsus_功受永久故障影响的用户 数，化为电网全部用户数，化=SNi; 62) 系统供电可靠率指标RS:式中，N功节点i的用户数，ULD3d_功节点i的年平均停电时间（小时/年）。7. 根据权利要求1步骤2)所述的基于隶属度函数的模糊综合评价法计算电网的脆弱性 评分，其特征在于，脆弱性评分的计算步骤为： 71 )Ξ标度法求取指标权重 通过Ν个指标综合评价电力系统的脆弱性时，需要求取每个指标的权重;采用Ξ标度 法求取指标的权重向量/= (W1，W2,…，WN); 72) 确定指标集与评语集 m个待评价方案，N个评价指标;其中第i个方案的指标集为化={uii，Ui2，…，ιην}， m，其中uu(l《j《N)为第i个方案的第j个指标数值； 对指标有1种评语，评语集V={vl，V2，…，vl}，例如，l = 5，V={优，良，中，合格，差}; 73) 评价矩阵标准化 由评价指标集可得m个评价方案的评价矩阵U:其中， 评价方案的评价指标分为效益型和成本型两类，将U标准化后得到标准评价矩阵C:其中，效益型指标按上述公式第1式处理，成本型指标按上述公式第2式处理； 74) 构造隶属度函数模糊评价矩阵 对于m个待评价方案，第i个方案标准化后的评价指标集为Cl = kil，Ci2，…，Cin}，1《i《 m，其第j个指标在评价集V上的模糊子集可通过其对V的隶属函数来计算;采用等腰Ξ角形 隶属度函数、等腰梯形隶属度函数和高斯型隶属度函数，分别得到第i个方案的模糊评价矩 阵R;"、Ri(2)、Ri(3);其中ru(vk)为第i个方案的第j个指标相对于评语vk(l《k《l)的隶属度； 75) 求取综合评价模糊子集 第i个方案的综合评价集Bi为V上的模糊子集，通过下式计算： 技。R 二斯，心 式中，其中wc为综合权重集，Ri为第i个方案的模糊评价矩阵，算子。采用M( ·，+ )模型， 则有：76) 利用模糊综合评价结果对方案进行评分 给评价结果的每一个评语赋予一个分值而将其量化，可得评分集V =巧Λ，··，w; 方案i的评价结果向量 Zi = [zil ,Zi2 , , Zil] %=4? 式中，Zik表示方案i相对于评语Vk的得分;那么，方案i的综合得分文m个方案的综合评价矩阵Z= 口1，Z2，· · ·，Zm]T; 等腰Ξ角形隶属度函数、等腰梯形隶属度函数和高斯型隶属度函数的综合评价矩阵分 别关77)组合系数加权求取Ξ种函数的综合评价结果向量ZE 考虑到各评价向量之间的相似度的不同，相似度大的两个评价向量在计算综合评价向 量时应予W更多考虑，采用组合系数加权的方法进行计算； 基于等腰Ξ角形-等腰梯形-高斯隶属度函数的综合评价结果向量.ZE.计算公式如下： 式中：ZW、ZW、ZW表示；种不同隶属度函数下的评价结果向量，ru、ri2、ri3为组合系 数。8.根据权利要求7步骤71)所述的采用Ξ标度法求取指标权重向量/，其特征在于，迭代 求取指标权重的计算步骤为： 81) 利用Ξ标度法建立互补型模糊判断矩阵F=(fu)NXN:乂"、./7分别表示第i、j个指标的重要程度； 82) 将互补型模糊判断矩阵F= (fu)NXN转变为模糊一致型判断矩阵R= (ru)NXN:83) 根据和行归一化方法求解得到排序向量wW:84) 将模糊一致性判断矩阵R=(町)NX娘变成互反型矩阵E= (eu)NXN:85) 选取排序向量wW作为迭代初始值，确定误差ε; 86) 求取特征向量的初值wW的无穷范数I |wW| μ，置k=l:88)进行判断 若I Iw^iicxH μ|<Ε，所得到的向量就是目标的权重向量，迭代结束;否则 置k = k+l，转步骤7)，继续计算； 迭代结束后所得到的权重向量即为最终的权重向量：
【文档编号】G06Q50/06GK105976257SQ201510958505
【公开日】2016年9月28日
【申请日】2015年12月17日
【发明人】魏文辉, 林昌年, 夏继红, 马强, 金逸, 费益军, 杨选怀, 高峰, 赵云军, 张辉, 秦可源, 马群, 史立勤, 林春龙
【申请人】国家电网公司, 北京科东电力控制系统有限责任公司, 国网山东省电力公司, 江苏省电力公司