CN104931806B - A kind of disturbing source theoretical based on empirical modal automates localization method - Google Patents
A kind of disturbing source theoretical based on empirical modal automates localization method Download PDFInfo
- Publication number
- CN104931806B CN104931806B CN201410466980.4A CN201410466980A CN104931806B CN 104931806 B CN104931806 B CN 104931806B CN 201410466980 A CN201410466980 A CN 201410466980A CN 104931806 B CN104931806 B CN 104931806B
- Authority
- CN
- China
- Prior art keywords
- disturbing source
- source
- empirical modal
- disturbing
- branch road
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000004807 localization Effects 0.000 title claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 3
- 238000000205 computational method Methods 0.000 claims 1
- 230000010355 oscillation Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 239000000284 extract Substances 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000003534 oscillatory effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- BYACHAOCSIPLCM-UHFFFAOYSA-N 2-[2-[bis(2-hydroxyethyl)amino]ethyl-(2-hydroxyethyl)amino]ethanol Chemical group OCCN(CCO)CCN(CCO)CCO BYACHAOCSIPLCM-UHFFFAOYSA-N 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Feedback Control In General (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
The present invention relates to a kind of disturbing source theoretical based on empirical modal to automate localization method, comprises the steps:(1) electrical power and angular frequency of each branch road in power network are gathered;(2) dominant component of disturbing source characteristic can be reacted by calculating;(3) empirical modal energy stream is calculated(4) rule of thumb mode energy streamMaximum and minimum value fit coenvelope line V respectively1And lower envelope line V (t)2(t);(5) empirical modal energy trend function V is calculated0(t);(6) calculation perturbation identifing source index Ab;(7) according to disturbing source distinguishing indexes Ab, automation positioning is carried out to disturbing source.Method provided by the invention can improve the accuracy of disturbance identifing source, and and can realizes automatic identification disturbing source, has stronger engineering practicability.
Description
Technical field
The present invention relates to field of power, and in particular to a kind of disturbing source theoretical based on empirical modal automates positioning
Method.
Background technology
With the development of extensive interconnected electric power system, Forced Power Oscillation in Power System problem increasingly attracts attention.
Forced oscillation have starting of oscillation is fast, can continue constant amplitude synchronized oscillation and oscillation-damped source after starting of oscillation after vibrate the spy to decay quickly
Point, if therefore the position of disturbing source can be found exactly, stablize for oscillation-damped, recovery system significant.
Some scholars have carried out the research of disturbance source locating method at present, and energy function method is shaken currently used for forcing
Swing one of main method of disturbance source locating.The electrical quantity gathered in real system includes the component that disturbing source determines, also includes
The component that non-disturbing source determines, when the component effect that non-disturbing source determines is larger, energy function method calculation error will be larger,
Sometimes it can not realize and be accurately positioned.
The present invention is directed to the deficiency of conventional energy function method, proposes a kind of disturbing source automation theoretical based on empirical modal
Localization method.This method application experience mode spatio-temporal filtering function extracts the dominant component of electrical quantity, utilizes empirical modal energy
Amount stream is extracted disturbing source feature, realizes disturbance according to the disturbing source distinguishing indexes obtained based on empirical modal energy trend function
The automatic identification in source.The present invention improves the accuracy of disturbance source locating, and provides reference for on-line automaticization of disturbing source.
The content of the invention
In view of the shortcomings of the prior art, the present invention provide it is a kind of can improve disturbance identifing source accuracy based on Empirical Mode
The theoretical disturbing source automation localization method of state.This method application experience mode spatio-temporal filtering function extracts the leading of electrical quantity
Component, extract using empirical modal energy stream disturbing source feature, according to being disturbed based on what empirical modal energy trend function obtained
Dynamic identifing source index realizes the automatic identification to disturbing source, has stronger engineering practicability.
In order to realize foregoing invention purpose, the present invention is realized using following technical proposals:
A kind of disturbing source theoretical based on empirical modal automates localization method, and it is theed improvement is that, this method includes
Following step:
(1) electrical power and angular frequency of each branch road in power network are gathered;
(2) dominant component of disturbing source characteristic can be reacted by calculating;
(3) empirical modal energy stream is calculated
(4) rule of thumb mode energy streamMaximum and minimum value fit coenvelope line V respectively1(t) and under
Envelope V2(t);
(5) empirical modal energy trend function V is calculated0(t);
(6) calculation perturbation identifing source index Ab;
(7) according to disturbing source distinguishing indexes Ab, automation positioning is carried out to disturbing source.
The dominant component can include branch road LijThe dominant component of electrical power variable quantityWith branch road LijAngular frequency becomes
The dominant component V of change amounti 1;
In the step (2), pass through empirical mode decomposition method (Empirical Mode Decomposition, abbreviation
EMD) can be in the hope of meeting the branch road L of following constraintsijThe dominant component of electrical power variable quantityWith branch road LijAngular frequency becomes
The dominant component V of change amounti 1:
①
In formula:PijFor branch road LijElectrical power;ωiFor the angular frequency of node i side;WithFor electrical power variable quantity
With the non-dominant component of angular frequency variable quantity, it is asked for by empirical mode decomposition method.
In the step (3), empirical modal energy streamIt can be calculated by following formula:
②
In formula, ω0=2 π f0, ω0On the basis of angular frequency, f0On the basis of frequency, t disturbance the time;
In the step (5), empirical modal energy trend function V0(t) can be calculated by following formula:
③
In the step (6), disturbing source distinguishing indexes AbIt can be tried to achieve by following methods:
A, repeated by least square method to empirical modal energy trend function V0(t) function of first order fitting is carried out, every time
Coefficient A, B are obtained after fitting, by all coefficient A composition sequences A;
B, by the maximum A in sequence AmaxAs a reference value, sequence A is normalized, obtains disturbing identifing source
Index Ab。
In the step (7), according to index AbThe method for disturbing source automate positioning is as follows:
A, A is worked asbWhen=1, represent to be disturbing source at this;
B, A is worked asb<=0, represent to be non-disturbing source at this;
C, when 0<Ab<When 1, AbNumerical value it is bigger, represent nearer from disturbing source at this.
Compared with prior art, the beneficial effects of the present invention are:
1st, the present invention has gathered the spatio-temporal filtering function of empirical modal and the concept of empirical modal energy stream, will can disturb
Dynamic source feature extraction comes out, and improves the accuracy of disturbance source locating;
2nd, the present invention judges disturbance using the disturbing source identification quantizating index obtained based on empirical modal energy trend function
The method of source position, that is, the accuracy of disturbance identifing source is improved, automatic identification disturbing source is realized again, there is stronger engineering
Practicality;
3rd, the vibration example analysis results of actual electric network demonstrate effectiveness of the invention, and it is more traditional to show the present invention
The advantage of energy function method, accuracy and practicality all increase.
Brief description of the drawings
Fig. 1 is the flow chart of the disturbing source automation localization method theoretical based on empirical modal of the present invention;
Fig. 2 is certain actual electric network vibration example provided by the invention;
Fig. 3 is empirical modal energy stream result of calculation provided by the invention.
Embodiment
The embodiment of the present invention is described in further detail below in conjunction with the accompanying drawings.
Referring to Fig. 1, the disturbing source automation localization method theoretical based on empirical modal of the invention comprises the steps:
(1) electrical power and angular frequency of each branch road in power network are gathered;
(2) the branch road L of disturbing source characteristic can be reacted by being calculated by empirical mode decomposition methodijElectrical power variable quantity is dominated
ComponentAnd the branch road L of disturbing source characteristic can be reactedijThe dominant component V of angular frequency variable quantityi 1;
Dominant component is calculated using equation below:
①
(3) it is based onAnd Vi 1Calculate empirical modal energy stream
Empirical modal energy stream is calculated using equation below
②
(4) rule of thumb mode energy streamMaximal point and minimal point fit coenvelope line V respectively1(t) and under
Envelope V2(t);
(5) empirical modal energy trend function V is calculated0(t);
Empirical modal energy trend function V is calculated using equation below0(t):
③
(6) it is based on empirical modal energy trend function V0(t) calculation perturbation identifing source index Ab;
Calculation perturbation identifing source index A with the following methodb:
A, to empirical modal energy trend function V0(t) least square fitting is carried out, you can obtain coefficient A, B, there is V0(t)
≈At+B;
B, sequence A maximum A is asked formax, and as a reference value, A sequences are normalized, disturbed
Dynamic identifing source index Ab。
(7) according to disturbing source distinguishing indexes Ab, realize the automatic identification of disturbing source.
According to index AbCharacteristic carry out disturbing source automatic identification, AbCharacteristic is as follows:
A, A is worked asbWhen=1, represent to be disturbing source at this;
B, A is worked asb<=0, represent to be non-disturbing source at this;
C, when 0<Ab<When 1, AbNumerical value it is bigger, represent nearer from disturbing source at this.
Embodiment
Below to the further detailed description of the present invention by taking the vibration example of certain actual electric network shown in Fig. 2 as an example, but this
Invention is not limited to given example.Under certain power network direct current sending decoupled mode, there is abnormal power and frequency in system
Rate oscillatory occurences, the PMU recordings of oscillatory process are as shown in Figure 2.Crash analysis and field demonstration, which work, to be shown, oscillation of power
Reason is the systems force oscillatory occurences that GD power plant governing system problem triggers.
Disturbance source locating analysis is carried out to vibration data at initial stage using method provided by the invention, step is as follows:
Step 1:The electrical power and angular frequency of each branch road in power network are gathered, as shown in Figure 2;
Step 2:Calculate the dominant component of branch road electrical power variable quantity and the branch road angular frequency of disturbing source characteristic can be reacted
The dominant component of variable quantity;
Step 3:Calculate empirical modal energy streamResult of calculation is as shown in Figure 3;
Step 4:Rule of thumb mode energy streamMaximal point and minimal point fit coenvelope line V respectively1(t)
And lower envelope line V2(t);
Step 5:Calculate empirical modal energy trend function V0(t);
Step 6:To empirical modal energy trend function V0(t) carry out least square fitting, can obtain every time coefficient A,
B, due to being to carry out function of first order fitting, therefore there is V0(t)≈At+B;Obtained coefficient A composition sequence A will be fitted every time, ask
Take sequence A maximum Amax, and as a reference value, A sequences are normalized, obtain disturbing source distinguishing indexes
Ab;
Step 7:According to disturbing source distinguishing indexes Ab, the automatic identification of disturbing source is realized, by being based on experience in this example
The analysis result that the disturbing source automation localization method of modal theory carries out automatic identification to disturbing source is as shown in table 1:
Table 1
It to be not difficult to find out from table 1, GD stations index is 1, is system disturbance source, and JP stations index is less than 0, is non-disturbing source, with
System actual conditions are consistent, so as to demonstrate validity of the method provided by the invention in real system analysis.
Finally it should be noted that:The above embodiments are merely illustrative of the technical scheme of the present invention and are not intended to be limiting thereof, to the greatest extent
The present invention is described in detail with reference to above-described embodiment for pipe, those of ordinary skills in the art should understand that:Still
The embodiment of the present invention can be modified or equivalent substitution, and without departing from any of spirit and scope of the invention
Modification or equivalent substitution, it all should cover among scope of the presently claimed invention.
Claims (1)
1. a kind of disturbing source theoretical based on empirical modal automates localization method, it is characterised in that methods described includes following
Step:
(1) electrical power and angular frequency of each branch road in power network are gathered;
(2) dominant component of disturbing source characteristic can be reacted by calculating;
(3) empirical modal energy stream is calculated
(4) rule of thumb mode energy streamMaximum and minimum value fit coenvelope line V respectively1And lower envelope line (t)
V2(t);
(5) empirical modal energy trend function V is calculated0(t);
(6) calculation perturbation identifing source index Ab;
(7) according to disturbing source distinguishing indexes Ab, automation positioning is carried out to disturbing source;
The dominant component includes branch road LijThe dominant component of electrical power variable quantityWith branch road LijAngular frequency variable quantity is dominated
Component
In the step (2), the branch road L for meeting following constraints is tried to achieve by empirical mode decomposition methodijElectrical power changes
The dominant component of amountWith branch road LijThe dominant component of angular frequency variable quantity
In formula:PijFor branch road LijElectrical power;ωiFor the angular frequency of node i side;WithFor electrical power variable quantity and angular frequency
The non-dominant component of rate variable quantity, it is asked for by empirical mode decomposition method;
In the step (3), empirical modal energy streamCalculation formula is as follows:
In formula, ω0=2 π f0, ω0On the basis of angular frequency, f0On the basis of frequency, t disturbance the time;
In the step (5), empirical modal energy trend function V0(t) calculation formula is as follows:
In the step (6), disturbing source distinguishing indexes AbComputational methods are as follows:
A, repeated by least square method to empirical modal energy trend function V0(t) function of first order fitting is carried out, every time after fitting
Coefficient A, B are obtained, by all coefficient A composition sequences A;
B, by the maximum A in sequence AmaxAs a reference value, sequence A is normalized, obtains disturbing source distinguishing indexes
Ab;
In the step (7), according to index AbThe method for disturbing source automate positioning is as follows:
A, A is worked asbWhen=1, represent to be disturbing source at this;
B, A is worked asb<=0, represent to be non-disturbing source at this;
C, when 0<Ab<When 1, AbNumerical value it is bigger, represent nearer from disturbing source at this.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410466980.4A CN104931806B (en) | 2014-09-13 | 2014-09-13 | A kind of disturbing source theoretical based on empirical modal automates localization method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410466980.4A CN104931806B (en) | 2014-09-13 | 2014-09-13 | A kind of disturbing source theoretical based on empirical modal automates localization method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104931806A CN104931806A (en) | 2015-09-23 |
CN104931806B true CN104931806B (en) | 2018-01-19 |
Family
ID=54119065
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410466980.4A Active CN104931806B (en) | 2014-09-13 | 2014-09-13 | A kind of disturbing source theoretical based on empirical modal automates localization method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104931806B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108195326B (en) * | 2017-12-14 | 2019-08-20 | 大连理工大学 | A kind of shock localization method based on empirical mode decomposition |
CN110890753B (en) * | 2019-12-03 | 2021-08-03 | 国网湖南省电力有限公司 | Generator set disturbance source positioning method based on VMD algorithm |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6738734B1 (en) * | 1996-08-12 | 2004-05-18 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Empirical mode decomposition apparatus, method and article of manufacture for analyzing biological signals and performing curve fitting |
CN101738551A (en) * | 2009-12-15 | 2010-06-16 | 西南交通大学 | Method for intelligent analysis of transient power quality disturbance based on networking |
CN103208808A (en) * | 2013-03-07 | 2013-07-17 | 武汉大学 | Power system sub-synchronous oscillation mode identification method |
US8660848B1 (en) * | 2010-08-20 | 2014-02-25 | Worcester Polytechnic Institute | Methods and systems for detection from and analysis of physical signals |
-
2014
- 2014-09-13 CN CN201410466980.4A patent/CN104931806B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6738734B1 (en) * | 1996-08-12 | 2004-05-18 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Empirical mode decomposition apparatus, method and article of manufacture for analyzing biological signals and performing curve fitting |
CN101738551A (en) * | 2009-12-15 | 2010-06-16 | 西南交通大学 | Method for intelligent analysis of transient power quality disturbance based on networking |
US8660848B1 (en) * | 2010-08-20 | 2014-02-25 | Worcester Polytechnic Institute | Methods and systems for detection from and analysis of physical signals |
CN103208808A (en) * | 2013-03-07 | 2013-07-17 | 武汉大学 | Power system sub-synchronous oscillation mode identification method |
Non-Patent Citations (2)
Title |
---|
基于EEMD自适应阈值去噪的电能质量;韩刚 等;《电测与仪表》;20140125;第51卷(第2期);第45-49,57页 * |
基于经验模态分解的电能质量扰动信号定位方法;黄奂 等;《电网技术》;20100531;第34卷(第5期);第41-45页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104931806A (en) | 2015-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103970801B (en) | Microblogging advertisement blog article recognition methods and device | |
CN106371801A (en) | Voice mouse system based on voice recognition technology | |
CN104931806B (en) | A kind of disturbing source theoretical based on empirical modal automates localization method | |
CN108200111A (en) | Resource allocation information update method, device and resource interface equipment | |
CN105528432B (en) | A kind of digital resource hot spot generation method and device | |
CN103984795A (en) | IEC (international electro-technical commission) 61850-based on-line partial discharge monitoring equipment communication simulation method | |
CN106126279A (en) | Automatically the method and system of interpolation BIM model race fileinfo | |
CN105488293B (en) | A kind of method and system configuring virtual terminator line in SCD | |
CN107134776A (en) | The centralized positioning of multi frequency system asynchronous oscillation and displacement path tracing system and method | |
CN104239643B (en) | WAMS applications algorithm checking system and its verification method based on fault recorder data | |
CN103955561A (en) | IEC61850-based on-line lightning arrester monitoring equipment communication simulation method | |
CN105653777A (en) | Intelligent mapping system of general survey of underground power line | |
CN110209834A (en) | A kind of hypergraph construction method for manufacturing industry process equipment Information Atlas | |
Abubakar et al. | Transient faults in wind energy conversion systems: Analysis, modelling methodologies and remedies | |
CN104331512B (en) | A kind of BBS pages automatic acquiring method | |
CN110955808B (en) | Method for matching and associating four remote monitoring information of main plant station | |
CN111723297B (en) | Dual-semantic similarity judging method for grid society situation research and judgment | |
CN207867479U (en) | A kind of man-machine interactive system based on NFC | |
CN106325852A (en) | Application dynamic interface customizing method | |
CN103455521A (en) | Alternating-current contactor type selection processing method and system | |
CN205590626U (en) | Long -range autonomous system of three -phase separator | |
CN105095461A (en) | Family karaoke ordering method and system | |
CN205334564U (en) | Safe inspection navigation management system | |
CN103838882A (en) | Method for dispatching, searching for and supplementing telemetry data of energy management system | |
CN103279899A (en) | Method and system for expanding electric network operation panoramic model |
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 |