CN105866585B - A kind of harmonic wave identifing source and division of responsibiltiy method based on distortion power - Google Patents
A kind of harmonic wave identifing source and division of responsibiltiy method based on distortion power Download PDFInfo
- Publication number
- CN105866585B CN105866585B CN201610248382.9A CN201610248382A CN105866585B CN 105866585 B CN105866585 B CN 105866585B CN 201610248382 A CN201610248382 A CN 201610248382A CN 105866585 B CN105866585 B CN 105866585B
- Authority
- CN
- China
- Prior art keywords
- harmonic
- power
- user
- pcc point
- wave
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000002307 isotope ratio mass spectrometry Methods 0.000 claims description 7
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000013139 quantization Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 8
- 238000005259 measurement Methods 0.000 description 8
- 230000005611 electricity Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 239000004020 conductor Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 230000002265 prevention Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000004378 air conditioning Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000005057 refrigeration Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 101100322243 Caenorhabditis elegans deg-3 gene Proteins 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010891 electric arc Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 239000003792 electrolyte Substances 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/001—Measuring interference from external sources to, or emission from, the device under test, e.g. EMC, EMI, EMP or ESD testing
Abstract
The present invention relates to a kind of harmonic wave identifing source and division of responsibiltiy method based on distortion power, specifically includes the following steps: step S1: measuring public interface, i.e. there is the related physical quantity of user in system and the place at PCC point；Step S2: the harmonic distortion power DH at PCC point is calculated；Step S3: judge the pollution level at PCC point；Step S4: the harmonic distortion power of each user at PCC point is calculated；Step S5: the responsibility that user should undertake is determined after repeatedly calculating.The present invention is the currently the only proprietary technology using distortion power quantization harmonic contributions, avoids the estimation of harmonic impedance needed for existing method, accuracy is higher, and operability is stronger, and in conjunction with the intelligent electric meter largely used in smart grid, is more suitable for engineer application.
Description
Technical field
The present invention relates to electrical network fields, more particularly to a kind of harmonic wave identifing source based on distortion power and division of responsibiltiy side
Method.
Background technique
With the development of social economy and smart grid, a large amount of power electronics, electric arc furnaces, impact, intermittent device with
Load access system more and more, the electric current and voltage waveform distortion in electric system are increasingly severe, thereby result in huge
Harmonic pollution, jeopardized the safety and economical operation of power network and electrical equipment itself.For this purpose, the analysis of harmonic problem and comprehensive
Close the project administered and have become domestic and international extensive concern.
Harmonic source is directed to the electrical equipment of harmonic electric current and generation harmonic voltage in power grid, typically refers to some non-
Linear device.The harmonic measure value of (PCC point) is that have user's harmonic emissions coefficient by the place at public interface
As a result, punished using economic means harmonic-producing load by the harmonic contribution degree for quantifying each load, to harmonic wave by
Evil person compensates, and guidance user takes the harmonic current of approach control injected system appropriate, is effectively arranging for Harmonic Prevention
It applies.Therefore, it correctly identifies harmonic source, measures harmonic contribution, carry out harmonic contributions and share, there is most important theories value and engineering meaning
Justice.
Existing harmonic wave identifing source is mainly based upon the detection of harmonic wave active power direction PH, and this method is simple and convenient,
But it is only just set up when nonlinear load is sufficiently large or other nonlinear loads decouple, precision deficiency.Therefore, have
Scholar suggests monitoring harmonic wave reactive power direction QH.Both methods complements one another, which kind of method to depend on the electricity of load using
The balance of resistance and reactance.But the determination of impedance characteristic is difficult in many cases,.
Existing harmonic contributions methodology is broadly divided into two classes: intervening method and non-intervention method.Intervention method is by artificially producing
Raw disturbance, non-intervention method utilize the harmonic source and measurable parameter of system itself, both by estimating harmonic impedance, into
And harmonic contribution degree calculating is carried out, specify the harmonic contributions of each load.Although the computational accuracy of harmonic impedance passes through various sides
Method is constantly improving, but on Practical Project, harmonic impedance value is time-varying, and it is very tired to obtain accurate harmonic impedance
Difficulty, it is not exclusively accurate so as to cause harmonic contribution degree calculated value, and it is not easy to Project Realization.
" intervene formula " method mainly generates disturbance by artificial mode, such as to system harmonic electric current, m-Acetyl chlorophosphonazo electric current,
Or the measurement of system side harmonic impedance is carried out by a certain branch of the system of cut-offfing.By this method, it accurate can count
Harmonic impedance is calculated, and thus carries out harmonic contributions division, but this disturbance may operate normally electric system and generate not
The influence of benefit, therefore such method cannot be used widely.
Non- intervention formula method estimates harmonic impedance and harmonic wave using the harmonic source of system itself and measurable parameter etc.
Voltage, and thus carry out harmonic contributions division.Existing method mainly has: Fluctuation Method, based on measured voltage undulate quantity to electric current
The theoretical research estimation method of the symbolic feature of undulate quantity ratio, the accuracy requirement that this method measures harmonic parameters is higher, together
When also need measured value to have sufficiently large fluctuation；Linear regression method passes through measured value (plural number) in Thevenin's equivalence circuit
Real part, imaginary part equationof structure estimate harmonic impedance, this method needs system relatively stable, but due to load, electrical network parameter and
The continuous variation of system operation mode, calculation method be not still overripened.
In practical application, the direction of harmonic wave active power, size harmonic source and quantization harmonic pollution for identification are commonly used.Recognize
For if it is linear load, PH > 0；If loading nonlinear and generating harmonic energy, PH < 0.In fact, this side
Method is only just set up when nonlinear load is sufficiently large or other nonlinear loads decouple, and is had scholar's emulation and is proved, several
Under kind actual conditions, this method is not 100% accurate.
The shortcomings that using above method, is as follows: (1) being based primarily upon around harmonic current and harmonic impedance, but to obtain accurate
Harmonic impedance it is extremely difficult；(2) " intervention method " progress harmonic contributions are shared, and the disturbance artificially generated may be to electric system
Normal operation has adverse effect on；(3) comparison for calculation methods of estimation harmonic impedance is complicated in " non-intervention method ", and harmonic wave hinders
Anti- estimated accuracy is not high；(4) prior art considers distortion power, is more not based on electricity not from the angle of electromagnetic power stream
Acquisition or metering device, engineer application are difficult.
Harmonic wave identifing source, harmonic contributions methodology are all based only on always harmonic voltage, harmonic current, harmonic impedance, and
The essence of electric system is the conversion and transmission of electric flux, and electric flux therein is mainly electrical power and electrical power to the time
Integral (electricity), so theory based on power carries out harmonic contribution metrization, harmonic contributions are shared with correctness and can
Row.
Therefore, the definition based on distortion power and its each component is incorporated in widely used intelligent electric meter in smart grid,
The distortion power of each user of public interface (PCC) is calculated, and then calculates its harmonic contribution degree, carries out harmonic contributions point
Booth, is the core content and key technology of this patent.
Summary of the invention
In view of this, the harmonic wave identifing source that the object of the present invention is to provide a kind of based on distortion power and division of responsibiltiy side
Method, the substance essence based on electromagnetic power stream, using can measure obtained distortion power, in conjunction with existing electric quantity collector and intelligence
The energy achievable function of ammeter, the distortion power concept defined based on IEEE Std.1459-2010.It avoids needed for existing method
Harmonic impedance estimation, overcomes existing method to depend on the defect of accurate harmonic impedance, and from the angle of power flow, scientific explarnation is each
The responsibility of harmonic source.
The present invention is realized using following scheme: a kind of harmonic wave identifing source and division of responsibiltiy method based on distortion power, tool
Body the following steps are included:
Step S1: measuring public interface, i.e. there is the related physical quantity of user in system and the place at PCC point；
Step S2: the harmonic distortion power D at PCC point is calculatedH；
Step S3: judge the pollution level at PCC point；
Step S4: the harmonic distortion power of each user at PCC point is calculated；
Step S5: the responsibility that user should undertake is determined after repeatedly calculating.
Further, in the step S1, required physical quantity, the physics are measured using intelligent electric meter and filter
Amount includes current effective value IRMS, voltage effective value URMS, apparent energy S, all-wave active-power P, fundamental active power P1, fundamental wave
Reactive power Q1, fundamental wave apparent energy S1。
Further, in the step S2, using the distortion power D of IEEE Std 1459-2014 propositionHDefinition:
In formula:
Current distortion rate THDIAre as follows:
Voltage distortion rate THDVAre as follows:
Further, in the step S3, judge DHSize relation between η；If the D at PCC pointHLess than η, then should
Harmonic pollution is lighter at PCC point, without investigating the harmonic pollution responsibility of each user, terminates to calculate；If the D at PCC pointHIt is greater than
η, then harmonic pollution is more serious at the PCC point, quantifies the harmonic contribution degree of each user.
Further, in the step S4, if the D at PCC pointHGreater than η, further each user of user side is calculated
Harmonic distortion power uses formula to each loadCalculate harmonic power DHi；By DHiGreatly
It is small to judge whether load is harmonic source, if DHiVery little, then load i anharmonic wave source；Otherwise it is harmonic source, harmonic wave is carried out to it
Division of responsibiltiy.
Further, in the step S5, when carrying out division of responsibiltiy to each user, by the harmonic contribution of m-th of user
Degree is defined as:
Compared with prior art, beneficial effects of the present invention are as follows:
(1) it proposes to share for research harmonic contributions and harmonic problems is waited to provide new approaches: is existing humorous from distortion power
Wave division of responsibiltiy method is all based on harmonic current and harmonic impedance, and " intervention " method and " non-intervention method " are for measurement harmonic wave
Impedance, most of researchs are dedicated to constantly improving the precision of measurement harmonic impedance, but accurately to measure harmonic impedance according to
It is old extremely difficult.The invention patent avoids measurement harmonic impedance, is set out based on distortion power, goes out from the inherent mechanism of power flow
Hair, explicit physical meaning, can harmonic contribution degree objective, that accurately reflect each disturbing source, share etc. for research harmonic contributions and ask
Topic provides new approaches.
(2) foundation is provided for electricity consumption both sides contradiction, raising power supply quality for coordination: if distinguished each in harmonic pollution
The responsibility of user reasonably carries out harmonic wave and shares, the contradiction between power supply company and user as caused by harmonic problem
It is readily solved, and is conducive to improve the power quality of user and the service level of power supply company.
(3) according to harmonic contribution degree, rewards and punishments is carried out to load, are conducive to Harmonic Prevention: at public interface (PCC point)
Harmonic measure value be to have user's harmonic emissions coefficient as a result, harmonic contribution by quantifying each load by the place
Degree, punishes harmonic-producing load using economic means, compensates to the victim of harmonic wave, and guidance user takes appropriate
Approach controls the harmonic current of injected system, is the effective measures of Harmonic Prevention.The method of patent according to the present invention, can be accurate
Quantization PCC point at each user harmonic contribution degree, provide foundation for harmonic wave control, Harmonic Prevention.
(4) method proposed can be in conjunction with intelligent electric meter, strong operability, mentions for metering and the update of detection device
For new direction, improve its economic benefit: the method for proposition is not needed to estimate grid side and load side harmonic impedance, can be grasped
The property made is stronger, can be more suitable for engineer application in conjunction with the intelligent electric meter largely used in smart grid.Only need to it is existing metering and
Detection device carries out small update, can realize harmonic wave identifing source harmony wave divisions of responsibility, improves the warp of the devices such as intelligent electric meter
Ji benefit.
Detailed description of the invention
Fig. 1 is the Poynting's vector spirogram of coaxial cable of the invention.
Fig. 2 is that the power grid of the invention containing background harmonics gives linear load electrical schematic.
Fig. 3 is linear load Poynting vector instantaneous power flow diagram of the invention.
Fig. 4 (a) is the basic circuit diagram of background harmonic voltage power supply nonlinear-load of the invention.Fig. 4 (b) is the present invention
Fig. 4 (a) equivalent circuit diagram.
Fig. 5 is nonlinear-load Poynting vector instantaneous power flow diagram of the invention.
Fig. 6 is harmonic wave identifing source and division of responsibiltiy flow chart of the invention.
Fig. 7 is Multi-harmonic Sources artificial circuit figure of the invention.
Fig. 8 is Load harmonic distortion power variation diagram of the invention.
Specific embodiment
The present invention will be further described with reference to the accompanying drawings and embodiments.
The present embodiment provides a kind of harmonic wave identifing source and division of responsibiltiy method based on distortion power, as shown in fig. 6, specifically
The following steps are included:
Step S1: measuring public interface, i.e. there is the related physical quantity of user in system and the place at PCC point；
Step S2: the harmonic distortion power D at PCC point is calculatedH；
Step S3: judge the pollution level at PCC point；
Step S4: the harmonic distortion power of each user at PCC point is calculated；
Step S5: the responsibility that user should undertake is determined after repeatedly calculating.
In the present embodiment, in the step S1, required physical quantity is measured using intelligent electric meter and filter, it is described
Physical quantity includes current effective value IRMS, voltage effective value URMS, apparent energy S, all-wave active-power P, fundamental active power P1、
First Harmonic Reactive Power Q1, fundamental wave apparent energy S1。
In the present embodiment, in the step S2, using the distortion power D of IEEE Std 1459-2014 propositionHDefinition:
In formula:
Current distortion rate THDIAre as follows:
Voltage distortion rate THDVAre as follows:
In the present embodiment, in the step S3, judge DHSize relation between η；If the D at PCC pointHLess than η,
Then harmonic pollution is lighter at the PCC point, without investigating the harmonic pollution responsibility of each user, terminates to calculate；If the D at PCC pointH
Greater than η, then harmonic pollution is more serious at the PCC point, quantifies the harmonic contribution degree of each user.
In the present embodiment, in the step S4, if the D at PCC pointHGreater than η, further to each user of user side
Harmonic distortion power is calculated, formula is used to each loadCalculate harmonic power DHi；By
DHiSize judges whether load is harmonic source, if DHiVery little, then load i anharmonic wave source；Otherwise be harmonic source, to its into
Row harmonic contributions are shared.
In the present embodiment, in the step S5, when carrying out division of responsibiltiy to each user, by the harmonic wave of m-th of user
Contribution degree is defined as:
In the present embodiment, the existing electric quantity collector or novel intelligent electric meter used in the step S1 is based on intelligence
Can power grid and exist, be made of electronic integrated circuit, can be realized the intelligent terminal of remote interaction, have digital sample,
CPU analysis, data storage, the features such as security level is high, stability is strong, it has also become the representative of smart grid intelligent terminal.Root
According to " GDW 1354_2013 intelligent electric energy meter functional specification ", intelligent electric meter can be with metering current virtual value IRMS, voltage effective value
URMS, apparent energy S, all-wave active-power P, all-wave active power Q etc., need to only add filtering to intelligent electric meter signal collected
Fundamental wave amount can be calculated in device, comprising: fundamental current I1, fundamental voltage U1, apparent energy S1, active-power P1, all-wave
Active power Q1.According to " DL_T_645-2007 multifunctional electric energy meter communication protocol ", intelligent electric meter can be incited somebody to action by communication channel
Data are passed at administration of power networks, and grid management systems can be calculated according to above-mentioned data according to formula (1), (2), (3), (4)
To voltage distortion rate THDI, current distortion rate THDV, distortion power DH, harmonic contribution degree HRm, to obtain the humorous of each user
Wave contribution degree.Therefore, the harmonic wave identifing source and division of responsibiltiy method, can be with intelligent electric meter widely used in smart grid and electricity
Amount acquisition device combines, and fully demonstrates its operability, has engineering significance.
In the present embodiment, the distortion power of harmonic source and anharmonic wave source has significant difference in power grid.By linear load
Harmonic current be this subharmonic voltage linear function, the harmonic current for flowing through nonlinear-load is answering for each harmonic voltage
Miscellaneous function.When power supply network includes the harmonic wave that has powerful connections, and load is linear load, harmonic voltage and curent change very little, distort function
Rate DHChange also very little；When load is nonlinear-load, load side harmonic wave (especially harmonic current) variation acutely, necessarily causes
Distortion power DHSignificant change occurs.Therefore, harmonic source is identified, distortion power is one of key factor.
In the present embodiment, intracorporal power transfer is led with Poynting vector analysis, clearly physical interpretation can be provided.Slope
Print pavilion vector definition is the turnover rate for describing space any point energy density, characterizes unit area, flows into the unit time
Or the electromagnetic energy of outflow.Flow density (power density) when electromagnetism energy flux is propagated along conductor can be quantified with Poynting vector.
Electric current generates electric field E and magnetic field H in conductor, and power is transmitted to load side along conductor from grid side.By taking coaxial cable as an example, such as
Shown in Fig. 1, both end voltage v, electric current is i, radius a, b, c in conductor, leads the slope that cylinder electrolyte electric and magnetic fields generate
Print court of a feudal ruler vector are as follows:
Poynting vector passes through conductors cross π (b2-a2) electromagnetic power are as follows:
As it can be seen that leading intracorporal power transfer with Poynting vector analysis, clearly physical interpretation can be provided.
In the present embodiment, as shown in Fig. 2, Fig. 3, Fig. 4, Fig. 5, power grid is analyzed respectively and contains background harmonics, to linear load
With power flow when charging to nonlinear-load, the representation in components of distortion power is obtained.
Wherein, it when Fig. 2 and Fig. 3 indicates that the power grid containing background harmonics is powered to linear load, is analyzed using Poynting vector
Instantaneous power.When grid side contains background harmonics, and load is linear load RL, grid side voltage such as formula (5), principle such as Fig. 2
It is shown.
In Fig. 2, the electric current of linear RL load is flowed through are as follows:
Grid side voltage vSIt can be analyzed to fundametal compoment vS1With harmonic component vSH, and be further broken into current in phase and
Two orthogonal components:
vS1=vsp1+vsq1
At this point, Poynting vector instantaneous power stream such as Fig. 3, wherein P1And vq1i1Corresponding to fundamental wave apparent energy S1, and
vq1iH、vqHi1、vqHiHCorrespond to voltage, electric current and distortion power DI、DVAnd DHInstantaneous power.
It can be obtained by Kirchhoff's second law:
Harmonic current iHCaused instantaneous power are as follows:
RiminIt is defined as the intrinsic power of second order (intrinsic power ofthe second-order), mean value zero,
Always exist with harmonic wave active power, power grid power loss will not be caused.
The distortion power D that RL absorbsHTransient component are as follows:
When Fig. 4 and Fig. 5 is that the power grid containing background harmonics is powered to nonlinear-load, instantaneous function is analyzed using Poynting vector
Rate.When grid side contains background harmonics, and load is nonlinear-load NL, basic principle and equivalent circuit such as Fig. 4 (a) and Fig. 4
(b), the Poynting vector instantaneous power stream of NL is as shown in Figure 5.
RSAnd LSFor the equivalent parameters of power grid, grid side voltage such as formula (5).
At this point, flowing through the electric current of NL are as follows:
Voltage on NL are as follows:
Similar to formula (16), vS、vNLIt is decomposable as follows:
v1=vp1+vq1
It is obtained by Kirchhoff's second law:
Instantaneous power caused by harmonic current are as follows:
(vSH-vH)i'H=vSHi'H-vHi'H (17)
Wherein, the Section 2 on the right side of formula (24) is instantaneous power caused by NL harmonic voltage and electric current, i.e. harmonic wave view is in function
Rate:
In formula, first item vsHWith i'HCause harmonic wave apparent energy SHTransient component；Section 2 be harmonic wave active power and
A part of the intrinsic power of mn second order；Last is LSWith vSHAnd vHBetween oscillation of power.
The distortion power D of NLHTransient component are as follows:
In formula, first item portrays the interaction of background harmonic voltage and harmonic current；Section 2 is power grid equivalent parameters
The harmonic distortion power of absorption.
The harmonic current generated due to nonlinear-load much larger than harmonic current caused by linear load, comparison expression (12) and
Formula (19), due to nonlinear-load harmonic current i'hGreater than linear load harmonic current ih, therefore, nonlinear-load distortion power
pDHNLGreater than linear load distortion pDHRL, this method proposition measures the technology of the harmonic contribution of each harmonic source with distortion power.
When each disturbing source Load harmonic contribution degree quantifies in actual electric network, without loss of generality, harmonic voltage and
The measurement of harmonic current meets IEEE Std 1159-2009 standard, (generally with a production cycle of each disturbing source load
24 hours), can specifically be determined according to actual load property or objectives of examination demand) it is examination period, harmonic voltage and harmonic wave electricity
Stream takes 95% probability value by existing standard, then calculates the harmonic contribution degree of harmonic source.
In the present embodiment, detailed process and the side of harmonic wave identifing source harmony wave contribution metrization are illustrated using instance analysis
The accuracy of method is proved in terms of theoretical model, physical model two respectively.
Theoretical model:
Simulation analysis is carried out with an electric network model with Multi-harmonic Sources, the nonlinear-load of common access power grid is main
Have: electronic control unit, Thyristor Controlled etc., these loads are usually typical harmonic source.For without loss of generality, to figure
Power grid shown in 7 is emulated, wherein L1 is purely load；L2 is compensating reactive power and the passive filter for filtering out harmonic wave in power grid
Wave device；L3 is that thyristor controls reactance (TCR) type static passive compensation device (SVC)；L4 is bridge rectifier load；L5 is PWM
Rectify load；L6 is that thyristor controls DC speed regulator load.
The background harmonic voltage of grid side such as table 1, fundamental wave voltage rating 220V, fundamental frequency 50Hz.
1 grid side background harmonic voltage of table
| Overtone order | Containing ratio/% | Phase/deg |
| 3 | 3.43 | 8.87 |
| 5 | 3.44 | 16.21 |
| 7 | 1.16 | 25.11 |
| 9 | 1.68 | 26.45 |
| 11 | 0.58 | 31.51 |
(1) each physical quantity of simulation calculation at PCC point, comprising: current effective value IRMS, voltage effective value URMS, view exist
Power S, fundamental active power P1, all-wave active-power P, reactive power Q, the measured value such as table 2 in this example；
(2) in PCC point survey calculation harmonic wave harmonic distortion power DH, according to formula 1 and table 3, the D that PCC is pointed out is calculatedH
=205.85, analysis show that harmonic pollution herein is more serious, needs to carry out it harmonic contribution degree for quantifying each load
Harmonic contributions are investigated；
(3) further to each user's survey calculation harmonic distortion power of user side, each load is calculated according to formula 1 humorous
Wave power DHi, calculated result such as table 3,
(4) division of responsibiltiy is carried out to each user, the harmonic contribution degree of each user, calculated result such as table is calculated by formula 4
3。
2 simulation result of table
| PCC | L1 | L2 | L3 | L4 | L5 | L6 | |
| IRMS/A | 38.05 | 7.34 | 6.74 | 10.31 | 11.77 | 12.07 | 20.41 |
| VRMS/V | 220.21 | 220.21 | 220.21 | 220.21 | 220.21 | 220.21 | 220.21 |
| S/VA | 8380.25 | 1615.56 | 1483.19 | 2270.09 | 2591.61 | 2658.27 | 4494.39 |
| P1/W | 5955.16 | 1610.99 | -526.88 | 1045.59 | 1713.62 | -1919.78 | 3992.43 |
| Q1/var | 5378.02 | 0 | 1384.14 | 1762.68 | -498.68 | 313.67 | -1922.01 |
| S1/VA | 8024.15 | 1610.99 | 1481.03 | 2049.47 | 1784.71 | 1945.24 | 4430.99 |
| P/W | 5894.92 | 1615.45 | -526.67 | 1040.37 | 1735.3 | -1866.21 | 3994.68 |
Each branch harmonic power of table 3 and harmonic contribution degree
| Load | Harmonic distortion power DH | Harmonic contribution degree HR/var |
| PCC | 205.85 | 205.85 |
| L1 | 1.01 | Anharmonic wave source |
| L2 | 0.69 | Anharmonic wave source |
| L3 | 51.37 | 27.86 |
| L4 | 97.55 | 52.91 |
| L5 | 190.97 | 103.57 |
| L6 | 37.98 | 20.60 |
Seen from table 3, L2 is passive filter, can inhibit harmonic wave, in the case where apparent energy is not much different, L2's
Yin Shi distortion power DHIt is significantly less than L1, can be acted on the harmonics restraint of reasonable dismissal filter, correctly reflect actual physics
Mechanism, and can effectively identify L3, L4, L5, the harmonic contribution degree of the loads such as L6.These results and theory analysis are completely the same, card
The correctness for the method for being illustrated.
Physical model:
In laboratory environments, make nonlinear-load using water heater, computer, refrigeration air-conditioner, cooling/heating air conditioner, carry out physics
Model verifying.The measurement of harmonic voltage and harmonic current follows IEEE Std 1159-2009, and measurement result result is through MATLAB
The harmonic contribution degree of measurand is calculated after processing, the distortion power of each measurand changes with time such as Fig. 8, measurement period
Interior distortion power average value such as table 4.
4 measured data result of table
| Water heater | Computer | Air-conditioning (refrigeration) | Air-conditioning (heating) | |
| IRMS/A | 6.26 | 0.51 | 7.05 | 13.36 |
| VRMS/V | 222.27 | 234.01 | 238.84 | 239.26 |
| S/VA | 1393.05 | 150.59 | 1711.49 | 3213.01 |
| P1/W | 1391.92 | 48.96 | 1683.45 | 2343.41 |
| Q1/var | -15.09 | -108.68 | 34.21 | 2175.39 |
| S1/VA | 1392.00 | 119.21 | 1683.79 | 3197.48 |
| PH/W | -0.92 | 0.78 | -1.06 | -0.76 |
| P/W | 1390.99 | 49.74 | 1682.39 | 2342.65 |
| THDV/ % | 2.65 | 2.79 | 2.72 | 2.67 |
| THDI/ % | 2.66 | 77.11 | 18.00 | 9.50 |
| DI/var | 38.36 | 91.92 | 303.09 | 303.71 |
| DV/var | 38.31 | 3.33 | 45.78 | 85.34 |
| DH/var | 0.51 | 2.45 | 8.17 | 8.07 |
By Fig. 8 and table 4 as it can be seen that measured result is consistent with actual physics phenomenon.
The foregoing is merely presently preferred embodiments of the present invention, all equivalent changes done according to scope of the present invention patent with
Modification, is all covered by the present invention.
Claims (5)
1. a kind of harmonic wave identifing source and division of responsibiltiy method based on distortion power, it is characterised in that: specifically includes the following steps:
Step S1: measuring public interface, i.e. there is the related physical quantity of user in system and the place at PCC point；
Step S2: the harmonic distortion power D at PCC point is calculatedH；
Step S3: judge the pollution level at PCC point；
Step S4: the harmonic distortion power of each user at PCC point is calculated；
Step S5: the responsibility that user should undertake is determined after repeatedly calculating；
In the step S5, when carrying out division of responsibiltiy to each user, by the harmonic contribution degree of m-th of user is defined as:
2. a kind of harmonic wave identifing source and division of responsibiltiy method, feature based on distortion power according to claim 1 exists
In: in the step S1, required physical quantity is measured using intelligent electric meter and filter, the physical quantity includes that electric current is effective
Value IRMS, voltage effective value URMS, apparent energy S, all-wave active-power P, fundamental active power P1, First Harmonic Reactive Power Q1, fundamental wave
Apparent energy S1。
3. a kind of harmonic wave identifing source and division of responsibiltiy method, feature based on distortion power according to claim 1 exists
In: in the step S2, using the distortion power D of IEEE Std 1459-2014 propositionHDefinition:
In formula:
Current distortion rate THDIAre as follows:
Voltage distortion rate THDVAre as follows:
4. a kind of harmonic wave identifing source and division of responsibiltiy method, feature based on distortion power according to claim 1 exists
In: in the step S3, judge DHSize relation between η；If the D at PCC pointHLess than η, then harmonic pollution at the PCC point
It is relatively light, without investigating the harmonic pollution responsibility of each user, terminate to calculate；If the D at PCC pointHIt is greater than η, then humorous at the PCC point
Wave pollution is more serious, quantifies the harmonic contribution degree of each user.
5. a kind of harmonic wave identifing source and division of responsibiltiy method, feature based on distortion power according to claim 1 exists
In: in the step S4, if the D at PCC pointHGreater than η, harmonic distortion power further is calculated to each user of user side,
Formula is used to each loadCalculate harmonic power DHi；By DHiSize judges that load is
No is harmonic source, if DHiVery little, then load i anharmonic wave source；Otherwise it is harmonic source, harmonic contributions is carried out to it and are shared.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610248382.9A CN105866585B (en) | 2016-04-20 | 2016-04-20 | A kind of harmonic wave identifing source and division of responsibiltiy method based on distortion power |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610248382.9A CN105866585B (en) | 2016-04-20 | 2016-04-20 | A kind of harmonic wave identifing source and division of responsibiltiy method based on distortion power |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105866585A CN105866585A (en) | 2016-08-17 |
| CN105866585B true CN105866585B (en) | 2018-12-11 |
Family
ID=56633571
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610248382.9A Active CN105866585B (en) | 2016-04-20 | 2016-04-20 | A kind of harmonic wave identifing source and division of responsibiltiy method based on distortion power |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105866585B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2752765C1 (en) * | 2020-12-10 | 2021-08-03 | федеральное государственное бюджетное образовательное учреждение высшего образования «Санкт-Петербургский горный университет» | Method for estimating contribution of nonlinear consumers to voltage distortion at point of common connection |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106771552B (en) * | 2016-12-02 | 2019-08-27 | 国网四川省电力公司电力科学研究院 | A kind of distortion power metering method |
| CN110320407A (en) * | 2019-07-31 | 2019-10-11 | 国家电网有限公司 | A kind of spot measurement power grid main harmonic source localization method based on apparent energy |
| CN115598435A (en) * | 2022-09-23 | 2023-01-13 | 海南电网有限责任公司（Cn） | Harmonic responsibility apportionment quantitative evaluation method |
Citations (3)
| 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 |
| CN103336265A (en) * | 2013-06-26 | 2013-10-02 | 武汉大学 | Metering error quantitative analysis method for electric energy meter under harmonic wave condition |
| CN103838959A (en) * | 2013-12-18 | 2014-06-04 | 国网上海市电力公司 | Method for applying partial least squares regression to power distribution network harmonic source positioning and detecting |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001004686A (en) * | 1999-06-22 | 2001-01-12 | Chubu Electric Power Co Inc | Method for measuring harmonic characteristics of electric power system |
| CN203071576U (en) * | 2012-11-30 | 2013-07-17 | 中国电力科学研究院 | Multi-voltage level harmonic disturbance source generation system |
| CN104833886B (en) * | 2015-05-20 | 2017-09-29 | 国网上海市电力公司 | A kind of direct current drop point harmonic wave evaluation method conducted from hiigh pressure stage to low-pressure stage |
-
2016
- 2016-04-20 CN CN201610248382.9A patent/CN105866585B/en active Active
Patent Citations (3)
| 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 |
| CN103336265A (en) * | 2013-06-26 | 2013-10-02 | 武汉大学 | Metering error quantitative analysis method for electric energy meter under harmonic wave condition |
| CN103838959A (en) * | 2013-12-18 | 2014-06-04 | 国网上海市电力公司 | Method for applying partial least squares regression to power distribution network harmonic source positioning and detecting |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2752765C1 (en) * | 2020-12-10 | 2021-08-03 | федеральное государственное бюджетное образовательное учреждение высшего образования «Санкт-Петербургский горный университет» | Method for estimating contribution of nonlinear consumers to voltage distortion at point of common connection |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105866585A (en) | 2016-08-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| He et al. | Incorporating non-intrusive load monitoring into building level demand response | |
| CN105866585B (en) | A kind of harmonic wave identifing source and division of responsibiltiy method based on distortion power | |
| Pfajfar et al. | Harmonic contributions evaluation with the harmonic current vector method | |
| CN108490288B (en) | A kind of stealing detection method and system | |
| Dell'Aquila et al. | New power-quality assessment criteria for supply systems under unbalanced and nonsinusoidal conditions | |
| CN107121609A (en) | A kind of electric network impedance on-line identification method and device that injection is disturbed based on PRBS | |
| CN101441231B (en) | Harmonic electric energy metering error analytical apparatus | |
| CN106093570B (en) | It is a kind of to judge system side for the harmonic wave source tracing method of main harmonic source | |
| Li et al. | Research of harmonic distortion power for harmonic source detection | |
| CN109061300A (en) | A kind of characteristic harmonics source localization method of PCC point | |
| CN105044649B (en) | A kind of current transformer error horizontal dynamic detection method and system | |
| Lee et al. | New power quality index in a distribution power system by using RMP model | |
| CN109829560A (en) | A kind of power distribution network renewable energy power generation cluster access planing method | |
| CN110048409A (en) | A kind of electric load structural recognition method based on load starting transient characterisitics | |
| Lin et al. | A method for harmonic sources detection based on harmonic distortion power rate | |
| CN111009893A (en) | Household user short-term load prediction method based on load decomposition technology | |
| CN109146265A (en) | Three-phase imbalance source source tracing method based on Dynamic Programming Time Series Similarity | |
| CN104375026A (en) | Method for identifying types of harmonic sources on basis of characteristic quantity analysis | |
| CN110571806B (en) | Feature extraction and identification method for load category of power distribution network | |
| CN105486945A (en) | Determination method for line loss abnormity of 10kV line | |
| CN207472983U (en) | A kind of electric network impedance on-line identification device based on PRBS disturbance injections | |
| CN109494743A (en) | The estimation method of electric system mains by harmonics state based on matrix reconstruction | |
| CN109038538A (en) | A method of static voltage stability and angle stability are evaluated with broad sense resistance nargin | |
| CN104483570B (en) | A kind of Distribution Network Harmonics improvement method for loss reduction effect reality testing | |
| Jo et al. | Kalman-filter-based multilevel analysis to estimate electric load composition |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |