CN117874497A - 一种负荷特征提取方法、装置、设备及介质 - Google Patents
一种负荷特征提取方法、装置、设备及介质 Download PDFInfo
- Publication number
- CN117874497A CN117874497A CN202410273266.7A CN202410273266A CN117874497A CN 117874497 A CN117874497 A CN 117874497A CN 202410273266 A CN202410273266 A CN 202410273266A CN 117874497 A CN117874497 A CN 117874497A
- Authority
- CN
- China
- Prior art keywords
- load
- frequency
- data
- clustering
- processing
- 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
Links
- 238000000605 extraction Methods 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 40
- 238000012545 processing Methods 0.000 claims abstract description 35
- 238000010586 diagram Methods 0.000 claims abstract description 34
- 238000012216 screening Methods 0.000 claims abstract description 19
- 238000007637 random forest analysis Methods 0.000 claims abstract description 11
- 230000009466 transformation Effects 0.000 claims abstract description 5
- 238000004590 computer program Methods 0.000 claims description 17
- 238000010606 normalization Methods 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 230000002776 aggregation Effects 0.000 claims description 4
- 238000004220 aggregation Methods 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000007405 data analysis Methods 0.000 claims description 3
- 239000000470 constituent Substances 0.000 claims description 2
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 238000009826 distribution Methods 0.000 abstract description 5
- 238000005457 optimization Methods 0.000 abstract description 2
- 238000004458 analytical method Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 238000010248 power generation Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 230000014759 maintenance of location Effects 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/24323—Tree-organised classifiers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
- G06N20/20—Ensemble learning
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Medical Informatics (AREA)
- Image Analysis (AREA)
Abstract
本发明属于配电网负荷优化技术领域,具体公开了一种负荷特征提取方法、装置、设备及介质。包括:获取预设区域内的负荷数据;用快速傅里叶变换生成频率‑幅值图;由频率‑幅值图获取主要构成频率,由主要构成频率和频率‑幅值图生成候选特征集;采用层次聚类法处理候选特征集得到第一聚类结果;采用随机森林分类器处理第一聚类结果得到分类结果,从而设重要限度值,由重要限度值筛选得到特征数据集;采用层次聚类法处理特征数据集得到第二聚类结果;比较第一聚类结果与第二聚类结果,吻合度≥预设值输出特征数据集;小于预设值则调整重要限度值。通过快速傅里叶变换将负荷数据转换至频域,通过不同频率下的幅值和相位特征,实现负荷特征量化描述。
Description
技术领域
本发明属于配电网负荷优化技术领域,具体涉及一种负荷特征提取方法、装置、设备及介质。
背景技术
随着新型电力系统建设,产生新型发用电负荷特征。在发电侧,随着分布式新能源发电渗透率不断提高,电源由确定可控的连续电源向不确定的随机波动电源转变;在用电侧,从机械电磁系统和高转动惯量系统向电力电子器件及弱转动惯量系统转变,对配电网规划和运行形成挑战。因此,针对新型电力负荷特征,有必要开发精准可靠的负荷特征提取和分析方法,熟悉电网负荷发用电特点,进行有效的分类,对负荷预测、需求响应管理、配网规划投资等意义重大。
目前的负荷特征分析,多局限于峰谷值、平均值、增长率等统计学指数,无法满足新型电力系统条件下,负荷预测、需求响应管理、配网规划等对高精细度负荷特征分析需求。同时,对于负荷的分类,也多局限于营销系统中对用户按照居民工商业等不同行业进行划分,无法体现不同规模、不同电压等级、不同区域、不同用户不同时间等因素对用电负荷的影响。
因此,亟需构建负荷特征分析的系统性框架,能够对新型电力负荷发用电特征进行精细化提取和分析。
发明内容
本发明的目的在于提供一种负荷特征提取方法、装置、设备及介质,以解决现有负荷特征提取多针对统计学指数,难以量化,导致提取结果准确的技术问题。
为实现上述目的,本发明采用如下技术方案予以实现:
第一方面,本发明提供一种负荷特征提取方法,包括以下步骤:
S1、获取预设区域内的若干组负荷数据;
S2、采用快速傅里叶变换处理若干组负荷数据生成频率-幅值图;
S3、根据频率-幅值图获取各组负荷数据的主要构成频率,根据主要构成频率和频率-幅值图生成候选特征集;
S4、采用层次聚类法处理候选特征集得到第一聚类结果;
S5、采用随机森林分类器处理第一聚类结果得到分类结果,根据分类结果设置重要限度值,根据重要限度值筛选分类结果得到特征数据集;
S6、采用层次聚类法处理特征数据集得到第二聚类结果;
S7、将第一聚类结果与第二聚类结果进行比较,若吻合度大于等于预设值则输出特征数据集;
若吻合度小于预设值则返回S5,调整重要限度值。
本发明的进一步改进在于:所述获取预设区域内的若干组负荷数据的步骤中,还包括对负荷数据进行预处理:
所述预处理包括:
对负荷数据进行数据治理,删除无效数据,统一负荷数据的时间跨度和时间颗粒度;
对负荷数据进行最大值归一化。
本发明的进一步改进在于:所述快速傅里叶变换的公式:
Xk= FFT(xn(t));
式中,Xk= (ak, θk)为快速傅里叶变换后的复数结果,k为频率,ak表示k频率下的振幅,θk表示k频率下的相角。
本发明的进一步改进在于:所述候选特征集中还包括预设频率的幅值和相角。
本发明的进一步改进在于:所述调整重要限度值的步骤中,重要限度值逐渐减小。
本发明的进一步改进在于:所述预设值为80%。
本发明的进一步改进在于:所述主要构成频率为频率-幅值图中的各峰值。
第二方面,本发明提供一种负荷特征提取装置,包括:
数据获取模块:用于获取预设区域内的若干组负荷数据;
数据处理模块:采用快速傅里叶变换处理若干组负荷数据生成频率-幅值图;
数据分析模块:用于根据频率-幅值图获取各组负荷数据的主要构成频率,根据主要构成频率和频率-幅值图生成候选特征集;
第一聚类模块:采用层次聚类法处理候选特征集得到第一聚类结果;
筛选模块:用于采用随机森林分类器处理第一聚类结果得到分类结果,根据分类结果设置重要限度值,根据重要限度值筛选分类结果得到特征数据集;
第二聚类模块:采用层次聚类法处理特征数据集得到第二聚类结果;
迭代输出模块:用于将第一聚类结果与第二聚类结果进行比较,若吻合度大于等于预设值则输出特征数据集;
若吻合度小于预设值则返回筛选模块,调整重要限度值。
第三方面,本发明提供一种计算机设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述的一种负荷特征提取方法。
第四方面,本发明提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现上述的一种负荷特征提取方法。
与现有技术相比,本发明至少包括以下有益效果:
1、本发明通过快速傅里叶变换将负荷数据转换至频域,通过不同频率下的幅值和相位特征,实现负荷特征量化描述;
2、本发明通过随机森林分类器对特征进行筛选,并将筛选前后的数据进行聚类比较,从而剔除无效数据,提高负荷特征提取的准确性;
3、本发明通过设置重要限度值进行特征筛选,并在后续迭代过程中对重要限度值进行调整,避免数据过度删减导致特征提取不准确。
附图说明
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为本发明一种负荷特征提取方法的流程图;
图2为本发明一种负荷特征提取装置的结构框图;
图3为本发明一种负荷特征提取方法中数据归一化后的时间-归一化负荷图;
图4为本发明一种负荷特征提取方法中频率-幅值图;
图5为本发明一种负荷特征提取方法中第一次聚类结果图;
图6为本发明一种负荷特征提取方法中随机森林分类器分类结果图;
图7为本发明一种负荷特征提取方法中第二次聚类结果图。
具体实施方式
下面将参考附图并结合实施例来详细说明本发明。需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。
以下详细说明均是示例性的说明,旨在对本发明提供进一步的详细说明。除非另有指明,本发明所采用的所有技术术语与本发明所属领域的一般技术人员的通常理解的含义相同。本发明所使用的术语仅是为了描述具体实施方式,而并非意图限制根据本发明的示例性实施方式。
实施例1
一种负荷特征提取方法,如图1所示,包括以下步骤:
S1、获取预设区域内的若干组负荷数据;
具体的,每组负荷数据可以是不同地区、不同行业、不同电压等级或不同用户;
作为本实施例的一种优选实施方式,获取预设区域内的负荷数据的步骤中,负荷数据的采集周期为5分钟一次数据采集或15分钟一次数据采集,采集周期不宜超过1小时一次,负荷数据的时间跨度大于等于1年,从而保证获取的负荷数据准确性和全面性,提高最终负荷提取的准确性。
具体的,在S1中还包括:
S11、对负荷数据进行数据治理,删除无效数据,统一负荷数据的时间跨度和时间颗粒度;
S12、对负荷数据进行最大值归一化;
Xn(t)=xm(t)/xpeak;
其中xm(t)为负荷数据中的值,xpeak为该组负荷数据中的最大值,xn(t)为归一化结果。
具体的,最大值归一化是指将找到每组负荷数据中的最大值,并每个负荷数据与本组最大负荷数据相除,得到最大值均为1的负荷数据;
根据归一化后的数据可以得到归一化结果和时间折线图,如图4所示;
通过两步预处理,有效清除无效数据,提高数据准确性,简化后续计算,提高准确性。
S2、采用快速傅里叶变换处理若干组负荷数据生成频率-幅值图;
具体的,快速傅里叶变换处理根据公式:
Xk=FFT(xn(t));
式中,Xk=(ak,θk)为快速傅里叶变换后的复数结果,k为频率,ak表示k频率下的振幅,θk表示k频率下的相角。
具体的,频率-幅值图也可称为周期-能量图,如图3所示,纵坐标为不同周期成分的能量,图3是由若干不同周期、不同强度的正弦波构成,某周期(频率)对应的能量越高,则该周期(频率)正弦波参与构成曲线的比例越高。
S3、根据频率-幅值图获取各组负荷数据的主要构成频率,根据主要构成频率和频率-幅值图生成候选特征集;
具体的,图3中的峰值所对应的频率,即主要构成频率。
作为本实施例的一种具体算例,在图3中主要构成频率为8736、6553、4368、3276、2184、1092、156、12、6和3等;
具体的,候选特征集中的数据为主要构成频率对应的幅值和相角;
具体的,根据实际需求可以在候选特征集中手动增加某段频率对应的幅值和相角,即预设频率的幅值和相角;
比如,希望比较居民用户用电负荷与工商业用电负荷特征,可定义幅值a12(日内2次峰值)与a168(每周1次峰值)之比作为候选特征。可预想,居民用户a12较高,而工商业用户a168较高,因此a12/a168对于两类负荷曲线将呈现较大区别,可作为区分两类用户的候选特征。此处举例仅用于说明用户可构造候选特征,但本实施例中并未涉及。
本实施例中,选取8736、6553、4368、3276、2184、1092、156、12、6、3等周期的归一化幅值和相位作为候选特征,幅值记为“m”,相位记为“ph”,候选特征记为[‘8736m’,‘8736ph’,‘6553m’,‘6553ph’,……];
S4、采用层次聚类法处理候选特征集得到第一聚类结果;
具体的,聚类接通过如图5所示,
以0/1/2/3分别代表不同类,则p1~p24分属类型为:
[0,1,1,2,0,2,1,3,1,1,2,0,0,1,2,1,3,0,1,0,1,1,2,0];
S5、采用随机森林分类器处理第一聚类结果得到分类结果,根据分类结果设置重要限度值,根据重要限度值筛选分类结果得到特征数据集;
具体的,通过随机森林分类器得到的结果如图6所示;
根据图6预设重要度下限为0.025;筛除重要度低于该值的特征,将['2184m','8736m','6553m','3276ph','1092ph','6m','3ph','12ph','6ph','4368ph']作为保留特征。
具体的,重要度限值非固定值,也是无法预先计算的。只有根据上述特征重要度分析结果,进行初步设定,然后在后续迭代过程中修正。
S6、采用层次聚类法处理特征数据集得到第二聚类结果;
具体的,第二次聚类结果如图7所示;
以0/1/2/3分别代表不同类,则p1~p24分属类型为:
[0,1,1,0,0,2,1,3,1,1,0,0,0,1,0,1,3,0,1,0,1,1,1,0];
S7、将第一聚类结果与第二聚类结果进行比较,若吻合度大于等于预设值则输出特征数据集;若吻合度小于预设值则返回步骤S5降低重要限度值;
作为本实施例的一种优选实施方式,预设值为80%,考虑到原候选特征分类可能存在“过拟合”情况,且经过筛选后能够获得更为清晰的聚类结果,因此对筛选后保留的特征予以保留,作为该系列负荷数据的频域特征。
实施例2
一种负荷特征提取装置,如图2所示,包括:
数据获取模块:用于获取预设区域内的若干组负荷数据;
具体的,每组负荷数据可以是不同地区、不同行业、不同电压等级或不同用户;
作为本实施例的一种优选实施方式,获取预设区域内的负荷数据的步骤中,负荷数据的采集周期为5分钟一次数据采集或15分钟一次数据采集,采集周期不宜超过1小时一次,负荷数据的时间跨度大于等于1年,从而保证获取的负荷数据准确性和全面性,提高最终负荷提取的准确性。
具体的,在数据获取模块中还包括:
对负荷数据进行数据治理,删除无效数据,统一负荷数据的时间跨度和时间颗粒度;
对负荷数据进行最大值归一化;
Xn(t)=xm(t)/xpeak;
其中xm(t)为负荷数据中的值,xpeak为该组负荷数据中的最大值,xn(t)为归一化结果。
具体的,最大值归一化是指将找到每组负荷数据中的最大值,并每个负荷数据与本组最大负荷数据相除,得到最大值均为1的负荷数据;
根据归一化后的数据可以得到归一化结果和时间折线图,如图4所示;
通过两步预处理,有效清除无效数据,提高数据准确性,简化后续计算,提高准确性。
数据处理模块:采用快速傅里叶变换处理若干组负荷数据生成频率-幅值图;
具体的,快速傅里叶变换处理根据公式:
Xk=FFT(xn(t));
式中,Xk=(ak,θk)为快速傅里叶变换后的复数结果,k为频率,ak表示k频率下的振幅,θk表示k频率下的相角。
具体的,频率-幅值图也可称为周期-能量图,如图3所示,纵坐标为不同周期成分的能量,图3是由若干不同周期、不同强度的正弦波构成,某周期(频率)对应的能量越高,则该周期(频率)正弦波参与构成曲线的比例越高。
数据分析模块:用于根据频率-幅值图获取各组负荷数据的主要构成频率,根据主要构成频率和频率-幅值图生成候选特征集;
具体的,图3中的峰值所对应的频率,即主要构成频率。
作为本实施例的一种具体算例,在图3中主要构成频率为8736、6553、4368、3276、2184、1092、156、12、6和3等;
具体的,候选特征集中的数据为主要构成频率对应的幅值和相角;
具体的,根据实际需求可以在候选特征集中手动增加某段频率的幅值比;
比如,希望比较居民用户用电负荷与工商业用电负荷特征,可定义幅值a12(日内2次峰值)与a168(每周1次峰值)之比作为候选特征。可预想,居民用户a12较高,而工商业用户a168较高,因此a12/a168对于两类负荷曲线将呈现较大区别,可作为区分两类用户的候选特征。此处举例仅用于说明用户可构造候选特征,但本实施例中并未涉及。
本实施例中,选取8736、6553、4368、3276、2184、1092、156、12、6、3等周期的归一化幅值和相位作为候选特征,幅值记为“m”,相位记为“ph”,候选特征记为[‘8736m’,‘8736ph’,‘6553m’,‘6553ph’,……];
第一聚类模块:采用层次聚类法处理候选特征集得到第一聚类结果;
具体的,聚类接通过如图5所示,
以0/1/2/3分别代表不同类,则p1~p24分属类型为:
[0,1,1,2,0,2,1,3,1,1,2,0,0,1,2,1,3,0,1,0,1,1,2,0];
筛选模块:用于采用随机森林分类器处理第一聚类结果得到分类结果,根据分类结果设置重要限度值,根据重要限度值筛选分类结果得到特征数据集;
具体的,通过随机森林分类器得到的结果如图6所示;
根据图6预设重要度下限为0.025;筛除重要度低于该值的特征,将['2184m','8736m','6553m','3276ph','1092ph','6m','3ph','12ph','6ph','4368ph']作为保留特征。
具体的,重要度限值非固定值,也是无法预先计算的。只有根据上述特征重要度分析结果,进行初步设定,然后在后续迭代过程中修正。
第二聚类模块:采用层次聚类法处理特征数据集得到第二聚类结果;
具体的,第二次聚类结果如图7所示;
以0/1/2/3分别代表不同类,则p1~p24分属类型为:
[0,1,1,0,0,2,1,3,1,1,0,0,0,1,0,1,3,0,1,0,1,1,1,0];
第二聚类模块:采用层次聚类法处理特征数据集得到第二聚类结果;
迭代输出模块:用于将第一聚类结果与第二聚类结果进行比较,若吻合度大于等于预设值则输出特征数据集;
若吻合度小于预设值则返回筛选模块,降低重要限度值。
作为本实施例的一种优选实施方式,预设值为80%,考虑到原候选特征分类可能存在“过拟合”情况,且经过筛选后能够获得更为清晰的聚类结果,因此对筛选后保留的特征予以保留,作为该系列负荷数据的频域特征。
实施例3
一种计算机设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现实施例1所述的一种负荷特征提取方法。
实施例4
一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现实施例1所述的一种负荷特征提取方法。
由技术常识可知,本发明可以通过其它的不脱离其精神实质或必要特征的实施方案来实现。因此,上述公开的实施方案,就各方面而言,都只是举例说明,并不是仅有的。所有在本发明范围内或在等同于本发明的范围内的改变均被本发明包含。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。
Claims (10)
1.一种负荷特征提取方法,其特征在于,包括以下步骤:
S1、获取预设区域内的若干组负荷数据;
S2、采用快速傅里叶变换处理若干组负荷数据生成频率-幅值图;
S3、根据频率-幅值图获取各组负荷数据的主要构成频率,根据主要构成频率和频率-幅值图生成候选特征集;
S4、采用层次聚类法处理候选特征集得到第一聚类结果;
S5、采用随机森林分类器处理第一聚类结果得到分类结果,根据分类结果设置重要限度值,根据重要限度值筛选分类结果得到特征数据集;
S6、采用层次聚类法处理特征数据集得到第二聚类结果;
S7、将第一聚类结果与第二聚类结果进行比较,若吻合度大于等于预设值则输出特征数据集;
若吻合度小于预设值则返回S5,调整重要限度值。
2.根据权利要求1所述的一种负荷特征提取方法,其特征在于,所述获取预设区域内的若干组负荷数据的步骤中,还包括对负荷数据进行预处理:
所述预处理包括:
对负荷数据进行数据治理,删除无效数据,统一负荷数据的时间跨度和时间颗粒度;
对负荷数据进行最大值归一化。
3.根据权利要求2所述的一种负荷特征提取方法,其特征在于,所述快速傅里叶变换的公式:
Xk = FFT(xn(t));
式中,Xk = (ak, θk)为快速傅里叶变换后的复数结果,xn(t)为负荷数据进行最大值归一化后的值;k为频率,ak表示k频率下的振幅,θk表示k频率下的相角。
4.根据权利要求1所述的一种负荷特征提取方法,其特征在于,所述候选特征集中还包括预设频率的幅值和相角。
5.根据权利要求1所述的一种负荷特征提取方法,其特征在于,所述调整重要限度值的步骤中,重要限度值逐渐减小。
6.根据权利要求1所述的一种负荷特征提取方法,其特征在于,所述预设值为80%。
7.根据权利要求1所述的一种负荷特征提取方法,其特征在于,所述主要构成频率为频率-幅值图中的各峰值。
8.一种负荷特征提取装置,其特征在于,包括:
数据获取模块:用于获取预设区域内的若干组负荷数据;
数据处理模块:采用快速傅里叶变换处理若干组负荷数据生成频率-幅值图;
数据分析模块:用于根据频率-幅值图获取各组负荷数据的主要构成频率,根据主要构成频率和频率-幅值图生成候选特征集;
第一聚类模块:采用层次聚类法处理候选特征集得到第一聚类结果;
筛选模块:用于采用随机森林分类器处理第一聚类结果得到分类结果,根据分类结果设置重要限度值,根据重要限度值筛选分类结果得到特征数据集;
第二聚类模块:采用层次聚类法处理特征数据集得到第二聚类结果;
迭代输出模块:用于将第一聚类结果与第二聚类结果进行比较,若吻合度大于等于预设值则输出特征数据集;
若吻合度小于预设值则返回筛选模块,调整重要限度值。
9.一种计算机设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1-7任一项所述的一种负荷特征提取方法。
10.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1-7中任一项所述的一种负荷特征提取方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410273266.7A CN117874497B (zh) | 2024-03-11 | 2024-03-11 | 一种负荷特征提取方法、装置、设备及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410273266.7A CN117874497B (zh) | 2024-03-11 | 2024-03-11 | 一种负荷特征提取方法、装置、设备及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117874497A true CN117874497A (zh) | 2024-04-12 |
CN117874497B CN117874497B (zh) | 2024-05-07 |
Family
ID=90588798
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410273266.7A Active CN117874497B (zh) | 2024-03-11 | 2024-03-11 | 一种负荷特征提取方法、装置、设备及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117874497B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095884A (zh) * | 2015-08-31 | 2015-11-25 | 桂林电子科技大学 | 一种基于随机森林支持向量机的行人识别系统及处理方法 |
CN108780661A (zh) * | 2016-03-16 | 2018-11-09 | 皇家飞利浦有限公司 | 用于改善将具有相似简档的患者聚类在一起的聚类模型的性能的相关性反馈 |
CN111914900A (zh) * | 2020-07-03 | 2020-11-10 | 深圳供电局有限公司 | 一种用户用电模式分类方法 |
CN113743977A (zh) * | 2021-06-28 | 2021-12-03 | 国网上海市电力公司 | 一种基于用户行为的用电数据特征提取方法及系统 |
US20220253641A1 (en) * | 2021-02-09 | 2022-08-11 | Samsung Sds Co., Ltd. | Method and apparatus for clustering images |
CN116933962A (zh) * | 2023-03-17 | 2023-10-24 | 安徽理工大学 | 一种面向电力用户负荷特性和可调节潜力的聚类方法 |
CN117113114A (zh) * | 2023-08-30 | 2023-11-24 | 无锡军工智能电气股份有限公司 | 基于信息熵的aco-fcm与特征选择的非侵入式负荷监测方法 |
CN117609818A (zh) * | 2023-11-24 | 2024-02-27 | 哈尔滨工业大学 | 基于聚类与信息熵的电网关联关系发现方法 |
-
2024
- 2024-03-11 CN CN202410273266.7A patent/CN117874497B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095884A (zh) * | 2015-08-31 | 2015-11-25 | 桂林电子科技大学 | 一种基于随机森林支持向量机的行人识别系统及处理方法 |
CN108780661A (zh) * | 2016-03-16 | 2018-11-09 | 皇家飞利浦有限公司 | 用于改善将具有相似简档的患者聚类在一起的聚类模型的性能的相关性反馈 |
US20190079938A1 (en) * | 2016-03-16 | 2019-03-14 | Koninklijke Philips N.V. | Relevance feedback to improve the performance of clustering model that clusters patients with similar profiles together |
CN111914900A (zh) * | 2020-07-03 | 2020-11-10 | 深圳供电局有限公司 | 一种用户用电模式分类方法 |
US20220253641A1 (en) * | 2021-02-09 | 2022-08-11 | Samsung Sds Co., Ltd. | Method and apparatus for clustering images |
CN113743977A (zh) * | 2021-06-28 | 2021-12-03 | 国网上海市电力公司 | 一种基于用户行为的用电数据特征提取方法及系统 |
CN116933962A (zh) * | 2023-03-17 | 2023-10-24 | 安徽理工大学 | 一种面向电力用户负荷特性和可调节潜力的聚类方法 |
CN117113114A (zh) * | 2023-08-30 | 2023-11-24 | 无锡军工智能电气股份有限公司 | 基于信息熵的aco-fcm与特征选择的非侵入式负荷监测方法 |
CN117609818A (zh) * | 2023-11-24 | 2024-02-27 | 哈尔滨工业大学 | 基于聚类与信息熵的电网关联关系发现方法 |
Non-Patent Citations (1)
Title |
---|
冉博路等: "基于二次聚类的工业用户负荷可调节潜力特征 提取与综合评估方法", 《电力系统保护与控制》, vol. 51, no. 18, 16 September 2023 (2023-09-16) * |
Also Published As
Publication number | Publication date |
---|---|
CN117874497B (zh) | 2024-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Qiu et al. | Power quality disturbances recognition using modified s transform and parallel stack sparse auto-encoder | |
CN108335010B (zh) | 一种风电出力时间序列建模方法及系统 | |
Cheng et al. | Enhanced state estimation and bad data identification in active power distribution networks using photovoltaic power forecasting | |
CN110909941B (zh) | 基于lstm神经网络的电力负荷预测方法、装置及系统 | |
Xue et al. | Adaptive ultra-short-term wind power prediction based on risk assessment | |
CN113705929B (zh) | 一种基于负荷特征曲线及典型特征值融合的春节假期负荷预测方法 | |
CN117874497B (zh) | 一种负荷特征提取方法、装置、设备及介质 | |
Yang et al. | Short-term wind power prediction based on CEEMDAN and parallel CNN-LSTM | |
Kamisan et al. | Forecasting Wind Speed Data by Using a Combination of ARIMA Model with Single Exponential Smoothing. | |
CN116578855A (zh) | 一种交流微电网的电能计量方法、系统、设备和介质 | |
Mashlakov et al. | Multi-timescale forecasting of battery energy storage state-of-charge under frequency containment reserve for normal operation | |
CN116365511A (zh) | 有源配电网模型构建方法、装置、终端及存储介质 | |
Jian-Hui et al. | Prediction of gold price based on WT-SVR and EMD-SVR model | |
Zhu et al. | Prediction of air quality index based on wavelet transform combination model | |
CN113949065B (zh) | 一种确定电力系统暂态电压稳定薄弱节点的方法及系统 | |
CN113627041A (zh) | 基于多源数据融合的时变动态负荷在线建模方法及装置 | |
CN114638171A (zh) | 电网项目投资预测方法、装置、存储介质及设备 | |
Liao et al. | A method for disturbance identification in power quality based on cross-attention fusion of temporal and spatial features | |
Kahawala et al. | Comparative Evaluation of Gradient Boosting with Active Thresholding and Model Explainability for Peak Demand Forecasting | |
CN115840078A (zh) | 基于特征优选与波形匹配的非侵入式负荷监测方法及系统 | |
CN117996838B (zh) | 一种基于非侵入负荷监测的分布式光伏识别装置 | |
CN117391343B (zh) | 一种基于城市规划的用电预测方法、装置和存储介质 | |
CN113807387B (zh) | 基于svm分类的特征指标风电出力时间序列构建方法和装置 | |
CN111626471A (zh) | 一种基于分解与聚合的风电场风速组合预测方法及装置 | |
Nayak et al. | Capacity value estimation of wind power incorporating hourly wind speed |
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 |