CN104376402B - 一种基于频域指标的负荷特性分类与综合方法 - Google Patents

一种基于频域指标的负荷特性分类与综合方法 Download PDF

Info

Publication number
CN104376402B
CN104376402B CN201410602513.XA CN201410602513A CN104376402B CN 104376402 B CN104376402 B CN 104376402B CN 201410602513 A CN201410602513 A CN 201410602513A CN 104376402 B CN104376402 B CN 104376402B
Authority
CN
China
Prior art keywords
distribution transformer
frequency domain
load
transformer area
load curve
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.)
Active
Application number
CN201410602513.XA
Other languages
English (en)
Other versions
CN104376402A (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.)
Hunan University
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan Electric Power Co Ltd
Original Assignee
Hunan University
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan 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 Hunan University, State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd, State Grid Hunan Electric Power Co Ltd filed Critical Hunan University
Priority to CN201410602513.XA priority Critical patent/CN104376402B/zh
Publication of CN104376402A publication Critical patent/CN104376402A/zh
Application granted granted Critical
Publication of CN104376402B publication Critical patent/CN104376402B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • 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

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种基于频域指标的负荷分类与综合方法，首先从监测控制和数据采集系统获取待分类的配变台区，对数据进行预处理，再对数据进行归一化处理，用标幺值表示，通过一周或几周数据取平均值，获得典型日负荷曲线；后基于定义的频域指标描述体系，采用指定的聚类方法，得到频域特征参数相似的配变台区类，最后得到各分类用户的典型负荷曲线及其频域指标描述集。此方案相比传统负荷分类与综合方法，能有效减少数据维数，降低计算开销和存储容量，负荷分类效果与传统方法相近。

Description

一种基于频域指标的负荷特性分类与综合方法
技术领域
本发明属于电力系统领域，具体涉及一种基于频域指标的负荷特性分类与综合方法。
背景技术
负荷分类与综合是电力系统规划、错峰管理、分时电价、负荷预测的基础，好的负荷分类与综合方法能给系统规划、错峰管理等提供正确的依据和指导。
电力系统用户类型繁多，各种类型的用户呈现不同的负荷特性，当前缺乏符合电网实际的、科学的负荷分类方法及完善的负荷特性指标体系。目前，负荷特性指标主要基于用户日负荷曲线(48点或96点数据)进行提取；指标体系主要采用时域指标，如日负荷率、日最小负荷率、日峰谷差，日峰谷差率、峰期负载率、平期负载率、谷期负载率等。经相关研究发现，仅通过这些指标都无法精确详细描述负荷特性。因此通常需要通过48点或96点日负荷曲线数据才能有效进行负荷特性聚类分析，随着电网的日益发展，需要统计和存储的数据量逐渐增大，维护成本较高。
以聚类为基础的数据挖掘算法是近年来新兴的智能负荷分析与处理手段，主要是采取各种聚类方法，如K-means(k均值聚类算法)、FCM(模糊聚类算法)、SOM(自组织映射算法)等对区域电力系统变压器或者用户的负荷数据进行聚类分析。已有专利提出了一种以K-means划分聚类为基础的典型日负荷曲线提取方法，其需要预先输入聚类数目K，且聚类以48点或96点日负荷曲线为依据，所需计算开销较大，数据存储量大。
发明内容
针对现有技术的不足，本发明提出一种基于频域指标的负荷特性分类与综合描述方法，其目的在于克服现有技术中进行负荷特性分类时，计算量大且数据存储量大的问题。
一种基于频域指标的负荷分类与综合方法，包括如下步骤：
步骤1：从监测控制和数据采集系统获取待分类的配变台区同属相似工作日的5-10组日负荷曲线数据；
所述相似工作日是指，对日类型相同的配变台区日负荷数据，依据气象数据根据灰色关联度分析得到关联度大于或等于0.6；即日负荷曲线具有相似性；
步骤2：对配变台区的每组日负荷曲线数据按下式进行归一化处理；
其中，Pn为配变台区日负荷曲线在第n个数据采集时刻的有功功率值，Pmax为同属相似工作日的日负荷曲线中该配变台区日负荷有功功率的最大值，xn为配变台区在第n个数据采集时刻的配变台区有功功率归一化数值；
步骤3：对步骤2得到的配变台区的多组已归一化处理的日负荷曲线数据在对应的数据采集点分别获取功率平均值，以每个数据采集点的功率平均值作为配变台区在对应相似工作日的典型日负荷曲线，典型日负荷曲线上每个数据采集点的功率值记为
步骤4：对步骤3得到的配变台区的典型日负荷曲线进行离散傅里叶变换得到第k次谐波幅值Ak与相位值αk；
其中，n＝0,1,2,...,95；
步骤5：计算每个配变台区典型日负荷曲线的频域指标：
幅值归一化指标：
其中，Aik为第i个配变台区的典型日负荷曲线离散傅里叶变换得到的第k次谐波幅值；
相位值归一化指标：
其中，αik为第i个配变台区的典型日离散傅里叶变换得到的第k次谐波相位值；
第i个配变台区的基波与前四次谐波构成的频域指标特性集：fi＝{(Ψ'ik,Ψ”ik),k＝0,1,2,3,4}；
步骤6：按照步骤1-步骤5所述的方法计算得到的每个待分类的配变台区的频域指标特性集，采用聚类算法对频域指标特性集进行聚类，将频域指标特性集中频域指标中的幅值归一化指标和相位归一化指标分别累加后得到两个累加值，同一类负荷特性的配变台区中任意两个配变台区的两个相应的累加值之差的百分比均在±10％范围内，实现负荷特性分类。
如被归类为具有同一负荷特性的配变台区A和配变台区B，配变台区A的幅值归一化指标累加值为c，相位归一化指标累加值为d，配变台区B的幅值归一化指标累加值为e，相位归一化指标累加值为f；则与两者均在±10％范围内；
所述聚类算法包括K均值或模糊C均值算法。
对某一类负荷特性的配变台区的频域指标进行统计，得到该类负荷相应的频域指标特性集，实现同一类负荷特性的配变台区的频域指标的综合。
将综合得到的同一类负荷特性频域指标用于描述同一类负荷特性特征。
有益效果
本发明提供了一种基于频域指标的负荷分类与综合方法，通过对配变台区日负荷曲线进行离散傅立叶变换，得到基波及各次谐波的幅值和频率，低次谐波主要描述曲线的轮廓，而高次谐波描述曲线的细节部分。仅通过保存基波和一部分低次谐波数据，就可以表征负荷特性的一些重要特征，还能重构日负荷曲线。将基波及低次谐波幅值和相位数据经过一定处理，作为负荷特性的频域指标，同时将高次谐波数据舍弃，其包含的较少数据可以基本重建日负荷曲线，表征负荷曲线的信息，达到较准确描述负荷特性的目标，能有效减少数据维数，以前为48点或96点数据，而本发明中只需10个频域指标数据，降低计算开销和存储容量，负荷分类效果与传统方法相近。
附图说明
图1为本发明所述方法的流程示意图；
图2为实施例中按本发明所述方法得到的配变台区1的典型日负荷曲线示意图；
图3为实施例中按本发明所述方法得到的配变台区2的典型日负荷曲线示意图；
图4为对应图2的频域指标示意图，其中，(a)为幅值归一化指标，(b)为相位归一化指标；
图5为对应图3的频域指标示意图，其中，(a)为幅值归一化指标，(b)为相位归一化指标。
具体实施方式
下面将结合附图和实施例对本发明做进一步的说明。
如图1所示，为本发明所述方法的流程示意图，一种基于频域指标的负荷分类与综合方法，包括如下步骤：
步骤1：从监测控制和数据采集系统获取待分类的配变台区同属相似工作日的5-10组日负荷曲线数据；
采集得到的负荷数据可能会出现丢失或失真，需进行数据预处理，即对失真数据采用数据平滑或邻近数据加权方法进行补全；
所述相似工作日是指，对日类型相同的配变台区日负荷数据，依据气象数据根据灰色关联度分析得到关联度大于或等于0.6；即日负荷曲线具有相似性；
步骤2：对配变台区的每组日负荷曲线数据按下式进行归一化处理；
其中，Pn为配变台区日负荷曲线在第n个数据采集时刻的有功功率值，Pmax为同属相似工作日的日负荷曲线中该配变台区日负荷有功功率的最大值，xn为配变台区在第n个数据采集时刻的配变台区有功功率归一化数值；
步骤3：对步骤2得到的配变台区的多组已归一化处理的日负荷曲线数据在对应的数据采集点分别获取功率平均值，以每个数据采集点的功率平均值作为配变台区在对应相似工作日的典型日负荷曲线，典型日负荷曲线上每个数据采集点的功率值记为
得到如图2和图3所示的配变台区1和配变台区2的典型日负荷曲线；
步骤4：对步骤3得到的配变台区的典型日负荷曲线进行离散傅里叶变换得到第k次谐波幅值Ak与相位值αk；
其中，n＝0,1,2,...,95；
步骤5：计算每个配变台区典型日负荷曲线的频域指标：
幅值归一化指标：
其中，Aik为第i个配变台区的典型日负荷曲线离散傅里叶变换得到的第k次谐波幅值；
相位值归一化指标：
其中，αik为第i个配变台区的典型日离散傅里叶变换得到的第k次谐波相位值；
第i个配变台区的基波与前四次谐波构成的频域指标特性集：fi＝{(Ψ'ik,Ψ”ik),k＝0,1,2,3,4}；
以上两式将幅值、相位归一化，将其幅值限制在(0,1]范围内；
得到图4和图5，图4为对应图2的频域指标示意图，其中，(a)为幅值归一化指标，(b)为相位归一化指标；图5为对应图3的频域指标示意图，其中，(a)为幅值归一化指标，(b)为相位归一化指标。
步骤6：按照步骤1-步骤5所述的方法计算得到的每个待分类的配变台区的频域指标特性集，采用聚类算法对频域指标特性集进行聚类，将频域指标特性集中频域指标中的幅值归一化指标和相位归一化指标分别累加后得到两个累加值，同一类负荷特性的配变台区中任意两个配变台区的两个相应的累加值之差的百分比均在±10％范围内，实现负荷特性分类。
表1
由图4和图5得到表1，利用表1求得这两个配变台区的频域指标的累加值的对应差值；
1号配变台区的幅值指标累加值与相位指标累加值分别为：
2号配变台区的幅值指标累加值与相位指标累加值分别为：
两配变台区的幅值累加值之差的百分比为：
两配变台区的相位累加值之差的百分比为：
由步骤6所述的评判方法，可以判定，配变台区1的负荷特性与配变台区2的负荷特性是相同的。
所述聚类算法包括K均值或模糊C均值算法。
对某一类负荷特性的配变台区的频域指标进行统计，得到该类负荷相应的频域指标特性集，实现同一类负荷特性的配变台区的频域指标的综合。
将综合得到的同一类负荷特性频域指标用于描述同一类负荷特性特征。

Claims (3)

1.一种基于频域指标的负荷分类与综合方法，其特征在于，包括如下步骤：
步骤1：从监测控制和数据采集系统获取待分类的配变台区同属相似工作日的5-10组日负荷曲线数据；
步骤2：对配变台区的每组日负荷曲线数据按下式进行归一化处理；
x n = P n P max , n = 0 , 1 , 2 , ... , 95
其中，Pn为配变台区日负荷曲线在第n个数据采集时刻的有功功率值，Pmax为同属相似工作日的日负荷曲线中该配变台区日负荷有功功率的最大值，xn为配变台区在第n个数据采集时刻的配变台区有功功率归一化数值；
步骤3：对步骤2得到的配变台区的多组已归一化处理的日负荷曲线数据在对应的数据采集点分别获取功率平均值，以每个数据采集点的功率平均值作为配变台区在对应相似工作日的典型日负荷曲线，典型日负荷曲线上每个数据采集点的功率值记为xn；
步骤4：对步骤3得到的配变台区的典型日负荷曲线进行离散傅里叶变换得到第k次谐波幅值Ak与相位值αk；
A k = real 2 ( X k ) + img 2 ( X k )
α k = a r c t a n [ i m g ( X k ) r e a l ( X k ) ]
其中，
步骤5：计算每个配变台区典型日负荷曲线的频域指标：
幅值归一化指标：
其中，Aik为第i个配变台区的典型日负荷曲线离散傅里叶变换得到的第k次谐波幅值；
相位值归一化指标：
其中，αik为第i个配变台区的典型日离散傅里叶变换得到的第k次谐波相位值；
第i个配变台区的基波与前四次谐波构成的频域指标特性集：fi＝{(Ψ'ik,Ψ”ik),k＝0,1,2,3,4}；
步骤6：按照步骤1-步骤5所述的方法计算得到的每个待分类的配变台区的频域指标特性集，采用聚类算法对频域指标特性集进行聚类，将频域指标特性集中频域指标中的幅值归一化指标和相位归一化指标分别累加后得到两个累加值，同一类负荷特性的配变台区中任意两个配变台区的两个相应的累加值之差的百分比均在±10％范围内，实现负荷特性分类；
其中，任意两个配变台区的幅值的累加值之差的百分比为：
任意两个配变台区的相位的累加值之差的百分比为：
2.根据权利要求1所述的基于频域指标的负荷分类与综合方法，其特征在于，所述聚类算法包括K均值或模糊C均值算法。
3.根据权利要求1所述基于频域指标的负荷分类与综合方法，其特征在于，对某一类负荷特性的配变台区的频域指标进行统计，得到该类负荷相应的频域指标特性集，实现同一类负荷特性的配变台区的频域指标的综合。
CN201410602513.XA 2014-10-31 2014-10-31 一种基于频域指标的负荷特性分类与综合方法 Active CN104376402B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410602513.XA CN104376402B (zh) 2014-10-31 2014-10-31 一种基于频域指标的负荷特性分类与综合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410602513.XA CN104376402B (zh) 2014-10-31 2014-10-31 一种基于频域指标的负荷特性分类与综合方法

Publications (2)

Publication Number Publication Date
CN104376402A CN104376402A (zh) 2015-02-25
CN104376402B true CN104376402B (zh) 2017-06-13

Family

ID=52555297

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410602513.XA Active CN104376402B (zh) 2014-10-31 2014-10-31 一种基于频域指标的负荷特性分类与综合方法

Country Status (1)

Country Link
CN (1) CN104376402B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105048450B (zh) * 2015-06-26 2017-03-08 国家电网公司 考虑220kV以下电源的电网稳定计算数据构建方法
CN107453346B (zh) * 2016-05-30 2022-03-18 中国电力科学研究院 一种适用于配电网重构的负荷曲线时段划分方法
CN107656898B (zh) * 2017-10-16 2020-11-10 国电南瑞科技股份有限公司 一种需求响应资源聚类方法
CN110334916A (zh) * 2019-06-17 2019-10-15 华北电力大学 一种基于频域变换的用户可调度负荷潜能评估方法
CN112686283A (zh) * 2020-12-15 2021-04-20 国网江苏省电力有限公司 一种基于谐波特性的电力负荷分类方法
CN112734261B (zh) * 2021-01-18 2023-05-16 国网山东省电力公司菏泽供电公司 一种配电网运行指标序列关联分析方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102156814A (zh) * 2011-04-06 2011-08-17 广东省电力设计研究院 基于聚类的典型日负荷曲线选取方法及装置
CN103177188A (zh) * 2013-04-02 2013-06-26 东南大学 一种基于特征映射的电力系统负荷动特性分类方法
CN103793788A (zh) * 2014-01-27 2014-05-14 国家电网公司 一种有序用电管理方法
CN104063480A (zh) * 2014-07-02 2014-09-24 国家电网公司 一种基于电力大数据的负荷曲线并行聚类方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102156814A (zh) * 2011-04-06 2011-08-17 广东省电力设计研究院 基于聚类的典型日负荷曲线选取方法及装置
CN103177188A (zh) * 2013-04-02 2013-06-26 东南大学 一种基于特征映射的电力系统负荷动特性分类方法
CN103793788A (zh) * 2014-01-27 2014-05-14 国家电网公司 一种有序用电管理方法
CN104063480A (zh) * 2014-07-02 2014-09-24 国家电网公司 一种基于电力大数据的负荷曲线并行聚类方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于用户日负荷曲线的用电行业分类与综合方法;李欣然;《电力系统自动化》;20100525;第34卷(第10期);第56-61页 *
电力系统的谐波分析算法及负荷谐波建模;任志莲;《中国优秀硕士学位论文全文数据库》;20100215(第2期);第1-93页 *

Also Published As

Publication number Publication date
CN104376402A (zh) 2015-02-25

Similar Documents

Publication Publication Date Title
CN104376402B (zh) 一种基于频域指标的负荷特性分类与综合方法
An et al. Short-term prediction of wind power using EMD and chaotic theory
US11043808B2 (en) Method for identifying pattern of load cycle
Qiu et al. Power quality disturbances recognition using modified s transform and parallel stack sparse auto-encoder
Hu et al. A hybrid forecasting approach applied to wind speed time series
CN102542262B (zh) 一种基于高铁运行特征工况波形库的波形识别方法
CN103177188B (zh) 一种基于特征映射的电力系统负荷动特性分类方法
CN102608441B (zh) 基于s变换的sdpqd信号识别方法、设备及系统
CN105572501A (zh) 一种基于sst变换和ls-svm的电能质量扰动识别方法
CN106447206A (zh) 基于用电信息采集数据的用电分析方法
CN104537433A (zh) 基于存量容量与业扩报装特征的售电量预测方法
CN105354866A (zh) 一种多边形轮廓相似度检测方法
Li et al. Life grade recognition method based on supervised uncorrelated orthogonal locality preserving projection and K-nearest neighbor classifier
CN102512158A (zh) 一种基于高维模糊识别的心电信号质量评估方法和装置
CN103323702A (zh) 复合电能质量扰动信号识别方法
CN110717619A (zh) 一种自下而上配电网多尺度时空负荷预测方法及系统
CN103400144A (zh) 一种基于k近邻的支持向量机主动学习的方法
CN104951787A (zh) 一种src框架下判别字典学习的电能质量扰动识别方法
CN103400039A (zh) 一种基于大风气象分类的风功率爬坡预测模型切换方法
Li et al. Adaptive S transform for feature extraction in voltage sags
Gao et al. Multi-step wind speed prediction based on LSSVM combined with ESMD and fractional-order beetle swarm optimization
CN105550450B (zh) 一种电能质量干扰源特征谐波建模方法
CN109726740A (zh) 一种基于聚类分析的行业用电行为分析方法
CN112819649A (zh) 确定台区户变关系的方法及装置
CN113128396B (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
GR01 Patent grant
GR01 Patent grant