CN113533831A - Rapid extraction method for flicker - Google Patents

Rapid extraction method for flicker Download PDF

Info

Publication number
CN113533831A
CN113533831A CN202110678396.5A CN202110678396A CN113533831A CN 113533831 A CN113533831 A CN 113533831A CN 202110678396 A CN202110678396 A CN 202110678396A CN 113533831 A CN113533831 A CN 113533831A
Authority
CN
China
Prior art keywords
flicker
value
voltage
fluctuation
time
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
CN202110678396.5A
Other languages
Chinese (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.)
Electric Power Research Institute of State Grid Anhui Electric Power Co Ltd
Original Assignee
Electric Power Research Institute of State Grid Anhui 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 Electric Power Research Institute of State Grid Anhui Electric Power Co Ltd filed Critical Electric Power Research Institute of State Grid Anhui Electric Power Co Ltd
Priority to CN202110678396.5A priority Critical patent/CN113533831A/en
Publication of CN113533831A publication Critical patent/CN113533831A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/12—Measuring rate of change
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Eye Examination Apparatus (AREA)

Abstract

The invention relates to a flicker rapid extraction method, which comprises the following steps: measuring a fluctuation component; subtracting the carrier effective value from the square root value of the cycle refreshing voltage, dividing the carrier effective value by the nominal voltage to obtain a flicker signal, and performing FFT (fast Fourier transform) on the flicker signal to obtain voltage fluctuation values of different frequency components in the fluctuation signal; calculating a visual sensitivity frequency characteristic coefficient; obtaining a fluctuation value of any frequency in a flicker perception frequency range according to the visual sensation frequency characteristic coefficient; the sine wave voltage fluctuation value of any frequency is equivalent to the voltage fluctuation of 8.8 Hz; the voltage fluctuation component is equivalent to a voltage fluctuation value of 8.8Hz to calculate a voltage flicker value; calculating a short-time flicker statistic value within a preset time according to the voltage flicker value; the statistical value of the long-time flicker is obtained by calculating the short-time flicker value contained in the measuring time period. The flicker extraction method is simple, the flicker value can be obtained only by calculating the fluctuation component, and the calculation amount is greatly reduced.

Description

Rapid extraction method for flicker
Technical Field
The invention relates to a rapid extraction method of flicker.
Background
The voltage fluctuation can cause obvious flicker of incandescent light, and the incandescent light is intolerable to human eyes when the voltage fluctuation is serious, but the electric light source can not generate the voltage fluctuation, and the lamp light can feel the flicker of the lamp light.
Disclosure of Invention
The technical problem to be solved by the invention is to provide a rapid extraction method for flicker, so as to overcome the defects in the prior art.
The technical scheme for solving the technical problems is as follows: a flash rapid extraction method comprises the following steps:
s100, measuring a fluctuation component;
s200, subtracting a carrier effective value from a square root value of a cycle refreshing voltage, dividing by a nominal voltage to obtain a flicker signal, and performing FFT (fast Fourier transform) on the flicker signal to obtain voltage fluctuation values of different frequency components in the fluctuation signal;
s300, calculating a visual sensitivity frequency characteristic coefficient;
s400, obtaining a fluctuation value of any frequency in a flicker perception frequency range according to the visual sensation frequency characteristic coefficient;
s500, enabling the sine wave voltage fluctuation value of any frequency to be equivalent to the voltage fluctuation of 8.8 Hz;
s600, enabling the voltage fluctuation component to be equivalent to a voltage fluctuation value of 8.8Hz, and calculating a voltage flicker value;
s700, calculating a short-time flicker statistic value within a preset time according to the voltage flicker value;
and S800, calculating the statistical value of the long-time flicker from the short-time flicker value contained in the measurement time period.
On the basis of the technical scheme, the invention can be further improved as follows.
Further, the square root of the cycle refresh voltage:
Figure BDA0003121769550000021
in the formula: n is the number of sampling points in each period; u (i) is the instantaneous value of the voltage waveform sampled the ith time;
k is the calculated window number (k is 1, 2, 3 …), i.e. the first value is obtained in one cycle, and the next value is calculated from sample 0.5N +1 to sample 0.5N + N.
Further, the time window has a width of 6144 × 2, corresponding to a duration of 122.88s, fundamental frequency f1＝0.00814Hz。
Further, the calculation formula of the visual sensation frequency characteristic coefficient is as follows:
Figure BDA0003121769550000022
further, a fluctuation value of an arbitrary frequency in the flicker perception frequency range at s (t) of 1 can be obtained from the visibility coefficient.
Further, s (t) is proportional to the voltage variation:
Figure BDA0003121769550000023
Figure BDA0003121769550000031
when d iss,8.8When 0.25, s (t) is 1;
Figure BDA0003121769550000032
Figure BDA0003121769550000033
for periodically stable rectangular voltage fluctuations:
P0.1＝P1＝P3＝P10＝P50＝S(t)
Figure BDA0003121769550000034
the obtained voltage fluctuation component is equivalent to a voltage fluctuation value of 8.8Hz to calculate a voltage flicker value:
Pst＝2.856ds,8.8
according to voltage flashVariable value PstCalculating the statistical value P of short-time flicker every 10minst.10；
Figure BDA0003121769550000035
Statistical value P of long-time flickerstShort-time flicker value P from 10min contained in the measurement periodst.10And (3) calculating to obtain:
Figure BDA0003121769550000036
wherein, Pst.10kIs the kth short-time flicker value in the measurement time.
The invention has the beneficial effects that: the flicker extraction method is simple, the flicker value can be obtained only by calculating the fluctuation component, the design of some filters is reduced, and the calculation amount is greatly reduced.
Drawings
FIG. 1 is a graph of a perceptibility factor frequency characteristic;
fig. 2 is a simulated sensitivity coefficient frequency characteristic curve.
Detailed Description
The principles and features of this invention are described below in conjunction with the following drawings, which are set forth by way of illustration only and are not intended to limit the scope of the invention.
A flash rapid extraction method comprises the following steps:
s100, measuring a fluctuation component;
s200, subtracting a carrier effective value from a square root value of a cycle refreshing voltage, dividing by a nominal voltage to obtain a flicker signal, and performing FFT (fast Fourier transform) on the flicker signal to obtain voltage fluctuation values of different frequency components in the fluctuation signal;
s300, calculating a visual sensitivity frequency characteristic coefficient;
s400, obtaining a fluctuation value of any frequency in a flicker perception frequency range according to the visual sensation frequency characteristic coefficient;
s500, enabling the sine wave voltage fluctuation value of any frequency to be equivalent to the voltage fluctuation of 8.8 Hz;
s600, enabling the voltage fluctuation component to be equivalent to a voltage fluctuation value of 8.8Hz, and calculating a voltage flicker value;
s700, calculating a short-time flicker statistic value within a preset time according to the voltage flicker value;
and S800, calculating the statistical value of the long-time flicker from the short-time flicker value contained in the measurement time period.
Square root of cycle refresh voltage:
Figure BDA0003121769550000041
in the formula: n is the number of sampling points in each period; u (i) is the instantaneous value of the voltage waveform sampled the ith time;
k is the calculated window number (k is 1, 2, 3 …), i.e. the first value is obtained in one cycle, and the next value is calculated from sample 0.5N +1 to sample 0.5N + N.
The time window has a width of 6144 × 2 points, corresponding to a duration of 122.88s, and a fundamental frequency of f1＝0.00814Hz。
The calculation formula of the visual sensation frequency characteristic coefficient is as follows:
Figure BDA0003121769550000051
from the visual sensation coefficient, a fluctuation value of an arbitrary frequency in the flicker perception frequency range at s (t) 1 can be obtained.
S (t) is proportional to the voltage variation:
Figure BDA0003121769550000052
Figure BDA0003121769550000053
when d iss,8.8When 0.25, s (t) is 1;
Figure BDA0003121769550000054
Figure BDA0003121769550000055
for periodically stable rectangular voltage fluctuations:
P0.1＝P1＝P3＝P10＝P50＝S(t)
Figure BDA0003121769550000056
the obtained voltage fluctuation component is equivalent to a voltage fluctuation value of 8.8Hz to calculate a voltage flicker value:
Pst＝2.856ds,8.8
according to the voltage flash value PstCalculating the statistical value P of short-time flicker every 10minst.10；
Figure BDA0003121769550000057
Statistical value P of long-time flickerstShort-time flicker value P from 10min contained in the measurement periodst.10And (3) calculating to obtain:
Figure BDA0003121769550000058
wherein, Pst.10kIs the kth short-time flicker value in the measurement time.
The flicker calculation method is simple, and the flicker value can be obtained only by calculating the fluctuation component. The calculation process shows that the method reduces the design of some filters and greatly reduces the calculation amount. In order to verify the accuracy of the method, the next section performs simulation verification on the method.
The flicker improvement algorithm uses a sensitivity frequency characteristic, and fig. 1 is a sensitivity frequency characteristic curve obtained from experimental data under a sinusoidal voltage fluctuation condition. As can be seen from fig. 1, when s (t) is 1, there is a minimum value of voltage fluctuation d (%) corresponding to the maximum sensitive frequency of flicker of 8.8Hz (see table 1, 8.8Hz corresponds to a minimum value of voltage fluctuation d (%) of 0.250). So that K (f) is less than or equal to 1;
from theoretical analysis of the mechanism of the flickering visual sensation, a strict mathematical expression method is necessary. The visual sensitivity frequency characteristic curve is approximated by 5 straight lines and asymptotes, namely the sum of logarithmic amplitude-frequency characteristics of 5 typical control links is used for representing, and the result is as follows:
Figure BDA0003121769550000061
further, the expression of the transfer function of the lamp-eye-brain link can be derived:
Figure BDA0003121769550000062
K＝1.74802,λ＝2π4.05981,ω1＝2π9.15494,ω2＝2π2.27979,ω3＝2π1.22535,ω4＝2π21.9
the visual sensation coefficient curve obtained by simulation is shown in fig. 2;
the voltage fluctuation values for a sine wave with apparent sensitivity s (t) of 1 perceived unit are shown in table 1. The partial visual perception frequency characteristic coefficients in the general range of flicker perceptibility are listed:
TABLE 1 Unit apparent sensitivity response of standard flicker instrument to sinusoidal wave voltage fluctuation
Figure BDA0003121769550000063
Figure BDA0003121769550000071
Although embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present invention, and that variations, modifications, substitutions and alterations can be made to the above embodiments by those of ordinary skill in the art within the scope of the present invention.

Claims (6)

1. A flicker rapid extraction method is characterized by comprising the following steps:
s100, measuring a fluctuation component;
s200, subtracting a carrier effective value from a square root value of a cycle refreshing voltage, dividing by a nominal voltage to obtain a flicker signal, and performing FFT (fast Fourier transform) on the flicker signal to obtain voltage fluctuation values of different frequency components in the fluctuation signal;
s300, calculating a visual sensitivity frequency characteristic coefficient;
s400, obtaining a fluctuation value of any frequency in a flicker perception frequency range according to the visual sensation frequency characteristic coefficient;
s500, enabling the sine wave voltage fluctuation value of any frequency to be equivalent to the voltage fluctuation of 8.8 Hz;
s600, enabling the voltage fluctuation component to be equivalent to a voltage fluctuation value of 8.8Hz, and calculating a voltage flicker value;
s700, calculating a short-time flicker statistic value within a preset time according to the voltage flicker value;
and S800, calculating the statistical value of the long-time flicker from the short-time flicker value contained in the measurement time period.
2. The flicker rapid extraction method according to claim 1, wherein: square root of cycle refresh voltage:
Figure FDA0003121769540000011
in the formula: n is the number of sampling points in each period; u (i) is the instantaneous value of the voltage waveform sampled the ith time;
k is the calculated window number (k is 1, 2, 3 …), i.e. the first value is obtained in one cycle, and the next value is calculated from sample 0.5N +1 to sample 0.5N + N.
3. The flicker rapid extraction method according to claim 2, wherein: the time window has a width of 6144 × 2 points, corresponding to a duration of 122.88s, and a fundamental frequency of f1＝0.00814Hz。
4. A flicker rapid extraction method according to claim 3, wherein: the calculation formula of the visual sensation frequency characteristic coefficient is as follows:
Figure FDA0003121769540000021
5. the flicker rapid extraction method according to claim 4, wherein: from the visual sensation coefficient, a fluctuation value of an arbitrary frequency in the flicker perception frequency range at s (t) 1 can be obtained.
6. The flash rapid extraction method according to claim 5, characterized in that: s (t) is proportional to the voltage variation:
Figure FDA0003121769540000022
Figure FDA0003121769540000023
when d iss,8.8When 0.25, s (t) is 1;
Figure FDA0003121769540000024
Figure FDA0003121769540000025
for periodically stable rectangular voltage fluctuations:
P0.1＝P1＝P3＝P10＝P50＝S(t)
Figure FDA0003121769540000026
the obtained voltage fluctuation component is equivalent to a voltage fluctuation value of 8.8Hz to calculate a voltage flicker value:
Pst＝2.856ds,8.8
according to the voltage flash value PstCalculating the statistical value P of short-time flicker every 10minst.10；
Figure FDA0003121769540000027
Statistical value P of long-time flickerstShort-time flicker value P from 10min contained in the measurement periodst.10And (3) calculating to obtain:
Figure FDA0003121769540000031
wherein, Pst.10kIs the kth short-time flicker value in the measurement time.
CN202110678396.5A 2021-06-18 2021-06-18 Rapid extraction method for flicker Pending CN113533831A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110678396.5A CN113533831A (en) 2021-06-18 2021-06-18 Rapid extraction method for flicker

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110678396.5A CN113533831A (en) 2021-06-18 2021-06-18 Rapid extraction method for flicker

Publications (1)

Publication Number Publication Date
CN113533831A true CN113533831A (en) 2021-10-22

Family

ID=78125232

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110678396.5A Pending CN113533831A (en) 2021-06-18 2021-06-18 Rapid extraction method for flicker

Country Status (1)

Country Link
CN (1) CN113533831A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114325022A (en) * 2021-11-24 2022-04-12 浙江中控技术股份有限公司 Method, system, equipment and medium for monitoring AO sinusoidal signal jump

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288807A (en) * 2011-05-14 2011-12-21 苏州大学 Method for measuring electric network voltage flicker
CN107271753A (en) * 2017-07-27 2017-10-20 华北电力大学 Voltage flicker detection and device
CN107402362A (en) * 2016-05-20 2017-11-28 西门子电力自动化有限公司 The method of testing and device of flickering

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288807A (en) * 2011-05-14 2011-12-21 苏州大学 Method for measuring electric network voltage flicker
CN107402362A (en) * 2016-05-20 2017-11-28 西门子电力自动化有限公司 The method of testing and device of flickering
CN107271753A (en) * 2017-07-27 2017-10-20 华北电力大学 Voltage flicker detection and device

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李洪池;赵不贿;: "基于DSP的电压闪变监测系统设计与实现", 低压电器, no. 19, pages 34 - 37 *
陈亮: "基于DSP的电能质量在线监测装置设计", 《硕士论文 工程科技II辑 》, pages 10 - 14 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114325022A (en) * 2021-11-24 2022-04-12 浙江中控技术股份有限公司 Method, system, equipment and medium for monitoring AO sinusoidal signal jump
CN114325022B (en) * 2021-11-24 2024-04-12 浙江中控技术股份有限公司 Method, system, equipment and medium for monitoring AO sinusoidal signal jump

Similar Documents

Publication Publication Date Title
KR100708036B1 (en) Method of characterizing a device under test
Ritter et al. Guideline for an assessment of electrochemical noise measurement devices
JP4783902B2 (en) Glare evaluation device and glare evaluation program
CN113533831A (en) Rapid extraction method for flicker
White et al. A discrete Matlab–Simulink flickermeter model for power quality studies
CN110940877A (en) Voltage flicker detection method based on windowed FFT interpolation algorithm
Gallo et al. Design and calibration of an objective flickermeter
CN109752614B (en) Lightning arrester parameter measurement method and system
CN104155542B (en) Detection method suitable for flickering generated by high-frequency inter-harmonics
CN112180177B (en) Power frequency electromagnetic field evaluation method and system fusing measured data
CN113094983A (en) Online simulation method for multi-dimensional time-varying characteristics of direct-current fault electric arc of photovoltaic system
Gutierrez et al. An alternative strategy to improve the flicker severity measurement
Gallo et al. On the Measurement of Power Electronic Devices' Frequency Coupling Admittance
Harris The Level Crossing Method applied to mean wind speeds from “mixed” climates
Leger et al. Automated system for determining frequency dependent parameter model of transmission line in a laboratory environment
EP3098614B1 (en) A method and device for flicker measurements performed on electrical appliances
CN107402362B (en) Flicker testing method and device
Macii et al. Fast detection of rapid voltage change events through dynamic RMS voltage tracking
Kukačka et al. Confidence intervals for luminous flicker measurements: Comparison of various approaches
Matthews et al. Evaluation of flicker measurement uncertainties by a Monte Carlo method
Mihai Statistical tools for analysis of the performance of a high resistance measurement bridge
KR100782450B1 (en) Weighting filters for 220voltage which using of flickermeter
Sansheng et al. Research on flicker measurement algorithm based on FFT
Chen et al. A Calibration test prototype of flickermeter for wind turbines
Velicka et al. The influence of fluctuation in the luminous flux of light sources and its quantification

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