CN102879681B - 一种闪变实时计算方法 - Google Patents

一种闪变实时计算方法 Download PDF

Info

Publication number
CN102879681B
CN102879681B CN201210363581.6A CN201210363581A CN102879681B CN 102879681 B CN102879681 B CN 102879681B CN 201210363581 A CN201210363581 A CN 201210363581A CN 102879681 B CN102879681 B CN 102879681B
Authority
CN
China
Prior art keywords
array
flickering
measured data
integers
real
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
CN201210363581.6A
Other languages
English (en)
Other versions
CN102879681A (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.)
Guodian Nanjing Automation Co Ltd
Original Assignee
Guodian Nanjing Automation 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 Guodian Nanjing Automation Co Ltd filed Critical Guodian Nanjing Automation Co Ltd
Priority to CN201210363581.6A priority Critical patent/CN102879681B/zh
Publication of CN102879681A publication Critical patent/CN102879681A/zh
Application granted granted Critical
Publication of CN102879681B publication Critical patent/CN102879681B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Complex Calculations (AREA)
  • Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)

Abstract

本发明公开了一种闪变实时计算方法，用于电能质量在线监测装置的闪变测量计算中，包括以下步骤：测量闪变瞬时视感度数据，将这些数据设为大于零的浮点数，称为被测数据浮点数集；将该被测数据浮点数集中的每个浮点数乘以一个系数，得到一组被测数据整数集；设定一组统计数组；将统计数组中的每一个元素清零，而后遍历被测数据整数集，每得到一个整数值，则统计数组中的相应变量加1；遍历过程结束，得到闪变瞬时视感度数据的概率。本发明提供了一种快速的闪变计算方法，使得闪变测量的实时计算得以完成，可满足对闪变指标的实时掌握，满足了对电网供电质量的实时质量监测。

Description

一种闪变实时计算方法
技术领域
本发明涉及一种电能质量在线监测装置闪变计算方法。
背景技术
闪变是电能质量监测中的一项重要的技术指标，其是电压波动引起的一种现象，可反映电压波动情况，因此，对于闪变的测量、计算的研究越来越多。对于闪变计算中的瞬时视感度的概率分布计算，目前常用的方法是分段线性插值法和统计排序法。按照IEC推荐的闪变测量方法，每10分钟计算一次瞬时闪变，如果按照800Hz的采样频率，10分钟内累积的瞬时视感度数据将达到480000点，使用以上的常规方法是很难在短时间内完成计算的，无法满足对闪变指标的实时掌握，不能满足对电网供电质量的实时质量监测。
发明内容
为了解决现有闪变计算技术中瞬时视感度数据的计算的实时性问题，本发明提供了一种既满足了精确性，又满足了实时性的闪变计算方法。
为了达到上述目的，本发明所采取的技术方案是：
一种闪变实时计算方法，其特征在于包括以下步骤：
（1）、测量闪变瞬时视感度数据，将这些数据均为大于零的浮点数称之为被测数据浮点数集。
（2）、将该被测数据浮点数集中的每个浮点数乘以一个系数，得到一组整数集称之为被测数据整数集，对被测数据浮点数集的概率统计实际上就等同于对被测数据整数集的概率统计。
（3）、设定一组整数数组称之为统计数组。
（4）、对被测数据整数集进行概率统计时，先将统计数组中的每一个变量清零，而后遍历被测数据整数集，每得到一个整数值，则统计数组中的相应变量加1。
（5）、遍历过程结束，则统计数组的每一个变量的数值就代表了被测数据整数集中相应整数在整个被测数据整数集中出现的次数，得到闪变瞬时视感度数据的概率。
前述的一种闪变实时计算方法，其特征在于包括以下步骤：
（1）、设需要统计的数据均为大于零的浮点数{Xn}称之为被测数据浮点数集，其中n={0，1，2，。。。，N},N为需要统计数据的数据点个数，确定{Xn}的测量精度即需要保留的小数点位数M。
（2）、将被测数据浮点数集{Xn}乘以10M转变成整数{Yn}称之为被测数据整数集，对{Xn}的概率统计实际上就是对{Yn}的概率统计。
（3）、在对{Yn}进行概率统计时，先估算{Yn}的最大值Ymax，根据该最大值Ymax，设定一整数数组sort[count]称为统计数组，统计数组的大小count>Ymax。
（4）、每次统计开始，先将统计数组sort[count]清零，而后从n=0到n=N，遍历{Yn}，每得到一个Yn，则相应的统计数组中的变量sort[Yn]加1。
（5）、该过程结束，则统计数组sort[count]中的每一个变量sort[y](0≤y≤Ymax)就代表了数值y在{Yn}中出现的次数m，m/N就为数值y在{Yn}中出现的概率，如y在{Xn}中对应得数值为x，则m/N也就是x在{Xn}中出现的概率。
本发明的有益效果是：
本发明可在很短的时间内完成闪变数据的计算，并且不受统计数据元素个数的影响，能实现大容量闪变测量数据的实时统计计算。可满足对闪变指标的实时掌握，满足了对电网供电质量的实时质量监测。
具体实施方式
以下对本发明的实施例作具体的介绍：
（1）、设x(n)＝f(n)，0≤n≤N，0≤x(n)≤Xmax。
x(n)为一离散浮点数数列，代表检测仪器检测过程中的中间计算结果，每得到N个点要对x(n)作一次概率统计计算。
假定该序列的最大值Xmax＜99，需要保留的小数点位数为2。
（2）、设y(n)＝x(n)*102，0≤n≤N。y(n)定义为整数集。
则0≤y(n)≤9999。于是就得到了与x(n)对应得整数集y(n)。
（3）、设定一组整数数组sort[count]，count=10000。
（4）、实际计算中，装置每采样一次就计算得到一个点x(n)。
当n=0时，sort[count]数组清零。
y(0)＝x(0)*102，sort[y(0)]=sort[y(0)]+1；
以后每得到一个x(n)，就得到一个y(n)，相应的sort[y(n)]＝sort[y(n)]+1。
当n=N-1时，进行概率统计。
如要得到浮点数x在x(n)中出现的概率，先将x变成整数，得y＝x*102。
则该点概率值为：sort[y]/N。
同样，可对其他数值进行概率统计。
以上显示和描述了本发明的基本原理、主要特征及优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界。

Claims (2)

1. 一种闪变实时计算方法，其特征在于，包括以下步骤：
（1）、测量闪变瞬时视感度数据，将这些数据设为大于零的浮点数，称为被测数据浮点数集；
（2）、将该被测数据浮点数集中的每个浮点数乘以一个系数，得到一组整数集，称之为被测数据整数集；
（3）、设定一组整数数组称之为统计数组；
（4）、将统计数组中的每一个元素清零，而后遍历被测数据整数集，每得到一个整数值，则统计数组中的相应变量加1；
（5）、上述遍历过程结束，则统计数组的每一个变量的数值就代表了被测数据整数集中相应整数在整个被测数据整数集中出现的次数，得到闪变瞬时视感度数据的概率。
2.根据权利要求1所述的一种闪变实时计算方法，其特征在于，包括以下步骤：
（1）、设闪变瞬时视感度数据均为大于零的浮点数{Xn}称之为被测数据浮点数集，其中n=0，1，2，…，N,N为需要统计数据的数据点个数，确定{Xn}的测量精度即需要保留的小数点位数M；
（2）、将被测数据浮点数集{Xn}乘以10M 转变成整数{Yn}称之为被测数据整数集；
（3）、在对{Yn}进行概率统计时，先估算{Yn}的最大值Ymax，根据该最大值Ymax，设定一组整数数组sort[count]，称为统计数组，统计数组的大小count>Ymax；
（4）、每次统计开始，先将统计数组sort[count]清零，而后从n=0到n=N，遍历{ Yn}，每得到一个Yn，则相应的统计数组中的变量sort[Yn]加1；
（5）、上述过程结束，则统计数组sort[count]中的每一个变量sort[y]就代表了数值y在{Yn}中出现的次数m，0≤y≤Ymax，m/N就为数值y在{Yn}中出现的概率，如y在{Xn}中对应得数值为x，则m/N也就是x在{Xn}中出现的概率。
CN201210363581.6A 2012-09-26 2012-09-26 一种闪变实时计算方法 Active CN102879681B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210363581.6A CN102879681B (zh) 2012-09-26 2012-09-26 一种闪变实时计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210363581.6A CN102879681B (zh) 2012-09-26 2012-09-26 一种闪变实时计算方法

Publications (2)

Publication Number Publication Date
CN102879681A CN102879681A (zh) 2013-01-16
CN102879681B true CN102879681B (zh) 2015-08-26

Family

ID=47481077

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210363581.6A Active CN102879681B (zh) 2012-09-26 2012-09-26 一种闪变实时计算方法

Country Status (1)

Country Link
CN (1) CN102879681B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103543317B (zh) * 2013-09-26 2016-06-01 光一科技股份有限公司 用于测量iec标准电压闪变的方法
CN103713223B (zh) * 2014-01-03 2016-01-27 烟台东方威思顿电气股份有限公司 一种自适应量程的低存储数据量闪变测量方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288807A (zh) * 2011-05-14 2011-12-21 苏州大学 一种测量电网电压闪变的方法
CN102565879A (zh) * 2011-08-04 2012-07-11 南京信息工程大学 一种风资料的处理方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6580659B1 (en) * 2000-08-25 2003-06-17 Micron Technology, Inc. Burst read addressing in a non-volatile memory device

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288807A (zh) * 2011-05-14 2011-12-21 苏州大学 一种测量电网电压闪变的方法
CN102565879A (zh) * 2011-08-04 2012-07-11 南京信息工程大学 一种风资料的处理方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
交流电弧炉的电压闪变及计算;刘诚等;《工业加热》;20120430;第41卷(第04期);第1-4页 *
基于FFT的电压波动与闪变测量算法;史三省等;《电力***及其自动化学报》;20101231;第22卷(第06期);第109-112页 *
电压闪变的评价与工程计算方法研究;余健明等;《西安理工大学学报》;19970131;第13卷(第01期);第46-50页 *
郭晓丽等.短时间闪变严重度Pst的计算方法的研究.《南通工学院学报(自然科学版)》.2004,第03卷(第03期), *

Also Published As

Publication number Publication date
CN102879681A (zh) 2013-01-16

Similar Documents

Publication Publication Date Title
CN105258664B (zh) 输电线路等值覆冰厚度获得方法
CN110031680A (zh) 一种***侧谐波阻抗估计方法和***
CN102879639A (zh) 一种电力***中频率的实时测量方法
CN103996084B (zh) 基于纵向时刻Markov链模型的风功率概率预测方法
CN109284933B (zh) 一种基于数理统计的电子式互感器状态评估***及方法
CN105098763A (zh) 一种基于wams和ems的暂态电压稳定在线评估方法
WO2015061227A1 (en) Performance analysis of power grid monitors
US8680842B2 (en) Method for measurement of total harmonic distortion
CN106383501A (zh) 对碳排放过程的监测控制方法
CN103424734A (zh) 电子式电能表在线校准方法
CN106526312A (zh) 基于r‑v (ⅲ)窗fft双峰插值的电能计量方法
CN105093166B (zh) 一种电子式电能表的现场检验方法
CN105469156A (zh) Moa健康管理和故障预测方法及***
CN108490257A (zh) 一种基于频谱曲线拟合的短时窗间谐波测量方法
CN106405280B (zh) 一种智能变电站在线监测参数趋势预警方法
CN104931541B (zh) 一种混凝土中gfrp筋的长期性能监测方法
Artale et al. Measurement of simplified single-and three-phase parameters for harmonic emission assessment based on IEEE 1459-2010
CN103472373B (zh) 一种评价自然污秽与人工污秽绝缘子外绝缘等价性的方法
CN102879681B (zh) 一种闪变实时计算方法
KR100576511B1 (ko) 시계열 데이터를 이용한 전력계통 실시간 전압안정도위험지수 계산 시스템 및 방법
CN105093163A (zh) 一种电子式电能表的误差调试方法
Jiang et al. Assessment method of harmonic emission level based on the improved weighted support vector machine regression
CN108415880A (zh) 一种基于样本熵与小波变换的线损特性分析方法
CN105158555A (zh) 一种配电线路升压运行节约电力电量的测量和验证方法
CN103576120A (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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20130116

Assignee: NANJING GUODIAN NANZI RURAL POWER GRID AUTOMATION ENGINEERING CO., LTD.

Assignor: Nanjing Automation Co., Ltd., China Electronics Corp.

Contract record no.: 2013320000276

Denomination of invention: Flickering real-time calculation method

License type: Exclusive License

Record date: 20130401

LICC Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model
C14 Grant of patent or utility model
GR01 Patent grant
EC01 Cancellation of recordation of patent licensing contract

Assignee: NANJING GUODIAN NANZI RURAL POWER GRID AUTOMATION ENGINEERING CO., LTD.

Assignor: Nanjing Automation Co., Ltd., China Electronics Corp.

Contract record no.: 2013320000276

Date of cancellation: 20181016

EC01 Cancellation of recordation of patent licensing contract