CN110441601A - Multi-harmonic Sources responsibility quantization method based on independent component analysis - Google Patents
Multi-harmonic Sources responsibility quantization method based on independent component analysis Download PDFInfo
- Publication number
- CN110441601A CN110441601A CN201910956177.1A CN201910956177A CN110441601A CN 110441601 A CN110441601 A CN 110441601A CN 201910956177 A CN201910956177 A CN 201910956177A CN 110441601 A CN110441601 A CN 110441601A
- Authority
- CN
- China
- Prior art keywords
- harmonic
- suspicious
- source
- bus
- voltage
- 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
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/16—Spectrum analysis; Fourier analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2134—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on separation criteria, e.g. independent component analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
The invention discloses the Multi-harmonic Sources responsibility quantization methods based on independent component analysis, the following steps are included: S1: implementing synchro measure to the connected bus of suspicious harmonic source and concern bus and obtain voltage, current measurement data, carry out Fourier decomposition to measurement data and obtain each harmonic voltage, current data;S2: according to step S1 using certain subharmonic voltage, current data as observation signal, source signal and hybrid matrix are solved with independent component analysis method;S3: repeating step S2, calculates each suspicious harmonic source to the harmonic voltage transmission level of concern bus and quantifies responsibility.The present invention utilizes the certainty of hybrid matrix and suspicious harmonic wave source signal estimated value product in independent component analysis method, each suspicious harmonic source responsibility of direct solution, each suspicious harmonic source need not be calculated to the transfger impedance of concern bus, it reduces and calculates error, and no matter there are several suspicious harmonic sources, it is only solved with an independent component analysis method, reduce calculation amount and avoids error accumulation effect.
Description
Technical field
The present invention relates to Multi-harmonic Sources responsibility quantification technique fields, and in particular to the Multi-harmonic Sources based on independent component analysis
Responsibility quantization method.
Background technique
Harmonic contributions quantitative evaluation is the premise of harmonic wave control.To certain points of common connection (a Point of Common
Coupling, PCC) harmonic contributions quantification problem, many solutions have been proposed in a large amount of scholars.However, in recent years,
With the large-scale grid connection of high permeability renewable energy, multi-infeed HVDC drop point gradually increase and high-speed rail, city rail
The non-linear customers such as traffic, electric automobile charging pile are increasing, and harmonic source quantity, distribution and harmonic content all exist in power grid
It is continuously increased, harmonic pollution and harm are on the rise, and the power grid accident of the initiations such as harmonic resonance and amplification happens occasionally.It is existing
Based on the harmonic contributions quantification technique of single PCC point measurement, only consider that PCC point two sides (system side and user side) is respective all humorous
The overall responsibility of wave source cannot specify the percentage contribution for the different harmonic sources that two sides are respectively contained, this is multiple for determination suspicious
The analysis of a possibility that respective responsibilities and network harmonic current distribution, harmonic resonance of harmonic source is inadequate.When certain pays close attention to bus
Harmonic excess is when being caused by multiple suspicious harmonic source collective effects, then to need clearly each suspicious harmonic source in the concern bus
Upper respective harmonic contribution degree.Study as a result, multiple PCC measurement point Multi-harmonic Sources responsibility quantification problem be it is very necessary,
The problem is also becoming a new research hotspot in frequency analysis field.
Summary of the invention
The technical problems to be solved by the present invention are: when being applied to Practical Project scene, existing distributing Multi-harmonic Sources
Its calculated result of responsibility quantization method is affected by background harmonics, there is coupling between each suspicious harmonic source responsibility calculated
Estimated result bigger error is caused, the present invention provides the Multi-harmonic Sources responsibilities based on independent component analysis to solve the above problems
Quantization method.
The present invention is achieved through the following technical solutions:
Multi-harmonic Sources responsibility quantization method based on independent component analysis, method includes the following steps:
S1: in power grid to be assessed the connected bus of suspicious harmonic source and concern bus implement synchro measure and obtain voltage,
Current measurement data carries out Fourier decomposition to measurement data and obtains each harmonic voltage, current data;
S2: it according to step S1 using certain subharmonic voltage, current data as observation signal, is solved with independent component analysis method
Source signal and hybrid matrix;
S3: repeating step S2, calculates each suspicious harmonic source to the harmonic voltage transmission level of concern bus and quantifies responsibility.
The present invention is in application, assume that there are multiple suspicious harmonic sources in power grid, and only where the suspicious harmonic source in part
Harmonic wave measuring equipment is installed on bus and certain concern bus, each metric data is synchronous.The present invention is concerned with equipped amount
The suspicious harmonic source and background harmonics for surveying device bus are to the harmonic contributions partition problem for paying close attention to bus.
Further, step S1 includes following sub-step:
S11: it is assumed that the harmonic contributions of 3 suspicious harmonic sources in power grid need to be assessed, to the connected bus of 3 suspicious harmonic sources and pass
Note bus implements synchro measure and simultaneously obtains voltage, current measurement data, to measurement data carry out Fourier decomposition obtain each time it is humorous
Wave voltage, current data;
S12: choosing certain subharmonic voltage, current data, establishes data observation matrix, the data observation matrixXXSuch as following formula institute
Show:
(1)
In formula,RespectivelyhThe port current measured value of secondary suspicious harmonic source A, B, C,It is female for concern
LineXThe of place's measurementhSubharmonic voltage,NFor sampled data total length.
Further, step S2 includes following sub-step:
S21: by data observation matrixXXAs observation signal, the mathematical model solved based on independent component analysis method is constructed,
Model is as follows:
(2)
In formula,RespectivelyhThe port current measured value of secondary suspicious harmonic source A, B, C,It is female for concern
LineXThe of place's measurementhSubharmonic voltage;H suspicious harmonic source of respectively suspicious harmonic source A, B, C transmitting
Electric current,For equivalent background harmonics electric current;Respectively suspicious harmonic wave ource electric currentTo concern bus harmonic voltageTransfger impedance;XXFor observation signal,XXIt is expressed as,SSFor source signal,SSIt is expressed as,AAFor hybrid matrix,AAIt is expressed as;AACoefficient in matrixIndicate the mixing of source signal to observation signal
Relationship;
S22: to data observing matrixXXZero averaging and albefaction are carried out, the correlation between each observation signal, zero-mean are weakened
Change formula are as follows:
(3)
In formula,It indicatesXX iColumn vector;
Albefaction formula are as follows:
(4)
In formula,QFor whitening matrix,,ForCharacteristic value constitute diagonal matrix,
ForEigenvectors matrix;
S23: initial random weight vector is set, iteration updates according to the following formula, and to updatedBe orthogonalized and
Normalization;
(5)
In formula,,It is led for the single order of G;ForSingle order lead;, subscript H table
Showing conjugate transposition, subscript T indicates transposition,ForXXValue after going equalization and albefaction;
It is decoupled by harmonic current of the step S23 to each suspicious harmonic source transmitting, harmonic emission level is absolutely blamed
Appointing is not influenced by other suspicious harmonic source variations, guarantees the subsequent reality that can more accurately estimate each suspicious harmonic source in this way
Border responsibility;
S24: judgementWhether restrain, if not restraining, repeats step S23;
S25: repeating step S23 ~ S24, until acquiring separation matrix,kFor isolated source signal
Number, it is rightWIt inverts to obtain the unknown hybrid matrix of order and scale, while acquire it is separated go out do not restore ruler
Degree suspicious harmonic wave ource electric current be;
S26: according to correlation analysis to source signalAnd hybrid matrixOrder restored, with the value after recovery
Update hybrid matrixAnd source signal。
Since harmonic wave source impedance is typically much deeper than system side harmonic impedance, so port current measured valueWith the suspicious harmonic wave ource electric current for not restoring scale having been estimated thatBetween have stronger correlation,
WithpearsonRelated coefficient indicates to the degree of correlation between two signal of source signal and hybrid matrix, by following formula filter out withThe maximum source signal of degree of correlation;
(6)
In formula,i=1,2,3,covFor covariance,For standard deviation,It indicatesXX iRow,It indicates j
Row;
Wheni=1,2,3When, successively calculate, and hybrid matrix is exchanged simultaneously i、Column and'si、Row;Most
After obtain hybrid matrixAAAnd source signalSSEstimated value, be shown below:
(7)
In formula,、Middle each element is respectively to correspond toAA、SSThe estimated value of middle each element,It is expressed as,It is expressed as;
There is the uncertainty of scale between the source signal and hybrid matrix recovered based on above-mentioned independent component analysis method, with
Relationship between true value is a uncertain proportionality coefficient, which indicates are as follows:
(8)
In formula,XXIt is expressed as,AAIt is expressed as,SSTable
It is shown as;
(8) formula equivalent representation are as follows:
(9)
In formula,For proportionality coefficient (plural number),KFor proportionality coefficient matrix,。
Further, step S3 includes following sub-step:
S31: hybrid matrix is obtained by independent component analysis methodAAAnd source signalSSEstimated value after, concern bus harmonic wave electricity
Pressure is decomposed into suspicious harmonic wave ource electric currentTo concern bus effect and, be shown below:
(10)
In formula,Suspicious harmonic source A, B, C and the suspicious harmonic source of background are respectively indicated to pass
Infuse the vector contribution of bus h subharmonic voltage;
S32: calculating separately suspicious harmonic source A, B, C and background harmonics to the h subharmonic voltage scalar contribution degree of target bus,
In, h subharmonic voltage scalar contribution degree of the suspicious harmonic source A to target bus, suspicious harmonic source B is to the h of target bus
Subharmonic voltage scalar contribution degree, suspicious harmonic source C is to the h subharmonic voltage scalar contribution degree of target bus, background
H subharmonic voltage scalar contribution degree of the harmonic wave to target busIt respectively indicates as follows;
(11)
(12)
(13)
(14)
(11) ~ (14) in formula, t indicates the time, and N indicates signal length,It indicatesWithBetween angle,It indicatesWithBetween angle,It indicatesWithBetween angle,It indicatesWithBetween angle;
S33: (i.e. according to the harmonic contributions of suspicious harmonic source A, B, C calculated in step S32 and background harmonics、、、), judge suspicious harmonic source A, B, C and background harmonics to busX'shSubharmonic voltage role, wherein
IfGreater than zero, indicate suspicious harmonic source A to busX'shSubharmonic voltage, which serves, helps increasing;IfLess than zero, indicate
Suspicious harmonic source A is to busX'shSubharmonic voltage plays inhibiting effect;
IfGreater than zero, indicate suspicious harmonic source B to busX'shSubharmonic voltage, which serves, helps increasing;IfLess than zero,
Indicate suspicious harmonic source B to busX'shSubharmonic voltage plays inhibiting effect;
IfGreater than zero, indicate suspicious harmonic source C to busX'shSubharmonic voltage, which serves, helps increasing;IfLess than zero,
Indicate suspicious harmonic source C to busX'shSubharmonic voltage plays inhibiting effect;
IfGreater than zero, indicate background harmonics to busX'shSubharmonic voltage, which serves, helps increasing;IfLess than zero, table
Show background harmonics to busX'shSubharmonic voltage plays inhibiting effect.
The present invention has the advantage that and the utility model has the advantages that
1, the present invention is based on Norton equivalent harmonic source electric current quantization Multi-harmonic Sources responsibility compared to conventional method the method for the present invention, to each
The harmonic current of suspicious harmonic source transmitting is decoupled, and harmonic emission level absolute responsibility is not become by other suspicious harmonic sources
The influence of change can more accurately estimate the practical responsibility of each harmonic source;
2, the present invention utilizes the certainty of hybrid matrix and suspicious harmonic wave source signal estimated value product in independent component analysis method,
Each suspicious harmonic source responsibility of direct solution, it is not necessary to calculate each suspicious harmonic source to the transfger impedance of concern bus, reduce to calculate and miss
Difference, and no matter have several suspicious harmonic sources, it is only solved with an independent component analysis method, reduces calculation amount and avoid mistake
Poor cumulative effect;
3, the present invention can obtain effective assessment result when background harmonics change greatly, and method error level is lower, application
Scene is extensive.
Detailed description of the invention
Attached drawing described herein is used to provide to further understand the embodiment of the present invention, constitutes one of the application
Point, do not constitute the restriction to the embodiment of the present invention.In the accompanying drawings:
Fig. 1 is that the present invention is based on the calculating flow chart of steps of the Multi-harmonic Sources responsibility quantization method of independent component analysis.
Fig. 2 is that the present invention is based on the analogue system figures of the Multi-harmonic Sources responsibility quantization method of independent component analysis.
Fig. 3 is that the present invention is based on the calculated result comparison diagrams of the Multi-harmonic Sources responsibility quantization method of independent component analysis.
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention clearer, below with reference to embodiment and attached drawing, to this
Invention is described in further detail, and exemplary embodiment of the invention and its explanation for explaining only the invention, are not made
For limitation of the invention.
Embodiment
As shown in Figure 1 to Figure 3, the Multi-harmonic Sources responsibility quantization method based on independent component analysis, this method includes following
Step:
S1: in power grid to be assessed the connected bus of suspicious harmonic source and concern bus implement synchro measure and obtain voltage,
Current measurement data carries out Fourier decomposition to measurement data and obtains each harmonic voltage, current data;
S2: it according to step S1 using certain subharmonic voltage, current data as observation signal, is solved with independent component analysis method
Source signal and hybrid matrix;
S3: repeating step S2, calculates each suspicious harmonic source to the harmonic voltage transmission level of concern bus and quantifies responsibility.
The present embodiment implement when, it is assumed that the suspicious harmonic source of different measurement points be it is independent or weak relevant, be based on this, mentioned
Method decomposites the harmonic current and background harmonics equivalent current of each suspicious harmonic source transmitting using independent component analysis method, together
In Shi Liyong independent component analysis method between hybrid matrix and suspicious harmonic source current estimation value product certainty, can be direct
Solve each suspicious harmonic source responsibility, it is not necessary to solve suspicious harmonic source and reduce to the transfger impedance of concern bus and calculate error, and nothing
It by there is several suspicious harmonic sources, is only solved with an independent component analysis method, reduces calculation amount and avoid error accumulation
Effect can obtain effective assessment result when background harmonics change greatly.The present invention passes through setting above-mentioned steps, it is assumed that
It is in weak correlation between each suspicious harmonic source current signal, can effectively solves the harmonic contributions of each suspicious harmonic source.
In order to further illustrate the effect that the present embodiment generates, tested using IEEE14 node system as shown in Figure 2
Card, which includes 2 generators, 3 synchronous capacitors, (1-14 label indicates 14 mothers to 14 buses in Fig. 2
Line), 15 transmission lines of electricity, 3 transformers.In calculating process, it is subtranient reactance by generator and synchronous capacitor equivalence, becomes
Depressor equivalence is impedance, and transmission line of electricity uses π type equivalent circuit, simulates transport properties of the harmonic wave in power grid.
Suspicious harmonic source is respectively arranged at bus 3,4,12,14, and harmonic wave measuring equipment is housed on median generatrix 4,12,14.
Bus 3 is not charged with harmonic wave measuring equipment, and the harmonic wave generated with the connected harmonic source of bus 3 is set as to the background harmonics of system.Assuming that
Bus 4,12,14 corresponds to suspicious harmonic source A, B, C, seeks suspicious harmonic source A, B, C respectively and blames to the harmonic wave of concern bus 11
Appoint.
5 subharmonic currents that suspicious harmonic wave source node 3,4,12,14 injects be set to mean value be 0.002,0.004,
0.008,0.005, variance is 0.003 and obeys the stochastic variable (per unit value) of laplacian distribution.Meanwhile considering measurement error
And network parameter variation influence, all voltage and current measurement values be added signal-to-noise ratio (signal to noise ratio,
SNR) 10% random fluctuation is added in the random noise for being 20dB, all load impedances.Use the method for the present invention above-mentioned formula
(11) ~ (14) solve each suspicious harmonic source harmonic contributions, and estimated result accuracy and other methods are compared, other
Method includes that true value method, original ICA method, covariance method, IGG method and data segment choose method, as shown in Figure 3.
Shown in table 1 under the variation of different background harmonic fluctuations, using the calculated each suspicious harmonic source A of the method for the present invention,
B, the harmonic contributions of C and background harmonics are indicated with " estimated value " in the table;" true value " in table is used as reference pair ratio
Reference value.
The calculated result contrast table of the present invention of table 1
Can be explicit from Fig. 3 and table 1 see, the present embodiment method is compared with the prior art, calculated result " estimated value "
Closer to " true value ", calculated result is more accurate, and is influenced by background harmonics small, and method is fast and effective, error level compared with
Low, application scenarios are extensive.And the harmonic contributions of calculated each suspicious harmonic source A, B, C and background harmonics are equal in the present embodiment
Greater than zero, illustrate suspicious harmonic source A, B, C and background harmonics to busX'shSubharmonic voltage plays the role of helping increasing.
Above-described specific embodiment has carried out further the purpose of the present invention, technical scheme and beneficial effects
It is described in detail, it should be understood that being not intended to limit the present invention the foregoing is merely a specific embodiment of the invention
Protection scope, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should all include
Within protection scope of the present invention.
Claims (5)
1. the Multi-harmonic Sources responsibility quantization method based on independent component analysis, it is characterised in that: method includes the following steps:
S1: in power grid to be assessed the connected bus of suspicious harmonic source and concern bus implement synchro measure and obtain voltage,
Current measurement data carries out Fourier decomposition to measurement data and obtains each harmonic voltage, current data;
S2: it according to step S1 using certain subharmonic voltage, current data as observation signal, is solved with independent component analysis method
Source signal and hybrid matrix;
S3: repeating step S2, calculates each suspicious harmonic source to the harmonic voltage transmission level of concern bus and quantifies responsibility.
2. the Multi-harmonic Sources responsibility quantization method according to claim 1 based on independent component analysis, it is characterised in that: step
Rapid S1 includes following sub-step:
S11: it is assumed that the harmonic contributions of 3 suspicious harmonic sources in power grid need to be assessed, to the connected bus of 3 suspicious harmonic sources and pass
Note bus implements synchro measure and simultaneously obtains voltage, current measurement data, to measurement data carry out Fourier decomposition obtain each time it is humorous
Wave voltage, current data;
S12: choosing certain subharmonic voltage, current data, establishes data observation matrix, the data observation matrixXXSuch as following formula institute
Show:
In formula,RespectivelyhThe port current measured value of secondary suspicious harmonic source A, B, C,To pay close attention to busX
The of place's measurementhSubharmonic voltage,NFor sampled data total length.
3. the Multi-harmonic Sources responsibility quantization method according to claim 1 or 2 based on independent component analysis, feature exist
In: step S2 includes following sub-step:
S21: certain subharmonic voltage, current data are established by data observation matrix according to step S1XXAs observation signal, building
Based on the mathematical model that independent component analysis method solves, model is as follows:
In formula,RespectivelyhThe port current measured value of secondary suspicious harmonic source A, B, C,To pay close attention to busXThe of place's measurementhSubharmonic voltage;H suspicious harmonic source electricity of respectively suspicious harmonic source A, B, C transmitting
Stream,For equivalent background harmonics electric current;Respectively suspicious harmonic wave ource electric currentTo concern bus harmonic voltageTransfger impedance;XXFor observation signal,XXIt is expressed as,SSFor source signal,SSIt is expressed as,AAFor hybrid matrix,AAIt is expressed as;AACoefficient in matrixMixed relationship of the expression source signal to observation signal;
S22: to data observing matrixXXZero averaging and albefaction are carried out, the correlation between each observation signal, zero-mean are weakened
Change formula are as follows:
In formula,It indicatesXX iColumn vector;
Albefaction formula are as follows:
In formula,QFor whitening matrix,,ForCharacteristic value constitute diagonal matrix,ForEigenvectors matrix;
S23: initial random weight vector is set, iteration updates according to the following formula, and to updatedBe orthogonalized and
Normalization;
In formula,,It is led for the single order of G;ForSingle order lead;, subscript H expression is altogether
Yoke transposition, subscript T indicate transposition,ForXXValue after going equalization and albefaction;
S24: judgementWhether restrain, if not restraining, repeats step S23;
S25: repeating step S23 ~ S24, until acquiring separation matrix,kFor isolated source signal
Number, it is rightWIt inverts to obtain the unknown hybrid matrix of order and scale, while acquire it is separated go out do not restore scale
Suspicious harmonic wave ource electric current be;
S26: according to correlation analysis to source signalAnd hybrid matrixOrder restored, more with the value after recovery
New hybrid matrixAnd source signal。
4. the Multi-harmonic Sources responsibility quantization method according to claim 3 based on independent component analysis, it is characterised in that: step
Rapid S26 includes following sub-step:
WithpearsonRelated coefficient indicates to filter out the degree of correlation between two signal of source signal and hybrid matrix by following formula
WithThe maximum source signal of degree of correlation;
In formula,i=1,2,3,covFor covariance,For standard deviation,It indicatesXX iRow,It indicates jRow;
Wheni=1,2,3When, successively calculate, and hybrid matrix is exchanged simultaneously i、Column and'si、Row;Finally
Obtain hybrid matrixAAAnd source signalSSEstimated value, be shown below:
In formula,、Middle each element is respectively to correspond toAA、SSThe estimated value of middle each element,It is expressed as,It is expressed as;
There is the uncertainty of scale between the source signal and hybrid matrix recovered based on above-mentioned independent component analysis method, with
Relationship between true value is a uncertain proportionality coefficient, which indicates are as follows:
In formula,XXIt is expressed as,AAIt is expressed as,SSIt indicates
For;
Equivalent representation are as follows:
In formula,For proportionality coefficient,KFor proportionality coefficient matrix,。
5. the Multi-harmonic Sources responsibility quantization method according to claim 4 based on independent component analysis, it is characterised in that: step
Rapid S3 includes following sub-step:
S31: hybrid matrix is obtained by independent component analysis methodAAAnd source signalSSEstimated value after, concern bus harmonic wave electricity
Pressure is decomposed into suspicious harmonic wave ource electric currentTo concern bus effect and, be shown below:
In formula,Suspicious harmonic source A, B, C and the suspicious harmonic source of background are respectively indicated to concern bus
The vector of h subharmonic voltage is contributed;
S32: calculating separately suspicious harmonic source A, B, C and background harmonics to the h subharmonic voltage scalar contribution degree of target bus,
In, h subharmonic voltage scalar contribution degree of the suspicious harmonic source A to target busIt is expressed as follows, is acquired using same procedure
The harmonic contributions of remaining suspicious harmonic source B, C and background harmonics;
In formula, t indicates the time, and N indicates signal length,It indicatesWithBetween angle;
S33: judge suspicious harmonic source A, B, C and background harmonics to busX'shSubharmonic voltage role, wherein ifGreater than zero, indicate suspicious harmonic source A to busX'shSubharmonic voltage, which serves, helps increasing;IfLess than zero, indicate suspicious
Harmonic source A is to busX'shSubharmonic voltage plays inhibiting effect;Remaining suspicious harmonic source B, C and background are judged using same procedure
Harmonic wave is to busX'shSubharmonic voltage role.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910956177.1A CN110441601B (en) | 2019-10-10 | 2019-10-10 | Independent component analysis-based multi-harmonic source responsibility quantification method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910956177.1A CN110441601B (en) | 2019-10-10 | 2019-10-10 | Independent component analysis-based multi-harmonic source responsibility quantification method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110441601A true CN110441601A (en) | 2019-11-12 |
CN110441601B CN110441601B (en) | 2020-01-07 |
Family
ID=68440635
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910956177.1A Active CN110441601B (en) | 2019-10-10 | 2019-10-10 | Independent component analysis-based multi-harmonic source responsibility quantification method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110441601B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112834821A (en) * | 2020-12-31 | 2021-05-25 | 国网上海市电力公司 | Multi-harmonic source harmonic responsibility quantification method for multi-direct-current feed-in full cable power supply area |
CN112858782A (en) * | 2021-01-07 | 2021-05-28 | 国网河南省电力公司电力科学研究院 | Harmonic responsibility quantification method and system under influence of multi-user interaction of power system |
CN114660362A (en) * | 2022-05-19 | 2022-06-24 | 四川大学 | System side harmonic impedance estimation method and device based on mutual information data optimization |
CN114709821A (en) * | 2022-04-20 | 2022-07-05 | 山东大学 | Power distribution network partition harmonic responsibility evaluation method and system based on limited measurement |
CN117269604A (en) * | 2023-09-18 | 2023-12-22 | 四川大学 | Multi-harmonic source responsibility quantification method and system considering impedance change |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003009430A (en) * | 2001-06-19 | 2003-01-10 | Central Res Inst Of Electric Power Ind | Method and apparatus for remotely monitoring electrical device and method and apparatus for estimating power consumption utilizing the same |
CN102323494A (en) * | 2011-05-24 | 2012-01-18 | 山东电力研究院 | Method for distinguishing multiple harmonic sources |
CN104865467A (en) * | 2015-05-11 | 2015-08-26 | 国家电网公司 | Harmonic wave source positioning and harmonic wave responsibility dividing method applied to active power distribution network |
CN108535547A (en) * | 2018-04-19 | 2018-09-14 | 西南交通大学 | A kind of adaptive independent component analysis method of harmonic impedance estimation |
CN108802491A (en) * | 2018-07-06 | 2018-11-13 | 四川大学 | Harmonic emission level estimation based on harmonic source independence solves new method |
CN109541305A (en) * | 2018-11-29 | 2019-03-29 | 广西电网有限责任公司电力科学研究院 | A kind of harmonic contributions partitioning model and harmonic contributions calculation method |
CN110309545A (en) * | 2019-06-04 | 2019-10-08 | 上海电力学院 | Based on the harmonic compensation method for improving particle swarm algorithm acquisition system harmonic impedance |
-
2019
- 2019-10-10 CN CN201910956177.1A patent/CN110441601B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003009430A (en) * | 2001-06-19 | 2003-01-10 | Central Res Inst Of Electric Power Ind | Method and apparatus for remotely monitoring electrical device and method and apparatus for estimating power consumption utilizing the same |
CN102323494A (en) * | 2011-05-24 | 2012-01-18 | 山东电力研究院 | Method for distinguishing multiple harmonic sources |
CN104865467A (en) * | 2015-05-11 | 2015-08-26 | 国家电网公司 | Harmonic wave source positioning and harmonic wave responsibility dividing method applied to active power distribution network |
CN108535547A (en) * | 2018-04-19 | 2018-09-14 | 西南交通大学 | A kind of adaptive independent component analysis method of harmonic impedance estimation |
CN108802491A (en) * | 2018-07-06 | 2018-11-13 | 四川大学 | Harmonic emission level estimation based on harmonic source independence solves new method |
CN109541305A (en) * | 2018-11-29 | 2019-03-29 | 广西电网有限责任公司电力科学研究院 | A kind of harmonic contributions partitioning model and harmonic contributions calculation method |
CN110309545A (en) * | 2019-06-04 | 2019-10-08 | 上海电力学院 | Based on the harmonic compensation method for improving particle swarm algorithm acquisition system harmonic impedance |
Non-Patent Citations (1)
Title |
---|
杨雪萍 等: "基于核独立分量分析的谐波责任划分", 《电力系统及其自动化学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112834821A (en) * | 2020-12-31 | 2021-05-25 | 国网上海市电力公司 | Multi-harmonic source harmonic responsibility quantification method for multi-direct-current feed-in full cable power supply area |
CN112834821B (en) * | 2020-12-31 | 2023-08-11 | 国网上海市电力公司 | Multi-DC feed-in full-cable power supply area multi-harmonic source harmonic responsibility quantification method |
CN112858782A (en) * | 2021-01-07 | 2021-05-28 | 国网河南省电力公司电力科学研究院 | Harmonic responsibility quantification method and system under influence of multi-user interaction of power system |
CN112858782B (en) * | 2021-01-07 | 2022-04-26 | 国网河南省电力公司电力科学研究院 | Harmonic responsibility quantification method and system under influence of multi-user interaction of power system |
CN114709821A (en) * | 2022-04-20 | 2022-07-05 | 山东大学 | Power distribution network partition harmonic responsibility evaluation method and system based on limited measurement |
CN114709821B (en) * | 2022-04-20 | 2022-10-18 | 山东大学 | Power distribution network partition harmonic responsibility evaluation method and system based on limited measurement |
CN114660362A (en) * | 2022-05-19 | 2022-06-24 | 四川大学 | System side harmonic impedance estimation method and device based on mutual information data optimization |
CN117269604A (en) * | 2023-09-18 | 2023-12-22 | 四川大学 | Multi-harmonic source responsibility quantification method and system considering impedance change |
CN117269604B (en) * | 2023-09-18 | 2024-04-30 | 四川大学 | Multi-harmonic source responsibility quantification method and system considering impedance change |
Also Published As
Publication number | Publication date |
---|---|
CN110441601B (en) | 2020-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110441601A (en) | Multi-harmonic Sources responsibility quantization method based on independent component analysis | |
Apostolopoulos et al. | Accurate fault location algorithm for double-circuit series compensated lines using a limited number of two-end synchronized measurements | |
CN103412242B (en) | A kind of harmonic source location method based on Fast Independent Component Analysis and mutual information | |
Špelko et al. | CIGRE/CIRED JWG C4. 42: Overview of common methods for assessment of harmonic contribution from customer installation | |
JP6331465B2 (en) | Transformer connection phase determination device, transformer connection phase determination method, and transformer connection phase determination program | |
CN105911334A (en) | Harmonic current estimation method under condition of unknown harmonic impedance | |
CN108802491B (en) | Harmonic emission level estimation based on harmonic source independence solves new method | |
CN104793082A (en) | Harmonic correlation analysis based electricity system harmonic source recognition device | |
CN111141995A (en) | Line double-end steady-state distance measuring method and system based on amplitude comparison principle | |
CN108614194A (en) | A kind of double line down distance measuring method and device based on distributed parameter model | |
Vallés et al. | A new distributed measurement index for the identification of harmonic distortion and/or unbalance sources based on the IEEE Std. 1459 framework | |
CN109541516A (en) | A kind of voltage transformer wideband error measurement method | |
Han et al. | A new assessment for the total harmonic contributions at the point of common coupling | |
CN104122452A (en) | Large grounded screen power-frequency grounding impedance measurement method based on multiple-synchronous different-frequency method | |
Cataliotti et al. | A DAQ-based sampling wattmeter for IEEE Std. 1459-2010 powers measurements. Uncertainty evaluation in nonsinusoidal conditions | |
CN109858061A (en) | Distribution network equivalence and simplifying method towards voltage power sensitivity estimation | |
CN102004191A (en) | Method for calculating line capacitance in arc suppression coil parallel running system | |
US11327105B2 (en) | Fault location in multi-terminal tapped lines | |
Cataliotti et al. | IEEE Std. 1459 power quantities ratio approaches for simplified harmonic emissions assessment | |
Zhang et al. | Research on power system harmonic state estimation | |
US20170023627A1 (en) | Method and Apparatus for Precision Phasor Measurements Through a Medium-voltage Distribution Transformer | |
CN108152782B (en) | Method for testing correction coefficient of high-supply high-count electric energy meter | |
CN109001589A (en) | Line parameter circuit value time-domain calculation method based on open-phase operation fault recorder data | |
CN109116136A (en) | The electrification recognition methods of low-voltage cable route | |
Shcherbakova et al. | Revisiting the active power direction method |
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 |