CN115271070A - 基于改进Relief-F算法的电压稳定裕度评估方法、设备及介质 - Google Patents
基于改进Relief-F算法的电压稳定裕度评估方法、设备及介质 Download PDFInfo
- Publication number
- CN115271070A CN115271070A CN202210814511.1A CN202210814511A CN115271070A CN 115271070 A CN115271070 A CN 115271070A CN 202210814511 A CN202210814511 A CN 202210814511A CN 115271070 A CN115271070 A CN 115271070A
- Authority
- CN
- China
- Prior art keywords
- voltage stability
- algorithm
- stability margin
- decision tree
- relief
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 76
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 66
- 238000012549 training Methods 0.000 claims abstract description 82
- 238000003066 decision tree Methods 0.000 claims abstract description 64
- 238000000034 method Methods 0.000 claims abstract description 34
- 239000013598 vector Substances 0.000 claims abstract description 13
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 230000007547 defect Effects 0.000 claims abstract description 9
- 230000008569 process Effects 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000005192 partition Methods 0.000 claims description 7
- 230000009467 reduction Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 238000013211 curve analysis Methods 0.000 claims description 2
- 230000003068 static effect Effects 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 13
- 238000004364 calculation method Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000013210 evaluation model Methods 0.000 description 2
- 238000000556 factor analysis Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010187 selection method Methods 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 108010074864 Factor XI Proteins 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000001149 cognitive effect Effects 0.000 description 1
- 230000001687 destabilization Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000013097 stability assessment Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Power Engineering (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本申请公开了一种基于改进Relief‑F算法的电压稳定裕度评估方法、设备及介质,以解决现有方法在评估***电压稳定裕度状态时易发生误判,进而造成***电压失稳的技术问题。该方法包括:根据预训练的决策树模型的路径结构提取电压稳定裕度评估规则;其中,用于构建所述决策树模型的训练样本集,通过引入了惩罚因子的Relief‑F算法进行划分,惩罚因子用于识别异类标签以弥补Relief‑F算法固有缺陷,增强场景适应性;将向量测量单元采集的实测数据与由决策树模型提取的所述电压稳定裕度评估规则进行匹配,得出***当前电压稳定状态。本申请实现了对不同类别标签差异的辨识,切实提高了Relief‑F算法在电力***静态电压稳定评估场景中的适用性,提升了电压稳定裕度评估的准确性。
Description
技术领域
本申请涉及电力***电压稳定技术领域,尤其涉及一种基于改进Relief-F算法的电压稳定裕度评估方法、设备及介质。
背景技术
电力的安全、可靠、稳定供应是国家及社会发展的基础,我国现阶段采取的电网调度运行模式通过一系列运行规则,将***限制在一定的安全稳定运行状态中。运行规则在保证我国电网安全稳定,防止大规模停电事故中发挥了决定性作用。然而,随着我国电网规模的不断扩大及新能源高比例并网,现有基于运行规则的运行模式弊端突显:首先,传统运行规则基于离线数据分析并依靠半人工方式总结、提炼,其可靠性和准确性取决于运行方式人员的个人经验,存在一定主观性,同时受制于人的精力及计算能力,所提取的运行规则对在线运行方式变化的适应性不强;其次,随着网架规模的不断扩大和以风电光伏为代表的新能源大规模接入,电网运行特性日益复杂,离线数据分析计算规模呈指数级增长。因此,现有依靠人工提炼***运行规则的工作模式已难以满足新环境下电网安全运行的需要,极易导致***状态的误判、漏判,降低电网运行的安全性、稳定性以及经济性。
现有基于人工神经网络、支持向量机等人工智能技术的电压稳定评估方法虽然能获得可靠的评估结果,但受制于其技术自身的局限性,无法展示评估的机理过程,不能为运行人员提供具有可解释性的运行经验,实用性有待进一步提升。
发明内容
本申请实施例提供了一种基于改进Relief-F算法的电压稳定裕度评估方法、设备及介质,用以解决现有的电压稳定裕度评估方法可解释性不强,并容易发生***电压失稳状态漏判的技术问题。
一方面,本申请实施例提供了一种基于改进Relief-F算法的电压稳定裕度评估方法,所述方法包括:根据预训练的决策树模型的路径提取电压稳定裕度评估规则;其中,所述决策树模型的训练样本集,利用引入了惩罚因子的Relief-F算法划分,所述惩罚因子用于识别异类标签以增强场景适应性;将向量测量单元采集的实测数据与提取出来的所述电压稳定裕度评估规则进行匹配,得出***当前电压稳定状态。
在本申请的一种实现方式中,所述决策树模型的训练过程为:基于母线参与因子、支路参与因子、发电机参与因子以及各自的第一预设阈值,筛选初始特征变量;在Relief-F算法中引入惩罚因子,并基于所述Relief-F算法计算所述初始特征变量的权重值;基于所述权重值确定关键特征变量;将训练样本集输入到决策树根节点,在决策树各级节点逐个计算所述关键特征变量的信息增益率;选取信息增益率最大的所述关键特征变量作为对应节点最终划分属性;利用所述最终划分属性对所述训练样本集进行逐级划分,以得到训练样本子集;将所述训练样本子集传输到所述决策树模型的次级节点,直至节点中训练样本子集的初始熵为零,则决策树模型训练完成。
在本申请的一种实现方式中,训练样本集的获取,具体为:在***初始运行点进行P-V曲线分析,并确定所述P-V曲线最大功率传输点;基于所述最大功率传输点,计算电压稳定裕度;基于所述电压稳定裕度,划分电压稳定状态区间;其中,所述电压稳定状态区间包括:正常状态区间、预警状态区间、紧急状态区间;在每个电压稳定状态区间内随机抽取等量运行点作为训练样本,并基于全部所述训练样本构建训练样本集。
在本申请的一种实现方式中,所述基于所述各母线参与因子、支路参与因子、发电机参与因子以及各自的第一预设阈值,确定初始特征变量,具体为:基于***运行状态变量构建潮流方程;其中,所述***运行状态变量包括:母线电压的幅值和相角、发电机无功功率以及支路流过的无功功率;对所述潮流方程依次进行线性化、矩阵降阶以及特征值分析,计算得到对应最小非零特征值;基于所述最小非零特征值,计算各条母线的参与因子、各条支路的参与因子以及各台发电机的参与因子数值;对各参与因子分别降序排列,设定各自的第一预设阈值以筛选所述初始特征变量。
在本申请的一种实现方式中,计算各母线参与因子,具体为:对所述潮流方程进行线性化、矩阵降阶以及特征值计算,以得到所述潮流方程降阶雅克比矩阵对应的特征值;计算最小特征值的左特征向量以及右特征向量;将所述左特征向量以及所述右特征向量进行相乘运算,得到母线参与因子。
在本申请的一种实现方式中,计算支路参与因子以及发电机参与因子,具体为:确定待测支路的无功损耗以及所有支路中无功损耗的最大值;将所述待测支路的无功损耗与所有支路中无功损耗的最大值做相除运算,以得到所述待测支路的参与因子;确定待测发电机输出的无功功率以及所有发电机输出的无功功率的最大值;将所述待测发电机输出的无功功率与所述所有发电机输出的无功功率的最大值做相除运算,以得到所述待测发电机的参与因子。
在本申请的一种实现方式中,所述根据所述关键特征变量对训练样本集进行划分,具体为:利用改进的Relief-F算法计算各初始特征变量权重值,并对所述权重值降序排列;设定第二预设阈值,确定相应关键特征变量;基于所述关键特征变量对所述训练样本集进行划分。
在本申请的一种实现方式中,所述将训练样本集输入决策树模型进行训练,具体为:将训练样本集输入决策树模型的根节点,并计算所述训练样本集的初始熵;计算所有所述关键特征变量划分所述训练样本集得到的属性熵;基于所述初始熵以及所述属性熵,依次计算信息增益以及信息增益率;根据所述信息增益率确定最终划分属性;利用所述最终划分属性对所述训练样本集进行划分,以得到训练样本子集;将所述训练样本子集传输到所述决策树模型的次级节点;重复上述过程,直至所述训练样本子集的初始熵为零,则决策树模型训练完成。
另一方面,本申请实施例还提供了一种基于改进Relief-F算法的电压稳定裕度评估设备,其特征在于,所述设备包括:至少一个处理器;以及,与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够:根据预训练的决策树模型的路径结构提取电压稳定裕度评估规则;其中,用于构建所述决策树模型的训练样本集,通过引入了惩罚因子的Relief-F算法划分,所述惩罚因子用于识别异类标签以弥补Relief-F算法固有缺陷,增强场景适应性;将向量测量单元采集的实测数据与由决策树模型提取的所述电压稳定裕度评估规则进行匹配,得出***当前电压稳定状态。
此外,本申请实施例还提供一种基于改进Relief-F算法的电压稳定裕度评估的非易失性计算机存储介质,存储有计算机可执行指令,其特征在于,所述计算机可执行指令设置为:根据预训练的决策树模型的路径结构提取电压稳定裕度评估规则;其中,用于构建所述决策树模型的训练样本集,通过引入了惩罚因子的Relief-F算法划分,所述惩罚因子用于识别异类标签以弥补Relief-F算法固有缺陷,增强场景适应性;将向量测量单元采集的实测数据与由决策树模型提取的所述电压稳定裕度评估规则进行匹配,得出***当前电压稳定状态。
本申请实施例提供的基于改进Relief-F算法的电压稳定裕度评估评估方法,通过在Relief-F算法中引入惩罚因子,有效的解决了其无法辨识不同标签差异的固有问题,避免了***电压稳定状态判断中易发生漏警、误警的情况,切实提高了Relief-F算法的场景适用性。融合参与因子分析和改进Relief-F算法的特征变量选取方法,从特征变量的物理意义和分类能力两个角度对其进行了优选,既可以剔除冗余特征变量,有效提高电压稳定裕度的评估速度,还可以去除负效力特征变量,切实提升评估准确性。通过决策树模型来评估***电压稳定状态,不仅能准确获知***当前状态,并且可以为***运行人员展示***状态评估的整个流程及机理,为其提供辅助决策依据,以便在***状态失稳时针对越限状态变量采取相应控制措施避免事故进一步扩大。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1为本申请实施例提供的基于改进Relief-F算法的电压稳定裕度评估评估方法流程图;
图2为本申请实施例提供的决策树模型的训练过程图;
图3为本申请实施例提供的基于改进Relief-F算法的电压稳定裕度评估评估设备示意图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,下面将结合本申请具体实施例及相应的附图对本申请技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
我国现行的电网调度运行模式主要依靠已有的安全运行规则,通过将各运行状态量与既定安全运行规则进行校核以判断***当前运行状态,并采取相应控制措施以维持电力***的安全稳定。但随着“碳达峰、碳中和”目标的提出及新型电力***构建的要求,现有电网调度运行模式已难以满足新形势下电网运行的需求。一方面,现有***安全运行规则主要依靠半人工模式对电网历史运行数据进行梳理、归纳及总结,不同运行人员对相同问题的认知深度不同,个人经验导致所形成的规则存在差异性;另一方面,随着电网规模的快速扩张及新能源渗透率日渐提高,电网运行特性愈加复杂,受制于人的精力,***运行人员无法及时、有效处理和分析***复杂运行方式下产生的海量数据及知识。
人工智能技术的出现为解决电力***稳定评估问题提供了新方向,但现有以人工神经网络为代表的求解模型绝大多数均为黑箱模型,即仅能通过模型输出获得***当前稳定状态结果,但无法提取支撑该结果的稳定知识,供***运行人员对***状态进行校验,因此工程实际中对此类方法认可度不高,实用性有待提升。同时,现有稳定评估技术的准确性有待进一步优化,一旦发生状态漏判或误判等情况,如果不采取紧急控制措施极易造成***电压失稳,进而导致大规模停电事故的发生。
本申请实施例提供了基于改进Relief-F算法的电压稳定裕度评估方法、设备及介质,用以解决现有的电压稳定裕度评估方法容易发生误判的技术问题。
下面通过附图对本申请实施例提出的技术方案进行详细的说明。
图1为本申请实施例提供的基于改进Relief-F算法的电压稳定裕度评估评估方法流程图。如图1所示,该方法主要包括以下步骤:
步骤101、根据预训练的决策树模型的路径结构提取电压稳定裕度评估规则。
本申请实施例,采用预先训练的决策树模型来实现对电压稳定裕度的评估,是由于决策树模型具有可以兼顾评估结果准确性以及可解释性的优点。因此,事先要对决策树模型进行预训练,然后,根据训练完成的决策树模型的路径提取电压稳定裕度评估规则。
步骤102、将向量测量单元采集的实测数据与由决策树模型提取出来的所述电压稳定裕度评估规则进行匹配,得出***当前电压稳定状态。
本申请实施例中,要得出电压稳定裕度,需要根据由决策树模型提取出来的电压稳定裕度评估规则来确定,因此需要预先根据决策树模型进行电压稳定裕度规则的提取。
本申请实施例中,决策树模型的训练过程如图2所示,图2为本申请实施例提供的决策树模型的训练过程图,主要包括如下步骤:
步骤201、创建训练样本集。
本申请实施例中,在构建模型之前,首先需要获取相当数量的训练样本数据。本申请中,在电力***初始运行点上逐步增加负荷,即可得到P-V曲线,然后将曲线的拐点定义为***的最大功率传输点,即在此点增加负荷会导致***逐步逼近失稳状态,因此可以把该拐点看作是***的电压稳定极限水平。考虑到实际***运行中需要保证运行状态远离稳定极限,本申请实施例中将***实时运行点到最大功率传输点的距离,定义为***的电压稳定裕度水平,以量化***当前运行状态距稳定极限的程度。电压稳定裕度的计算公式如下:
其中,Pmax为***的最大功率传输点,Pi为***的第i个运行点的功率,Mi为运行点Pi对应的电压稳定裕度。
为保证***运行人员在***即将失稳的状态下有充足的时间采取控制措施保证***稳定性,本申请中根据电压稳定裕度计算结果,将***的电压稳定状态划分为三类区间,分别为正常状态区间(5%<Mi≤8%)、预警状态区间(1%<Mi≤5%)、紧急状态区间(0%<Mi≤1%)。
***在正常状态区间运行时,***电压状态稳定,无需采取控制措施。***运行于预警状态区间时,***有电压失稳的倾向和可能,***运行人员此时应密切关注***运行状态,并采取一定的控制措施,避免***状态进一步恶化。而***运行于紧急状态区间时,则表示此时***已逼近失稳状态,如不采取紧急措施将会导致电压失稳,进而发生停电事故。
本申请实施例中,兼顾静态电压稳定裕度在线评估的速度和准确度,分别在每个状态区间取等量且充裕的***运行点,并将抽取的全部***运行点创建训练样本集。
步骤202、基于母线参与因子、支路参与因子、发电机参与因子以及各自的第一预设阈值,筛选初始特征变量。
在本申请实施例中,由于影响电力***电压稳定的主要原因是无功功率支撑不足,因此选取与无功功率密切相关的发电机、母线、支路等***设备元件的相应状态变量作为原始特征变量,即各条母线的电压幅值和相角,各台发电机所能提供的无功功率以及各条支路损耗的无功功率。在本申请实施例中,各个原始特征变量的确定,需要经过相应的计算操作来完成,根据***运行状态变量构建潮流方程,***运行状态变量包括:母线电压的幅值和相角、发电机无功功率以及支路流过的无功功率等。
进一步地,对所述潮流方程进行线性化、矩阵降阶、特征值计算等操作,可以得到潮流方程雅克比矩阵对应的一系列特征值λi,λi可以表示为:
其中,ΔV电压幅值增量,ΔQ为无功功率增量。可以看出当λi很小时,无功功率的微小变化都会导致***电压的剧烈波动。因此λi越小表示***越接近电压失稳状态。基于此在λi中选取非零最小值进行母线参与因子计算,即量化各母线对***电压稳定的贡献度,也就是说通过特征值λi可以得出特征值λi的左特征向量ξi以及右特征向量ηi。
将前述左特征向量ξi以及右特征向量ηi进行相乘操作,即可得到母线参与因子,公式如下:
Pki=ξkiηki
其中,Pki为λi对应的***状态下第k条母线的参与因子,ξki和ηki分别为λi的左、右特征向量的第k个元素。
各条支路对***电压稳定的影响由其无功功率损耗确定,计算公式如下:
其中,Pji为λi对应的***状态下第j条支路的参与因子,ΔQloss.j为支路j的无功损耗,ΔQloss.max.i为所有支路中无功功率损耗的最大值。
各台发电机对***电压稳定的影响由其输出的无功功率确定,计算公式如下:
其中,ΔQg.i为发电机i提供的无功功率;ΔQgen.max.i为所有发电机的无功功率最大值。
在本申请实施例中,计算出各个设备的参与因子之后,分别对其进行降序排列,然后设定第一预设阈值对母线电压的幅值和相角、支路无功功率以及发电机无功功率等各状态变量进行筛选,并作为模型评估的初始特征变量。
步骤203、基于改进的Relief-F算法计算初始特征变量的权重值,并确定关键特征变量。
本申请实施例中,考虑到决策树模型一旦选取到冗余特征变量或负效力特征变量会直接影响静态电压稳定评估结果的准确性和速度,因此利用特征选择算法进一步降维。同时考虑到现有特征选择算法无法衡量异类标签的差异,将其直接应用于电压稳定评估场景中易引发漏警或误警情况,场景适应性不高,因此本发明专利在Relief-F算法中引入惩罚因子,提出了一种改进的Relief-F算法:
其中,R为任一训练样本,Hj为与R同属一类的样本,Mj和Nj为与R异类的样本,ωA为任一初始特征变量A的权重值,初始值为零,P(class(*))表示样本集中属于类别class(Ri)的样本的比例,k1和k2分别为异类标签的惩罚因子,n为迭代计算的次数,m表示选取m个距离最近的样本,diffA(*)为样本距离函数:
其中,vA(I1)和vA(I2)分别表示样本I1和I2在特征变量A下的取值;vA(Ii)max和vA(Ii)min则分别表示所有样本在特征变量A下的最大值和最小值。
本申请实施例中,惩罚因子的取值取决于样本R类别与其异类样本M、N类别的相对关系。当异类样本相对于样本R为漏警情况时,样本的误判会导致***电压进一步失稳,因此为避免误警情况的发生应增加公式中该异类样本的迭代比重,即将惩罚因子设定为正值,进而使具备区分漏警能力的特征变量的权重值增大;而当异类样本相对于样本R为误警情况时,样本的误判并不会对***电压稳定造成严重后果,因此无需设定惩罚因子以增加特征变量权重值在迭代计算中的权重。惩罚因子具体设置情况见表1:
表1改进Relief-F算法的惩罚因子设置
进一步地,利用改进的Relief-F算法公式计算所有初始特征变量的权重值,然后对各个初始特征变量的权重值进行降序排列,通过设定阈值选取关键特征变量以构建静态电压稳定在线评估模型。
步骤204、利用关键特征变量在决策树的根节点和每个子节点对训练样本集进行逐级划分,直到决策树到达叶节点。
本申请实施例中,为了兼顾***电压稳定状态评估结果的准确性和可解释性,采用决策树模型来构建电压稳定裕度评估模型。
本申请实施例中,所述训练样本集存储于决策树模型的根节点中,然后选取任一所述关键特征变量对其进行划分,所得训练样本子集分别存储于下一级两个子节点中,重复此过程直到某个节点中全部训练样本同属一类,并将该节点作为叶节点。从根节点到任意一个叶节点的路径,清楚地展现了含有不同类别的训练样本的集合逐步细化为多个仅包含同类训练样本的子集的全流程。
进一步地,首先将完整的训练样本集输入到根节点中,然后计算初始熵,计算过程如下:
其中,S为由4.2节得到的训练样本集,m=3对应***状态个数,pi为训练样本集S中属于类别i(i=1,2,3)的训练样本占训练样本总数的比例。
在此基础上,计算利用不同关键特征变量作为划分属性对训练样本集进行划分后的属性熵:
其中,A是任一前述关键特征变量,SL和SR是依据关键特征变量A将训练样本集划分后的两个训练样本子集,|SL|、|SR|以及|S|表示对应训练样本子集中样本的个数。
本申请实施例中,将初始熵与属性熵做差即可量化各关键特征变量降低训练样本集混乱程度的能力,并将其定义为信息增益,信息增益的计算过程如下:
Gain(A)=Entropy(S)-EntropyA(S)
然后,引入***信息,即可量化各关键特征变量的增益率,其中***信息的公式如下:
式中,k为划分属性能将训练样本划分成的类个数,Si是由这个划分属性划分后的训练样本数。
增益率的计算过程为:
分别计算所有关键特征变量的增益率,然后将最大增益率对应的关键特征变量作为根节点处训练样本的最终划分属性,然后利用最终划分属性对训练样本集进行划分,以得到最终训练样本子集。
进一步地,将最终训练样本子集分别传递给次级节点,再次计算各关键特征变量的增益率,并选取各节点处的划分属性,重复此过程直到节点中训练样本同属一类,即初始熵为0,并将此节点作为叶节点,当全部路径都达到即可认为决策树模型构建完成。
本申请实施例中,决策树模型构建完成之后,可以根据决策树的路径提取电压稳定评估规则。***运行人员可以将向量测量单元采集的实测数据与所提取的规则进行匹配,以评估***处于何种电压稳定状态。
本申请实施例提供的一种基于改进Relief-F算法的电压稳定裕度评估评估方法,通过在Relief-F算法中引入惩罚因子,然后通过改进后的Relief-F算法确定关键特征变量,有效的解决了Relief-F算法无法辨识不同标签差异的固有问题,避免了漏警、误警情况的发生,切实提高了Relief-F算法在静态电压稳定评估场景中的适用性。融合参与因子分析和改进Relief-F算法的关键特征变量选取方法,从变量的物理意义和分类能力两个角度对特征变量进行了优选,既可以剔除冗余特征变量,有效提高电压稳定裕度的评估速度,还可以去除负效力特征变量,切实提升评估结果的准确性。利用决策树模型对电力***静态电压稳定状态进行评估,不仅可以得到准确的评估结果,并且可以为***运行人员展现***的失稳机理,为其提供辅助决策依据,以便针对越限状态变量采取相应的控制措施。
以上是本申请实施例提供的基于改进Relief-F算法的电压稳定裕度评估方法,基于同样的发明构思,本申请实施例还提供了基于改进Relief-F算法的电压稳定裕度评估设备,图3为本申请实施例提供的基于改进Relief-F算法的电压稳定裕度评估评估设备示意图,如图3所示,该设备主要包括:至少一个处理器301;以及,与至少一个处理器通信连接的存储器302;其中,存储器302存储有可被至少一个处理器301执行的指令,指令被至少一个处理器301执行,以使至少一个处理器301能够完成:根据预训练的决策树模型的路径结构提取电压稳定裕度评估规则;其中,用于构建所述决策树模型的训练样本集,通过引入了惩罚因子的Relief-F算法划分,所述惩罚因子用于识别异类标签以弥补Relief-F算法固有缺陷,增强场景适应性;
将向量测量单元采集的实测数据与由决策树模型提取的所述电压稳定裕度评估规则进行匹配,得出***当前电压稳定状态。
除此之外,本申请实施例还提供了基于改进Relief-F算法的电压稳定裕度评估的非易失性计算机存储介质,存储有计算机可执行指令,其特征在于,所述计算机可执行指令设置为:根据预训练的决策树模型的路径结构提取电压稳定裕度评估规则;其中,用于构建所述决策树模型的训练样本集,通过引入了惩罚因子的Relief-F算法划分,所述惩罚因子用于识别异类标签以弥补Relief-F算法固有缺陷,增强场景适应性;
将向量测量单元采集的实测数据与由决策树模型提取的所述电压稳定裕度评估规则进行匹配,得出***当前电压稳定状态。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
本申请中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (10)
1.基于改进Relief-F算法的电压稳定裕度评估方法,其特征在于,所述方法包括:
根据预训练的决策树模型的路径结构提取电压稳定裕度评估规则;其中,用于构建所述决策树模型的训练样本集,通过引入了惩罚因子的Relief-F算法划分,所述惩罚因子用于识别异类标签以弥补Relief-F算法固有缺陷,增强场景适应性;
将向量测量单元采集的实测数据与由决策树模型提取的所述电压稳定裕度评估规则进行匹配,得出***当前电压稳定状态。
2.根据权利要求1所述的基于改进Relief-F算法的电压稳定裕度评估方法,其特征在于,所述决策树模型的训练过程为:
基于母线参与因子、支路参与因子、发电机参与因子以及各自的第一预设阈值,筛选初始特征变量;
在Relief-F算法中引入惩罚因子,并基于所述改进Relief-F算法计算所述初始特征变量的权重值;
基于所述权重值确定关键特征变量;
根据所述关键特征变量在决策树模型各级节点对训练样本集进行逐级划分,并将划分后得到的训练样本子集输入决策树模型的次级节点进行训练,直到节点中训练样本同属一类。
3.根据权利要求2所述的基于改进Relief-F算法的电压稳定裕度评估方法,其特征在于,获取训练样本集的过程为:
在***初始运行点进行P-V曲线分析,确定所述P-V曲线的最大功率传输点;
基于所述最大功率传输点,计算电压稳定裕度;
基于所述电压稳定裕度,划分电压稳定状态区间;其中,所述电压稳定状态区间包括:正常状态区间、预警状态区间、紧急状态区间;
在每个电压稳定状态区间内随机抽取等量运行点作为训练样本。
4.根据权利要求2所述的基于改进Relief-F算法的电压稳定裕度评估方法,其特征在于,所述基于所述母线参与因子、支路参与因子、发电机参与因子以及各自的第一预设阈值,筛选初始特征变量,具体为:
基于***运行状态变量构建潮流方程;其中,所述***运行状态变量包括:母线电压的幅值和相角、发电机无功功率以及支路流过的无功功率;
对所述潮流方程依次进行线性化、矩阵降阶以及特征值分析,计算得到对应最小非零特征值;
基于所述最小非零特征值,计算各条母线的参与因子、各条支路的参与因子以及各台发电机的参与因子数值;
对各参与因子分别降序排列,设定各自的第一预设阈值以筛选所述初始特征变量。
5.根据权利要求4所述的基于改进Relief-F算法的电压稳定裕度评估方法,其特征在于,计算各条母线的参与因子,具体为:
对所述潮流方程进行线性化、矩阵降阶以及特征值分析,以得到所述潮流方程降阶雅克比矩阵对应的最小非零特征值;
计算所述最小非零特征值的左特征向量以及右特征向量;
将所述左特征向量以及所述右特征向量进行相乘运算,以得到母线参与因子。
6.根据权利要求4所述的基于改进Relief-F算法的电压稳定裕度评估方法,其特征在于,计算各条支路的参与因子以及各台发电机的参与因子,具体为:
确定待测支路的无功损耗以及所有支路中无功损耗的最大值;
将所述待测支路的无功损耗以及所有支路中无功损耗的最大值做相除运算,以得到所述待测支路的支路参与因子;
确定待测发电机输出的无功功率以及所有发电机输出的无功功率最大值;
将所述待测发电机输出的无功功率以及所述所有发电机输出的无功功率最大值做相除运算,以得到所述待测发电机的参与因子。
7.根据权利要求2所述的基于改进Relief-F算法的电压稳定裕度评估方法,其特征在于,所述根据所述关键特征变量在决策树模型各级节点对训练样本集进行逐级划分,具体为:
利用改进的Relief-F算法计算各初始特征变量权重值;
对所述各初始特征变量按照权重值进行降序排列;
设定第二预设阈值,确定相应初始特征变量作为关键特征变量;
基于所述关键特征变量对所述训练样本集进行划分。
8.根据权利要求2所述的基于改进Relief-F算法的电压稳定裕度评估方法,其特征在于,所述将划分后得到的训练样本子集输入决策树模型的次级节点进行训练,具体为:
将训练样本集输入决策树模型的根节点,并计算所述训练样本集的初始熵;
计算所有所述关键特征变量划分所述训练样本集得到的属性熵;
基于所述初始熵以及所述属性熵,依次计算信息增益以及信息增益率;
根据所述信息增益率确定最终划分属性;
利用所述最终划分属性对所述训练样本集进行划分,以得到训练样本子集;
将所述训练样本子集传输到所述决策树模型的次级节点;
重复上述过程,直至所述训练样本子集的初始熵为零,则决策树模型训练完成。
9.基于改进Relief-F算法的电压稳定裕度评估设备,其特征在于,所述设备包括:至少一个处理器;以及,
与所述至少一个处理器通信连接的存储器;其中,
所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够:
根据预训练的决策树模型的路径结构提取电压稳定裕度评估规则;其中,用于构建所述决策树模型的训练样本集,通过引入了惩罚因子的Relief-F算法划分,所述惩罚因子用于识别异类标签以弥补Relief-F算法固有缺陷,增强场景适应性;
将向量测量单元采集的实测数据与由决策树模型提取的所述电压稳定裕度评估规则进行匹配,得出***当前电压稳定状态。
10.基于改进Relief-F算法的电压稳定裕度评估的非易失性计算机存储介质,存储有计算机可执行指令,其特征在于,所述计算机可执行指令设置为:
根据预训练的决策树模型的路径结构提取电压稳定裕度评估规则;其中,用于构建所述决策树模型的训练样本集,通过引入了惩罚因子的Relief-F算法划分,所述惩罚因子用于识别异类标签以弥补Relief-F算法固有缺陷,增强场景适应性;
将向量测量单元采集的实测数据与由决策树模型提取的所述电压稳定裕度评估规则进行匹配,得出***当前电压稳定状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210814511.1A CN115271070A (zh) | 2022-07-12 | 2022-07-12 | 基于改进Relief-F算法的电压稳定裕度评估方法、设备及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210814511.1A CN115271070A (zh) | 2022-07-12 | 2022-07-12 | 基于改进Relief-F算法的电压稳定裕度评估方法、设备及介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115271070A true CN115271070A (zh) | 2022-11-01 |
Family
ID=83766354
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210814511.1A Pending CN115271070A (zh) | 2022-07-12 | 2022-07-12 | 基于改进Relief-F算法的电压稳定裕度评估方法、设备及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115271070A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117215205A (zh) * | 2023-11-09 | 2023-12-12 | 国网经济技术研究院有限公司 | 基于决策树和iss理论的直流***控制参数分析方法 |
-
2022
- 2022-07-12 CN CN202210814511.1A patent/CN115271070A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117215205A (zh) * | 2023-11-09 | 2023-12-12 | 国网经济技术研究院有限公司 | 基于决策树和iss理论的直流***控制参数分析方法 |
CN117215205B (zh) * | 2023-11-09 | 2024-02-06 | 国网经济技术研究院有限公司 | 基于决策树和iss理论的直流***控制参数分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110311376B (zh) | 一种电力***动态安全评估综合模型及时空可视化方法 | |
CN110417011B (zh) | 一种基于互信息与迭代随机森林的在线动态安全评估方法 | |
CN109214708B (zh) | 基于交叉熵理论优化支持向量机的电力***风险评估方法 | |
CN108053128A (zh) | 一种基于elm和tf的电网暂态稳定快速评估方法 | |
CN112200694B (zh) | 基于图神经网络的主导失稳模式识别模型构建及应用方法 | |
CN107292502B (zh) | 一种配电网可靠性评估方法 | |
CN107067143A (zh) | 一种设备安全等级分类方法 | |
CN111401749A (zh) | 一种基于随机森林与极限学习回归的动态安全评估方法 | |
CN110417005B (zh) | 结合深度学习和仿真计算的暂态稳定严重故障筛选方法 | |
CN106649479A (zh) | 一种基于概率图的变压器状态关联规则挖掘方法 | |
CN112821424B (zh) | 一种基于数据-模型融合驱动的电力***频率响应分析方法 | |
Zheng et al. | Real-time transient stability assessment based on deep recurrent neural network | |
CN104578053A (zh) | 基于受扰电压轨迹簇特征的电力***暂态稳定性预测方法 | |
CN112069727A (zh) | 具备高可信度的电力***暂态稳定智能化评估***及方法 | |
Yang et al. | Analysis of power system transient stability characteristics with the application of massive transient stability simulation data | |
CN111401476A (zh) | 基于边界区域重要性采样与核向量机的暂态安全评估方法 | |
CN110705831A (zh) | 电力***故障后功角失稳模式预判模型构建方法及其应用 | |
CN115271070A (zh) | 基于改进Relief-F算法的电压稳定裕度评估方法、设备及介质 | |
CN115526258A (zh) | 基于Spearman相关系数特征提取的电力***暂稳评估方法 | |
CN113128130B (zh) | 一种判断直流配电***稳定性的实时监测方法及装置 | |
Tiwary et al. | Multi-dimensional ANN application for active power flow state classification on a utility system | |
CN116011982A (zh) | 一种磨煤机磨辊断裂在线监测方法及*** | |
Huang et al. | Sequential detection of microgrid bad data via a data-driven approach combining online machine learning with statistical analysis | |
Tomin et al. | Development of automatic intelligent system for on-line voltage security control of power systems | |
Hu et al. | Power grid's Intelligent Stability Analysis based on big data technology |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |