CN103177188A - 一种基于特征映射的电力系统负荷动特性分类方法 - Google Patents

一种基于特征映射的电力系统负荷动特性分类方法 Download PDF

Info

Publication number
CN103177188A
CN103177188A CN2013101120504A CN201310112050A CN103177188A CN 103177188 A CN103177188 A CN 103177188A CN 2013101120504 A CN2013101120504 A CN 2013101120504A CN 201310112050 A CN201310112050 A CN 201310112050A CN 103177188 A CN103177188 A CN 103177188A
Authority
CN
China
Prior art keywords
load
data
modeling
mapping
group
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.)
Granted
Application number
CN2013101120504A
Other languages
English (en)
Other versions
CN103177188B (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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN201310112050.4A priority Critical patent/CN103177188B/zh
Publication of CN103177188A publication Critical patent/CN103177188A/zh
Application granted granted Critical
Publication of CN103177188B publication Critical patent/CN103177188B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种基于特征映射的电力系统负荷动特性分类方法,包括下述步骤:1)对扰动数据进行负荷建模;2)确定负荷动特性分类特征量的映射;3)计算建模组数据和训练组数据之间的关联度;4)估计训练组中每条数据的负荷模型参数;5)指标映射的自适应修正;6)负荷动特性分类;该分类方法利用负荷建模平台中已有的电压扰动数据,形成分类特征量映射数据表,以便对负荷动特性进行快速分类,有利于解决负荷建模过程中的时变性问题。

Description

一种基于特征映射的电力系统负荷动特性分类方法
技术领域
本发明属于电力系统负荷建模领域,涉及一种负荷特性分类方法,具体是涉及一种基于特征量映射的电力系统负荷动特性分类方法。
背景技术
负荷模型是决定电力系统仿真结果可信程度的关键因素之一,因此,负荷建模是电力系统重要的基础性课题。目前常用的负荷建模方法主要有总体测辨法和统计综合法两种。总体测辨法以现场大量的实测数据为基础,将负荷群看成一个整体,然后通过系统辨识理论确定负荷模型结构和参数。总体测辨法无需知道各个用户的负荷组成及参数,不依赖于用户统计资料,节约了大量的人力物力,在现阶段得到了广泛地使用。
就总体测辨法负荷建模而言,其最大的困难在于负荷的时变性。负荷时变性的本质在于负荷组成随着时间的不同在不断地变化。对现场实测的每组数据建立一个负荷模型虽然能够解决负荷时变性问题,但是每个变电站下面就会有许多种不同的模型,在系统仿真中是没有实际应用价值的。理想的做法是对实测建模数据首先进行分类,将不同时间采集的负荷扰动数据中负荷组成相近的划分为一类,然后在此基础上对每类数据建立一个通用的负荷模型。分类和综合方法能在一定程度上解决负荷的时变性问题,是负荷建模走向实用化的关键。
特征向量的选取是负荷动特性分类的关键因素之一,特征向量的选择必须能够反映样本的本质特征。现阶段特征向量的选取主要有实测响应(某一电压扰动下负荷实测的功率响应)、模型参数(根据模型的电压扰动输入和功率输出辨识得到的模型参数)、标准电压激励下的模型响应(对辨识得到的模型参数施加同一个电压扰动后得到的功率响应)、实测数据的采集时间(日类型、负荷水平等)。上述这些特征向量的选择都存在不足。由于电压扰动的不同,实测的功率响应很难具有可比性;模型参数和标准电压激励下的模型响应都是以参数的辨识为基础,辨识过程的加入不仅延长了分类时间,而且不可避免地引入了辨识误差;而以采集时间作为特征向量仅仅是一种粗略的分类方法,过多地依靠人为经验。
发明内容
技术问题:本发明提供一种能够有效地解决负荷建模过程中负荷的时变性问题的基于特征映射的电力系统负荷动特性分类方法。
技术方案:本发明采用的基于特征映射的电力系统负荷动特性分类方法,该分类方法包括下述步骤:
1)对扰动数据进行负荷建模,具体流程为:
11)生成用于负荷建模的电压和功率波形:首先利用电能质量监测装置采集各个典型变电站的扰动数据,扰动数据为扰动的电压瞬时值波形和电流瞬时值波形,然后对采集的扰动数据提取基波正序分量,计算三相有功功率、无功功率波形,接着对电压基波正序分量波形和计算得到的三相有功功率、无功功率波形进行降维处理,最后提取降维后的电压基波正序分量波形和三相有功功率、无功功率波形中的扰动波段,作为负荷建模数据;
12)将得到的负荷建模数据随机分为两组,一组为建模组,另一组为训练组,对建模组的负荷建模数据进行负荷建模,将负荷等值为电动机与静态负荷模型并联构成的一个综合负荷模型,采用遗传算法对综合负荷模型进行模型参数辨识,从而得到建模组中每条数据的负荷模型参数;
2)确定负荷动特性分类特征量的映射:负荷动特性分类的特征量包含扰动数据的采集时间、日类型、季节、负荷水平、温度指标,将上述特征量映射到[0,1]的区间内,对每一条扰动数据形成一个对应的特征向量,特征向量由各个特征量映射组成;
3)计算建模组数据和训练组数据之间的关联度:根据每一条扰动数据所对应的特征向量,求解建模组数据与训练组数据之间的关联度,将关联度按照
Figure BDA00003000127100021
进行归一化,得到一个n1×n2的关联度矩阵,其中n1表示建模组数据个数,n2表示训练组数据个数;
4)估计训练组中每条数据的负荷模型参数:根据关联度矩阵和建模组中每条数据的负荷模型参数,计算训练组中每条数据的负荷模型参数;
5)指标映射的自适应修正:根据训练组中每条数据的负荷模型参数,计算训练组中每条数据在对应电压波形的扰动波段下的模型有功、无功输出,然后将模型有功、无功输出同三相功率波形中的扰动波段进行比较,具体方法为:计算它们之间的欧式距离,得到目标函数 J = { 1 n Σ k = 1 n [ ( p ^ ( k ) - p ( k ) ) 2 + ( q ^ ( k ) - q ( k ) ) 2 ] } 1 / 2 , 其中p(k)、q(k)为模型有功、无功输出,
Figure BDA00003000127100032
Figure BDA00003000127100033
为步骤11)中得到的三相功率波形中的有功、无功扰动波段,k=1,2…,n为每条数据的采样点个数;
采用优化算法修正步骤2)中确定的所有负荷动特性分类特征量的映射,使得目标函数J最小,得到各个负荷动特性分类特征量修正后的映射;
6)负荷动特性分类:根据步骤5)得到的各个指标修正后的映射,形成新的特征向量,对新的特征向量采用模糊聚类方法进行负荷动特性分类。
本发明的步骤11)中,对采集的扰动数据提取基波正序分量的方法为:利用电能质量监测装置采集典型变电站扰动波形,对实测的电压电流瞬时值,采用傅里叶分解的方法提取各相基波分量,然后采用对称分量法,求取基波的正序分量,最后计算三相有功功率和无功功率;
提取降维后的电压基波正序分量波形和有功功率、无功功率波形中的扰动波段的方法为:对降维后的电压基波正序分量波形和三相有功功率、无功功率波形,首先依照式(1)进行异常数据波段的检测,如果式(1)成立,则以第n个采样点作为起始点,否则由用户自己选择起始点;然后再依照式(2)进行截止点的检测,如果式(2)成立,则以第m个采样点作为截止点,否则由用户自己选择截止点;
|y(n+i)-y(n)|>ε i=1,2        式(1)
|y(m)-y(m-i)|>ε i=1,2,…,5         式(2)
式中y为降维后的电压基波正序分量波形和三相有功功率、无功功率波形,n=1,2…,N,m=N,N-1…,1,其中N为采样点数;ε为根据扰动的大小设定的阀值。
本发明的步骤2)的具体方法为:将负荷动特性分类的特征量分为两类,一类为定量指标,包括温度指标和负荷水平指标,负荷水平指标包括扰动发生前稳态时刻的电压U0、有功P0、无功Q0,另一类为定性指标,包括数据采集时间、日类型、季节,数据采集时间分为工作时间段和休息时间段,日类型包括周一到周日;
对于定量指标,采用分段映射的方法确定其映射,对于分段点之间的值,采用线性插值的方法确定其映射;对于定性指标,将其具体指标映射为[0,1]之间的数值。
本发明的步骤3)中,采用灰色关联度的方法求解建模组数据与训练组数据之间的关联度。
本发明的步骤4)的具体方法为:以步骤3)求得的建模组数据与训练组数据之间的关联度作为权重,利用步骤12)得到的建模组中每条数据的负荷模型参数,通过加权平均的方法得到训练组中每条数据的负荷模型参数
Figure BDA00003000127100041
其中θ为建模组数据负荷模型参数,θ'为训练组数据负荷模型参数,i=1,2…,n1为建模组的数据个数,j=1,2…,n2为训练组的数据个数,t为模型参数个数。
本发明的步骤5)中使用遗传优化算法,修正步骤2)中确定的所有负荷动特性分类特征量的映射。
本发明的步骤6)采用类平均距离法进行负荷动特性分类。
有益效果:与现有技术相比,本发明具有以下优点:
(1)确定了负荷动特性分类特征量的映射。负荷动特性分类特征量映射的确定有利于负荷动特性分类问题的更简便解决,随着变电站扰动数据的不断积累,各个分类特征量的映射将不断完善,可以为将来新增数据的分类提供依据。同时,由于不同变电站所带负荷类型不同,影响负荷特性的因素也不尽相同,不同变电站可以根据实际情况灵活选择分类特征量。
(2)有利于提高负荷动特性分类的效率。在本分类方法中,特征向量是根据分类特征量映射得到的,不需要进行参数辨识,不但避免了模型结构选择和模型参数辨识所带来的误差,而且大大减少了特征向量的计算时间,提高了分类的效率。这样,用户可以对负荷特性进行实时分类,将极大地推动在线负荷建模,特别是负荷综合建模的发展。
本发明方法把不同量纲的值映射到一个特定的区间,通过自适应训练算法修正这些映射,从而使得各个量之间具有数值上的可比性,然后以特征量之间的相似度作为分类的标准,将负荷组成相近的划分为一类,能够有效地解决负荷建模过程中负荷的时变性问题。
附图说明
图1是本发明的流程框图。
图2是本发明中步骤12)中所用的负荷模型结构。
图3是本发明中步骤12)和步骤5)中所用的遗传算法流程图。
具体实施方式
下面结合附图,对本发明的技术方案进行具体介绍。
如图1所示,本发明提供了一种基于特征映射的电力系统负荷动特性分类方法。负荷动特性分类指的是将同一综合负荷不同时间采集的动态负荷扰动数据中负荷组成相近的划分为一类,以便为每一类数据建立一个综合负荷模型。该分类方法包括下述步骤:
1)对扰动数据进行负荷建模,具体流程为:
11)生成用于负荷建模的电压和功率波形。
该步骤11)通过以下步骤111)、112)和113)实现。
111)对实测的录波波形提取基波正序分量:利用电能质量监测装置采集典型变电站扰动波形,对实测的电压电流瞬时值,采用傅里叶分解的方法提取各相基波分量。然后采用对称分量法,求取基波的正序分量。最后计算三相有功无功功率。
112)数据降维、平滑:对步骤111)得到的基波的正序分量进行降维,以减少计算量。例如,原先的采样频率为1kHz,1秒钟数据1000个,将采样频率降低为100Hz,1秒钟数据100个。然后采用五点二次平滑方法滤去信号中的高频噪声,具体方法如式(1)所示:
y(n)=[-3x(n-2)+12x(n-1)+17x(n)+12x(n+1)-3x(n+2)]/35]  式(1)
式中y(n)为平滑后的信号,x(n)为平滑前信号,n=2,…,N-3为时间采样点,N为采样点数。
113)提取扰动波段:对步骤112)得到的降维、平滑后的电压基波正序分量波形和三相有功功率、无功功率波形,首先依照式(1)进行异常数据波段的检测,如果式(1)成立,则以第n个采样点作为起始点,否则由用户自己选择起始点;然后再依照式(2)进行截止点的检测,如果式(2)成立,则以第m个采样点作为截止点,否则由用户自己选择截止点;
|y(n+i)-y(n)|>ε i=1,2      式(2)
|y(m)-y(m-i)|>ε i=1,2,…,5     式(3)
式中y为降维后的电压基波正序分量波形和三相有功功率、无功功率波形,n=1,2…,N,m=N,N-1…,1,其中N为采样点数;ε为根据扰动的大小设定的阀值;
12)将步骤11)中得到的负荷建模数据随机分为两组,一组为建模组,另一组为训练组,对建模组的负荷建模数据进行负荷建模,将负荷等值为电动机与静态负荷模型并联构成的一个综合负荷模型,采用遗传算法对综合负荷模型进行模型参数辨识,从而得到建模组中每条数据的负荷模型参数;
该步骤12)中采用三阶感应电动机并联ZIP(恒阻抗、恒电流、恒功率)静态模型的综合负荷模型,其结构如图2所示,状态方程为:
de d dt = 1 T d [ - e d + ( x - x ′ ) I q + 314 e q * s * T d ]
de q dt = 1 T d [ - e q + ( x - x ′ ) I d + 314 e d * s * T d ]    式(4)
ds dt = 1 T j { K L [ a + ( 1 - a ) ( 1 - s ) b ] - ( e d I d + e q I q ) }
I d = 1 r 1 2 + x ′ 2 ( r 1 * ( e d - u d ) + x ′ * ( e q - u q ) )    式(5)
I q = 1 r 1 2 + x ′ 2 ( r 1 * ( e q - u q ) + x ′ * ( e d - u d ) )
输出方程为:
PL=-UdId+PS0[PZ(u/u0)2+PI(u/u0)+Pp]        式(6)
QL=UdIq+QS0[QZ(u/u0)2+QI(u/u0)+Qp]
其中xm=ωs*r2*Td-x2,为定转子互感抗;为转子暂态电抗;x=x1+xm,为转子稳态电抗。
动态负荷模型的状态向量为[Ed Eq s],输入向量为[U],输出向量为[P Q],待辨识的参数为[r1 x1 Td x2 s0 r2 Tj a b PZ Pp QZ Qp Pm]。试图寻找一组最优的负荷模型参数,使目标函数的值最小。目标函数定义为
W ( k ) = 0.4 | U ( k ) - U ‾ | Σ k = 1 n | U ( k ) - U ‾ | + 0.6 n    式(6)
J = { 1 n Σ k = 1 n W ( k ) [ ( p ^ ( k ) - p ( k ) ) 2 + ( q ^ ( k ) - q ( k ) ) 2 ] } 1 / 2
式中:J为目标函数,
Figure BDA00003000127100073
Figure BDA00003000127100074
分别为负荷模型计算得到的第k个采样点有功功率、无功功率,p(k)、q(k)分别为步骤113)确定的三相有功、无功扰动波段中的第k个采样点的功率值,W(k)为误差权重。
该步骤12)采用遗传算法进行负荷模型参数的辨识,具体步骤是:
(1)设定待辨识参数的取值区间,选择参数编码策略,随机初始化生成群体P;
(2)定义适应度函数f(x),计算个体的适应度;
(3)进行选择操作,选择适应度较高的个体复制到下一代。
(4)进行交叉操作,采用两点交叉的方法,使两个个体相互交换部分基因,以扩大解的空间。
(5)进行变异操作,在个体中引入新的基因,防止个体陷入局部最优解,最终形成下一代群体。
(6)判断新的群体是否满足某一指标,或者达到最大迭代次数,若不满足则返回步骤(3)。
2)确定负荷动特性分类特征量的映射。该步骤与步骤1)同步进行。
负荷动特性分类的特征量包含扰动数据的采集时间、日类型、季节、负荷水平、温度指标,将上述特征量映射到[0,1]的区间内,对每一条扰动数据形成一个对应的特征向量,特征向量由各个特征量映射组成;
对于定量指标,采用分段映射的方法,这里分段点可以灵活选取,可以分散程度较大,也可以比较密集。对于分段点之间的值,可以采用线性插值的方法。例如各个变电站扰动发生前稳态时刻的有功在1.2~0.8(标幺值)之间,其中在0.9~1.1的分布尤为密集,那么可以设定0.8、0.9、1.1、1.2四个阈值,在0.8~0.9、0.9~1.1、1.1~1.2分别采用一组线性映射;对于定性指标,将其具体指标映射为[0,1]之间的数值。例如将工作时间段映射为0.2,将休息时间段映射为0.6,以体现不同时间段的区别。
3)计算建模组数据和训练组数据之间的关联度:根据步骤2)得到的每一条扰动数据所对应的特征向量,求解建模组数据与训练组数据之间的关联度,将关联度按照
Figure BDA00003000127100081
进行归一化,得到一个n1×n2的关联度矩阵,其中n1表示建模组数据个数,n2表示训练组数据个数;
该步骤3)中,采用灰色关联度法形成关联度矩阵,具体过程是:
(1)建立数据矩阵。选择任意一条训练组数据和全部的n1条建模组数据,将他们的p个特征向量形成数据矩阵X。
X = x 11 x 12 · · · x 1 p x 21 x 22 · · · x 2 p · · · · · · · · · · · · x n 1 x n 2 · · · x np
矩阵中第1行为训练组数据的特征向量,第2至n行为建模组数据的特征向量。
(2)求特征向量差值矩阵。
Δij=|xi+1,j-x1j|
其中i=1,2,...,n1,j=1,2,...,p,求得的差值矩阵为:
Δ = Δ 11 Δ 12 · · · Δ 1 p Δ 21 Δ 22 · · · Δ 2 p · · · · · · · · · · · · Δ n 1 1 Δ n 1 2 · · · Δ n 1 p
(3)求特征量之间的关联系数。
ξ i ( j ) = Δ min + ρΔ max Δ ij + ρΔ max
ρ为分辨率系数,一般取ρ为0.5。
(4)求训练数据与建模数据之间的关联度
r i = 1 n 1 Σ j = 1 p ξ i ( j ) , i = 1,2 , . . . , n 1
(5)依次选择全部的n2条训练组数据,计算它们与建模组数据的关联度,形成关联度矩阵Y。
Y = r 11 r 12 · · · r 1 n 2 r 21 r 22 · · · r 2 n 2 · · · · · · · · · · · · r n 1 1 r n 1 2 · · · r n 1 n 2
(6)将关联度按照进行归一化,其中n1为建模数据的个数,得到归一化后的关联度矩阵。
4)估计训练组中每条数据的负荷模型参数:根据步骤12)辨识得到的建模组中每条数据的负荷模型参数和步骤3)得到的关联度矩阵,计算训练组中每条数据的负荷模型参数;
该步骤4)中,训练组数据负荷模型参数的估计方法是:以步骤3)求得的建模组数据与训练组数据之间的关联度作为权重,利用步骤12)得到的建模组中每条数据的负荷模型参数,通过加权平均的方法得到训练组中每条数据的负荷模型参数其中θ为建模组数据负荷模型参数,θ'为训练组数据负荷模型参数,i=1,2…,n1为建模组的数据个数,j=1,2…,n2为训练组的数据个数,t为模型参数个数。
5)指标映射的自适应修正:根据步骤4)得到的训练组中每条数据的负荷模型参数,计算训练组中每条数据在对应电压波形的扰动波段下的模型有功、无功输出,然后将模型有功、无功输出同步骤11)中得到的三相功率波形中的扰动波段进行比较,具体方法为:计算它们之间的欧式距离,得到目标函数 J = { 1 n Σ k = 1 n [ ( p ^ ( k ) - p ( k ) ) 2 + ( q ^ ( k ) - q ( k ) ) 2 ] } 1 / 2 , 其中p(k)、q(k)为模型有功、无功输出,
Figure BDA00003000127100095
Figure BDA00003000127100096
为步骤11)中得到的三相功率波形中的有功、无功扰动波段,k=1,2…,n为每条数据的采样点个数;
采用遗传优化算法修正步骤2)中确定的所有负荷动特性分类特征量的映射,使得目标函数J最小,遗传算法的流程与步骤12)相同。最后得到各个负荷动特性分类特征量修正后的映射;
6)负荷动特性分类:根据步骤5)得到的各个指标修正后的映射,形成新的特征向量,对新的特征向量采用模糊聚类方法进行负荷动特性分类。
该步骤中负荷动特性分类采用类平均距离法,它是系统聚类法的一种。系统聚类法的基本思路为:设有N个样本,初始这N个样本自成一类,然后计算样品之间的距离,将距离最小的类并为一新类,再计算并类后的新类与其它类的距离,又将距离最小的两类合并为一新类,这样每次减少一些类,直到将这N个样本合并成一类为止。类平均距离法中类与类之间的距离定义为其中np为第p类中的样本个数;nq为第q类中的样本个数;dij表示p类中的样本与q类中的样本之间的距离,这里采用欧氏距离。

Claims (7)

1.一种基于特征映射的电力系统负荷动特性分类方法,其特征在于,该分类方法包括下述步骤:
1)对扰动数据进行负荷建模,具体流程为:
11)生成用于负荷建模的电压和功率波形:首先利用电能质量监测装置采集各个典型变电站的扰动数据,所述扰动数据为扰动的电压瞬时值波形和电流瞬时值波形,然后对采集的扰动数据提取基波正序分量,计算三相有功功率、无功功率波形,接着对电压基波正序分量波形和计算得到的三相有功功率、无功功率波形进行降维处理,最后提取降维后的电压基波正序分量波形和三相有功功率、无功功率波形中的扰动波段,作为负荷建模数据;
12)将所述步骤11)中得到的负荷建模数据随机分为两组,一组为建模组,另一组为训练组,对所述建模组的负荷建模数据进行负荷建模,将负荷等值为电动机与静态负荷模型并联构成的一个综合负荷模型,采用遗传算法对所述综合负荷模型进行模型参数辨识,从而得到建模组中每条数据的负荷模型参数;
2)确定负荷动特性分类特征量的映射:负荷动特性分类的特征量包含扰动数据的采集时间、日类型、季节、负荷水平、温度指标,将上述特征量映射到[0,1]的区间内,对每一条扰动数据形成一个对应的特征向量,所述特征向量由各个特征量映射组成;
3)计算建模组数据和训练组数据之间的关联度:根据所述步骤2)得到的每一条扰动数据所对应的特征向量,求解建模组数据与训练组数据之间的关联度,将所述关联度按照
Figure FDA00003000127000011
进行归一化,得到一个n1×n2的关联度矩阵,其中n1表示建模组数据个数,n2表示训练组数据个数;
4)估计训练组中每条数据的负荷模型参数:根据步骤12)辨识得到的建模组中每条数据的负荷模型参数和步骤3)得到的关联度矩阵,计算训练组中每条数据的负荷模型参数;
5)指标映射的自适应修正:根据所述步骤4)得到的训练组中每条数据的负荷模型参数,计算训练组中每条数据在对应电压波形的扰动波段下的模型有功、无功输出,然后将所述模型有功、无功输出同所述步骤11)中得到的三相功率波形中的扰动波段进行比较,具体方法为:计算它们之间的欧式距离,得到目标函数 J = { 1 n Σ k = 1 n [ ( p ^ ( k ) - p ( k ) ) 2 + ( q ^ ( k ) - q ( k ) ) 2 ] } 1 / 2 , 其中p(k)、q(k)为模型有功、无功输出,
Figure FDA00003000127000022
为步骤11)中得到的三相功率波形中的有功、无功扰动波段,k=1,2…,n为每条数据的采样点个数;
采用优化算法修正步骤2)中确定的所有负荷动特性分类特征量的映射,使得目标函数J最小,得到各个负荷动特性分类特征量修正后的映射;
6)负荷动特性分类:根据步骤5)得到的各个指标修正后的映射,形成新的特征向量,对所述新的特征向量采用模糊聚类方法进行负荷动特性分类。
2.根据权利要求1所述的基于特征映射的电力系统负荷动特性分类方法,其特征在于,所述的步骤11)中,对采集的扰动数据提取基波正序分量的方法为:利用电能质量监测装置采集典型变电站扰动波形,对实测的电压电流瞬时值,采用傅里叶分解的方法提取各相基波分量,然后采用对称分量法,求取基波的正序分量,最后计算三相有功功率和无功功率;
提取降维后的电压基波正序分量波形和有功功率、无功功率波形中的扰动波段的方法为:对降维后的电压基波正序分量波形和三相有功功率、无功功率波形,首先依照式(1)进行异常数据波段的检测,如果式(1)成立,则以第n个采样点作为起始点,否则由用户自己选择起始点;然后再依照式(2)进行截止点的检测,如果式(2)成立,则以第m个采样点作为截止点,否则由用户自己选择截止点;
|y(n+i)-y(n)|>ε i=1,2          式(1)
|y(m)-y(m-i)|>ε i=1,2,…,5           式(2)
式中y为降维后的电压基波正序分量波形和三相有功功率、无功功率波形,n=1,2…,N,m=N,N-1…,1其中N为采样点数;ε为根据扰动的大小设定的阀值。
3.根据权利要求1所述的基于特征映射的电力系统负荷动特性分类方法,其特征在于,所述步骤2)的具体方法为:将负荷动特性分类的特征量分为两类,一类为定量指标,包括温度指标和负荷水平指标,所述负荷水平指标包括扰动发生前稳态时刻的电压U0、有功P0、无功Q0,另一类为定性指标,包括数据采集时间、日类型、季节,所述数据采集时间分为工作时间段和休息时间段,所述日类型包括周一到周日;
对于所述定量指标,采用分段映射的方法确定其映射,对于分段点之间的值,采用线性插值的方法确定其映射;对于定性指标,将其具体指标映射为[0,1]之间的数值。
4.根据权利要求1所述的基于特征映射的电力系统负荷动特性分类方法,其特征在于,所述的步骤3)中,采用灰色关联度的方法求解建模组数据与训练组数据之间的关联度。
5.根据权利要求1所述的基于特征映射的电力系统负荷动特性分类方法,其特征在于,所述的步骤4)的具体方法为:以步骤3)求得的建模组数据与训练组数据之间的关联度作为权重,利用步骤12)得到的建模组中每条数据的负荷模型参数,通过加权平均的方法得到训练组中每条数据的负荷模型参数
Figure FDA00003000127000031
其中θ为建模组数据负荷模型参数,θ'为训练组数据负荷模型参数,i=1,2…,n1为建模组的数据个数,j=1,2…,n2为训练组的数据个数,t为模型参数个数。
6.根据权利要求1所述的基于特征映射的电力系统负荷动特性分类方法,其特征在于,所述的步骤5)中使用遗传优化算法,修正步骤2)中确定的所有负荷动特性分类特征量的映射。
7.根据权利要求1所述的基于特征映射的电力系统负荷动特性分类方法,其特征在于,所述的步骤6)采用类平均距离法进行负荷动特性分类。
CN201310112050.4A 2013-04-02 2013-04-02 一种基于特征映射的电力系统负荷动特性分类方法 Expired - Fee Related CN103177188B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310112050.4A CN103177188B (zh) 2013-04-02 2013-04-02 一种基于特征映射的电力系统负荷动特性分类方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310112050.4A CN103177188B (zh) 2013-04-02 2013-04-02 一种基于特征映射的电力系统负荷动特性分类方法

Publications (2)

Publication Number Publication Date
CN103177188A true CN103177188A (zh) 2013-06-26
CN103177188B CN103177188B (zh) 2015-12-09

Family

ID=48637038

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310112050.4A Expired - Fee Related CN103177188B (zh) 2013-04-02 2013-04-02 一种基于特征映射的电力系统负荷动特性分类方法

Country Status (1)

Country Link
CN (1) CN103177188B (zh)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103617371A (zh) * 2013-12-10 2014-03-05 国家电网公司 一种基于灰色理论的电能质量综合评估方法
CN103915841A (zh) * 2014-04-16 2014-07-09 华北电力大学 一种用于电力系统负荷特性仿真的建模方法
CN104052058A (zh) * 2014-06-13 2014-09-17 华北电力大学 一种基于马尔可夫链蒙特卡罗法的系统谐波概率评估方法
CN104200032A (zh) * 2014-09-05 2014-12-10 山东大学 广义负荷建模中基于时段性的横向时间轴聚类方法
CN104376402A (zh) * 2014-10-31 2015-02-25 国家电网公司 一种基于频域指标的负荷特性分类与综合方法
CN104809337A (zh) * 2015-04-15 2015-07-29 华南理工大学 一种变电站主要站用负荷用电实用模型建模方法
CN105391062A (zh) * 2015-12-07 2016-03-09 国网浙江省电力公司宁波供电公司 一种基于直流潮流模型的有功不良数据辨识方法
CN107958019A (zh) * 2017-10-23 2018-04-24 中国南方电网有限责任公司 基于电网运行管理系统的图形化应用方法
CN109213028A (zh) * 2018-07-25 2019-01-15 华南理工大学 一种安装简便的智慧用电系统及方法
CN110163304A (zh) * 2019-06-14 2019-08-23 福州大学 一种采用线性关系聚类的谐波源耦合参数辨识方法
CN110571806A (zh) * 2019-09-25 2019-12-13 北方民族大学 一种用于配电网负荷类别的特征提取与辨识方法
CN110866840A (zh) * 2019-11-12 2020-03-06 南昌工程学院 基于知识图谱的电力负荷特征量训练的数据库建模方法
CN112001550A (zh) * 2020-08-25 2020-11-27 杭州电力设备制造有限公司 一种变电站多点电能质量监测系统
CN112288598A (zh) * 2020-12-24 2021-01-29 中国电力科学研究院有限公司 一种确定变电站的负荷元件的构成的方法和系统
CN113627041A (zh) * 2021-10-12 2021-11-09 国网江西省电力有限公司电力科学研究院 基于多源数据融合的时变动态负荷在线建模方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0784513A (ja) * 1993-09-14 1995-03-31 Toshiba Corp 電力系統運用訓練シミュレータ
CN101625381A (zh) * 2009-08-06 2010-01-13 江苏省电力试验研究院有限公司 电压暂降模糊评估方法
CN102646155A (zh) * 2011-10-17 2012-08-22 中国电力科学研究院 一种适用于短路电流计算的负荷建模方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0784513A (ja) * 1993-09-14 1995-03-31 Toshiba Corp 電力系統運用訓練シミュレータ
CN101625381A (zh) * 2009-08-06 2010-01-13 江苏省电力试验研究院有限公司 电压暂降模糊评估方法
CN102646155A (zh) * 2011-10-17 2012-08-22 中国电力科学研究院 一种适用于短路电流计算的负荷建模方法

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103617371A (zh) * 2013-12-10 2014-03-05 国家电网公司 一种基于灰色理论的电能质量综合评估方法
CN103915841A (zh) * 2014-04-16 2014-07-09 华北电力大学 一种用于电力系统负荷特性仿真的建模方法
CN103915841B (zh) * 2014-04-16 2015-12-09 华北电力大学 一种用于电力系统负荷特性仿真的建模方法
CN104052058A (zh) * 2014-06-13 2014-09-17 华北电力大学 一种基于马尔可夫链蒙特卡罗法的系统谐波概率评估方法
CN104200032B (zh) * 2014-09-05 2017-04-05 山东大学 广义负荷建模中基于时段性的横向时间轴聚类方法
CN104200032A (zh) * 2014-09-05 2014-12-10 山东大学 广义负荷建模中基于时段性的横向时间轴聚类方法
CN104376402A (zh) * 2014-10-31 2015-02-25 国家电网公司 一种基于频域指标的负荷特性分类与综合方法
CN104376402B (zh) * 2014-10-31 2017-06-13 国家电网公司 一种基于频域指标的负荷特性分类与综合方法
CN104809337B (zh) * 2015-04-15 2018-01-05 华南理工大学 一种变电站主要站用负荷用电实用模型建模方法
CN104809337A (zh) * 2015-04-15 2015-07-29 华南理工大学 一种变电站主要站用负荷用电实用模型建模方法
CN105391062A (zh) * 2015-12-07 2016-03-09 国网浙江省电力公司宁波供电公司 一种基于直流潮流模型的有功不良数据辨识方法
CN107958019A (zh) * 2017-10-23 2018-04-24 中国南方电网有限责任公司 基于电网运行管理系统的图形化应用方法
CN107958019B (zh) * 2017-10-23 2022-02-25 中国南方电网有限责任公司 基于电网运行管理系统的图形化应用方法
CN109213028A (zh) * 2018-07-25 2019-01-15 华南理工大学 一种安装简便的智慧用电系统及方法
CN110163304A (zh) * 2019-06-14 2019-08-23 福州大学 一种采用线性关系聚类的谐波源耦合参数辨识方法
CN110571806A (zh) * 2019-09-25 2019-12-13 北方民族大学 一种用于配电网负荷类别的特征提取与辨识方法
CN110571806B (zh) * 2019-09-25 2023-05-12 北方民族大学 一种用于配电网负荷类别的特征提取与辨识方法
CN110866840A (zh) * 2019-11-12 2020-03-06 南昌工程学院 基于知识图谱的电力负荷特征量训练的数据库建模方法
CN112001550A (zh) * 2020-08-25 2020-11-27 杭州电力设备制造有限公司 一种变电站多点电能质量监测系统
CN112288598A (zh) * 2020-12-24 2021-01-29 中国电力科学研究院有限公司 一种确定变电站的负荷元件的构成的方法和系统
CN113627041A (zh) * 2021-10-12 2021-11-09 国网江西省电力有限公司电力科学研究院 基于多源数据融合的时变动态负荷在线建模方法及装置

Also Published As

Publication number Publication date
CN103177188B (zh) 2015-12-09

Similar Documents

Publication Publication Date Title
CN103177188B (zh) 一种基于特征映射的电力系统负荷动特性分类方法
CN106650767B (zh) 基于聚类分析和实时校正的洪水预报方法
CN102074955B (zh) 基于知识发现技术的电力系统稳定评估及控制方法
Ma et al. Ultra-short-term wind generation forecast based on multivariate empirical dynamic modeling
Khoshkhoo et al. On-line dynamic voltage instability prediction based on decision tree supported by a wide-area measurement system
CN105911476B (zh) 一种基于数据挖掘的电池储能系统soc预测方法
CN104732300B (zh) 一种基于模糊分区理论的神经网络风功率短期预测方法
CN103323755A (zh) Gis特高频局部放电信号识别方法及系统
CN110070282A (zh) 一种基于综合关联度的低压台区线损影响因素分析方法
CN105160437A (zh) 基于极限学习机的负荷模型预测方法
CN111525587B (zh) 一种基于无功负荷态势的电网无功电压控制方法及系统
CN103915841B (zh) 一种用于电力系统负荷特性仿真的建模方法
CN110659693A (zh) 基于k近邻分类的配电网快速拓扑识别方法、系统及可读存储介质
CN107565553A (zh) 一种基于ukf的配电网抗差动态状态估计方法
CN109871622A (zh) 一种基于深度学习的低压台区线损计算方法及系统
CN103279803A (zh) 基于综合信息论及现代内点理论的负荷建模方法和系统
CN109066819B (zh) 一种基于案件推理的配电网的无功优化方法
CN109066651B (zh) 风电-负荷场景的极限传输功率的计算方法
CN110738232A (zh) 一种基于数据挖掘技术的电网电压越限成因诊断方法
CN103324841A (zh) 基于故障录波器数据的在线动态负荷建模方法
CN105406461A (zh) 一种配电网停电事件的负荷自适应动态监测方法
CN103793556B (zh) 综合出线负荷的变电站负荷建模方法
CN112636328B (zh) 一种中压配电网设备负荷状态辨识方法
CN109061354A (zh) 一种基于改进pso和svm的电能质量扰动识别方法
CN106372440A (zh) 一种并行计算的配电网自适应抗差状态估计方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
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: 20151209

Termination date: 20200402