【精彩论文】电网发展诊断的改进模糊评价模型与评级
电网发展诊断的改进模糊评价模型与评级
艾欣1, 胡寰宇1, 任大鹏1, 彭冬2, 刘汇川3, 薛雅玮2, 张天琪2
(1. 华北电力大学 电气与电子工程学院,北京 102206; 2. 国网经济技术研究院有限公司,北京 102209; 3. 国网江苏省电力有限公司经济技术研究院,江苏 南京 210008)
引文信息
艾欣, 胡寰宇, 任大鹏, 等. 电网发展诊断的改进模糊评价模型与评级[J]. 中国电力, 2022, 55(5): 66-75, 165.
AI Xin, HU Huanyu, REN Dapeng, et al. Improved fuzzy evaluation model and assessment of power grid development diagnosis[J]. Electric Power, 2022, 55(5): 66-75, 165.
引言
鉴于此,本文延续传统综合评价的框架,对电网发展诊断的问题进行了研究。首先建立比较完备的动态评价指标体系,构建三层级体系结构以克服电网成长过程中造成的指标增减影响。其次,利用主成分法融合三级指标;然后通过分析二级指标集合,构建合适的隶属度函数,表达其关于各状态的趋势程度,并利用聚类分析方法获取其趋近类型;进而通过隐马尔可夫模型挖掘不同时刻下趋势类型背后的实际电网发展等级,对省级电网的发展进行动态综合评价。最后利用所提的改进模糊评价模型对电网实际时序数据进行计算和分析。
1 电网发展诊断动态评价指标体系构建
1.1 指标的选取
在国网经研院最新发布的《电网发展诊断分析指标体系(2018年版)》《电网发展诊断分析大纲》中,以电网发展协调、安全、质量、效率、经营与政策6个方面,共计百余项基础指标描述电网发展状态。该体系经过逐年完善,已具备详尽的指标构建模型且辐射面广泛,明确了电网发展与经济发展、环境影响间紧密的耦合关系,可适用于电网发展诊断的指标构建。本文聚焦于电网发展的诊断,拟对电网发展进行目标时序下的动态综合评价研究。而传统动态综合评价以静态评价方式为基础,其固定的评价体系往往忽略了电网成长不同时刻下存在的差异。因此,基于上述资料与分析,本文通过研究各省市电网发展的时空特性,并结合现有文献对电网发展诊断评价指标体系的研究成果[4-13,18-20],同时采纳专家评价意见与数据获取情况,综合考虑上述内容以构建三级指标体系,如图1所示。其中,三级指标的筛选[4]遵循一致性原则、独立性原则、可比性原则。由图1可知,指标层表现了三级指标数量随电网发展而产生的增长变化,反映了电网发展的时空差异性,包括最大数额等三级指标88项,在图1中以矩形的面积变化表示。关于电网评价的基础指标研究已较为充足,逐步规范了指标算式,为本文面向电网发展整体视角下的指标广泛选取提供了支持,现列举其中部分指标定义:以电网发展安全为例,“N−1通过率”为某一电压等级电网主要输变电元件(不考虑母线)满足N−1标准的比例;“同塔双回线N−2通过率”为某一电压等级同塔双回线中满足N−2的线路比例;“安全隐患”是指电网在N−2、N−1−1等特殊故障方式下可能发生特别重大事故、重大事故、较大事故、一般事故的隐患的数量。以电网发展质量为例,“变压器/线路停运率”为反映各电压等级变压器/线路的强迫停运统计比例;“自动化率”反映综合自动化变电站配置比例。控制层由电网发展的速度与规模、电网发展的安全与质量、电网发展的效率与效益、企业经营与政策4类二级指标构成,涵盖发电、输电、负荷的基本环节,同时考虑了经济、政策、环境的交互影响,表现了电网本身的发展潜力与多方协调性。目标层由预警、亚健康、健康3种诊断状态构成,各状态为考虑电网全时序下各二级指标在某评价时刻的综合发展水平体现,也是对综合评价结果的直观表示形式。预警状态表明电网明显处于较差的发展水平,电网各方面发展出现显著的不协调情况导致短板效应的出现。亚健康状态代表电网处于中庸的发展水平,各项二级指标值将集中于评价分数的中位数周围,值得引起调控规划人员的注意,以把握电网未来朝健康发展的趋势。健康状态反映了电网综合发展呈最优状态,电网各方面及整体协调性都处在最优水平。图1 三层动态指标体系结构
Fig.1 Architecture of three-level dynamic index system
1.2 三级指标标准化处理
电网发展诊断的综合评价所涉及的指标数量庞大,为消除多类型指标量纲不同造成的影响,需首先对三级指标集进行标准化处理[21]。随着电网各项技术的发展,其各静态时刻下三级指标数量不尽相同,应参与电网评级诊断工作的指标数将动态增长。进行标准化之前需对指标类型进行分析,采用适配的标准化公式处理。对于极大型指标,如发输配环节的电源装机增长率、售用电增长率及经营环节的人均输配电线路长度等,都属于指标值越大反映发展效应越好的指标,即式中:
对于综合线损率,环境类的碳排放、硫排放量及经营环节的单位输配电成本、资产运维费用等,都属于极小型指标,其物理意义同极大型恰好相反,即
而对于合格型指标,如重大事故、较大事故等安全隐患类指标,根据是否超出规定限额将其标准化更加合理。故本文引入0−1型处理,即
因为指标在各评价时间截面取值不同,所以利用
动态综合评价旨在量化评价目标的动态成长,但目前的研究实质为多次静态综合评价结果的加权集结,其任意时刻下评价体系不变。而本文研究的电网发展问题,将考虑电网发展造成的基础指标增减影响,这须打破传统评价体系的固定结构,相应影响静态指标的权重确定问题。
目前权重计算的方式可分为主观法、客观法、主客观结合法三大类,已有较成熟的研究基础。但动态指标体系的三级指标在不同时刻下的指标数量不一定相同,再通过对不同时刻分别进行权重计算的方式将出现评价计算时指标维度不同的问题,同时导致对整体评价方向的影响。因此,引入主成分分析法(principal component analysis,PCA)将不同时间截面上数量不同的三级指标进行分步融合,使其转换为数量统一的二级指标,以将不同时刻的三级指标集结问题置于同维度解决,同时尽可能考虑到电网发展带来的指标变化。
PCA是多元统计分析中经典方法之一[23]。该方式能在保留基础指标所反映关键信息的同时规避各时刻下指标集结对整体评价结果的影响。各主成分融合所取权重直接为对应主成分的方差贡献率。
2 改进模糊综合评价模型
2.1 模糊综合评价模型
模糊综合评价方法[24]是本文开展电网发展诊断研究工作的重要手段,其关键是通过模糊数学的隶属度理论来描述评价目标的多种受限因素。根据本文所建立的三级指标体系,融合后的电网发展二级指标将参与本节的模型计算。其具体数学思路简述如下。(1)利用某种模糊函数映射将含有Q个二级指标的评价目标2.2 隶属度矩阵特征分析
隶属度函数兼并了模糊概念的特性及不确定性的量化,可作为电网评价指标与评级计算间实现相应数学运算和处理的桥梁。其函数的选择、建立与评价目标模糊度的描述准确性息息相关。本文针对电网发展动态进行诊断评级,希望各指标值在不同等级间隶属程度明显,故选用岭型函数作为各等级间隶属度函数,如图2所示。
图2 岭型隶属度函数模糊分布
Fig.2 The fuzzy distribution of ridge function
式中:Ci为各观测值类型中心。依据上述工作将获得电网发展在任时刻下的隶属度矩阵特征,可由隐马尔可夫模型进一步挖掘各时刻下不同指标所反映的整体状态。
2.3 基于隐马尔可夫的评价模型隐马尔可夫模型可以推导出非平稳时间序列中存在的关系,并利用概率统计方式进行预测[25],目前已经广泛应用于电网如故障检测,电力系统预测[26-27]等问题中。隐马尔可夫问题可分为学习、评估、预测3类[28],在模型构建环节常采用前向、后向,Baum-Welch等算法计算模型参数。其中,Baum-Welch算法属于期望最大化算法,其目标是在已知观测序列情况下求得模型参数,以使得观测序列的条件概率最大[29]。隐马尔可夫理论能深入挖掘指标背后反映的评价目标整体状态,其模型结构也同综合评价结构类似。因此,若利用概率模型推演指标评价的过程,取代依靠权重集结的数理统计方式,能更好保留传统评价过程中忽略的不确定性,使评价结果更符合实际。
首先简要回顾一下隐马尔可夫模型基础参数及概念,搭建电网发展诊断评级模型为
式中:N表示电网发展等级种类;在模型中一般由 θt 表示各等级,研究电网发展等级全时序列可表示为 (θ1,θ2,⋯,θT) ,其中T代表时序评级长度,时刻t依据时序数据的时间单位决定。隐马尔可夫计算过程的每一步骤被称为状态,由于本文对电网发展动态进行评价,故马尔可夫过程的每一步骤同时刻相对应,由St表示当前时刻t的马尔可夫状态。隐马尔可夫模型观测值由观测序列 {O1,O2,⋯,OT} 表示,代表每时刻电网处于某发展等级下的观测结果。其中,A为状态转移概率矩阵,各元素分别代表 P(St+Δt=θj|St=θt),1⩽i,j⩽N ,对应电网不同时刻的等级变化趋势。B为观测概率矩阵,各元素分别对应电网处在不同等级下输出各类型观测值的概率。πN×1={πi} 代表初始状态矩阵,各元素分别代表 πi=P(S1=θt),1⩽i⩽N ,表示电网发展综合评价起始时刻电网所处各等级的估计概率。本文利用Baum-Welch算法迭代计算估计模型基础参数,其计算过程由式(6)和式(7)表示,以构建 λ(A,B,π) 隐马尔可夫模型。基于隐马尔可夫模型的综合评级过程如下。
(1)初始化模型各参数 A0,B0,π0→λ0 。
(2)利用式(6)(7)迭代计算并更新隐马尔可夫模型参数。其中 αt(i) 是总观测序列Ot中的前向序列部分 O={O1,O2,⋯,Ot} 时刻, βt(i) 是余下的后向序列部分 O={Ot+1,Ot+2,⋯,OT} ,状态 θt 处在t时刻。
(3)迭代计算 P(O|S) ,当满足停止条件时输出 Aiter,Biter,πiter ,完成隐马尔可夫模型构建。
(4)利用动态规划法求取二级指标隶属度矩阵所反映的电网发展诊断等级序列。
基于上述思路,本文提出基于电网发展诊断的改进模糊评价模型,其诊断工作的流程如图3所示。
图3 电网发展诊断整体框架
Fig.3 The overall framework for power grid development diagnosis
3 算例分析
本文以沿海省份的实际历史数据为研究对象,数据时间为2010—2017年。采用所提改进模糊评价模型对其中某省的电网发展进行动态综合评级分析。
3.1 指标整理
原始数据主要来源于国网经研院。参考本文提出的指标体系进行计算,三级指标分别隶属于电网发展速度与规模(SS)、电网发展安全与质量(SQ)、电网发展效率与效益(EB)、企业经营与政策(MP)4个二级指标下。实际参与算例的指标情况由表1统计,由各年的总指标统计值可见其逐年增长态势。通过前述分析,首先对原始指标标准化处理,以消除指标间量纲的不同,且将多元指标进行类型上的统一;再利用多次主成分法将三级指标融为对应各时刻下的二级指标值。通过三级指标间融合,可将数据包含的评价信息最大化折射至二级指标中。同时使评价体系由动态转为静态。
表1 指标数量
Table 1 The value of each index
图4描述了该省电网各二级指标值的时序变化,由任一曲线特征可以得到对应方面的时序发展特性,获得单时间截面下的静态评价结果。如该省电网除安全与质量方面,其余二级指标都于2016年达到全年发展的最大值,反映了该省于2015—2016年发展优异。而电网整体发展的状态由这4个方面耦合体现,需进一步分析二级指标集与评价目标间的映射关系。
图4 该省电网二级指标时序图
Fig.4 Sequence diagram of the provincial secondary indexes
3.2 隶属度特征计算
通过前文讨论及指标分析研究,确定岭型分布隶属度函数来描述二级指标与各等级的模糊关系。本文预设电网实际状态数量为N=3,分别代表电网发展等级处于“健康”“亚健康”“预警”。表2以该省2015年为例展示了电网发展二级指标隶属度计算结果。不难看出,不同二级指标隶属度向量表示对各状态不同的趋近程度,包含了状态间转移的趋势影响。其中,利用最大隶属度准则判断[29],可知效率与效益比速度与规模、经营与政策两方面趋近于更好状态。针对电网发展的整体状态,当前时间截面的实际状态由4个二级指标构成的隶属度矩阵耦合体现,故无法直接判别比较。
表2 2015年该省电网发展二级指标隶属度
Table 2 The secondary index membership of the provincial power grid development in 2015
考虑电网发展本身具有规律的特性,故利用聚类分析法挖掘全时序下电网处在每时刻截面的隶属度矩阵特性,以获得电网整体的实际状态趋势,将聚类结果以观测值Ⅰ、观测值Ⅱ、观测值Ⅲ、观测值Ⅳ表示。各观测值类型代表与实际状态的趋近类别。其中,属于观测值Ⅳ的二级指标隶属度矩阵具有趋近更好状态等级的最大概率,趋近于最差等级的最小概率,观测值Ⅲ、观测值Ⅱ在趋近程度上依次减小,观测值Ⅰ则具有趋近于最差状态的最大概率,同观测值Ⅳ相反。算例结果如表3所示。
表3 该省电网二级指标隶属度特征模式
Table 3 The membership characteristics of the provincial power grid secondary indexes
通过上述聚类分析,已将二级指标对各状态等级的趋近程度完成最大程度挖掘并表现,实现由传统的评价分值排序至评级划分的转变。各观测值类型包含了与不同状态等级的映射关系,也具有序列特性,可构成隐马尔可夫模型的观测序列。
3.3 隐马尔可夫模型计算
通过前述小节的隶属度矩阵特征划分,一定程度上揭示了隐藏在指标后的实际电网发展状态与二级指标集之间的关系,最后利用隐马尔可夫模型挖掘全时序下二级指标集所处的发展等级。将马尔可夫模型隐藏状态确定为 θ1、θ2、θ3 ,代表电网发展所处的预警、亚健康、健康等级。利用Baum-Welch算法,获得隐马尔可夫模型参数,如图5所示。经过分析可以确定隐马尔可夫模型的隐藏层所代表状态的实际含义。构建隐马尔可夫模型的状态转移概率和观测概率如图5所示。图5a)状态之间的变化具有以下特点。(1)该省电网发展状态由预警转移至亚健康有较高的转移概率,达到46.65%。(2)当该省电网处于亚健康时,未来维持该状态的概率为40.81%;而未来跌回至预警状态的概率略大于转移至健康的概率。(3)当该省电网实际处于健康时,具有较高的概率维持当前状态;未来转移至亚健康的概率是转移至预警等级的2倍。图5 转移概率矩阵与观测概率矩阵
Fig.5 The hidden state transition probabilities and observation probabilities
各状态之间的转移概率符合实际意义:(1)由亚健康状态往高或低转移的概率均衡,且由预警状态易恢复到亚健康状态,当处于健康状态时易保持当前状态,且下滑至亚健康状态的概率略低于维持当前状态概率。各状态变化转移明显具有平滑性、连贯性、合理性。(2)预警往亚健康及健康状态的高转移概率体现了电网公司每年对电网规划诊断的合理把控。(3)各状态间跨等级转移的概率较低,维持当前状态的概率较高,体现了电网发展的惯性特征。以上的分析表明,该参数有效地描述了电网发展特性,可以用来开展电网发展诊断评级工作。图5b)直观地展示了各状态与观测值之间的关系:(1)当电网处于预警状态时,输出观测值Ⅳ的概率接近于零。这表明电网处于预警发展状态时的特征明显易辨识。(2)当电网处于亚健康时易输出观测值Ⅱ,其概率为56.45%;而电网处于健康状态时,输出观测值Ⅳ与观测值Ⅲ的概率分别为12.67%与43.96%。不难看出观测值Ⅱ、观测值Ⅲ处于隶属度矩阵特征的中间类型,同电网各级状态都具有一定的映射关系,可随状态的转移而反映在输出观测值概率的变化上。随着电网状态由亚健康向健康状态的好转,可见输出相同类型观测值Ⅱ、观测值Ⅲ的概率发生了翻转,符合实际。(3)当电网处于健康状态时具有输出观测值Ⅳ较高的概率23.86%,同时该值为输出矩阵内各状态下输出观测值Ⅳ的最高概率,具有唯一性。观测概率矩阵值符合实际意义:当电网发展状态处于预警或健康的两极时分别具有与之相对应的高概率出现的观测值结果,且当电网发展处于中间状态时,概率分布反映了观测结果与状态转移间的相关性。输出概率矩阵有效地揭示了评价指标结果与电网实际状态划分之间的概率分布情况。
通过上述分析研究,构建的隐马尔可夫模型可以反映电网发展的实际情况、变化特征,可以利用其对电网发展进行诊断评级。根据上述隐马尔可夫模型,利用动态规划算法Baum-Welch可以求取电网发展最大概率状态序列,对电网实际发展进行诊断评级。同时,在本文动态评价体系的基础上采用G1-熵权法[30]的综合评价方法进行最后的评价计算,其结果如表4所示。不难看出,常规的综合评价优劣排序结果同本文判定的电网发展等级相似。但在2010年与2014年的结果出现了差异,综合评价分数较高的年份却处于预警等级,通过分析可知2014年的效率与效益指标处于往年最低值。综上所述,本文所提模型对电网发展诊断提供了更直观、可视化的评级结果,且能够兼顾指标优缺点对总体结果造成的影响,使评价结果更科学。
表4 该省电网发展诊断动态评级结果
Table 4 The dynamic rating results of the provincial grid development diagnosis
该省电网逐年发展趋于平稳。根据二级指标变化趋势对等级的动态变化展开分析,前5年电网发展各方面不均衡,指标差异性大;尤其体现在2011年电网发展速度与规模、安全与质量环节处于最低值,短板效应导致电网发展处于预警期。而近3年的健康发展等级明显体现在二级指标集体较好发展,得益于公司的优化管理,同时近年以来的电网诊断评估工作也起到了改善性作用。
4 结论
加速推进电网新基建发展的背景下,为了实现电网发展诊断涵盖全面、全时序,并能够合理评价发展状态,改进评价体系及方法,提出发展等级指标是本文的主要研究目标。本文的主要结论有:(1)电网的飞速发展使其结构、运行方式等都产生变化,同时也影响电网综合评价指标的确立和选取。在电网发展诊断工作中,电网发展的时空差异性反映在评价指标的增减上,利用动态评价指标体系可考虑指标增减对评价结果的影响,相比于固定指标结构考虑到更多信息,使评价结果更切合实际情况。本文构建的动态指标体系仍可适应未来评价研究。(2)目前综合评价方法的研究聚焦于指标间的权重赋值问题,优化结果局限于评价对象的评分排序。利用隐马尔可夫模型对综合评价结果进行判定,用概率计算的模式取代传统综合评价中依据权重以线性集结各指标的方式,为综合评价提供新的思路。同时,相比于传统评价下分值高低排序的方式也使电网状态层次反映更加直观。(3)本文的算例分析以某省历史数据为例,所得出的8年来实际发展诊断的等级结果符合实际,有助于探寻其发展态势变化规律并预测其未来变化趋势。所提的电网发展分级指标有利于直观表示评价结果,同时可用作新的指标参数辅助相关电网数字信息化技术的顺利实现。
本文所提出的方法同样可以适用于时间步长更小的分析,随着电网新基建发展进程中建设新兴数字智能产业,未来将产生庞大数据、信息流,可以实现对电网发展状态更加精确的评估和评级,实现对电网发展的态势感知。
(责任编辑 张重实)
作者介绍
艾欣(1964—),男,教授,博士生导师,从事新能源电力系统及微网研究,E-mail:aixin@ncepu.edu.cn;★
胡寰宇(1995—),男,通信作者,博士研究生,从事电力系统分析与控制研究,E-mail:huanyuncepu@163.com;
★
任大鹏(1987—),男,博士研究生,从事电力系统分析与控制研究,E-mail:rocperson@126.com.
往期回顾
◀审核:方彤
根据国家版权局最新规定,纸媒、网站、微博、微信公众号转载、摘编《中国电力》编辑部的作品,转载时要包含本微信号名称、二维码等关键信息,在文首注明《中国电力》原创。个人请按本微信原文转发、分享。欢迎大家转载分享。