CN107069708B - 一种基于极限学习机的输电网线路有功安全校正方法 - Google Patents

一种基于极限学习机的输电网线路有功安全校正方法 Download PDF

Info

Publication number
CN107069708B
CN107069708B CN201710151119.2A CN201710151119A CN107069708B CN 107069708 B CN107069708 B CN 107069708B CN 201710151119 A CN201710151119 A CN 201710151119A CN 107069708 B CN107069708 B CN 107069708B
Authority
CN
China
Prior art keywords
generator
power
load
sensitivity
active
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.)
Expired - Fee Related
Application number
CN201710151119.2A
Other languages
English (en)
Other versions
CN107069708A (zh
Inventor
李淼
周强明
周悦
鲁鸿毅
曾鹏
姜盛波
杨军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
State Grid Hubei Electric Power Co Ltd
Original Assignee
Wuhan University WHU
State Grid Hubei Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU, State Grid Hubei Electric Power Co Ltd filed Critical Wuhan University WHU
Priority to CN201710151119.2A priority Critical patent/CN107069708B/zh
Publication of CN107069708A publication Critical patent/CN107069708A/zh
Application granted granted Critical
Publication of CN107069708B publication Critical patent/CN107069708B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明属于电力***运行与控制领域,是涉及一种基于极限学习机的输电网线路有功安全校正方法。本发明在考虑现有研究基于灵敏度的启发式有功安全校正方法时,在调整策略制定的快速性和精确度方面往往顾此失彼提出一种基于极限学习机的输电线路有功安全校正方法,均衡有功安全校正方法的快速性和精确度。该方法首先基于极限学习机对不同状态下的电网运行数据进行分析学习,获得各节点注入功率对输电线路的灵敏度;然后基于节点注入功率对输电线路的灵敏度建立有功安全校正模型;最后利用粒子群智能算法对有功安全校正模型进行求解获得发电机和负荷调整方案。

Description

一种基于极限学习机的输电网线路有功安全校正方法
技术领域
本发明属于电力***运行与控制领域,是涉及一种基于极限学习机的输电网线路有功安全校正方法。
背景技术
智能电网和大电网互联技术的发展,在解决我国能源分布格局不均、提高电网坚强性方面起到重大作用,但同时造成我国电网互联程度及复杂性空前提高,对电网安全稳定运行提出更高要求。研究表明,连锁故障是造成大停电事故的重要原因,故障元件引起潮流转移会导致正常运行线路变为重载、过载运行状态,这是连锁故障蔓延的重要因素。因此,对处于重载、过载状态的运行输电线路进行有功安全校正,调度部门需要采取快速有效的发电机出力调整及切负荷措施,使输电线路恢复正常运行,这对于抑制连锁故障发展保证电力***安全稳定运行具有重要意义。
现阶段,在输电线路的有功安全校正领域,专家学者主要采用优化规划类算法和灵敏度算法。优化规划类算法通过建立有功安全校正数学模型,设定目标函数和各类约束条件,并利用数学解析方法求解数学模型获得校正措施,此类方法考虑全面,具有较好的安全性和经济性,但可能存在求解时间过长、调整设备过多、潮流收敛条件难以满足等问题;灵敏度方法通过求取发电机对过载线路的灵敏度,进而根据线路过载量获得校正措施,此类方法可以快速获得调整调整措施消除过载,不存在潮流收敛性问题,但是灵敏度的精确性将直接影响校正方案的有效性。
发明内容
本发明的上述技术问题主要是通过下述技术方案得以解决的:
一种基于极限学习机的输电线路有功安全校正方法,其特征在于,包括:
步骤1,采集发电机灵敏度样本集合以及负荷灵敏度样本集合,基于极限学习机对采集的电网样本集合进行学习训练,获得节点注入功率对输电线路的灵敏度与电网运行数据间的映射关系;
步骤2,基于步骤1的灵敏度分析,建立输电线路的有功安全校正模型,所述有功安全校正模型基于目标函数:
Figure GDA0001296865170000021
Figure GDA0001296865170000022
为减出力发电机调整量,
Figure GDA0001296865170000023
为加出力发电机调整量,ΔPG,z为平衡机调整量,
Figure GDA0001296865170000024
为负荷切除量,M为惩罚系数,为远大于1的正常数,保证消除线路过载时优先调整发电机出力,再进行切负荷,目标函数以发电机和负荷调整量最小为目标:
步骤3,采集待控制电网数据,基于步骤2的有功安全校正模型,求解输电线路有功安全校正模型,到发电机和负荷的调整量,确定控制方案。
在上述的一种基于极限学习机的输电线路有功安全校正方法,所述步骤1具体包括:
步骤1.1、生成极限学习机训练样本集合,包括发电机灵敏度样本集合和负荷灵敏度样本集合,在任意电网***中,首先确定电网的初始运行状态和网络拓扑结构,记录电网初始状态下的发电机有功出力PG 0、负荷有功
Figure GDA0001296865170000025
以及各传输线路的有功功率
Figure GDA0001296865170000026
具体是:
集合一,发电机灵敏度样本集合:针对任意一台发电机(除平衡节点外),随机改变该发电机的有功出力
Figure GDA0001296865170000031
随机改变出力可保证得到的样本数据能全面反应电网的不同运行状态。为保证***有功功率平衡,在改变发电机有功出力的同时,需要反向调整平衡机的有功出力
Figure GDA0001296865170000032
调整发电出力之后,更新电网的潮流信息,并记录各传输线路的有功功率
Figure GDA0001296865170000033
及有功变化量
Figure GDA0001296865170000034
根据所记录的信息可以计算被调整的发电机对输电线路的灵敏度
Figure GDA0001296865170000035
可得到一组样本数据(x1,t1),其中
Figure GDA0001296865170000036
该过程重复N次,每一次可以得到一组(xi,ti),其中
Figure GDA0001296865170000037
最终可得到发电机灵敏度训练样本集合
Figure GDA0001296865170000038
其中n=3,m=1。
集合二,负荷灵敏度样本集合:针对任意一个除平衡节点外的负荷节点,随机改变有功负荷
Figure GDA0001296865170000039
同时同向调整平衡机出力
Figure GDA00012968651700000310
更新电网潮流信息,记录各传输线路的有功功率
Figure GDA00012968651700000311
及有功变化量
Figure GDA00012968651700000312
计算被调整的负荷节点对输电线路的灵敏度
Figure GDA00012968651700000313
可得到一组样本数据(x1,t1),其中
Figure GDA00012968651700000314
该过程重复N次,每次可得一组(xi,ti),其中
Figure GDA00012968651700000315
最终可得到负荷灵敏度训练样本集合
Figure GDA00012968651700000316
其中n=3,m=1。
步骤1.2、基于极限学习机对样本集合进行训练:基于步骤1.1给定的灵敏度训练样本集合
Figure GDA00012968651700000317
选定隐含层节点个数为L,则将极限学习机的数学模型表示为:
Figure GDA00012968651700000318
式(1)中,gi为激励函数wi=[wi1,wi2,...,win]为输入层的权重向量,βi=[βi1i2,...,βim]为输出权重向量,即隐含层到输出层的连接权重向量,bi为第i个隐含层的偏置。
矩阵表达形式为:
Hβ=Y (2)
式(2)中,
Figure GDA0001296865170000046
Figure GDA0001296865170000041
H为神经网络的隐含层输出,第i列向量对应第i个隐含层节点的输出。
已知当激励函数g无限可微时,网络参数并不需要全部调整。因此随机设置输入权值wi和隐含层的偏置bi,这样对于给定的训练样本集合,隐含层的输出矩阵H可被唯一确定,模型中仅存在变量β。
构造极限学习机关于β的误差函数:
Figure GDA0001296865170000042
此时极限学习机对样本集合的训练过程可转化为求解输出权值矩阵的最小二乘解问题:
minβ||Hβ-T|| (4)
求解式(4)可得:
Figure GDA0001296865170000043
式(5)中
Figure GDA0001296865170000044
是矩阵H的Moore-Penrose广义逆。
最终得到灵敏度和电网运行信息之间的映射关系:
Figure GDA0001296865170000045
在上述的一种基于极限学习机的输电线路有功安全校正方法,所述步骤2基于定义:电网在某一运行状态下,线路Lm发生过载,过载量为
Figure GDA0001296865170000051
其中Pm为线路Lm的实际有功功率,
Figure GDA0001296865170000052
为线路Lm的有功传输极限,为消除过载并保留一定裕度,需要将线路Lm的有功功率降低
Figure GDA0001296865170000053
已知电网各发电机的有功出力PG和负荷有功大小PL,因此可确定式(6)的输入向量xG=[ΔPc,Pm,PG]和xL=[ΔPc,Pm,PL],求得各发电机和负荷对线路Lm的灵敏度;具体包括有约束条件的目标函数,其中约束条件包括:
约束条件a.***有功功率平衡约束,保证发电机调整量和负荷调整量相等。
Figure GDA0001296865170000054
约束条件b.消除线路过载约束,保证调整后过载现象消除。
Figure GDA0001296865170000055
式(9)中,
Figure GDA0001296865170000056
为减出力发电机对过载线路Lm的灵敏度,
Figure GDA0001296865170000057
为加出力发电机对过载线路Lm的灵敏度,
Figure GDA0001296865170000058
为负荷对过载线路Lm的灵敏度。
约束条件c.支路有功运行安全约束,保证调整后其他线路安全运行。
Figure GDA0001296865170000059
式(10)中,t为除过载线路以外的其它线路编号,
Figure GDA00012968651700000510
分别为线路t的初始有功功率、最小有功功率限值、最大有功功率限值。
约束条件d.发电机组出力调整量约束。
Figure GDA00012968651700000511
式(11)中,
Figure GDA00012968651700000512
分别为减出力发电机的初始有功出力和有功出力下限,
Figure GDA00012968651700000513
分别为加出力发电机的初始有功出力和有功出力上限。
约束条件e.负荷调整量约束。
Figure GDA0001296865170000061
式(12)中,
Figure GDA0001296865170000062
为初始有功负荷。
在上述的一种基于极限学习机的输电线路有功安全校正方法,所述步骤4基于定义电网中某条线路Lm过载时,待调整发电机有u台,待切除负荷有v个。基于粒子群智能算法,对2.2建立的有功安全校正模型进行求解,得到发电机和负荷的调整量,确定控制方案。具体步骤如下:
步骤4.1、设定最大迭代次数并初始化粒子群:设定最大迭代次数为MaxD,粒子群规模为N,粒子维数为u+v,并按照式(8)、式(11)和式(12)所示约束条件,对每个粒子初始位置xi和初始速度vi进行初始化,xi=[ΔPG,1,ΔPG,2,...,ΔPG,u,ΔPL,1,ΔPL,2,...ΔPL,v],vi=[vi1,vi2,...vi(u+v)]。
步骤4.2、计算每个粒子的初始适应度值,确定粒子个体极值pi和粒子群全局极值pg,适应度函数如式(13)所示。
Figure GDA0001296865170000063
式(13)中ΔPG,i为待调发电机的调整量,ΔPL,j为待切负荷的切除量,M为惩罚系数。
步骤4.3、更新每个粒子的速度vi和位置xi,保证xi满足式(8)、式(11)和式(12)所示约束条件,按照式(14)进行更新。
vi=w*vi+c1r1(pi-xi)+c2r2(pg-xi),xi=xi+vi (14)
式(14)中w为惯性权重,c1和c2为学习因子,r1和r2为[0,1]范围内的均匀随机数。
步骤4.4、计算更新位置后的每个粒子xi的适应度值F(xi),并与个体极值pi的适应度值F(pi)比较,若F(xi)<F(pi),则用xi代替pi;对于每个粒子xi,将其适应度值F(xi)与全局极值F(pg)比较,若F(xi)<F(pg),则用xi代替pg
步骤4.5、针对每个粒子xi,按照式(9)进行过载消除校验,按照式(10)进行线路安全运行约束校验。若满足校验,则保存局部最优pi和全局最优pg;对不满足校验的粒子按步骤1进行粒子更新。
步骤4.6、检验是否满足搜索终止条件,即迭代次数是否大于最大迭代次数。若满足终止条件,则输出全局最优pg,得到控制方案;否则返回步骤3。
在上述的一种基于极限学***衡,优先调整发电机出力再进行负荷调整,同时达到发电机和负荷调整量最小的目标;根据各发电机和负荷对过载线路Lm的灵敏度,对发电机和负荷进行分组:发电机灵敏度为正,则归为减出力机组G-;发电机灵敏度为负,则归为加出力机组G+;发电机灵敏度为零,则归为平衡机组G0;对于负荷节点,选取灵敏度为负的节点作为切负荷集合L-。然后按照绝对值大小分别对集合内的发电机和负荷进行降序排列。具体调整策略包括:
步骤5.1、从灵敏度为负的机组中按照绝对值大小顺序进行加出力调整,从灵敏度为负的机组中按照绝对值大小顺序进行减出力时调整,要保证一个机组加出力时存在另一个机组减出力,且调整量绝对值相等。
步骤5.2、一台机组可与多台机组进行顺序配对调整,直至其调整能力完全发挥。
步骤5.3、当集合G-中机组没有减出力的能力时,集合G0中的机组可以减出力;当集合G+中机组没有加出力的能力时,集合中的机组G0可以加出力。
步骤5.4、当发电机调整不能消除线路过载时,再进行切负荷处理,直到支线路过载消除。
因此,本发明具有如下优点:本发明以电力***分析理论为基础,综合考虑规划类算法和灵敏度算法的优缺点,建立基于灵敏度的启发式有功安全校正模型,同时考虑直流灵敏度存在的精确性不足问题,基于极限学习机方法,高效学习历史数据集获得精确灵敏度值,实现有功安全校正方案的快速生成的同时,保证有功安全校正方案的精确度。
附图说明
附图1极限学习机网络结构图。
附图2优化模型求解流程图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。
实施例:
下面具体介绍本发明的方法步骤,具体包括:
1.利用极限学习机对电网运行数据进行学习训练,获得节点注入功率对输电线路的灵敏度与电网运行数据间的映射关系。
1.1极限学习机训练样本集合的生成
给定某一电网***,首先确定电网的初始运行状态和网络拓扑结构,记录电网初始状态下的发电机有功出力PG 0、负荷有功
Figure GDA0001296865170000081
以及各传输线路的有功功率
Figure GDA0001296865170000082
(1)发电机灵敏度样本集合
针对任意一台发电机(除平衡节点外),随机改变该发电机的有功出力
Figure GDA0001296865170000091
随机改变出力可保证得到的样本数据能全面反应电网的不同运行状态。为保证***有功功率平衡,在改变发电机有功出力的同时,需要反向调整平衡机的有功出力
Figure GDA0001296865170000092
调整发电出力之后,更新电网的潮流信息,并记录各传输线路的有功功率
Figure GDA0001296865170000093
及有功变化量
Figure GDA0001296865170000094
根据所记录的信息可以计算被调整的发电机对输电线路的灵敏度
Figure GDA0001296865170000095
可得到一组样本数据(x1,t1),其中
Figure GDA0001296865170000096
该过程重复N次,每一次可以得到一组(xi,ti),其中
Figure GDA0001296865170000097
最终可得到发电机灵敏度训练样本集合
Figure GDA0001296865170000098
其中n=3,m=1。
(2)负荷灵敏度样本集合
类似于发电机灵敏度样本集合的生成,针对任意一个负荷节点(除平衡节点外),随机改变有功负荷
Figure GDA0001296865170000099
同时同向调整平衡机出力
Figure GDA00012968651700000910
更新电网潮流信息,记录各传输线路的有功功率
Figure GDA00012968651700000911
及有功变化量
Figure GDA00012968651700000912
计算被调整的负荷节点对输电线路的灵敏度
Figure GDA00012968651700000913
可得到一组样本数据(x1,t1),其中
Figure GDA00012968651700000914
该过程重复N次,每次可得一组(xi,ti),其中
Figure GDA00012968651700000915
最终可得到负荷灵敏度训练样本集合
Figure GDA00012968651700000916
其中n=3,m=1。
1.2基于极限学习机对样本集合进行训练
极限学习机是由黄广斌提出来的求解神经网络算法,网络结构图如图1所示。其最大的特点是相对于传统的神经网络,尤其是单隐层前馈神经网络,极限学习机的训练速度更快,训练过程不需要迭代调整输入权重和偏置,仅需求得输出权重β,非常适合在线计算和更新。
基于1.1节给定的灵敏度训练样本集合
Figure GDA00012968651700000917
选定隐含层节点个数为L,则可以将极限学习机的数学模型表示为:
Figure GDA0001296865170000101
式(1)中,gi为激励函数wi=[wi1,wi2,...,win]为输入层的权重向量,βi=[βi1i2,...,βim]为输出权重向量,即隐含层到输出层的连接权重向量,bi为第i个隐含层的偏置。
矩阵表达形式为:
Hβ=Y (2)
式(2)中,
Figure GDA0001296865170000102
Figure GDA0001296865170000103
H为神经网络的隐含层输出,第i列向量对应第i个隐含层节点的输出。
已知当激励函数g无限可微时,网络参数并不需要全部调整。因此随机设置输入权值wi和隐含层的偏置bi,这样对于给定的训练样本集合,隐含层的输出矩阵H可被唯一确定,模型中仅存在变量β。
构造极限学习机关于β的误差函数:
Figure GDA0001296865170000104
此时极限学习机对样本集合的训练过程可转化为求解输出权值矩阵的最小二乘解问题:
minβ||Hβ-T|| (4)
求解式(4)可得:
Figure GDA0001296865170000105
式(5)中
Figure GDA0001296865170000111
是矩阵H的Moore-Penrose广义逆。
最终得到灵敏度和电网运行信息之间的映射关系:
Figure GDA0001296865170000112
2.基于灵敏度分析,建立输电线路的有功安全校正模型。
假设电网在某一运行状态下,线路Lm发生过载,过载量为
Figure GDA0001296865170000113
其中Pm为线路Lm的实际有功功率,
Figure GDA0001296865170000114
为线路Lm的有功传输极限,为消除过载并保留一定裕度,需要将线路Lm的有功功率降低
Figure GDA0001296865170000115
已知电网各发电机的有功出力PG和负荷有功大小PL,因此可确定式(6)的输入向量xG=[ΔPc,Pm,PG]和xL=[ΔPc,Pm,PL],求得各发电机和负荷对线路Lm的灵敏度。
2.1基于等量反向配对原则的发电机、负荷调整策略
通过调整发电机出力和切负荷是常用的有功安全校正措施,在调整过程中,需要保证***的功率平衡,优先调整发电机出力再进行负荷调整,同时达到发电机和负荷调整量最小的目标。
根据各发电机和负荷对过载线路Lm的灵敏度,对发电机和负荷进行分组:发电机灵敏度为正,则归为减出力机组G-;发电机灵敏度为负,则归为加出力机组G+;发电机灵敏度为零,则归为平衡机组G0;对于负荷节点,选取灵敏度为负的节点作为切负荷集合L-。然后按照绝对值大小分别对集合内的发电机和负荷进行降序排列。具体调整策略如下:
a.从灵敏度为负的机组中按照绝对值大小顺序进行加出力调整,从灵敏度为负的机组中按照绝对值大小顺序进行减出力时调整,要保证一个机组加出力时存在另一个机组减出力,且调整量绝对值相等。
b.一台机组可与多台机组进行顺序配对调整,直至其调整能力完全发挥。
c.当集合G-中机组没有减出力的能力时,集合G0中的机组可以减出力;当集合G+中机组没有加出力的能力时,集合中的机组G0可以加出力。
d.当发电机调整不能消除线路过载时,再进行切负荷处理,直到支线路过载消除。
2.2建立输电线路有功安全校正模型
针对2.1介绍的基于等量反向配对原则的发电机、负荷调整策略,用数学公式进行表达,建立输电线路的有功安全校正模型。
(1)目标函数
目标函数以发电机和负荷调整量最小为目标:
Figure GDA0001296865170000121
式(7)中,
Figure GDA0001296865170000122
为减出力发电机调整量,
Figure GDA0001296865170000123
为加出力发电机调整量,ΔPG,z为平衡机调整量,
Figure GDA0001296865170000124
为负荷切除量,M为惩罚系数,为远大于1的正常数,保证消除线路过载时优先调整发电机出力,再进行切负荷。
(2)约束条件
a.***有功功率平衡约束,保证发电机调整量和负荷调整量相等。
Figure GDA0001296865170000125
b.消除线路过载约束,保证调整后过载现象消除。
Figure GDA0001296865170000126
式(9)中,
Figure GDA0001296865170000127
为减出力发电机对过载线路Lm的灵敏度,
Figure GDA0001296865170000128
为加出力发电机对过载线路Lm的灵敏度,
Figure GDA0001296865170000129
为负荷对过载线路Lm的灵敏度。
c.支路有功运行安全约束,保证调整后其他线路安全运行。
Figure GDA00012968651700001210
式(10)中,t为除过载线路以外的其它线路编号,
Figure GDA0001296865170000131
分别为线路t的初始有功功率、最小有功功率限值、最大有功功率限值。
d.发电机组出力调整量约束。
Figure GDA0001296865170000132
式(11)中,
Figure GDA0001296865170000133
分别为减出力发电机的初始有功出力和有功出力下限,
Figure GDA0001296865170000134
分别为加出力发电机的初始有功出力和有功出力上限。
e.负荷调整量约束。
Figure GDA0001296865170000135
式(12)中,
Figure GDA0001296865170000136
为初始有功负荷。
3.基于粒子群智能算法,求解输电线路有功安全校正模型。
粒子群算法通过粒子间的竞争和协作实现在复杂搜索空间中寻找全局最优点,具有易实现和全局搜索能力强的特点,在解决优化全局优化问题方面具有独特优势。
假设电网中某条线路Lm过载时,待调整发电机有u台,待切除负荷有v个。基于粒子群智能算法,对2.2建立的有功安全校正模型进行求解,得到发电机和负荷的调整量,确定控制方案。具体步骤如下,流程图如图2所示。
1.设定最大迭代次数并初始化粒子群:设定最大迭代次数为MaxD,粒子群规模为N,粒子维数为u+v,并按照式(8)、式(11)和式(12)所示约束条件,对每个粒子初始位置xi和初始速度vi进行初始化,
xi=[ΔPG,1,ΔPG,2,...,ΔPG,u,ΔPL,1,ΔPL,2,...ΔPL,v],vi=[vi1,vi2,...vi(u+v)]。
2.计算每个粒子的初始适应度值,确定粒子个体极值pi和粒子群全局极值pg,适应度函数如式(13)所示。
Figure GDA0001296865170000137
式(13)中ΔPG,i为待调发电机的调整量,ΔPL,j为待切负荷的切除量,M为惩罚系数。
3.更新每个粒子的速度vi和位置xi,保证xi满足式(8)、式(11)和式(12)所示约束条件,按照式(14)进行更新。
vi=w*vi+c1r1(pi-xi)+c2r2(pg-xi),xi=xi+vi (14)
式(14)中w为惯性权重,c1和c2为学习因子,r1和r2为[0,1]范围内的均匀随机数。
4.计算更新位置后的每个粒子xi的适应度值F(xi),并与个体极值pi的适应度值F(pi)比较,若F(xi)<F(pi),则用xi代替pi;对于每个粒子xi,将其适应度值F(xi)与全局极值F(pg)比较,若F(xi)<F(pg),则用xi代替pg
5.针对每个粒子xi,按照式(9)进行过载消除校验,按照式(10)进行线路安全运行约束校验。若满足校验,则保存局部最优pi和全局最优pg;对不满足校验的粒子按步骤1进行粒子更新。
6.检验是否满足搜索终止条件,即迭代次数是否大于最大迭代次数。若满足终止条件,则输出全局最优pg,得到控制方案;否则返回步骤3。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

Claims (4)

1.一种基于极限学习机的输电线路有功安全校正方法,其特征在于,包括:
步骤1,采集发电机灵敏度样本集合以及负荷灵敏度样本集合,基于极限学习机对采集的电网样本集合进行学习训练,获得节点注入功率对输电线路的灵敏度与电网运行数据间的映射关系;
步骤2,基于步骤1的灵敏度分析,建立输电线路的有功安全校正模型,所述有功安全校正模型基于目标函数:
Figure FDA0002731259490000011
Figure FDA0002731259490000012
为减出力发电机调整量,
Figure FDA0002731259490000013
为加出力发电机调整量,ΔPG,z为平衡机调整量,
Figure FDA0002731259490000014
为负荷切除量,M为惩罚系数,为远大于1的正常数,保证消除线路过载时优先调整发电机出力,再进行切负荷,目标函数以发电机和负荷调整量最小为目标,发电机灵敏度为正,则归为减出力机组G-;发电机灵敏度为负,则归为加出力机组G+;发电机灵敏度为零,则归为平衡机组G0;对于负荷节点,选取灵敏度为负的节点作为切负荷集合L-
步骤3,采集待控制电网数据,基于步骤2的有功安全校正模型,求解输电线路有功安全校正模型,到发电机和负荷的调整量,确定控制方案;
所述步骤1具体包括:
步骤1.1、生成极限学习机训练样本集合,包括发电机灵敏度样本集合和负荷灵敏度样本集合,在任意电网***中,首先确定电网的初始运行状态和网络拓扑结构,记录电网初始状态下的发电机有功出力PG 0、负荷有功
Figure FDA0002731259490000015
以及各传输线路的有功功率
Figure FDA0002731259490000016
具体是:
集合一,发电机灵敏度样本集合:针对任意一台发电机,随机改变该发电机的有功出力
Figure FDA0002731259490000017
随机改变出力可保证得到的样本数据能全面反应电网的不同运行状态;为保证***有功功率平衡,在改变发电机有功出力的同时,需要反向调整平衡机的有功出力
Figure FDA0002731259490000021
调整发电出力之后,更新电网的潮流信息,并记录各传输线路的有功功率
Figure FDA0002731259490000022
及有功变化量
Figure FDA0002731259490000023
根据所记录的信息可以计算被调整的发电机对输电线路的灵敏度
Figure FDA0002731259490000024
可得到一组样本数据(x1,t1),其中
Figure FDA0002731259490000025
该过程重复N次,每一次可以得到一组(xi,ti),其中
Figure FDA0002731259490000026
最终可得到发电机灵敏度训练样本集合
Figure FDA0002731259490000027
其中n=3,m=1;
集合二,负荷灵敏度样本集合:针对任意一个除平衡节点外的负荷节点,随机改变有功负荷ΔPL 1,同时同向调整平衡机出力
Figure FDA0002731259490000028
更新电网潮流信息,记录各传输线路的有功功率
Figure FDA0002731259490000029
及有功变化量
Figure FDA00027312594900000210
计算被调整的负荷节点对输电线路的灵敏度
Figure FDA00027312594900000211
可得到一组样本数据(x1,t1),其中
Figure FDA00027312594900000212
该过程重复N次,每次可得一组(xi,ti),其中
Figure FDA00027312594900000213
最终可得到负荷灵敏度训练样本集合
Figure FDA00027312594900000214
其中n=3,m=1;
步骤1.2、基于极限学习机对样本集合进行训练:基于步骤1.1给定的灵敏度训练样本集合
Figure FDA00027312594900000215
选定隐含层节点个数为L,则将极限学习机的数学模型表示为:
Figure FDA00027312594900000216
式(1)中,gi为激励函数wi=[wi1,wi2,...,win]为输入层的权重向量,βi=[βi1i2,...,βim]为输出权重向量,即隐含层到输出层的连接权重向量,bi为第i个隐含层的偏置;
矩阵表达形式为:
Hβ=Y(2)
式(2)中,
Figure FDA0002731259490000031
Figure FDA0002731259490000032
H为神经网络的隐含层输出,第i列向量对应第i个隐含层节点的输出;
已知当激励函数wi无限可微时,网络参数并不需要全部调整;因此随机设置输入权值wi和隐含层的偏置bi,这样对于给定的训练样本集合,隐含层的输出矩阵H可被唯一确定,模型中仅存在变量β;
构造极限学习机关于β的误差函数:
Figure FDA0002731259490000033
此时极限学习机对样本集合的训练过程可转化为求解输出权值矩阵的最小二乘解问题:
minβ||Hβ-T||(4)
求解式(4)可得:
Figure FDA0002731259490000034
式(5)中
Figure FDA0002731259490000035
是矩阵H的Moore-Penrose广义逆;
最终得到灵敏度和电网运行信息之间的映射关系:
Figure FDA0002731259490000036
2.根据权利要求1所述的一种基于极限学习机的输电线路有功安全校正方法,其特征在于,所述步骤2基于定义:电网在某一运行状态下,线路Lm发生过载,过载量为
Figure FDA0002731259490000037
其中Pm为线路Lm的实际有功功率,
Figure FDA0002731259490000041
为线路Lm的有功传输极限,为消除过载并保留一定裕度,需要将线路Lm的有功功率降低
Figure FDA0002731259490000042
已知电网各发电机的有功出力PG和负荷有功大小PL,因此可确定式(6)的输入向量xG=[ΔPc,Pm,PG]和xL=[ΔPc,Pm,PL],求得各发电机和负荷对线路Lm的灵敏度;具体包括有约束条件的目标函数,其中约束条件包括:
约束条件a.***有功功率平衡约束,保证发电机调整量和负荷调整量相等;
Figure FDA0002731259490000043
约束条件b.消除线路过载约束,保证调整后过载现象消除;
Figure FDA0002731259490000044
式(9)中,
Figure FDA0002731259490000045
为减出力发电机对过载线路Lm的灵敏度,
Figure FDA0002731259490000046
为加出力发电机对过载线路Lm的灵敏度,
Figure FDA0002731259490000047
为负荷对过载线路Lm的灵敏度;
约束条件c.支路有功运行安全约束,保证调整后其他线路安全运行;
Figure FDA0002731259490000048
式(10)中,t为除过载线路以外的其它线路编号,Pt o,Pt min,Pt max分别为线路t的初始有功功率、最小有功功率限值、最大有功功率限值;
约束条件d.发电机组出力调整量约束;
Figure FDA0002731259490000049
式(11)中,
Figure FDA00027312594900000410
分别为减出力发电机的初始有功出力和有功出力下限,
Figure FDA00027312594900000411
分别为加出力发电机的初始有功出力和有功出力上限;
约束条件e.负荷调整量约束;
Figure FDA00027312594900000412
式(12)中,
Figure FDA00027312594900000413
为初始有功负荷。
3.根据权利要求2所述的一种基于极限学习机的输电线路有功安全校正方法,其特征在于,步骤2基于定义电网中某条线路Lm过载时,待调整发电机有u台,待切除负荷有v个;基于粒子群智能算法,对建立的有功安全校正模型进行求解,得到发电机和负荷的调整量,确定控制方案;具体步骤如下:
步骤4.1、设定最大迭代次数并初始化粒子群:设定最大迭代次数为MaxD,粒子群规模为N,粒子维数为u+v,并按照式(8)、式(11)和式(12)所示约束条件,对每个粒子初始位置x′i和初始速度vi进行初始化,x′i=[ΔPG,1,ΔPG,2,...,ΔPG,u,ΔPL,1,ΔPL,2,...ΔPL,v],vi=[vi1,vi2,...vi(u+v)];
步骤4.2、计算每个粒子的初始适应度值,确定粒子个体极值pi和粒子群全局极值pg,适应度函数如式(13)所示;
Figure FDA0002731259490000051
式(13)中ΔPG,i为待调发电机的调整量,ΔPL,j为待切负荷的切除量,M为惩罚系数;
步骤4.3、更新每个粒子的速度vi和位置x′i,保证x′i满足式(8)、式(11)和式(12)所示约束条件,按照式(14)进行更新;
vi=w*vi+c1r1(pi-x′i)+c2r2(pg-x′i),x′i=x′i+vi (14)
式(14)中w为惯性权重,c1和c2为学习因子,r1和r2为[0,1]范围内的均匀随机数;
步骤4.4、计算更新位置后的每个粒子x′i的适应度值F(xi),并与个体极值pi的适应度值F(pi)比较,若F(xi)<F(pi),则用x′i代替pi;对于每个粒子x′i,将其适应度值F(xi)与全局极值F(pg)比较,若F(xi)<F(pg),则用x′i代替pg
步骤4.5、针对每个粒子x′i,按照式(9)进行过载消除校验,按照式(10)进行线路安全运行约束校验;若满足校验,则保存局部最优pi和全局最优pg;对不满足校验的粒子按步骤1进行粒子更新;
步骤4.6、检验是否满足搜索终止条件,即迭代次数是否大于最大迭代次数;若满足终止条件,则输出全局最优pg,得到控制方案;否则返回步骤4.3。
4.根据权利要求1所述的一种基于极限学***衡,优先调整发电机出力再进行负荷调整,同时达到发电机和负荷调整量最小的目标;根据各发电机和负荷对过载线路Lm的灵敏度,对发电机和负荷进行分组:然后按照绝对值大小分别对集合内的发电机和负荷进行降序排列;具体调整策略包括:
步骤5.1、从灵敏度为负的机组中按照绝对值大小顺序进行加出力调整,从灵敏度为负的机组中按照绝对值大小顺序进行减出力时调整,要保证一个机组加出力时存在另一个机组减出力,且调整量绝对值相等;
步骤5.2、一台机组可与多台机组进行顺序配对调整,直至其调整能力完全发挥;
步骤5.3、当减出力机组G-中机组没有减出力的能力时,平衡机组G0中的机组可以减出力;当加出力机组G+中机组没有加出力的能力时,集合中的平衡机组G0可以加出力;
步骤5.4、当发电机调整不能消除线路过载时,再进行切负荷处理,直到支线路过载消除。
CN201710151119.2A 2017-03-14 2017-03-14 一种基于极限学习机的输电网线路有功安全校正方法 Expired - Fee Related CN107069708B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710151119.2A CN107069708B (zh) 2017-03-14 2017-03-14 一种基于极限学习机的输电网线路有功安全校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710151119.2A CN107069708B (zh) 2017-03-14 2017-03-14 一种基于极限学习机的输电网线路有功安全校正方法

Publications (2)

Publication Number Publication Date
CN107069708A CN107069708A (zh) 2017-08-18
CN107069708B true CN107069708B (zh) 2021-01-15

Family

ID=59617849

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710151119.2A Expired - Fee Related CN107069708B (zh) 2017-03-14 2017-03-14 一种基于极限学习机的输电网线路有功安全校正方法

Country Status (1)

Country Link
CN (1) CN107069708B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107967497A (zh) * 2017-08-28 2018-04-27 苏州明逸智库信息科技有限公司 基于卷积神经网络和极限学习机的手写体识别方法
CN108879665B (zh) * 2018-07-03 2021-07-27 河海大学 以调整设备数目最少为目标的电力***安全校正优化方法
CN110676852B (zh) * 2019-08-26 2020-11-10 重庆大学 考虑潮流特征的改进极限学习机快速概率潮流计算方法
CN114243690A (zh) * 2021-12-15 2022-03-25 国网河北省电力有限公司 一种电网有功安全校正方法、装置、电子设备及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104821605A (zh) * 2015-04-13 2015-08-05 国家电网公司 一种基于改进粒子群混合优化算法的有功安全校正方法
CN106356856A (zh) * 2016-09-18 2017-01-25 国电南瑞科技股份有限公司 一种基于分区负荷控制的安全校正计算方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104821605A (zh) * 2015-04-13 2015-08-05 国家电网公司 一种基于改进粒子群混合优化算法的有功安全校正方法
CN106356856A (zh) * 2016-09-18 2017-01-25 国电南瑞科技股份有限公司 一种基于分区负荷控制的安全校正计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于广域信息量的潮流转移问题研究;姜臻;《CNKI中国优秀硕士学位论文全文库》;20120715;第I,23-33页 *
基于粒子群优化的核极限学习机模型的风电功率区间预测方法;杨锡运 等;《中国电机工程学报》;20150930;第35卷;第146-153页 *
姜臻.基于广域信息量的潮流转移问题研究.《CNKI中国优秀硕士学位论文全文库》.2012, *

Also Published As

Publication number Publication date
CN107069708A (zh) 2017-08-18

Similar Documents

Publication Publication Date Title
CN107069708B (zh) 一种基于极限学习机的输电网线路有功安全校正方法
CN108694467B (zh) 一种对配电网线损率进行预测的方法及***
Sivalingam et al. A modified whale optimization algorithm-based adaptive fuzzy logic PID controller for load frequency control of autonomous power generation systems
CN109861202B (zh) 一种柔性互联配电网动态优化调度方法及***
Raglend et al. Comparison of AI techniques to solve combined economic emission dispatch problem with line flow constraints
CN110932281B (zh) 一种基于电网准稳态灵敏度的多断面协同校正方法和***
CN114362196B (zh) 一种多时间尺度主动配电网电压控制方法
CN105048479B (zh) 一种光伏电站无功分组调节方法
CN110738344A (zh) 电力***负荷预测的分布式无功优化方法及装置
CN109638815B (zh) 一种电力***中长期电压安全稳定预防控制策略确定方法
CN113489015B (zh) 一种基于强化学习的配电网多时间尺度无功电压控制方法
CN102129259A (zh) 用于砂尘环境试验风洞的基于神经网络pi的智能温度控制***和方法
CN108551175B (zh) 配电网储能容量配置方法
CN113300380B (zh) 一种基于负荷曲线分段的配电网无功优化补偿方法
Abaci et al. Hybrid Artificial Neural Network by Using Differential Search Algorithm for Solving Power Flow Problem.
CN107436971A (zh) 适用于非正定型相关性控制的改进拉丁超立方抽样方法
CN107516892A (zh) 基于处理有功优化约束条件提高电能质量的方法
CN115085202A (zh) 电网多区域智能功率协同优化方法、装置、设备及介质
Sun et al. Global state estimation for whole transmission and distribution networks
CN113872213B (zh) 一种配电网电压自主优化控制方法及装置
CN111030089A (zh) 一种基于飞蛾扑火优化算法的pss参数优化方法及***
CN111900767B (zh) 一种用于控制电力***多断面有功潮流的方法及***
CN116995682A (zh) 一种可调节负荷参与有功潮流连续调整方法与***
CN109193662A (zh) 计及不平衡功率分摊的最危险负荷裕度计算方法及***
CN105576653A (zh) 一种220kV片区电网供电能力优化方法

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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210115