CN109409762A - Based on the harmonic contributions division methods for improving data relation analysis - Google Patents
Based on the harmonic contributions division methods for improving data relation analysis Download PDFInfo
- Publication number
- CN109409762A CN109409762A CN201811320306.XA CN201811320306A CN109409762A CN 109409762 A CN109409762 A CN 109409762A CN 201811320306 A CN201811320306 A CN 201811320306A CN 109409762 A CN109409762 A CN 109409762A
- Authority
- CN
- China
- Prior art keywords
- harmonic
- impedance
- voltage
- user
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 91
- 238000004458 analytical method Methods 0.000 title claims abstract description 24
- 239000013598 vector Substances 0.000 claims description 13
- 230000011218 segmentation Effects 0.000 claims description 6
- 230000005611 electricity Effects 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 238000001914 filtration Methods 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 8
- 230000008569 process Effects 0.000 abstract description 5
- 238000007796 conventional method Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 11
- 238000012098 association analyses Methods 0.000 description 6
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 238000012417 linear regression Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010891 electric arc Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001373 regressive effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- 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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06395—Quality analysis or management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
The invention discloses a kind of based on the harmonic contributions division methods for improving data relation analysis, comprising: step 1: in the case where considering harmonic wave phase angle, establishing the relationship between PCC point h subharmonic voltage and harmonic current；Step 2: the characteristic estimation power system impedance angle for being zero based on independent random variable covariance；Step 3: based on the relationship between power system impedance angle and harmonic voltage and harmonic current, calculating Harmonious Waves in Power Systems impedance；Step 4: being based on Harmonious Waves in Power Systems impedance, the harmonic voltage for obtaining user side and electric system side in PCC point is contributed；User's harmonic impedance is considered in calculating process, reduces conventional method due to ignoring user's harmonic impedance bring error；Simulation analysis and Practical Project example, which demonstrate this method than existing methods, can effectively inhibit system harmonics variation and user's harmonic impedance bring to influence, and obtain more accurate harmonic contributions division result.
Description
Technical field
The present invention relates to harmonic contributions to divide field, and in particular, to a kind of based on the harmonic wave for improving data relation analysis
Divisions of responsibility method.
Background technique
With the development of industry, various power electronic equipments and nonlinear-load are answered extensively in each department of national economy
With the harmonic pollution problems in electric system especially distribution are got worse.It, can be in public interface after harmonic current flows into power grid
(PCC point) generates harmonic voltage, causes voltage distortion, seriously affects power quality.Therefore, in PCC point to electric system and
The harmonic voltage responsibility of user is rationally assessed, and is the premise for realizing Power quality management and control.
The method for dividing harmonic contributions both at home and abroad at present mainly surrounds system and user's harmonic impedance is estimated to be unfolded, main
It wants that " intervening formula " [method and " non-intervention formula " method can be divided into." intervening formula " method mainly passes through artificially to system harmonic
Electric current carrys out measuring system side harmonic impedance, and this method, which may operate normally electric system, to have an adverse effect, and has
Biggish evaluated error." non-intervention formula " method estimates harmonic impedance by the measured value of PCC point harmonic voltage electric current, quantitatively to draw
Subharmonic responsibility, mainly has: 1, " undulate quantity method " is the symbol based on institute's measuring point magnitude of a voltage fluctuation to current wave momentum ratio
The estimation method of feature, the accuracy requirement that this method measures harmonic parameters is higher, while it is sufficiently large also to need measured value to have
Fluctuation；2, linear regression method: including robustness regression method, binary regression method, Partial Least Squares etc., essence is in accordance with Dai Wei
Southern equivalent circuit or promise equivalent circuit, construct corresponding regression equation, by the real and imaginary parts minute of PCC point voltage and current phasor
Regressive Solution harmonic impedance is opened, so not being the accurate least square solution of former problem.But current linear regression method exists as follows
Limitation: the constant term (i.e. background harmonic voltage) in regression model is really a fluctuation larger quantities, results in recurrence
Error is larger.Data relation analysis in existing part is nearest method, is to change computing system using user's harmonic current
Harmonic impedance, since this method has ignored harmonic wave phase angle, this method can equally generate large error.
Relationship and user's harmonic current under traditional association analysis and utilization same frequency between PCC point harmonic voltage electric current
Variation carry out computing system impedance.If Fig. 1 a is equivalent circuit, subscript ' s ' and ' c ' respectively refer to system side and user side.It is measuring
Measured value is obtained in period, and obtains PCC point h subharmonic voltage with Fourier transformationAnd harmonic currentIt drawsIt closes
InImage such as Fig. 1 b, and can be obtained with linear regression method:
A is obtained according to formula (1) and is equal to system impedance, and it is zero that b, which is equal to user's harmonic current,When PCC point harmonic voltage it is flat
Mean value, i.e. system harmonics voltage.Known by the condition of data relation analysis, system harmonics voltage (b) answers PCC point harmonic currentAnd system impedance (a) Xiang Duli.Ignore user's impedance, thenEqual to user's harmonic currentWith system harmonics voltage b
It is unrelated, meet conditions above, system impedance and system harmonics voltage can be obtained by this method.
Known by formula (1), when traditional association analytical calculation system impedance has ignored harmonic wave phase angle, and calculates user's harmonic voltage
User's impedance is had ignored when contribution, in this way when user's impedance is smaller, as a result can generate large error；In addition, system in practice
Harmonic voltage not instead of constant, a variable quantity.
Summary of the invention
The present invention provides a kind of based on the harmonic contributions division methods for improving data relation analysis, solves existing side
The deficiency of method considers user's harmonic impedance in calculating process, reduces conventional method due to ignoring user's harmonic impedance band
The error come.Simulation analysis and Practical Project example demonstrate this method effectively can inhibit system harmonics to change than existing methods
It is influenced with user's harmonic impedance bring, obtains more accurate harmonic contributions division result.
For achieving the above object, this application provides based on the harmonic contributions division side for improving data relation analysis
Method, which comprises
Step 1: in the case where considering harmonic wave phase angle, establishing the pass between PCC point h subharmonic voltage and harmonic current
System；
Step 2: the characteristic estimation power system impedance angle for being zero based on independent random variable covariance；
Step 3: based on the relationship between power system impedance angle and harmonic voltage and harmonic current, calculating electric system
Harmonic impedance；
Step 4: being based on Harmonious Waves in Power Systems impedance, obtain user side and electric system side in the harmonic voltage tribute of PCC point
It offers.
Further, it further comprises the steps of: to filter out using time series segmentation method before step 1 and meets analysis condition
Harmonic voltage current data.
This method is quantitatively to divide Harmonics source customer in the harmonic contributions of points of common connection (PCC), is proposed with improvement number
According to association analysis computing system harmonic impedance and the method for dividing harmonic contributions.It is filtered out first with time series segmentation method full
The harmonic voltage current data of sufficient analysis condition.System harmonic impedance is acquired using the characteristic that independent random variable covariance is zero
Phase angle, then harmonic wave phase angle is introduced in available data relation analysis model, in the fluctuation of user's harmonic current, utilize user's harmonic wave
The fluctuation of electric current acquires accurate system harmonic impedance modulus value.When user's harmonic current is steady, system harmonics voltage is utilized
Fluctuation quantization user's harmonic voltage responsibility.This method considers user's harmonic impedance in calculating process, reduces conventional method
Due to ignoring user's harmonic impedance bring error.Simulation analysis and Practical Project example demonstrate this method than existing methods
It can effectively inhibit system harmonics variation and user's harmonic impedance bring to influence, obtain more accurate harmonic contributions division result.
More accurate system harmonic impedance in order to obtain, this method have fully considered harmonic wave phase angle on the basis of existing.
In addition, this method proposes that a kind of new harmonic contributions quantify division methods, this method utilizes system harmonics voltage (background harmonics electricity
Pressure) variation quantitatively calculate user and contributed in the harmonic voltage of PCC point, and user's harmonic impedance is considered in calculating process.
Theoretical simulation and engineering verification are carried out to this method, the results showed that the more existing correlation fractal dimension of this method can effectively inhibit system humorous
Wave variation and influence of user's harmonic impedance to result, assessment result accuracy is higher, there is certain theoretical value and engineer application
Prospect.
Method in the application includes:
Computing system harmonic impedance:
Harmonic wave phase angle is not considered in formula (1), therefore respectively amount is scalar in formula, can generate mistake in actually calculating in this way
Difference.The error can be reduced using this method, formula (2) indicates the pass in Fig. 1 a between PCC point h subharmonic voltage and harmonic current
System:
WhereinFor system impedance,For system harmonics voltage.(2) polar plot of formula is as shown in Fig. 2, consider each vector phase
Angle, formula (2), which is rewritten into following scalar form, to be had:
γ in above formulahFor vectorWithBetween angle, λhFor vectorWithBetween angle, can by Fig. 2
Find out WithRespectivelyWith?On projection.According to γhDefinition can γh
=αh+βh, wherein αhFor vectorWithPhase angle difference, βhFor system impedance angle, polar plot such as Fig. 3.Known by definition: WithRespectively vectorWithPhase angle), but βhIt need to be estimated for unknown quantity.This
The characteristic estimation system impedance angle beta that method is zero based on independent random variable covarianceh。
The characteristic estimation system impedance angle for being zero based on independent random variable covariance:
It is N number of in PCC point samplingWithValue, in formula (2) respectively amount average, have:
It enablesFormula (2) subtracts formula (4)
Have:
Due to user's impedanceMuch larger than system impedanceIt is believed thatI.e.Mainly byContribution, withSolely
It is vertical.
According to probability theory, the covariance of two independent random vectors is zero, thereforeWithCovariance be approximately 0.That is:
Wherein * indicates the conjugation of vector.It can be approximately mean value by expectation, therefore formula (6) can when sample size is sufficiently large
It is written as:
By formula (5) both sides simultaneously multiplied byIt averages again:
Formula (8), which is split into real and imaginary parts two parts, to be obtained:
Each amount on the right of equation can be obtained by sample value.Obtain βhAfterwards willIt is rightRecurrence can obtain system
Impedance magnitudeSuch as Fig. 1 b.
Harmonic contributions divide
User's impedance is ignored in traditional harmonic contributions division, and large error can be brought when user's impedance is smaller.It is existing
Method is assumedFor steady state value, but in practiceIt is variation.The main thought of this method is before considering user's impedance
It puts, usesContribution of the variation assessment user to PCC point harmonic voltage.It is obtained by calculated above.User's harmonic voltageWith Between relationship are as follows:
WhereinFor user's impedance.ObviouslyMain component fromThereforeWithIt is not independent mutually.For
It is set to meet the condition of association analysis, it willWithSubstitution formula (10):
WhereinFor system harmonics electric current.It willFollowing formula can be obtained by moving on to the equation left side:
Have:
ByWithDefinition know,WithIt independently of each other, then can be rightWithMake association analysis.It is sharp first
It is calculated with Fig. 4 and known quantitySecondly component can be obtained with recurrenceAnd
The respectively slope and Y intercept of formula (13).The former is system and subscriber multiple impedance, and the latter is user's harmonic voltage in PCC
The contribution of point.It can with due regard to user's impedance when computing system and user's harmonic voltage are contributed using this methodBut it is not required to
It obtainsExact value:
And have:
WhereinRespectively user side and system side are contributed in the harmonic voltage of PCC point,It is humorous for user
Wave voltage,It is contributed for user's harmonic voltage in PCC point.
Data selection scheme based on time series segmentation
Computing system impedance mainly utilizesVariation, divide harmonic contributions then utilizeVariation, therefore in Practical Project
In suitable data segment need to be selected to be respectively calculated.This method uses the data selection scheme based on time series segmentation, i.e.,
Select each amount (Or) the fluctuation biggish period.A fluctuation threshold can be preset in practice (is in this method
10%), select undulate quantity be more than the threshold value period be concern the period.Fluctuation threshold is defined as: the electric current in measurement period
The ratio between the sum of electric current (voltage) absolute value of (voltage) value/in this time.
Such as Fig. 5, within the T1 period(ignore user's impedance i.e.) fluctuation is larger, then it is T1 sections optional in measured value meter
Calculate system impedance；In T2 sectionsSteadily (i.e.Substantially without fluctuation), andIt fluctuates larger, it is known that system harmonics electricity in the section
Press (system harmonics electric current) fluctuation is larger, then it is T2 section optional in data division harmonic contributions.
One or more technical solution provided by the present application, has at least the following technical effects or advantages:
Accurate computing system harmonic impedance is the premise realizing harmonic contributions and quantitatively dividing, and this method is for traditional association point
Ignore the deficiency of harmonic wave phase angle when analysis method solving system impedance, proposes to consider harmonic wave phase angle in equation, to obtain more smart
True system harmonic impedance.In addition, this method ignores user's impedance and vacation when dividing harmonic contributions for the conventional linear Return Law
If system harmonics voltage is the deficiency of definite value, proposes the fluctuation for utilizing system harmonics voltage, obtain user in conjunction with correlation fractal dimension
In the quantitative harmonic voltage contribution of PCC point.Simulation analysis and Practical Project example demonstrate this method and can obtain than existing methods
To more accurate system impedance, more satisfied result has been obtained to the division of harmonic contributions.
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 a is equivalent circuit schematic diagram when ignoring user's harmonic impedance；
Fig. 1 b isWithRelation schematic diagram；
Fig. 2 is equilibrium relationships polar plot；
Fig. 3 is harmonic wave phase angle polar plot；
Fig. 4 is equivalent circuit schematic diagram when considering user's harmonic impedance；
Fig. 5 is time series segmentation schematic diagram；
Fig. 6 is IEEE14 node standard test system schematic diagram；
Fig. 7 a isWithRelation schematic diagram；
Fig. 7 b isWithRelation schematic diagram；
Fig. 8 a-b is 5 subharmonic current schematic diagrames of harmonic source HS injection；
Fig. 9 is 5 subharmonic voltage of PCC pointWith system harmonics ource electric currentRelation schematic diagram；
Figure 10 is 35kV/10kV industrial system schematic equivalent circuit；
Figure 11 a-b is 5 subharmonic voltage electric current modulus value schematic diagram of PCC point；
Figure 12 is the actual measurement of PCC pointSchematic diagram.
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.
Emulation experiment
This method does example using the IEEE14 node standard test system in Fig. 6, and 1-14 refer to system loading in Fig. 6,
HL1 is the non-linear harmonic wave source load of PCC point concern, and L2 is linear load, and HS is system harmonics source load.Use Digsilent
Software emulation.No. 11 nodes are chosen as concern bus (PCC point), HL1 is the non-linear harmonic wave source load of PCC point concern, and L2 is
Linear load, HS are system harmonics source load.The other parts for removing load HL1, L2 are considered as system side, example will be by imitative
The true harmonic voltage contribution for calculating harmonic-producing load HL1 at No. 11 buses.By taking 5 subharmonic as an example.
Computing system harmonic impedance:
It is scanned to obtain 5 subharmonic system impedance actual value of PCC point with Digsilent software harmonic impedance first Ignore user's impedance in the portion.
The 5 subharmonic current amplitudes of harmonic source HS are set as 114.8A, phase angle is -76.56 °, ignores user's impedance.First
It specifies N number ofSample value carries out harmonic flow calculation to each sample value, acquires N number of harmonic voltage electric current and its phase in PCC point
Angle value, by N number ofly defined aboveValue, can obtain system harmonics phase angle estimated value using formula (9)With actual valueError Absolute Value be only 2.25 °, can be ignored.?After can obtain γ5Vector: It willIt is rightIt returns and acquires system impedance modulus value
Fig. 7 a and Fig. 7 b are the simulation result obtained respectively with method one, two, and method one, two respectively refers to traditional association analysis
And this method.Know that with the system impedance that method one obtains be 2.54 Ω by figure, obtaining system impedance with method two is 2.75 Ω.It is right
It is 4.18% with this method error, hence it is evident that lower than the error 11.49% of conventional method than 2.87 Ω of actual value.
The system impedance value that table 1 is obtained with method one, two respectively when being user's impedance difference under 5 subharmonicFrom table 1
It can be seen that, the precision of method two is apparently higher than method one, while can obtain the system impedance value obtained by method twoUser is hindered
It is anti-The susceptibility of variation is lower than method one.
The system impedance value that 1 user's impedance of table does not obtain simultaneously
Harmonic contributions quantitatively divide
5 subharmonic currents of system harmonics source load HS are using classical curve, and 1 minute correspondence, one sample point, one day altogether
There are 1440 points, as shown in Figure 8.
Assuming that user's harmonic current is definite value, size 10A, phase angle is -74.25 °.System impedance is 79.75 ° of 2.87 ∠
Ω, user's impedance are 78.2 ∠, 90 ° of Ω.The simulation result obtained with this method is as shown in Figure 9.User is in PCC point as seen from the figure
Harmonic voltage be 29.2V, it is focus that take 5 subharmonic voltage of PCC point, which be 50V, then the contribution of user's harmonic voltage is 58.4%,
It is very close with actual value 55.4%.In order to prove the accuracy of this method, harmonic contributions are differently calculated respectively: being ignored
The linear regression method (method 1) of user's impedance, traditional association analyze (method 2), this method (method 3).
User's harmonic voltage contribution rate that table 2 is obtained with three kinds of methods respectively when being user's impedance difference.The each column from table
It is found that this method accuracy is apparently higher than method 1,2, and as the reduction of user's harmonic impedance is more obvious.Because of method 1,2
User's impedance is had ignored in calculating process (to regard in equivalent circuit),Smaller, the error generated is got over
Greatly.Since this method considers user's impedance, even ifWhen smaller, as a result also within the error range.
5 subharmonic voltage of the user contribution that 2 user's impedance of table does not obtain simultaneously
Case Simulation:
Figure 10 is the circuit equivalent figure of certain 35kV/10kV industrial system, is contained at capacity of short circuit 10MVA, 35kV bus
2.5% 5 order harmonic components are connected to linear load and nonlinear-load electric arc furnaces at 10kV.Sample frequency 3200Hz, per minute
It is interior that each harmonic measured value can be obtained to sampled data progress Fast Fourier Transform (FFT).PCC point 5 times in surveyed continuous 300min one day
Harmonic voltageAnd electric currentComponent size is as shown in figure 11.Preceding 240min load operates normally, rear 60min load synthesis.
By context data selection scheme it is found that harmonic current undulate quantity is big during startup by user in T1 sections in figure,
Therefore data computing system impedance in T1 sections can be used, obtaining 5 subharmonic impedance angle of system by this method is 67.89 °, size 2.16
Ω。
Selection T2 segment data asks user's harmonic voltage to contribute, and table 3 is to takeWhen 95% probability value is 446.5V, use respectively
The system and user that traditional association analysis and this method obtain are contributed in 5 subharmonic voltage of PCC point.
Table 3WithCalculated result comparison
Using 5 subharmonic voltages in measured data when load synthesis as system harmonics voltage contributions, as shown in figure 12, obtain
It is 175.2V to its 95% probability value.3 data of contrast table, it is known that the calculated value of this method closest to this as a result, hence it is demonstrated that
User is 267.9V's the result is that believable in the harmonic voltage that PCC point generates.
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 (7)
1. based on the harmonic contributions division methods for improving data relation analysis, which is characterized in that the described method includes:
Step 1: in the case where considering harmonic wave phase angle, establishing the relationship between PCC point h subharmonic voltage and harmonic current；
Step 2: the characteristic estimation power system impedance angle for being zero based on independent random variable covariance；
Step 3: based on the relationship between power system impedance angle and harmonic voltage and harmonic current, calculating Harmonious Waves in Power Systems
Impedance；
Step 4: being based on Harmonious Waves in Power Systems impedance, the harmonic voltage for obtaining user side and electric system side in PCC point is contributed.
2. according to claim 1 based on the harmonic contributions division methods for improving data relation analysis, which is characterized in that
The harmonic voltage current data for filtering out using time series segmentation method and meeting analysis condition is further comprised the steps of: before step 1.
3. according to claim 2 based on the harmonic contributions division methods for improving data relation analysis, which is characterized in that
Harmonic current in certain section of periodGreater than fluctuation threshold, then the measured value in this section is selected to calculate power system impedance；In certain section
Interior harmonic voltageGreater than fluctuation threshold, then the data in this section is selected to divide harmonic contributions.
4. according to claim 1 based on the harmonic contributions division methods for improving data relation analysis, which is characterized in that step
Rapid 1 specifically includes:
Relationship between PCC point h subharmonic voltage and harmonic current:
WhereinFor power system impedance, Vs hFor Harmonious Waves in Power Systems voltage,For PCC point h subharmonic voltage,For PCC
Point h subharmonic currentFormula (2), which is rewritten into following scalar form, to be had:
γ in formula (3)hFor vectorWithBetween angle, λhFor vector Vs hWithBetween angle, γh=αh+βh,
Wherein αhFor vectorWithPhase angle difference, βhFor power system impedance angle, WithPoint
It Wei not vectorWithPhase angle.
5. according to claim 1 based on the harmonic contributions division methods for improving data relation analysis, which is characterized in that electricity
Force system impedance angle calculating formula are as follows:
6. according to claim 4 based on the harmonic contributions division methods for improving data relation analysis, which is characterized in that
It is N number of in PCC point samplingWithValue, in formula (2) respectively amount average, have:
It enables
Formula (2), which subtracts formula (4), to be had:
The characteristic for being zero based on independent random variable covariance obtains:
Wherein * indicates the conjugation of vector；Formula (6) is writeable are as follows:
By formula (5) both sides simultaneously multiplied byIt averages again:
Formula (8), which is split into real and imaginary parts two parts, to be obtained:
7. according to claim 1 based on the harmonic contributions division methods for improving data relation analysis, which is characterized in that
Wherein,Respectively user side and electric system side are contributed in the harmonic voltage of PCC point,For user's resistance
It is anti-,For user's harmonic voltage,It is contributed for user's harmonic voltage in PCC point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811320306.XA CN109409762A (en) | 2018-11-07 | 2018-11-07 | Based on the harmonic contributions division methods for improving data relation analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811320306.XA CN109409762A (en) | 2018-11-07 | 2018-11-07 | Based on the harmonic contributions division methods for improving data relation analysis |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109409762A true CN109409762A (en) | 2019-03-01 |
Family
ID=65472256
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811320306.XA Pending CN109409762A (en) | 2018-11-07 | 2018-11-07 | Based on the harmonic contributions division methods for improving data relation analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109409762A (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109839538A (en) * | 2019-03-29 | 2019-06-04 | 云南电网有限责任公司电力科学研究院 | A kind of harmonic source identification method and system based on correlation analysis |
| CN109870615A (en) * | 2019-03-12 | 2019-06-11 | 西安科技大学 | Harmonic contributions assessment method suitable for industrial and mining enterprises |
| CN110244121A (en) * | 2019-06-12 | 2019-09-17 | 国电南瑞南京控制系统有限公司 | A kind of harmonic contributions evaluation method based on power quality statistical data |
| CN110456159A (en) * | 2019-08-21 | 2019-11-15 | 四川大学 | System side harmonic impedance estimation method and system based on amendment independent random vector |
| CN110927491A (en) * | 2019-12-03 | 2020-03-27 | 南方电网科学研究院有限责任公司 | Multi-harmonic source responsibility division method and device based on phase-free data |
| CN111783035A (en) * | 2020-05-21 | 2020-10-16 | 南方电网科学研究院有限责任公司 | Harmonic responsibility calculation method and device adapting to different harmonic fluctuation conditions |
| CN114137444A (en) * | 2021-11-29 | 2022-03-04 | 国网山东省电力公司日照供电公司 | Transformer running state monitoring method and system based on acoustic signals |
| CN114325021A (en) * | 2021-12-29 | 2022-04-12 | 国网上海市电力公司 | Electricity charge calculation method considering harmonic characteristics of high-speed rail |
| CN116467556A (en) * | 2023-03-28 | 2023-07-21 | 四川大学 | Harmonic impedance calculation method considering background harmonic voltage fluctuation and impedance variation |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101726663A (en) * | 2008-10-30 | 2010-06-09 | 华北电力科学研究院有限责任公司 | Method and system for monitoring user-side harmonic pollution |
| CN102478605A (en) * | 2010-11-23 | 2012-05-30 | 上海市电力公司 | Harmonic calculation method of 256 sampling points for power quality |
| CN103278685A (en) * | 2013-05-17 | 2013-09-04 | 国家电网公司 | Harmonic duty quantitative allocation method based on statistical data correlation analysis |
| CN105759117A (en) * | 2015-10-15 | 2016-07-13 | 华北电力大学 | Harmonic responsibility quantification method based on harmonic wave analysis integrated equivalent circuit |
| CN107167659A (en) * | 2017-06-29 | 2017-09-15 | 西安科技大学 | Harmonic contributions differentiating method containing background harmonic voltage |
-
2018
- 2018-11-07 CN CN201811320306.XA patent/CN109409762A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101726663A (en) * | 2008-10-30 | 2010-06-09 | 华北电力科学研究院有限责任公司 | Method and system for monitoring user-side harmonic pollution |
| CN102478605A (en) * | 2010-11-23 | 2012-05-30 | 上海市电力公司 | Harmonic calculation method of 256 sampling points for power quality |
| CN103278685A (en) * | 2013-05-17 | 2013-09-04 | 国家电网公司 | Harmonic duty quantitative allocation method based on statistical data correlation analysis |
| CN105759117A (en) * | 2015-10-15 | 2016-07-13 | 华北电力大学 | Harmonic responsibility quantification method based on harmonic wave analysis integrated equivalent circuit |
| CN107167659A (en) * | 2017-06-29 | 2017-09-15 | 西安科技大学 | Harmonic contributions differentiating method containing background harmonic voltage |
Non-Patent Citations (1)
| Title |
|---|
| 余发山等: "一种谐波源定位方法的研究 ", 《测控技术》 * |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109870615A (en) * | 2019-03-12 | 2019-06-11 | 西安科技大学 | Harmonic contributions assessment method suitable for industrial and mining enterprises |
| CN109839538A (en) * | 2019-03-29 | 2019-06-04 | 云南电网有限责任公司电力科学研究院 | A kind of harmonic source identification method and system based on correlation analysis |
| CN110244121A (en) * | 2019-06-12 | 2019-09-17 | 国电南瑞南京控制系统有限公司 | A kind of harmonic contributions evaluation method based on power quality statistical data |
| CN110244121B (en) * | 2019-06-12 | 2021-06-01 | 国电南瑞南京控制系统有限公司 | Harmonic responsibility estimation method based on power quality statistical data |
| CN110456159A (en) * | 2019-08-21 | 2019-11-15 | 四川大学 | System side harmonic impedance estimation method and system based on amendment independent random vector |
| CN110456159B (en) * | 2019-08-21 | 2020-03-17 | 四川大学 | System side harmonic impedance estimation method and system based on corrected independent random vector |
| CN110927491B (en) * | 2019-12-03 | 2021-07-20 | 南方电网科学研究院有限责任公司 | Multi-harmonic source responsibility division method and device based on phase-free data |
| CN110927491A (en) * | 2019-12-03 | 2020-03-27 | 南方电网科学研究院有限责任公司 | Multi-harmonic source responsibility division method and device based on phase-free data |
| CN111783035A (en) * | 2020-05-21 | 2020-10-16 | 南方电网科学研究院有限责任公司 | Harmonic responsibility calculation method and device adapting to different harmonic fluctuation conditions |
| CN114137444A (en) * | 2021-11-29 | 2022-03-04 | 国网山东省电力公司日照供电公司 | Transformer running state monitoring method and system based on acoustic signals |
| CN114137444B (en) * | 2021-11-29 | 2024-04-02 | 国网山东省电力公司日照供电公司 | Transformer running state monitoring method and system based on acoustic signals |
| CN114325021A (en) * | 2021-12-29 | 2022-04-12 | 国网上海市电力公司 | Electricity charge calculation method considering harmonic characteristics of high-speed rail |
| 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 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109409762A (en) | Based on the harmonic contributions division methods for improving data relation analysis | |
| Karimzadeh et al. | Method for determining utility and consumer harmonic contributions based on complex independent component analysis | |
| Kontis et al. | Measurement-based dynamic load modeling using the vector fitting technique | |
| Zygarlicki et al. | A reduced Prony's method in power-quality analysis—parameters selection | |
| Xu et al. | Study on constraints for harmonic source determination using active power direction | |
| CN109936169A (en) | Uneven and harmonic power distribution control method and device between shunt chopper | |
| CN105981252B (en) | The system and method inhibited for load harmonic | |
| Morsi et al. | Reformulating three-phase power components definitions contained in the IEEE standard 1459–2000 using discrete wavelet transform | |
| CN109030941A (en) | Tri- spectral line interpolation harmonic analysis method of Hanning involution convolution window FFT | |
| CN103401238B (en) | A kind of power load modelling approach based on Measurement-based approach | |
| Štremfelj et al. | Nonparametric estimation of power quantities in the frequency domain using Rife-Vincent windows | |
| Serfontein et al. | Continuous event-based harmonic impedance assessment using online measurements | |
| Dehghani et al. | Nonlinear state space model identification of synchronous generators | |
| Ujile et al. | A dynamic approach to identification of multiple harmonic sources in power distribution systems | |
| CN109541305A (en) | A kind of harmonic contributions partitioning model and harmonic contributions calculation method | |
| CN105548736A (en) | Harmonic wave responsibility distribution quantitative evaluation method | |
| Thomas et al. | Practical evaluation of unbalance and harmonic distortion in power conditioning | |
| Yang et al. | A novel algorithm for accurate frequency measurement using transformed consecutive points of DFT | |
| CN103424621A (en) | Artificial neural network detecting method of harmonic current | |
| CN108802491A (en) | Harmonic emission level estimation based on harmonic source independence solves new method | |
| Zolfaghari et al. | Evaluation of windowed ESPRIT virtual instrument for estimating Power Quality Indices | |
| CN110007147A (en) | A kind of harmonic contributions calculation method of single busbar multi-user | |
| Menezes et al. | Moving average-based mitigation of exponentially decaying DC components | |
| Han et al. | A new assessment for the total harmonic contributions at the point of common coupling | |
| CN104076203B (en) | A kind of intrasonic harmonic detection method considering that negative frequency affects |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190301 |