CN110456159A - System side harmonic impedance estimation method and system based on amendment independent random vector - Google Patents

System side harmonic impedance estimation method and system based on amendment independent random vector Download PDF

Info

Publication number
CN110456159A
CN110456159A CN201910774203.9A CN201910774203A CN110456159A CN 110456159 A CN110456159 A CN 110456159A CN 201910774203 A CN201910774203 A CN 201910774203A CN 110456159 A CN110456159 A CN 110456159A
Authority
CN
China
Prior art keywords
harmonic
impedance
system side
time period
correlation
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
Application number
CN201910774203.9A
Other languages
Chinese (zh)
Other versions
CN110456159B (en
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.)
Sichuan University
Original Assignee
Sichuan University
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 Sichuan University filed Critical Sichuan University
Priority to CN201910774203.9A priority Critical patent/CN110456159B/en
Publication of CN110456159A publication Critical patent/CN110456159A/en
Application granted granted Critical
Publication of CN110456159B publication Critical patent/CN110456159B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R27/00Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
    • G01R27/02Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
    • G01R27/08Measuring resistance by measuring both voltage and current

Abstract

The invention discloses system side harmonic impedance estimation methods and system based on amendment independent random vector, comprising: the correlation analysis based on PCC point harmonic data determines weak relevant time period by the threshold value screening to normalizated correlation coefficient.In weak relevant time period, PCC point harmonic current has very weak correlation with system side harmonic voltage source, and then can derive the corresponding system harmonic impedance calculation formula of any time interval in weak relevant time period.Using the mean value of the corresponding impedance value of time intervals all in weak relevant time period as the system side harmonic impedance value of PCC point.The present invention is according to the correlation properties of PCC point harmonic data, the less quick computing system side harmonic impedance of sample point can be used, and it is suitable for the situation that background harmonics are unstable or two sides impedance ratio is close, reduce error caused by due to electric system changes the influence harmonic data correlation.

Description

System side harmonic impedance estimation method and system based on amendment independent random vector
Technical field
The present invention relates to the division fields of user side and system side harmonic contributions, and in particular, to based on amendment it is independent with The system side harmonic impedance estimation method and system of machine vector.
Background technique
As a large amount of power electronic equipments access power grid, the trend of China's electric system power electronics gradually shows.It is public Common-battery net nonlinear-load type increases, is distributed augmentation, and harmonic content increases, the frequency increases, and harmonic pollution problems are more serious, The accident that harmonic resonance and Harmonics amplification cause also happens occasionally, and mains by harmonics has jeopardized the safe and stable operation of power grid.Harmonic wave Administering becomes public electric wire net problem and task in the urgent need to address.Power grid points of common connection (point of common Coupling, PCC) system side and user side harmonic contributions it is effective distinguish be harmonic wave control premise and difficult point, this needs The accurate estimating system side harmonic impedance under the power grid environment of various complexity.
But existing harmonic impedance estimation method usually imply it is assumed hereinafter that: 1) system side harmonic fluctuations are put down substantially Surely;2) the user-side harmonic impedance of PCC point is much larger than system side.But in real system, various new energy equipment, high-speed rail, charging pile The access of equal non-linear customers is so that there are a large amount of harmonic source in power grid, background harmonics increase;User side installs filter additional, So that the relationship that user side and system side harmonic impedance are no longer much larger than.When assuming that being unsatisfactory for, existing method is calculated Harmonic impedance error is larger.
Summary of the invention
To solve the influence of background harmonics fluctuation and impedance ratio variation to conventional method, the present invention falls wind power plant, direct current A large amount of engineering measurement harmonic datas such as point and traditional nonlinear-load such as electric arc furnaces make correlation analysis, discovery PCC point harmonic wave electricity The auto-correlation and cross correlation of pressure and harmonic current are rapid decays, and apparent short range dependence is presented, this characteristic is humorous The calculating of wave impedance provides new approach.
For engineering, background harmonics fluctuation and two sides impedance ratio change the influence to existing method to the present invention in practice, propose A kind of method based on the weak dependence of points of common connection harmonic voltage and harmonic current estimation harmonic impedance.The method of the present invention Effectively inhibit the interference of background harmonics, and this method is equally applicable when two sides impedance ratio is close, greatly reduce because Electric system influence of the variation to harmonic data correlation and caused by error.It is said by emulation and the analysis of measured data verifying The validity of clear mentioned inventive method.
The present invention is achieved through the following technical solutions:
System side harmonic impedance estimation method based on amendment independent random vector, comprising the following steps:
A, the correlation analysis based on PCC point harmonic data, it is weak by the threshold value screening determination to normalizated correlation coefficient Relevant time period;
B, the equivalent circuit under certain subharmonic is established, circuit equation is established by Kirchhoff's second law, by equation Correlation function is asked to it after carrying out equivalent transformation;PCC point harmonic current and system side harmonic wave in the weak relevant time period determined in A Voltage source has very weak correlation, i.e. the correlation function of gained equation sees 0 as;
C, weak correlation can be obtained based on the harmonic impedance calculation formula derived based on correlation function equations obtained in B The corresponding impedance value of any time interval in period makees the mean value of the corresponding impedance value of time intervals all in weak relevant time period For PCC dot system side harmonic impedance value.
Further, in the system side harmonic impedance estimation method based on amendment independent random vector, the step A equation The foundation of group comprises the steps of:
A1, it is based on PCC point harmonic data, obtains the expression formula of the auto-correlation function of harmonic current and harmonic voltage are as follows:
And the expression formula of the cross-correlation function of harmonic voltage and harmonic current are as follows:
In formula,PCC point harmonic current and harmonic voltage are respectively indicated, N indicates that sample point sum, n represent letter Any moment of number sequence, k indicate the time interval of two signals, and E { } indicates mathematic expectaion, and conjugation is sought in * expression.
A2, it is screened by the threshold value to normalizated correlation coefficient, determines weak relevant time period.Harmonic current and harmonic voltage The expression formula of normalized autocorrelation coefficient are as follows:
The expression formula of the normalized-cross-correlation function of harmonic current and harmonic voltage are as follows:
In probability theory and statistics, related coefficient shows the intensity of linear relationship and direction between two stochastic variables. The attenuation degree for being spaced k at any time of correlation function is usually characterized with normalizated correlation coefficient.Normalizated correlation coefficient is just Negative sign only indicates relevant direction, the relevant degree of absolute value representation.It may be defined as when normalizated correlation coefficient is less than 0.3 weak It is related.
The present invention sets the threshold value of metric data related coefficient as 0.3, remembers the auto-correlation coefficient of harmonic current at PCC for the first time It is k1 less than time interval corresponding to threshold value, remembers that the cross-correlation coefficient of harmonic voltage and electric current is less than threshold value pair for the first time at PCC The time interval answered is k2, and maximum value therein is denoted as kmax, for N number of sample, define k value model corresponding to weak relevant time period It encloses for kmax~N-1.
Further, in the system side harmonic impedance estimation method based on amendment independent random vector, in the step B Formula equivalent transformation the following steps are included:
B1, equivalent circuit figure, column write circuit equation:
In formula, ZuIndicate system side harmonic impedance,Indicate the vector value of system side harmonic current source,Indicate system side The vector value of harmonic voltage source;
Formula (6) both ends are multiplied togetherCorrelation function is sought it afterwards, as a result as follows:
When metric data is that second-order statistic is ergodic, formula (7) is approximately:
In B2, the weak relevant time period determined in A, the cross correlation and harmonic current of harmonic voltage and harmonic current Autocorrelation is all very weak, it can thus be assumed that PCC point harmonic current has very weak correlation, formula with system side harmonic voltage source (8) 0 is seen as.
In system side harmonic impedance estimation method based on amendment independent random vector, system harmonics hinder in the step C The derivation of anti-calculation formula the following steps are included:
C1, formula (8) see 0 as, derive the table of system harmonic impedance corresponding to any time interval in weak relevant time period It is up to formula
In formula, the value range of k is kmax~N-1.
C2, to obtained (N-kmax) a impedance averages, obtain the harmonic impedance based on correction independent random vector Are as follows:
On the other hand, corresponding with method, the present invention also provides a kind of system side based on amendment independent random vector is humorous Wave impedance estimation system, the system comprises:
Weak relevant time period obtaining unit is analyzed for the correlation to PCC point harmonic data, obtains PCC point harmonic wave The normalized autocorrelation coefficient of harmonic current at the normalized-cross-correlation function and PCC of electric current and harmonic voltage, based on described Cross-correlation coefficient and the auto-correlation coefficient, given threshold screening obtain weak relevant time period;
Correlation function obtaining unit, the equivalent circuit for establishing under certain subharmonic are fixed based on Kirchoff s voltage Rule establishes circuit equation, and the circuit equation of foundation is carried out equivalent transformation, seeks the correlation function of circuit equation after equivalent transformation;
Harmonic impedance value computing unit, the correlation function for being obtained based on correlation function obtaining unit are derived described Power grid points of common connection system side harmonic impedance calculation formula corresponding to any time interval in weak relevant time period, based on deriving Harmonic impedance calculation formula, using the mean value of the corresponding impedance value of time intervals all in weak relevant time period as PCC dot system Side harmonic impedance value.
Compared with prior art, the present invention having the following advantages and benefits:
Existing method is usually intended to based on the assumption that 1) system side harmonic fluctuations are substantially steady;2) PCC point user side is humorous Wave impedance is much larger than system side.The present invention determines weak correlation on the basis of correlation analysis, through the screening to time interval k Period estimates harmonic impedance, only depends on the attribute of measured data itself, and reaction is system self-characteristic, weakens because of back The fluctuation of scape harmonic wave and influence of the variation to source data correlation of two sides impedance ratio and caused by error, the scope of application is wider.
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;
Fig. 1 is the equivalent circuit figure in the present invention;
Fig. 2 is the system side harmonic impedance estimating system composition schematic diagram based on amendment independent random vector in the present invention.
Specific embodiment
To better understand the objects, features and advantages of the present invention, with reference to the accompanying drawing and specific real Applying mode, the present invention is further described in detail.It should be noted that in the case where not conflicting mutually, the application's Feature in embodiment and embodiment can be combined with each other.
In the following description, numerous specific details are set forth in order to facilitate a full understanding of the present invention, still, the present invention may be used also Implemented with being different from the other modes being described herein in range using other, therefore, protection scope of the present invention is not by under The limitation of specific embodiment disclosed in face.
Embodiment one
System side harmonic impedance estimation method based on amendment independent random vector, comprising the following steps:
A, the correlation analysis based on PCC point harmonic data, it is weak by the threshold value screening determination to normalizated correlation coefficient Relevant time period;
B, the equivalent circuit under certain subharmonic is established, circuit equation is established by Kirchhoff's second law, by equation Correlation function is asked to it after carrying out equivalent transformation;PCC point harmonic current and system side harmonic wave in the weak relevant time period determined in A Voltage source has very weak correlation, i.e. the correlation function of gained equation sees 0 as;
C, based on correlation function equations obtained in B, system corresponding to any time interval in weak relevant time period is derived Harmonic impedance calculation formula, it is based on the harmonic impedance calculation formula derived, time intervals all in weak relevant time period are corresponding Impedance value system side harmonic impedance value of the mean value as the method for the present invention PCC point.
Further, in the system side harmonic impedance estimation method based on amendment independent random vector, the step A equation The foundation of group comprises the steps of:
A1, it is based on PCC point harmonic data, obtains the expression formula of the auto-correlation function of harmonic current and harmonic voltage are as follows:
And the expression formula of the cross-correlation function of harmonic voltage and harmonic current are as follows:
In formula,PCC point harmonic current and harmonic voltage are respectively indicated, N indicates that sample point sum, n represent letter Any moment of number sequence, k indicate the time interval of two signals, and E { } indicates mathematic expectaion, and conjugation is sought in * expression.
A2, it is screened by the threshold value to normalizated correlation coefficient, determines weak relevant time period.Harmonic current and harmonic voltage The expression formula of normalized autocorrelation coefficient are as follows:
The expression formula of the normalized-cross-correlation function of harmonic current and harmonic voltage are as follows:
In probability theory and statistics, related coefficient shows the intensity of linear relationship and direction between two stochastic variables. The attenuation degree for being spaced k at any time of correlation function is usually characterized with normalizated correlation coefficient.Normalizated correlation coefficient is just Negative sign only indicates relevant direction, the relevant degree of absolute value representation.It may be defined as when normalizated correlation coefficient is less than 0.3 weak It is related.
The present invention sets the threshold value of metric data related coefficient as 0.3, remembers the auto-correlation coefficient of harmonic current at PCC for the first time It is k1 less than time interval corresponding to threshold value, remembers that the cross-correlation coefficient of harmonic voltage and electric current is less than threshold value pair for the first time at PCC The time interval answered is k2, and maximum value therein is denoted as kmax, for N number of sample, define k value model corresponding to weak relevant time period It encloses for kmax~N-1.
Further, in the system side harmonic impedance estimation method based on amendment independent random vector, in the step B Formula equivalent transformation the following steps are included:
B1, equivalent circuit figure, column write circuit equation:
In formula, ZuIndicate system side harmonic impedance,Indicate the vector value of system side harmonic current source,Indicate system side The vector value of harmonic voltage source;
Formula (6) both ends are multiplied togetherCorrelation function is sought it afterwards, as a result as follows:
When metric data is that second-order statistic is ergodic, formula (7) is approximately:
In B2, the weak relevant time period determined in A, the cross correlation and harmonic current of harmonic voltage and harmonic current Autocorrelation is all very weak, it can thus be assumed that PCC point harmonic current has very weak correlation, formula with system side harmonic voltage source (8) 0 is seen as.
In system side harmonic impedance estimation method based on amendment independent random vector, system harmonics hinder in the step C The derivation of anti-calculation formula the following steps are included:
C1, formula (8) see 0 as, derive the table of system harmonic impedance corresponding to any time interval in weak relevant time period It is up to formula
In formula, the value range of k is kmax~N-1.
C2, to obtained (N-kmax) a impedance averages, obtain the harmonic impedance based on correction independent random vector Are as follows:
Principle based on the above embodiment, the present embodiment disclose specific embodiment:
Simulation model is built according to Fig. 1 equivalent circuit, design parameter setting is as follows.
1) user-side harmonic electric currentAmplitude be 100A,Amplitude is made within the entirely estimation period relative to initial value 20% sinusoidal fluctuation and ± 15% random fluctuation;System side harmonic wave ource electric currentAmplitude is user side electric currentThe n of amplitude Again (n takes 0.3,0.5,0.8,1.0,1.2), make the random wave relative to amplitude size ± 15% within the entirely estimation period It is dynamic.Phase angle initial value is -30 °,Phase angle initial value is 30 °, adds the random perturbation of relative phase angle size ± 10%.
2) system side impedance ZuPerseverance is (5+j10) Ω.Two kinds of impedance conditions, i.e. user side impedance Z are setcRespectively (5+ J10) Ω and (20+j40) Ω, influence of the comparative analysis impedance ratio to Algorithm Error.ZuAnd ZcAdded within the entirely estimation period 10% sinusoidal fluctuation.
100 data are emulated, every 40 data are a period to carry out recurrence calculation, and the mean value of 61 segment datas is as one The impedance magnitude of secondary operation and the error amount of phase angle, program are run 100 times, and taking the root-mean-square value of 100 operation results to be used as should The error result of method.Meanwhile comparing the calculated result of context of methods and other methods.(method 1 is undulate quantity method, and method 2 is Multiple ICA method, method 3 is independent vector covariance method, and method 4 is binary regression method, and method 5 is context of methods).
The present embodiment is divided into 2 scenes.In scene 1, system side impedance ZuFor (5+j10) Ω, user side impedance ZcFor (20+ J40) Ω, two sides impedance ratio are 4, and the RMS error of system harmonic impedance amplitude and phase angle is respectively such as table 1, shown in 2.
1 impedance magnitude error user side impedance of table/system surveys impedance=4
2 impedance angle error user side impedance of table/system surveys impedance=4
In scene 2, system side impedance ZuFor (5+j10) Ω, user side impedance ZcFor (5+j10) Ω, two sides impedance ratio is 1, The RMS error of system harmonic impedance amplitude and phase angle is respectively such as table 3, shown in 4.
3 impedance magnitude error user side impedance of table/system surveys impedance=1
4 impedance angle error user side impedance of table/system surveys impedance=1
Embodiment two
Referring to FIG. 2, the present invention also provides a kind of system side harmonic impedance estimations based on amendment independent random vector System, the system comprises:
Weak relevant time period obtaining unit is analyzed for the correlation to PCC point harmonic data, obtains PCC point harmonic wave The normalized autocorrelation coefficient of harmonic current at the normalized-cross-correlation function and PCC of electric current and harmonic voltage, based on described mutual Related coefficient and the auto-correlation coefficient obtain weak relevant time period by threshold value screening;
Correlation function obtaining unit, the equivalent circuit for establishing under certain subharmonic are fixed based on Kirchoff s voltage Rule establishes circuit equation, and the circuit equation of foundation is carried out equivalent transformation, seeks the correlation function of circuit equation after equivalent transformation;
Harmonic impedance value computing unit, the correlation function for being obtained based on correlation function obtaining unit are derived described Power grid points of common connection system side harmonic impedance calculation formula corresponding to any time interval in weak relevant time period, based on deriving Harmonic impedance calculation formula, using the mean value of the corresponding impedance value of time intervals all in weak relevant time period as PCC dot system Side harmonic impedance value.
Although preferred embodiments of the present invention have been described, it is created once a person skilled in the art knows basic Property concept, then additional changes and modifications may be made to these embodiments.So it includes excellent that the following claims are intended to be interpreted as It selects embodiment and falls into all change and modification of the scope of the invention.
Obviously, various changes and modifications can be made to the invention without departing from essence of the invention by those skilled in the art Mind and range.In this way, if these modifications and changes of the present invention belongs to the range of the claims in the present invention and its equivalent technologies Within, then the present invention is also intended to include these modifications and variations.

Claims (10)

1. the system side harmonic impedance estimation method based on amendment independent random vector, which is characterized in that the described method includes:
Step A: analyzing the correlation of PCC point harmonic data, obtains the normalization of PCC point harmonic current and harmonic voltage The normalized autocorrelation coefficient of cross-correlation coefficient and harmonic current is based on the cross-correlation coefficient and the auto-correlation coefficient, Given threshold screening obtains weak relevant time period;
Step B: establishing the equivalent circuit under certain subharmonic, establishes circuit equation based on Kirchhoff's second law, will build Vertical circuit equation carries out equivalent transformation, seeks the correlation function of circuit equation after equivalent transformation;
Step C: it based on the correlation function obtained in step B, derives in the weak relevant time period corresponding to any time interval PCC dot system side harmonic impedance calculation formula can be obtained in weak relevant time period based on the harmonic impedance calculation formula derived The corresponding impedance value of all time intervals, using the mean value of the corresponding impedance value of time intervals all in weak relevant time period as PCC Dot system side harmonic impedance value.
2. the system side harmonic impedance estimation method according to claim 1 based on amendment independent random vector, feature It is, in the weak relevant time period, the correlation function of circuit equation sees 0 as after equivalent transformation.
3. the system side harmonic impedance estimation method according to claim 1 based on amendment independent random vector, feature It is, the step A further includes that the cross-correlation coefficient and the auto-correlation coefficient are normalized.
4. the system side harmonic impedance estimation method according to claim 1 based on amendment independent random vector, feature It is, the step A is specifically included:
Step A1: it is based on PCC point harmonic data, obtains the expression formula of the auto-correlation function of harmonic current and harmonic voltage are as follows:
The expression formula of the cross-correlation function of harmonic voltage and harmonic current are as follows:
Wherein, PCC point harmonic current and harmonic voltage are respectively indicated, N indicates sample point sum, n representation signal sequence Any moment, k indicate two signals time interval, E { } indicate mathematic expectaion, * expression seek conjugation;
Step A2: the auto-correlation function of cross-correlation function and harmonic current based on harmonic voltage and harmonic current obtains phase Relationship number, given threshold screen related coefficient, determine weak relevant time period.
5. the system side harmonic impedance estimation method according to claim 4 based on amendment independent random vector, feature It is, sets the threshold value of metric data related coefficient;The auto-correlation coefficient of harmonic current is less than for the first time corresponding to threshold value at note PCC Time interval be k1, remember at PCC that harmonic voltage and the cross-correlation coefficient of electric current are less than the corresponding time interval of threshold value for the first time and are K2, maximum value therein are denoted as kmax;For N number of sample, defining k value range corresponding to weak relevant time period is kmax~N-1.
6. the system side harmonic impedance estimation method according to claim 1 based on amendment independent random vector, feature It is, step B includes:
Step B1: according to equivalent circuit figure, column write the circuit equation of equivalent circuit:
In formula, ZuIndicate system side harmonic impedance,Indicate the vector value of system side harmonic current source,Indicate system side harmonic wave The vector value of voltage source;
Formula (4) both ends are multiplied togetherCorrelation function is sought it afterwards, as a result are as follows:
When metric data is that second-order statistic is ergodic, formula (5) transformation are as follows:
Step B2: in the weak relevant time period determined in step, formula (6) sees 0 as.
7. the system side harmonic impedance estimation method according to claim 6 based on amendment independent random vector, feature It is, based on the correlation function obtained in step B, derives that power grid corresponding to any time interval is public in the weak relevant time period Tie point system side harmonic impedance calculation formula altogether, specifically includes:
Step C1: formula (6) sees 0 as, derives the table of system harmonic impedance corresponding to any time interval in weak relevant time period Up to formula are as follows:
Wherein, the value range of k is kmax~N-1;
Step C2: to obtained N-kmaxA impedance is averaged, and the harmonic impedance based on correction independent random vector is obtained are as follows:
8. a kind of system side harmonic impedance estimating system based on amendment independent random vector, which is characterized in that the system packet It includes:
Weak relevant time period obtaining unit, for PCC point harmonic data carry out correlation analysis, obtain PCC point harmonic current with The normalized autocorrelation coefficient of normalized-cross-correlation function and PCC the point harmonic current of harmonic voltage is based on the cross correlation The several and auto-correlation coefficient, given threshold screening obtain weak relevant time period;
Correlation function obtaining unit, the equivalent circuit for establishing under certain subharmonic, is built based on Kirchhoff's second law The circuit equation of foundation is carried out equivalent transformation, seeks the correlation function of circuit equation after equivalent transformation by vertical circuit equation;
Harmonic impedance value computing unit, the correlation function for being obtained based on correlation function obtaining unit derive the weak phase The calculation formula for closing PCC dot system side harmonic impedance corresponding to any time interval in the period, based on the harmonic wave resistance derived The corresponding impedance value of all time intervals in weak relevant time period can be obtained, by institute's having time in weak relevant time period in anti-calculation formula The mean value of corresponding impedance value is spaced as PCC dot system side harmonic impedance value.
9. the system side harmonic impedance estimating system according to claim 8 based on amendment independent random vector, feature It is, weak relevant time period obtaining unit is specifically used for:
Based on PCC point harmonic data, the expression formula of the auto-correlation function of harmonic current and harmonic voltage is obtained are as follows:
And the cross-correlation function of harmonic voltage and harmonic current:
Wherein,PCC point harmonic current and harmonic voltage are respectively indicated, N indicates sample point sum, n representation signal sequence Any moment of column, k indicate the time interval of two signals, and E { } indicates mathematic expectaion, and conjugation is sought in * expression;
The auto-correlation function of cross-correlation function and harmonic current based on harmonic voltage and harmonic current, obtains related coefficient, The threshold value of related coefficient is screened, determines weak relevant time period.
10. the system side harmonic impedance estimating system according to claim 9 based on amendment independent random vector, feature It is, sets the threshold value of metric data related coefficient, remembers that the auto-correlation coefficient of harmonic current at PCC is less than for the first time corresponding to threshold value Time interval be k1, remember at PCC that harmonic voltage and the cross-correlation coefficient of electric current are less than the corresponding time interval of threshold value for the first time and are K2, maximum value therein are denoted as kmax, for N number of sample, defining k value range corresponding to weak relevant time period is kmax~N-1.
CN201910774203.9A 2019-08-21 2019-08-21 System side harmonic impedance estimation method and system based on corrected independent random vector Active CN110456159B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910774203.9A CN110456159B (en) 2019-08-21 2019-08-21 System side harmonic impedance estimation method and system based on corrected independent random vector

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910774203.9A CN110456159B (en) 2019-08-21 2019-08-21 System side harmonic impedance estimation method and system based on corrected independent random vector

Publications (2)

Publication Number Publication Date
CN110456159A true CN110456159A (en) 2019-11-15
CN110456159B CN110456159B (en) 2020-03-17

Family

ID=68488289

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910774203.9A Active CN110456159B (en) 2019-08-21 2019-08-21 System side harmonic impedance estimation method and system based on corrected independent random vector

Country Status (1)

Country Link
CN (1) CN110456159B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112285421A (en) * 2020-10-16 2021-01-29 国网四川省电力公司电力科学研究院 System side harmonic impedance estimation method based on minimum norm criterion
CN114660362A (en) * 2022-05-19 2022-06-24 四川大学 System side harmonic impedance estimation method and device based on mutual information data optimization
CN116467556A (en) * 2023-03-28 2023-07-21 四川大学 Harmonic impedance calculation method considering background harmonic voltage fluctuation and impedance variation

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020180458A1 (en) * 2001-03-29 2002-12-05 Ho Eddy Ying Yin Method and system for automatic determination of inductance
CN102998535A (en) * 2012-12-10 2013-03-27 华北电力大学(保定) Method for computing harmonic impedance of system based on maximum likelihood estimation theory
US20170310348A1 (en) * 2016-04-22 2017-10-26 Uti Limited Partnership Method and System for Extraction of Noise Parameters of Nonlinear Devices
CN107576853A (en) * 2017-07-12 2018-01-12 国网上海市电力公司 Distribution Network Harmonics impedance computation method based on canonical correlation analysis
CN108089059A (en) * 2017-12-19 2018-05-29 西南交通大学 A kind of harmonic impedance evaluation method based on sparse independent component analysis
CN109409762A (en) * 2018-11-07 2019-03-01 国网四川省电力公司成都供电公司 Based on the harmonic contributions division methods for improving data relation analysis
CN109839538A (en) * 2019-03-29 2019-06-04 云南电网有限责任公司电力科学研究院 A kind of harmonic source identification method and system based on correlation analysis

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020180458A1 (en) * 2001-03-29 2002-12-05 Ho Eddy Ying Yin Method and system for automatic determination of inductance
CN102998535A (en) * 2012-12-10 2013-03-27 华北电力大学(保定) Method for computing harmonic impedance of system based on maximum likelihood estimation theory
US20170310348A1 (en) * 2016-04-22 2017-10-26 Uti Limited Partnership Method and System for Extraction of Noise Parameters of Nonlinear Devices
CN107576853A (en) * 2017-07-12 2018-01-12 国网上海市电力公司 Distribution Network Harmonics impedance computation method based on canonical correlation analysis
CN108089059A (en) * 2017-12-19 2018-05-29 西南交通大学 A kind of harmonic impedance evaluation method based on sparse independent component analysis
CN109409762A (en) * 2018-11-07 2019-03-01 国网四川省电力公司成都供电公司 Based on the harmonic contributions division methods for improving data relation analysis
CN109839538A (en) * 2019-03-29 2019-06-04 云南电网有限责任公司电力科学研究院 A kind of harmonic source identification method and system based on correlation analysis

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
JIN HUI: "Assessing Utility Harmonic Impedance Based on the Covariance Characteristic of Random Vectors", 《IEEE TRANSACTIONS ON POWER DELIVERY》 *
TOMAZ PFAJFAR: "Harmonic Contributions Evaluation With the Harmonic Current Vector Method", 《IEEE TRANSACTIONS ON POWER DELIVER》 *
ZHENXIN ZHONG: "Assessment Method of Wind Farm Harmonic Emission Level Based on Improved Covariance Characteristic of Random Vectors", 《2018 CHINA INTERTIONAL CONFERENCE ON ELECTRICITY DISTRIBUTION》 *
张龙 等: "谐波阻抗和背景谐波电压估算方法的对比研究", 《电气技术》 *
惠锦 等: "基于独立随机矢量协方差特性的谐波发射水平评估方法", 《电力系统自动化》 *
赵俊屹 等: "基于二阶盲辨识的系统谐波阻抗计算", 《电测与仪表》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112285421A (en) * 2020-10-16 2021-01-29 国网四川省电力公司电力科学研究院 System side harmonic impedance estimation method based on minimum norm criterion
CN112285421B (en) * 2020-10-16 2023-05-23 国网四川省电力公司电力科学研究院 System side harmonic impedance estimation method based on minimum norm criterion
CN114660362A (en) * 2022-05-19 2022-06-24 四川大学 System side harmonic impedance estimation method and device based on mutual information data optimization
CN116467556A (en) * 2023-03-28 2023-07-21 四川大学 Harmonic impedance calculation method considering background harmonic voltage fluctuation and impedance variation
CN116467556B (en) * 2023-03-28 2023-12-15 四川大学 Harmonic impedance calculation method considering background harmonic voltage fluctuation and impedance variation

Also Published As

Publication number Publication date
CN110456159B (en) 2020-03-17

Similar Documents

Publication Publication Date Title
CN110456159A (en) System side harmonic impedance estimation method and system based on amendment independent random vector
Suonan et al. An accurate fault location algorithm for transmission line based on R–L model parameter identification
CN109409762A (en) Based on the harmonic contributions division methods for improving data relation analysis
CN108918932A (en) Power signal adaptive filter method in load decomposition
Wang et al. Faulty feeder detection based on the integrated inner product under high impedance fault for small resistance to ground systems
Barros et al. A new method for measurement of harmonic groups in power systems using wavelet analysis in the IEC standard framework
CN108535613A (en) A kind of voltage flicker parameter detection method based on combination window function
CN108170885A (en) Multi-harmonic Sources recognition methods in a kind of power distribution network
Kumar et al. Efficient detection of false data injection attacks on AC state estimation in smart grids
Zhang et al. Estimation of interarea electromechanical modes during ambient operation of the power systems using the RDT–ITD method
CN106651150A (en) Method and device for assessing harmonic influence on power grid
CN107167658B (en) A kind of jamproof electric system fundamental frequency of high-precision and Method for Phase Difference Measurement
Samet et al. Improvement of reactive power calculation in electric arc furnaces utilising Kalman filter
Jiao et al. An approach for electrical harmonic analysis based on interpolation DFT
Hwang et al. Identification of interarea modes from an effectual impulse response of ringdown frequency data
CN110069883A (en) A kind of lower transformer core saturation unstability method of discrimination of direct current disturbance
Niewiadomski et al. Comparative analysis of deterministic and random modulations based on mathematical models of transmission errors in series communication
Redondo et al. A strategy for improving the accuracy of flicker emission measurement from wind turbines
Shu et al. Novel estimation method of utility harmonic impedance based on short‐term impedance minimum variance criterion
Zhuang et al. Four harmonic analysis and energy metering algorithms based on a new cosine window function
CN112285421B (en) System side harmonic impedance estimation method based on minimum norm criterion
Inamdar et al. A review of methods employed to identify flicker producing sources
Martins et al. A new approach for EPLL-based frequency estimation under severe disturbances
Ahmed et al. Application of modal transformations for power system harmonic impedance measurement
Rodríguez-Pajarón et al. Residential Harmonic Injection Models Based on Field Measurements

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