CN110930035A - 一种基于区间统计值的谐波危害评估方法 - Google Patents

一种基于区间统计值的谐波危害评估方法 Download PDF

Info

Publication number
CN110930035A
CN110930035A CN201911181655.2A CN201911181655A CN110930035A CN 110930035 A CN110930035 A CN 110930035A CN 201911181655 A CN201911181655 A CN 201911181655A CN 110930035 A CN110930035 A CN 110930035A
Authority
CN
China
Prior art keywords
harmonic
harmonic current
value
interval
monitoring period
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
Application number
CN201911181655.2A
Other languages
English (en)
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.)
State Grid Fujian Electric Power Co Ltd
Putian Power Supply Co of State Grid Fujian Electric Power Co Ltd
Original Assignee
State Grid Fujian Electric Power Co Ltd
Putian Power Supply Co of State Grid Fujian 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 State Grid Fujian Electric Power Co Ltd, Putian Power Supply Co of State Grid Fujian Electric Power Co Ltd filed Critical State Grid Fujian Electric Power Co Ltd
Priority to CN201911181655.2A priority Critical patent/CN110930035A/zh
Publication of CN110930035A publication Critical patent/CN110930035A/zh
Pending legal-status Critical Current

Links

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
    • 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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression 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
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Data Mining & Analysis (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Tourism & Hospitality (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Mathematical Physics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Health & Medical Sciences (AREA)
  • Computational Mathematics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Primary Health Care (AREA)
  • Probability & Statistics with Applications (AREA)
  • General Health & Medical Sciences (AREA)
  • Algebra (AREA)
  • Water Supply & Treatment (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Public Health (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

本发明提供了一种基于区间统计值的谐波危害评估方法，谐波监测装置设置于系统与用户之间，用于记录用户注入系统的谐波电流值，并生成区间谐波电流的统计值；谐波电流值计算法具体为：所述谐波监测装置以3min或10min为单位时间记录2至25次的谐波电流的最大值、最小值、平均值、95％最大值，取最大值和最小值构成一个区间用以表征该时段用户注入谐波电流的范围，取多个时段构成区间谐波电流的时序统计值X。应用本技术方案可实现预测用户未来注入系统的区间谐波电流，并提出相应的谐波危害评估方法从而实现用户的谐波危害预警。

Description

一种基于区间统计值的谐波危害评估方法
技术领域
本发明涉及电网维护领域，具体是指一种基于区间统计值的谐波危害评估方法。
背景技术
近年来，谐波用户的并网给电网带来了严重的谐波污染，极大地影响了电网和用户的安全、经济运行。目前，谐波监测算法及装置研制方面已有较好的研究进展，特别是随着硬件技术、信号处理技术和人工智能技术的进步，谐波监测装置基本上实现了自动化，大量网络化监测系统已经得到了工程应用。现有谐波监测主要为电网公司提供监管服务，为用户能否接入提供定性判断依据。然而电力公司通过对超标用户的管制来维护电网的高质量运行水平，却无法从工业用户的用电需求出发，指导用户合理安排运行方式，从而避免因谐波污染停电而造成较大的经济损失。因此，目前电网的谐波监管方式具有鲜明的定性评判、事后评价特征，尚未能够为用户的优质运行提供决策依据。
目前，电网公司监管部门主要采取相关检测手段监测用户注入系统的谐波电流数值，并根据系统供电容量、短路容量水平、用户用电容量等参数计算其许可限值。当用户注入系统的谐波电流低于限值要求时，方可允许用户接入运行。然而用户的谐波污染大小与用户参数、电网运行方式密切相关，电网的谐波监测数据只能反映当前电网运行方式和用户参数作用下用户注入的谐波电流。由于用户注入的谐波电流因系统运行方式改变而改变，而谐波监测数据只能反映用户当前的谐波危害，缺无法预测用户未来的谐波危害。因此，根据当前的谐波监测数据来评估用户的谐波污染具有一定的片面性。鉴于尚无专利基于区间的谐波监测数据来预测用户注入系统的区间谐波电流，因此本发明提出一种基于区间数据来预测区间谐波电流的方法，在未来的用户参数和电网运行方式变化后通过预测用户注入的谐波电流，从而评估用户的谐波危害，以实现用户的谐波危害预警。目前针对用户谐波危害预警存在几个技术问题：
(一)由于用户在一段监测时段内注入的谐波电流具有一定的波动性，谐波监测装置记录着海量的监测数据，为了直观地描述该监测时段内用户注入的谐波电流大小，谐波监测装置一般以3min或10min为一监测时段，并生成该监测时段所对应的谐波电流的最大、最小、平均和95大值的统计值，一般取该监测时段的最大和最小值以构成一个区间谐波电流来表征用户注入系统的谐波电流的水平。因此如何基于历史的区间谐波电流统计值去预测用户未来注入系统的区间谐波电流值是一个关键的技术难题。
(二)用户注入系统的谐波电流是多维的，并且预测的谐波电流是一个区间量，如何基于多维的区间量来评估用户未来的谐波污染，从而实现用户谐波危害预警也是一个急需解决的问题。
因此，实际工程上迫切需要一种基于谐波监测装置获得的区间谐波电流来预测用户未来注入系统的区间谐波电流，并提出相应的谐波危害评估方法从而实现用户的谐波危害预警。
发明内容
本发明的目的在于提供一种基于区间统计值的谐波危害评估方法，预测用户未来注入系统的区间谐波电流，并提出相应的谐波危害评估方法从而实现用户的谐波危害预警。
为了解决上述技术问题，本发明提供了一种基于区间统计值的谐波危害评估方法，谐波监测装置设置于系统与用户之间，用于记录用户注入系统的谐波电流值，并生成区间谐波电流的统计值；
谐波电流值计算法具体为：所述谐波监测装置以3min或10min为单位时间记录2至25次的谐波电流的最大值、最小值、平均值、95％最大值，取最大值和最小值构成一个区间用以表征该时段用户注入谐波电流的范围，取多个时段构成区间谐波电流的时序统计值X为：
Figure BDA0002291435110000031
其中，n＝1,2,3,4……,
Figure BDA0002291435110000032
表示2次谐波在第一个监测时段内谐波电流的最小值，
Figure BDA0002291435110000033
表示2次谐波在第一个监测时段内谐波电流的最大值，
Figure BDA0002291435110000034
表示2次谐波在第二个监测时段内谐波电流的最小值，
Figure BDA0002291435110000035
表示2次谐波在第二个监测时段内谐波电流的最大值，
Figure BDA0002291435110000036
表示2次谐波在第n个监测时段内谐波电流的最小值，
Figure BDA0002291435110000037
表示2次谐波在第n个监测时段内谐波电流的最大值；
Figure BDA0002291435110000038
表示3次谐波在第一个监测时段内谐波电流的最小值，
Figure BDA0002291435110000039
表示3次谐波在第一个监测时段内谐波电流的最大值，
Figure BDA00022914351100000310
表示3次谐波在第二个监测时段内谐波电流的最小值，
Figure BDA00022914351100000311
表示3次谐波在第二个监测时段内谐波电流的最大值，
Figure BDA0002291435110000041
表示3次谐波在第n个监测时段内谐波电流的最小值，
Figure BDA0002291435110000042
表示3次谐波在第n个监测时段内谐波电流的最大值；
Figure BDA0002291435110000043
表示25次谐波在第一个监测时段内谐波电流的最小值，
Figure BDA0002291435110000044
表示25次谐波在第一个监测时段内谐波电流的最大值，
Figure BDA0002291435110000045
表示25次谐波在第二个监测时段内谐波电流的最小值，
Figure BDA0002291435110000046
表示25次谐波在第二个监测时段内谐波电流的最大值，
Figure BDA0002291435110000047
表示25次谐波在第n个监测时段内谐波电流的最小值，
Figure BDA0002291435110000048
表示25次谐波在第n个监测时段内谐波电流的最大值；其余数据含义依据此规律类推；
根据区间谐波电流的时序统计值X搭建区间自回归模型，对于第t个监测时段内的区间2次谐波电流统计值
Figure BDA0002291435110000049
该第t个监测时段内的区间内的任一点ht,2表示为：
Figure BDA00022914351100000410
其中，
Figure BDA00022914351100000411
表示2次谐波在第t个监测时段内谐波电流的最小值，
Figure BDA00022914351100000412
表示2次谐波在第t个监测时段内谐波电流的最大值；γt为区间比例参数，0≤γt≤1，给定一个γt值就得到对应区间谐波电流统计值
Figure BDA00022914351100000413
内的一个值；
根据区间谐波电流时序统计值X，以2次谐波区间谐波电流为例，构建该监测时段区间2次谐波电流上下限的自回归模型为：
Figure BDA0002291435110000051
Figure BDA0002291435110000052
t＝p+1,…n (16)
其中，
Figure BDA0002291435110000053
为第t个监测时间段内2次谐波电流的最小值，
Figure BDA0002291435110000054
为下限常数系数，
Figure BDA0002291435110000055
为下限的回归系数，
Figure BDA0002291435110000056
为下限误差，
Figure BDA0002291435110000057
为第t个监测时间段内2次谐波电流的最大值，
Figure BDA0002291435110000058
为上限常数系数，
Figure BDA0002291435110000059
为最大值的回归系数，
Figure BDA00022914351100000510
为上限误差，p为样本个数即取第t个监测时间段前的p个样本数据；
将公式(2)代入公式(3)可得：
Figure BDA00022914351100000511
Figure BDA00022914351100000512
t＝p+1,…n (17)
分别计算上下限自回归模型参数的回归系数：
计算最小值的回归系数，定义系数αi min为下限的最小值回归系数，定义
Figure BDA00022914351100000513
为下限的最大值回归系数，使其满足：
Figure BDA00022914351100000514
Figure BDA00022914351100000515
则最小值的自回归方程表示为：
Figure BDA0002291435110000061
故表示为矩阵形式即：
Ymin＝Xminβmin+εmin (20)
其中
Figure BDA0002291435110000062
Figure BDA0002291435110000063
这里的T表示矩阵的转置，
Figure BDA0002291435110000064
利用最小乘法辨识最小值的回归系数
Figure BDA0002291435110000065
为:
βmin＝((Xmin)TXmin)-1(Xmin)TYmin (21)
则，
Figure BDA0002291435110000066
计算最大值的回归系数，定义系数
Figure BDA0002291435110000067
为上限的最小值回归系数，定义
Figure BDA0002291435110000068
为上限的最大值回归系数使其满足：
Figure BDA0002291435110000071
Figure BDA0002291435110000072
则最大值的自回归方程表示为：
Figure BDA0002291435110000073
故表示为矩阵形式即：
Ymax＝Xmaxβmax+εmax (25)
其中，
Figure BDA0002291435110000074
Figure BDA0002291435110000075
Figure BDA0002291435110000076
利用最小乘法辨识最小值的回归系数
Figure BDA0002291435110000077
为:
βmax＝((Xmax)TXmax)-1(Xmax)TYmax (26)
则，
Figure BDA0002291435110000078
根据上下限的自回归系数，预测t+1时段用户注入系统的区间2次谐波电流值：
Figure BDA0002291435110000081
Figure BDA0002291435110000082
同理依次计算第t+1时段用户注入系统的区间3至25次的谐波电流值；
设定一定长度的监测时间长度，取该时间长度2至25次谐波电流的监测数据的最大值和最小值，计算该时间长度区间谐波电流值；取各该时间长度的区间谐波电流值上述谐波电流值计算法来预测下一所述时间长度的区间谐波电流值，取区间2至25次谐波电流值的最大值与计算得到的谐波电流限值进行比较，判断并计算该用户的2至25维谐波电流超标的维度，结合谐波危害等级以判断该用户下一所述时间长度的谐波电流危害，从而进行用户的一定时间长度谐波危害预警。
2.根据权利要求1所述的一种基于区间统计值的谐波危害评估方法，其特征在于，所述时间长度具体为一年或一月或一日。
3.根据权利要求1所述的一种基于区间统计值的谐波危害评估方法，其特征在于，所述超标维度为0时的危害等级具体为无危害，所述超标维度为1至2时的危害等级具体为轻微危害，超标维度为2至5时的危害等级具体为轻度危害，超标维度为5至10时的危害等级具体为中度危害，超标维度为10以上时的危害等级具体为重度危害。
相较于现有技术，本发明的技术方案具备以下有益效果：
该方法根据谐波监测装置测量得到用户的区间谐波电流统计值，从而构建用户的区间自回归模型，并预测用户未来注入系统的区间谐波电流上下限，根据预测的区间谐波电流评估用户的谐波危害，以实现用户的谐波危害预警。因此，该方法能预测用户的谐波电流并预判用户未来的谐波危害，能够克服目前电网只针对当前电网的运行情况无视未来用户参数和系统运行方式的变化去评价用户的谐波污染的局限性，具有较强的工程实用性和推广应用价值。
附图说明
图1为本发明优选实施例中用户谐波危害预警评估方法步骤。
具体实施方式
下文结合附图和具体实施方式对本发明做进一步说明。
一种基于区间统计值的谐波危害评估方法，其特征在于，谐波监测装置设置于系统与用户之间，用于记录用户注入系统的谐波电流值，并生成区间谐波电流的统计值；
谐波电流值计算法具体为：所述谐波监测装置以3min或10min为单位时间记录2至25次的谐波电流的最大值、最小值、平均值、95％最大值，取最大值和最小值构成一个区间用以表征该时段用户注入谐波电流的范围，取多个时段构成区间谐波电流的时序统计值X为：
Figure BDA0002291435110000091
其中，n＝1,2,3,4……,
Figure BDA0002291435110000092
表示2次谐波在第一个监测时段内谐波电流的最小值，
Figure BDA0002291435110000101
表示2次谐波在第一个监测时段内谐波电流的最大值，
Figure BDA0002291435110000102
表示2次谐波在第二个监测时段内谐波电流的最小值，
Figure BDA0002291435110000103
表示2次谐波在第二个监测时段内谐波电流的最大值，
Figure BDA0002291435110000104
表示2次谐波在第n个监测时段内谐波电流的最小值，
Figure BDA0002291435110000105
表示2次谐波在第n个监测时段内谐波电流的最大值；
Figure BDA0002291435110000106
表示3次谐波在第一个监测时段内谐波电流的最小值，
Figure BDA0002291435110000107
表示3次谐波在第一个监测时段内谐波电流的最大值，
Figure BDA0002291435110000108
表示3次谐波在第二个监测时段内谐波电流的最小值，
Figure BDA0002291435110000109
表示3次谐波在第二个监测时段内谐波电流的最大值，
Figure BDA00022914351100001010
表示3次谐波在第n个监测时段内谐波电流的最小值，
Figure BDA00022914351100001011
表示3次谐波在第n个监测时段内谐波电流的最大值；
Figure BDA00022914351100001012
表示25次谐波在第一个监测时段内谐波电流的最小值，
Figure BDA00022914351100001013
表示25次谐波在第一个监测时段内谐波电流的最大值，
Figure BDA00022914351100001014
表示25次谐波在第二个监测时段内谐波电流的最小值，
Figure BDA00022914351100001015
表示25次谐波在第二个监测时段内谐波电流的最大值，
Figure BDA00022914351100001016
表示25次谐波在第n个监测时段内谐波电流的最小值，
Figure BDA00022914351100001017
表示25次谐波在第n个监测时段内谐波电流的最大值；其余数据含义依据此规律类推；
根据区间谐波电流的时序统计值X搭建区间自回归模型，对于第t个监测时段内的区间2次谐波电流统计值
Figure BDA0002291435110000111
该第t个监测时段内的区间内的任一点ht,2表示为：
Figure BDA0002291435110000112
其中，
Figure BDA0002291435110000113
表示2次谐波在第t个监测时段内谐波电流的最小值，
Figure BDA0002291435110000114
表示2次谐波在第t个监测时段内谐波电流的最大值；γt为区间比例参数，0≤γt≤1，给定一个γt值就得到对应区间谐波电流统计值
Figure BDA0002291435110000115
内的一个值；
根据区间谐波电流时序统计值X，以2次谐波区间谐波电流为例，构建该监测时段区间2次谐波电流上下限的自回归模型为：
Figure BDA0002291435110000116
Figure BDA0002291435110000117
t＝p+1,…n (30)
其中，
Figure BDA0002291435110000118
为第t个监测时间段内2次谐波电流的最小值，
Figure BDA0002291435110000119
为下限常数系数，
Figure BDA00022914351100001110
为下限的回归系数，
Figure BDA00022914351100001111
为下限误差，
Figure BDA00022914351100001112
为第t个监测时间段内2次谐波电流的最大值，
Figure BDA00022914351100001113
为上限常数系数，
Figure BDA00022914351100001114
为最大值的回归系数，
Figure BDA00022914351100001115
为上限误差，p为样本个数即取第t个监测时间段前的p个样本数据；
将公式(2)代入公式(3)可得：
Figure BDA0002291435110000121
Figure BDA0002291435110000122
t＝p+1,…n (31)
分别计算上下限自回归模型参数的回归系数：
计算最小值的回归系数，定义系数
Figure BDA0002291435110000123
为下限的最小值回归系数，定义
Figure BDA0002291435110000124
为下限的最大值回归系数，使其满足：
Figure BDA0002291435110000125
Figure BDA0002291435110000126
则最小值的自回归方程表示为：
Figure BDA0002291435110000127
故表示为矩阵形式即：
Ymin＝Xminβmin+εmin (34)
其中
Figure BDA0002291435110000128
Figure BDA0002291435110000129
这里的T表示矩阵的转置，
Figure BDA0002291435110000131
利用最小乘法辨识最小值的回归系数
Figure BDA0002291435110000132
为:
βmin＝((Xmin)TXmin)-1(Xmin)TYmin (35)
则，
Figure BDA0002291435110000133
计算最大值的回归系数，定义系数
Figure BDA0002291435110000134
为上限的最小值回归系数，定义
Figure BDA0002291435110000135
为上限的最大值回归系数使其满足：
Figure BDA0002291435110000136
Figure BDA0002291435110000137
则最大值的自回归方程表示为：
Figure BDA0002291435110000138
故表示为矩阵形式即：
Ymax＝Xmaxβmax+εmax (39)
其中，
Figure BDA0002291435110000139
Figure BDA00022914351100001310
Figure BDA0002291435110000141
利用最小乘法辨识最小值的回归系数
Figure BDA0002291435110000142
为:
βmax＝((Xmax)TXmax)-1(Xmax)TYmax (40)
则，
Figure BDA0002291435110000143
根据上下限的自回归系数，预测t+1时段用户注入系统的区间2次谐波电流值：
Figure BDA0002291435110000144
Figure BDA0002291435110000145
同理依次计算第t+1时段用户注入系统的区间3至25次的谐波电流值；
设定一定长度的监测时间长度，取该时间长度2至25次谐波电流的监测数据的最大值和最小值，计算该时间长度区间谐波电流值；取各该时间长度的区间谐波电流值上述谐波电流值计算法来预测下一所述时间长度的区间谐波电流值，取区间2至25次谐波电流值的最大值与计算得到的谐波电流限值进行比较，判断并计算该用户的2至25维谐波电流超标的维度，结合谐波危害等级以判断该用户下一所述时间长度的谐波电流危害，从而进行用户的一定时间长度谐波危害预警。在本实施例中，所述时间长度具体为一年或一月或一日，就可以预测年谐波危害预警、月谐波危害预警、日谐波危害预警。
所述超标维度为0时的危害等级具体为无危害，所述超标维度为1至2时的危害等级具体为轻微危害，超标维度为2至5时的危害等级具体为轻度危害，超标维度为5至10时的危害等级具体为中度危害，超标维度为10以上时的危害等级具体为重度危害。
以上所述，仅为本发明较佳的具体实施方式，但本发明的设计构思并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，利用此构思对本发明进行非实质性的改动，均属于侵犯本发明保护范围的行为。

Claims (3)

1.一种基于区间统计值的谐波危害评估方法，其特征在于，谐波监测装置设置于系统与用户之间，用于记录用户注入系统的谐波电流值，并生成区间谐波电流的统计值；
谐波电流值计算法具体为：所述谐波监测装置以3min或10min为单位时间记录2至25次的谐波电流的最大值、最小值、平均值、95％最大值，取最大值和最小值构成一个区间用以表征该时段用户注入谐波电流的范围，取多个时段构成区间谐波电流的时序统计值X为：
Figure FDA0002291435100000011
其中，
Figure FDA0002291435100000012
表示2次谐波在第一个监测时段内谐波电流的最小值，
Figure FDA0002291435100000013
表示2次谐波在第一个监测时段内谐波电流的最大值，
Figure FDA0002291435100000014
表示2次谐波在第二个监测时段内谐波电流的最小值，
Figure FDA0002291435100000015
表示2次谐波在第二个监测时段内谐波电流的最大值，
Figure FDA0002291435100000016
表示2次谐波在第n个监测时段内谐波电流的最小值，
Figure FDA0002291435100000017
表示2次谐波在第n个监测时段内谐波电流的最大值；
Figure FDA0002291435100000018
表示3次谐波在第一个监测时段内谐波电流的最小值，
Figure FDA0002291435100000019
表示3次谐波在第一个监测时段内谐波电流的最大值，
Figure FDA00022914351000000110
表示3次谐波在第二个监测时段内谐波电流的最小值，
Figure FDA0002291435100000021
表示3次谐波在第二个监测时段内谐波电流的最大值，
Figure FDA0002291435100000022
表示3次谐波在第n个监测时段内谐波电流的最小值，
Figure FDA0002291435100000023
表示3次谐波在第n个监测时段内谐波电流的最大值；
Figure FDA0002291435100000024
表示25次谐波在第一个监测时段内谐波电流的最小值，
Figure FDA0002291435100000025
表示25次谐波在第一个监测时段内谐波电流的最大值，
Figure FDA0002291435100000026
表示25次谐波在第二个监测时段内谐波电流的最小值，
Figure FDA0002291435100000027
表示25次谐波在第二个监测时段内谐波电流的最大值，
Figure FDA0002291435100000028
表示25次谐波在第n个监测时段内谐波电流的最小值，
Figure FDA0002291435100000029
表示25次谐波在第n个监测时段内谐波电流的最大值；其余数据含义依据此规律类推；
根据区间谐波电流的时序统计值X搭建区间自回归模型，对于第t个监测时段内的区间2次谐波电流统计值
Figure FDA00022914351000000210
该第t个监测时段内的区间内的任一点ht,2表示为：
Figure FDA00022914351000000211
其中，
Figure FDA00022914351000000212
表示2次谐波在第t个监测时段内谐波电流的最小值，
Figure FDA00022914351000000213
表示2次谐波在第t个监测时段内谐波电流的最大值；γt为区间比例参数，0≤γt≤1，给定一个γt值就得到对应区间谐波电流统计值
Figure FDA0002291435100000031
内的一个值；
根据区间谐波电流时序统计值X，以2次谐波区间谐波电流为例，构建该监测时段区间2次谐波电流上下限的自回归模型为：
Figure FDA0002291435100000032
Figure FDA0002291435100000033
其中，
Figure FDA0002291435100000034
为第t个监测时间段内2次谐波电流的最小值，
Figure FDA0002291435100000035
为下限常数系数，
Figure FDA0002291435100000036
为下限的回归系数，
Figure FDA0002291435100000037
为下限误差，
Figure FDA0002291435100000038
为第t个监测时间段内2次谐波电流的最大值，
Figure FDA0002291435100000039
为上限常数系数，
Figure FDA00022914351000000310
为最大值的回归系数，
Figure FDA00022914351000000311
为上限误差，p为样本个数即取第t个监测时间段前的p个样本数据；
将公式(2)代入公式(3)可得：
Figure FDA00022914351000000312
Figure FDA00022914351000000313
分别计算上下限自回归模型参数的回归系数：
计算最小值的回归系数，定义系数为下限的最小值回归系数，定义
Figure FDA0002291435100000041
为下限的最大值回归系数，使其满足：
Figure FDA0002291435100000042
Figure FDA0002291435100000043
则最小值的自回归方程表示为：
Figure FDA0002291435100000044
故表示为矩阵形式即：
Ymin＝Xminβmin+εmin (6)
其中
Figure FDA0002291435100000045
Figure FDA0002291435100000046
Figure FDA0002291435100000047
这里的T表示矩阵的转置，
Figure FDA0002291435100000048
利用最小乘法辨识最小值的回归系数
Figure FDA0002291435100000049
为:
βmin＝((Xmin)TXmin)-1(Xmin)TYmin (7)
则，
Figure FDA0002291435100000051
计算最大值的回归系数，定义系数
Figure FDA0002291435100000052
为上限的最小值回归系数，定义
Figure FDA0002291435100000053
为上限的最大值回归系数使其满足：
Figure FDA0002291435100000055
则最大值的自回归方程表示为：
Figure FDA0002291435100000056
故表示为矩阵形式即：
Ymax＝Xmaxβmax+εmax (11)
其中，
Figure FDA0002291435100000057
Figure FDA0002291435100000058
Figure FDA0002291435100000059
利用最小乘法辨识最小值的回归系数
Figure FDA00022914351000000510
为:
βmax＝((Xmax)TXmax)-1(Xmax)TYmax (12)
则，
Figure FDA0002291435100000061
根据上下限的自回归系数，预测t+1时段用户注入系统的区间2次谐波电流值：
Figure FDA0002291435100000062
Figure FDA0002291435100000063
同理依次计算第t+1时段用户注入系统的区间3至25次的谐波电流值；
设定一定长度的监测时间长度，取该时间长度2至25次谐波电流的监测数据的最大值和最小值，计算该时间长度区间谐波电流值；取各该时间长度的区间谐波电流值上述谐波电流值计算法来预测下一所述时间长度的区间谐波电流值，取区间2至25次谐波电流值的最大值与计算得到的谐波电流限值进行比较，判断并计算该用户的2至25维谐波电流超标的维度，结合谐波危害等级以判断该用户下一所述时间长度的谐波电流危害，从而进行用户的一定时间长度谐波危害预警。
2.根据权利要求1所述的一种基于区间统计值的谐波危害评估方法，其特征在于，所述时间长度具体为一年或一月或一日。
3.根据权利要求1所述的一种基于区间统计值的谐波危害评估方法，其特征在于，所述超标维度为0时的危害等级具体为无危害，所述超标维度为1至2时的危害等级具体为轻微危害，超标维度为2至5时的危害等级具体为轻度危害，超标维度为5至10时的危害等级具体为中度危害，超标维度为10以上时的危害等级具体为重度危害。
CN201911181655.2A 2019-11-27 2019-11-27 一种基于区间统计值的谐波危害评估方法 Pending CN110930035A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911181655.2A CN110930035A (zh) 2019-11-27 2019-11-27 一种基于区间统计值的谐波危害评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911181655.2A CN110930035A (zh) 2019-11-27 2019-11-27 一种基于区间统计值的谐波危害评估方法

Publications (1)

Publication Number Publication Date
CN110930035A true CN110930035A (zh) 2020-03-27

Family

ID=69852177

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911181655.2A Pending CN110930035A (zh) 2019-11-27 2019-11-27 一种基于区间统计值的谐波危害评估方法

Country Status (1)

Country Link
CN (1) CN110930035A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111999580A (zh) * 2020-08-14 2020-11-27 佳都新太科技股份有限公司 地铁站台闸门的检测方法、装置、计算机设备和存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2352038A1 (en) * 2008-10-27 2011-08-03 Tsinghua University Method for detecting single phase grounding fault based on harmonic component of residual current
CN106485089A (zh) * 2016-10-21 2017-03-08 福州大学 谐波用户典型工况的区间参数获取方法
CN110472195A (zh) * 2019-08-15 2019-11-19 国网福建省电力有限公司莆田供电公司 一种基于区间监测数据的用户谐波危害预警评估方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2352038A1 (en) * 2008-10-27 2011-08-03 Tsinghua University Method for detecting single phase grounding fault based on harmonic component of residual current
CN106485089A (zh) * 2016-10-21 2017-03-08 福州大学 谐波用户典型工况的区间参数获取方法
CN110472195A (zh) * 2019-08-15 2019-11-19 国网福建省电力有限公司莆田供电公司 一种基于区间监测数据的用户谐波危害预警评估方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111999580A (zh) * 2020-08-14 2020-11-27 佳都新太科技股份有限公司 地铁站台闸门的检测方法、装置、计算机设备和存储介质
CN111999580B (zh) * 2020-08-14 2023-12-01 佳都科技集团股份有限公司 地铁站台闸门的检测方法、装置、计算机设备和存储介质

Similar Documents

Publication Publication Date Title
CN110472195B (zh) 一种基于区间监测数据的用户谐波危害预警评估方法
CN108199795A (zh) 一种设备状态的监测方法和装置
KR20200128232A (ko) 전력 수요 예측 장치 및 그 방법
CN111984705B (zh) 一种电能质量监测装置的精度评估管理方法和装置
CN116308306B (zh) 一种基于5g的新能源场站智慧管理系统及方法
CN116805068A (zh) 一种基于大数据的故障监测预警系统及方法
CN111931849A (zh) 一种水电机组运行数据趋势预警方法
Qu et al. A parameter and resolution adaptive algorithm for rapid detection of ramp events in different timescale databases of the power system
CN117933974B (zh) 水质数据采集监测设备的远程维护方法及系统
CN117439256A (zh) 一种基于物联网的电站设备管理方法及系统
Ren et al. Online anomaly detection using machine learning and HPC for power system synchrophasor measurements
CN108376293A (zh) 一种基于模糊数学改进的层次分析法的zj17卷烟设备维修智能预测方法
CN110930035A (zh) 一种基于区间统计值的谐波危害评估方法
CN112949201B (zh) 一种风速预测方法、装置、电子设备及存储介质
CN114548493A (zh) 一种电能表电流过载预测方法与系统
CN116756505B (zh) 一种基于大数据的光伏设备智能管理系统及方法
CN113822533A (zh) 基于实时事件驱动的风险评估量化模型构建方法及系统
CN117494009A (zh) 基于绝缘材料热解分析的电气设备状态评估方法及云平台
CN107368979A (zh) 一种天气敏感型区域配电网扰动事件预测方法及装置
CN110674242A (zh) 一种实现配电网线损数据分析和评估降损系统的方法
CN116430272A (zh) 基于电气扰动的煤矿设备异常诊断方法、系统及电子设备
CN118193959B (zh) 一种电涌能量的监测方法及系统
CN118211837B (zh) 一种智能电表的性能评估方法及装置
CN117113267B (zh) 基于大数据的预测模型训练方法、光伏发电性能检测方法
CN118152970B (zh) 一种基于边缘计算算法的设备状态趋势感知方法

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200327