CN112986675A - Method for detecting frequency change rate - Google Patents
Method for detecting frequency change rate Download PDFInfo
- Publication number
- CN112986675A CN112986675A CN202110117906.1A CN202110117906A CN112986675A CN 112986675 A CN112986675 A CN 112986675A CN 202110117906 A CN202110117906 A CN 202110117906A CN 112986675 A CN112986675 A CN 112986675A
- Authority
- CN
- China
- Prior art keywords
- frequency
- value
- change rate
- delta
- 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.)
- Granted
Links
Images
Classifications
-
- 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)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
The invention discloses a method for detecting frequency change rate, which comprises the following steps: firstly, judging whether the detection method of the invention is operated for the first time, and initializing the frequency basic quantity f according to the first and non-first assignmentnfThe current frequency value f is detected and assigned as the frequency basis quantity fnf(ii) a Calculating a frequency change rate value delta f/delta t, and performing amplitude limiting treatment on the calculated value; calculating to obtain the frequency variation delta fnfBy the amount of frequency change Δ fnfFor frequency fundamental quantity fnfUpdating is carried out; judging whether the value of the frequency change rate delta f/delta t exceeds a limit value, if so, accumulating the duration time, otherwise, clearing the duration time; judging whether the duration time exceeds the error judgment window time, and if the judgment result is true, triggering the fault; otherwise, the detection is finished; the method for detecting the frequency change rate is simple and easy to implement, can be used as a general method in a passive island detection method, and can also be used in combination with other active island detection.
Description
Technical Field
The invention relates to the field of anti-islanding detection of distributed power generation systems, in particular to a frequency change rate detection method.
Background
When the distributed power generation system is in grid-connected operation with a power grid, the frequency is basically unchanged, when an island is formed, serious power unbalance may exist between a power supply and a load, and the frequency of the system can be changed, so that the island can be detected by testing the change rate of the frequency. The section, e.g. LoM in the EN50549-1 safety code, explicitly states that a distributed power generation system should include a frequency rate of change detection function, but is not limited to including this function only. At present, a general frequency change rate detection method is not described, but a detection method using different frequency change rates directly influences the detection of the functions of the inverter, and the inverter may be frequently tripped due to detection errors.
Therefore, in combination with the above-mentioned technical problems, there is a need to provide a new technical solution.
Disclosure of Invention
In order to solve the technical problems in the prior art, the invention provides a general frequency change rate detection method, in particular to a frequency change rate-based detection method in a passive island detection method, and the specific technical scheme is as follows:
s1: reading the value of an initial mark, and if the initial mark is 0, operating the detection method for the first time;
s2: initializing a frequency basis quantity fnfDetecting the current frequency value f and assigning the current frequency value f as the frequency basic quantity fnfWhile the initial flag is assigned a value of 1;
s3: the frequency rate of change value deltaf/deltat is calculated according to the following formula,
wherein Δ f ═ f1-f0＝f1-fnf， 
Wherein, the time interval Δ t is obtained by a sampling count quantity rocef _ counter and a sampling frequency rocef _ adc _ freq;
s4: setting a limit value for the frequency change rate calculation value Δ f/Δ t calculated in step S3;
s5: the frequency variation Δ f is calculated by the following equationnf，
Change amount of frequency delta fnfSubstituting into the frequency basis quantity formula fnf＝fnf+ΔfnfObtaining a frequency basis quantity fnfAnd for the frequency base quantity fnfUpdating is carried out;
s6: judging whether the value of the frequency change rate delta f/delta t calculated and obtained in the step S3 exceeds the limit value in the step S4, wherein the value exceeding the limit value is recorded as true, and the value not exceeding the limit value is recorded as false;
s7: if the judgment result of the step S6 is true, accumulating the duration time;
s8: judging whether the duration time in the step S7 exceeds the error window time, recording the time exceeding the error window as true, and recording the time not exceeding the error window as false;
s9: and if the judgment result of the step S8 is true, triggering the fault.
Preferably, in step S1, if the initial flag is 1, and the detection method is not operated for the first time, steps S3, S4, S5, S6, S7, S8, and S9 are sequentially performed.
Preferably, in step S7, the duration time is cleared if the determination result in step S6 is false.
Preferably, in step S9, if the determination result in step S8 is false, the detection ends.
Preferably, in step S3, the sampling count and the sampling frequency are sampled by a digital-to-analog converter adc.
Preferably, in step S3, the time interval Δ t is calculated by the adc sampling count obtained in the zero-crossing interrupt and the adc sampling frequency rocaf _ adc _ freq.
Preferably, in step S4, the limit value is a threshold value FreqChangeMax required in the EN safety specification.
Preferably, in step S4, the frequency change rate calculation value Δ f/Δ t is defined as follows:
preferably, in step S1, the mains frequency is sampled once per mains period, and the frequency change rate detection process is executed once.
The frequency change rate detection method has the following beneficial effects:
(1) the method for detecting the frequency change rate is simple in a passive island detection method, easy to implement, capable of serving as a general method and capable of being popularized and used;
(2) the frequency change rate detection method can be used in combination with other active anti-islanding detection.
Additional aspects and advantages of the invention will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the invention.
Drawings
In order to more clearly illustrate the technical solution of the present invention, the drawings needed for the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without creative efforts.
FIG. 1 is a schematic diagram of the frequency change rate detection method of the present invention;
fig. 2 is a flow chart of the frequency change rate detection method of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Examples
As shown in fig. 1, fig. 1 is a schematic diagram illustrating the principle of the frequency change rate detection method of the present invention, and the principle is as follows:
t1time relative to t0The frequency variation Δ f at a time and the time interval Δ t such asLower part
Δf＝f1-f0＝f1-fnf (1)
Δt＝t1-t0 (2)
t0And t1The interval is a commercial power grid period, and the time interval delta t is obtained by calculating an adc sampling count value Rocof _ counter and a sampling frequency Rocof _ adc _ freq which are acquired in zero-crossing interruption, and is as follows:
the frequency change rate Δ f/Δ t can be obtained from the formulas (1) and (3), and is represented by the following formula (4):
the calculated frequency change rate delta f/delta t is set with a limiting value FreqChangeMax, and the limiting rule is as follows as formulas (5-1), (5-2) and (5-3):
in an example, the limit value FreqChangeMax is a threshold value required in the EN safety specification.
The frequency change amount Δ f is calculated from the following equation (6)nf，
The frequency variation Δ f calculated by equation (6)nfSubstituting into equation (7) to obtain the frequency basis quantity fnfThe following are:
fnf＝fnf+Δfnf (7)
preferably, the commercial power frequency value f is sampled once every commercial power period, and a detection process of the frequency change rate is executed once, in an example, the commercial power frequency adopted is 50hz, so that 20ms obtains the primary commercial power frequency. Of course, the detection method provided by the invention is not limited to the mains frequency of 50hz, and is also applicable to other mains frequencies.
As shown in fig. 2, fig. 2 is a flowchart of the method for detecting a frequency change rate according to the present invention, and the flow steps are as follows:
s1: reading the value of an initial mark, and if the initial mark is 0, operating the detection method for the first time; (ii) a
In another embodiment, if the detection method is not run for the first time, step S2 is skipped, and steps S3, S4, S5, S6, S7, S8 and S9 are sequentially performed.
S2: initializing a frequency basis quantity fnfDetecting the current frequency value f and assigning the current frequency value f as the frequency basic quantity fnfWhile the initial flag is assigned a value of 1;
s3: calculating a frequency change rate delta f/delta t value according to the formula (4);
s4: setting a limit value for the frequency change rate Δ f/Δ t calculated in step S3; in the example, the limit value is a threshold value FreqChangeMax required in the EN safety specification;
s5: the frequency variation Δ f is calculated from equation (6)nfAmount of frequency change Δ f to be calculatednfCalculating the frequency basic quantity f by substituting formula (7)nfAnd for the frequency base quantity fnfUpdating is carried out;
s6: judging whether the value of the frequency change rate delta f/delta t calculated and obtained in the step S3 exceeds the limit value in the step S4, wherein the value exceeding the limit value is recorded as true, and the value not exceeding the limit value is recorded as false;
s7: if the judgment result of the step S6 is true, accumulating the duration time; in other embodiments, if the determination result of step S6 is false, the duration amount is cleared.
S8: judging whether the duration time of the step S7 exceeds the error window time, recording the time exceeding the error window as true, and recording the time not exceeding the error window as false;
s9: and if the judgment result of the step S8 is true, triggering the fault.
In other embodiments, in step S9, if the determination result in step S8 is false, the detection ends.
The method for detecting the frequency change rate has the advantages that:
the method for detecting the frequency change rate is simple and easy to implement, can be used as a general method in a passive island detection method, and can also be used in combination with other active island detection.
In the description of the specification, reference to the description of the term "one embodiment", "some embodiments", "an example", "a specific example" or "some examples", etc., means that a particular feature or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples described in this specification can be combined and combined by one skilled in the art.
While embodiments of the present invention have been shown and described above, it is to be understood that the above embodiments are exemplary and not to be construed as limiting the present invention, and that changes, modifications and variations may be made therein by those of ordinary skill in the art within the scope of the present invention.
Claims (9)
1. A method of detecting a rate of change of frequency, comprising the steps of:
s1: reading the value of an initial mark, and if the initial mark is 0, operating the detection method for the first time;
s2: initializing a frequency basis quantity fnfDetecting the current frequency value f and assigning the current frequency value f as the frequency basic quantity fnfWhile the initial flag is assigned a value of 1;
s3: the frequency rate of change value deltaf/deltat is calculated according to the following formula,
wherein Δ f ═ f1-f0＝f1-fnf，
Wherein, the time interval Δ t is obtained by a sampling count quantity rocef _ counter and a sampling frequency rocef _ adc _ freq;
s4: setting a limit value for the frequency change rate calculation value Δ f/Δ t calculated in step S3;
s5: the frequency variation Δ f is calculated by the following equationnf，
Change amount of frequency delta fnfSubstituting into the frequency basis quantity formula fnf＝fnf+ΔfnfObtaining a frequency basis quantity fnfAnd for the frequency base quantity fnfUpdating is carried out;
s6: judging whether the value of the frequency change rate delta f/delta t calculated and obtained in the step S3 exceeds the limit value in the step S4, wherein the value exceeding the limit value is recorded as true, and the value not exceeding the limit value is recorded as false;
s7: if the judgment result of the step S6 is true, accumulating the duration time;
s8: judging whether the duration time in the step S7 exceeds the error window time, recording the time exceeding the error window as true, and recording the time not exceeding the error window as false;
s9: and if the judgment result of the step S8 is true, triggering the fault.
2. The method of claim 1, wherein in step S1, if the initial flag is 1, and the detection method is not a first run, steps S3, S4, S5, S6, S7, S8, and S9 are sequentially performed.
3. The method according to claim 1, wherein in step S7, if the determination result in step S6 is false, the duration is cleared.
4. The method of claim 1, wherein in step S9, if the determination result in step S8 is false, the detecting ends.
5. The method according to claim 1, wherein in step S3, the sampling count and the sampling frequency are sampled by a digital-to-analog converter adc.
6. The method according to claim 5, wherein in step S3, the time interval Δ t is calculated by the adc sample count obtained in the zero-crossing interrupt and the adc sampling frequency Rocof _ adc _ freq.
7. The method of claim 1, wherein the limit value is a threshold value FreqChangeMax required in EN safety specification in step S4.
8. The method as claimed in claim 7, wherein in step S4, the criterion between the calculated frequency change rate Δ f/Δ t and the limiting value FreqChangeMax is as follows:
9. the method according to claim 1, wherein in step S1, the commercial power frequency is sampled once per commercial power period, and a frequency change rate detection procedure is performed once.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110117906.1A CN112986675B (en) | 2021-01-28 | 2021-01-28 | Method for detecting frequency change rate |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110117906.1A CN112986675B (en) | 2021-01-28 | 2021-01-28 | Method for detecting frequency change rate |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112986675A true CN112986675A (en) | 2021-06-18 |
| CN112986675B CN112986675B (en) | 2023-03-24 |
Family
ID=76345620
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110117906.1A Active CN112986675B (en) | 2021-01-28 | 2021-01-28 | Method for detecting frequency change rate |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112986675B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09243669A (en) * | 1996-03-07 | 1997-09-19 | Mitsubishi Electric Corp | Relay for detection of change rate |
| JPH09329630A (en) * | 1996-06-07 | 1997-12-22 | Fuji Electric Co Ltd | Frequency measuring method and frequency measuring circuit |
| EP0908735A2 (en) * | 1997-10-10 | 1999-04-14 | Daimler-Benz Aerospace AG | Method for determining the frequency of a signal |
| JP2008267917A (en) * | 2007-04-18 | 2008-11-06 | Mitsubishi Electric Corp | Frequency change part measuring device, frequency change rate measuring device and electric power system control protecting device |
| TW201133020A (en) * | 2010-03-18 | 2011-10-01 | O2Micro Inc | Method and device of inspecting carrier frequency of GPS satellite and GPS receiver |
| CN103383412A (en) * | 2013-07-10 | 2013-11-06 | 珠海许继芝电网自动化有限公司 | Adaptive software and hardware frequency tracking and sampling method |
| CN107422196A (en) * | 2017-02-27 | 2017-12-01 | 国电南瑞科技股份有限公司 | A kind of island detection method based on percent harmonic distortion and frequency |
| CN207502605U (en) * | 2017-09-07 | 2018-06-15 | 江汉大学 | A kind of frequency measuring system based on VCXO references |
-
2021
- 2021-01-28 CN CN202110117906.1A patent/CN112986675B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09243669A (en) * | 1996-03-07 | 1997-09-19 | Mitsubishi Electric Corp | Relay for detection of change rate |
| JPH09329630A (en) * | 1996-06-07 | 1997-12-22 | Fuji Electric Co Ltd | Frequency measuring method and frequency measuring circuit |
| EP0908735A2 (en) * | 1997-10-10 | 1999-04-14 | Daimler-Benz Aerospace AG | Method for determining the frequency of a signal |
| JP2008267917A (en) * | 2007-04-18 | 2008-11-06 | Mitsubishi Electric Corp | Frequency change part measuring device, frequency change rate measuring device and electric power system control protecting device |
| TW201133020A (en) * | 2010-03-18 | 2011-10-01 | O2Micro Inc | Method and device of inspecting carrier frequency of GPS satellite and GPS receiver |
| CN103383412A (en) * | 2013-07-10 | 2013-11-06 | 珠海许继芝电网自动化有限公司 | Adaptive software and hardware frequency tracking and sampling method |
| CN107422196A (en) * | 2017-02-27 | 2017-12-01 | 国电南瑞科技股份有限公司 | A kind of island detection method based on percent harmonic distortion and frequency |
| CN207502605U (en) * | 2017-09-07 | 2018-06-15 | 江汉大学 | A kind of frequency measuring system based on VCXO references |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112986675B (en) | 2023-03-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gupta et al. | Active ROCOF relay for islanding detection | |
| US5805395A (en) | Half-cycle summation V/Hz relay for generator and transformer over-excitation protection | |
| US6185482B1 (en) | System and method for rms overcurrent backup function | |
| US11828815B2 (en) | AC/DC leakage detection method | |
| Yin et al. | Islanding detection using proportional power spectral density | |
| CN111130074B (en) | Direct-current transmission line ultra-high-speed protection system and method based on voltage first extreme value time | |
| RU2765380C2 (en) | Method for identifying a failure event in a sector of an electrical distribution network | |
| Bastos et al. | Comparison of methods for determining inception and recovery points of voltage variation events | |
| CN109145825A (en) | A kind of coherent noise filtering method and system | |
| JP2019062730A (en) | System and method for detecting and evaluating oscillation in electrical power grid | |
| CN112986675B (en) | Method for detecting frequency change rate | |
| CN110247380A (en) | A kind of comprehensive criterion type boundary protection method and system for multiterminal element power grid | |
| US20020105818A1 (en) | Power conversion device with power source voltage drop detection unit | |
| CN104821598B (en) | A kind of control method of grid-connected inverter and control device | |
| CN113302835A (en) | Method for stabilizing a rotor of a wind turbine | |
| CN114583657B (en) | Relay protection starting method for improving reliability of single-CPU double AD sampling system | |
| CN110196354B (en) | Method and device for detecting switching event of load | |
| CN111865126B (en) | Neutral point potential balance control method, device, medium and motor of three-level inverter | |
| JP3405786B2 (en) | Frequency change rate detection method for grid connection of distributed power supply | |
| CN112842075A (en) | Liquid level detection method and device, storage medium and kettle | |
| JP2005045914A (en) | Inverter device | |
| JPH0968556A (en) | Insulation diagnostic device for electrical equipment | |
| Zayyani et al. | Frequency estimation of unbalanced three-phase power system using a new LMS algorithm | |
| CN116908537B (en) | Current voltage frequency calculation circuit and method | |
| JPH11136849A (en) | Overcorrect detector |
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 |