CN109193787A - A kind of harmonic wave path identification method of the power distribution network containing new energy - Google Patents
A kind of harmonic wave path identification method of the power distribution network containing new energy Download PDFInfo
- Publication number
- CN109193787A CN109193787A CN201811263740.9A CN201811263740A CN109193787A CN 109193787 A CN109193787 A CN 109193787A CN 201811263740 A CN201811263740 A CN 201811263740A CN 109193787 A CN109193787 A CN 109193787A
- Authority
- CN
- China
- Prior art keywords
- branch
- harmonic
- formula
- current
- entropy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 claims description 3
- 230000015556 catabolic process Effects 0.000 claims 1
- 239000004744 fabric Substances 0.000 claims 1
- 230000005611 electricity Effects 0.000 description 4
- 238000010248 power generation Methods 0.000 description 2
- 238000009827 uniform distribution Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
Classifications
-
- H02J3/382—
-
- 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
Abstract
The present invention is a kind of harmonic wave path identification method of power distribution network containing new energy, its main feature is that relationship including building harmonic current source and branch harmonic current, being solved to the probability density function of line harmonic electric current using the probability density function of function variable, the definition of comentropy being utilized to solve contacting between line harmonic current information entropy and harmonic source comentropy.Since the present invention is to analyze from network perspective branch harmonic current, abundant information is obtained to determine the distribution approach of harmonic current, its clear physics conception, with scientific and reasonable, the advantages that calculating is simple, identification path is accurate, be easily achieved and engineering value is high, can be improved harmonic wave Path Recognition efficiency.
Description
Technical field
The present invention relates to system for distribution network of power the safe and economic operation fields, are a kind of harmonic wave roads of power distribution network containing new energy
Diameter recognition methods, the harmonic wave Path Recognition technology applied to the power distribution network containing new energy.
Background technique
Wind-powered electricity generation, photovoltaic are at home and abroad widely applied as clean energy resource.Wind-powered electricity generation, photovoltaic installed capacity and power generation
Amount accounts for national total installation of generating capacity and the accounting of total power generation steps up.Generation of electricity by new energy is bringing the same of a large amount of clean energy resourcies
When, also a large amount of harmonic waves are introduced to power grid because of itself inverter structure.Power grid is humorous after how assessing and analyzing new energy access
The distribution problem of wave electric current is of great significance for new energy power distribution network the safe and economic operation.
Currently, it is main there are two in terms of about the harmonic study of new energy power distribution network, first is that new energy harmonic wave Influencing Mechanism
Research, such as by being modeled to new energy inverter, the relationship of the influence factors such as analysis new energy harmonic wave and irradiation, but such
Research haves the defects that certain, cannot analyze the influence caused by new energy power distribution network of each factor comprehensively；Second is that new energy is humorous
Wave inhibits control research to adjust the output impedance of inverter such as by the algorithm of optimization inverter, to avoid inverter and electricity
Resonance occurs for net, but this method control effect is limited and realizes considerably complicated.
Summary of the invention
The object of the present invention is to the shortcomings that overcoming the prior art, provide a kind of scientific and reasonable, calculates simple, identification path
Accurately, it is easily achieved and engineering value is high, can be improved the harmonic wave path of the power distribution network containing new energy of harmonic wave Path Recognition efficiency
Recognition methods.The purpose of the present invention is what is realized by following technical scheme: a kind of harmonic wave path knowledge of the power distribution network containing new energy
Other method, characterized in that it the following steps are included:
1) functional relation of harmonic current source and branch harmonic current is constructed；
1. the nodal voltage equation formula containing voltage source is (1) formula:
U=ZI (1)
Wherein: Z is nodal impedance matrix, and I is equivalent nodal current source vector, and U is node voltage column vector；
2. branch current IijRelationship between node voltage is (2) formula:
Iij=(Ui-Uj)/Zij (2)
Wherein: ZijFor branch impedance, Ui、UjThe respectively voltage of the side same branch i, the side j, IijFor branch current；
The linear relation that harmonic current source and branch harmonic current are obtained in conjunction with (1) formula and (2) formula is (3) formula:
Ibc=ZbcZI (3)
Wherein: ZbcFor branch impedance matrix, IbcFor branch harmonic current vector；
2) it is solved using probability density function of the probability density function of function variable to line harmonic electric current；
3. branch harmonic current obeys probability density function identical with current source, expression formula is (4) formula:
G(Ibc)=K F ((ZbcZ)-1I) (4)
Wherein: F ((ZbcZ)-1It I) is (ZbcZ)-1The probability density function of I, K are (ZbcZ)-1Determinant, G (Ibc) it is branch
The probability density function of road harmonic current；
3) connection between line harmonic current information entropy and harmonic source comentropy is solved using the definition of comentropy；
4. the expression formula of comentropy is (5) formula:
Wherein: x is stochastic variable, xiFor i-th of stochastic variable, m xiValue number, H (x) are information entropy function, P
(xi) represent x=xiProbability, loga(P(xi)) it is P (xi) logarithmic function, a be logarithm bottom；
5. the comentropy expression formula of harmonic current source is (6) formula:
Wherein: P (x1, x2。。。, xn) it is stochastic variable xiDensity function；i1For first variable x1Cyclic variable, in
For variable xnCyclic variable.N is the number of x variable, k1To dependent variable x1It is distributed value number, knFor xnValue in distribution
Number；
Comentropy is the Measure Indexes of information distribution, and for Harmonious Waves in Power Systems, when entropy maximum, information is obeyed uniform
Distribution, corresponding electric system is constantly in stable state, and when being transitioned into transient state from stable state, entropy is then from large to small；Current source entropy is greater than
The entropy of branch current, entropy flow in load process from current source, are gradually reduced, i.e., energy is in the flowing of power grid, and entropy is continuous
Loss.
A kind of harmonic wave path identification method of power distribution network containing new energy of the invention with it is existing humorous from inverter angle analysis
Wave method is compared, since the present invention is to analyze from network perspective branch harmonic current, to determine the distribution of harmonic current
Approach and obtain abundant information, clear physics conception has scientific and reasonable, calculates that simple, identification path is accurate, is easy to real
The advantages that existing and engineering value is high, can be improved harmonic wave Path Recognition efficiency.
Specific embodiment
A kind of harmonic wave path identification method of power distribution network containing new energy of the invention is described in detail below.
A kind of harmonic wave path identification method of power distribution network containing new energy of the invention, comprising the following steps:
1) functional relation of harmonic current source and branch harmonic current is constructed；
1. the nodal voltage equation formula containing voltage source is (1) formula:
U=ZI (1)
Wherein: Z is nodal impedance matrix, and I is equivalent nodal current source vector, and U is node voltage column vector；
2. branch current IijRelationship between node voltage is (2) formula:
Iij=(Ui-Uj)/Zij (2)
Wherein: ZijFor branch impedance, Ui、UjThe respectively voltage of the side same branch i, the side j, IijFor branch current；
The linear relation that harmonic current source and branch harmonic current are obtained in conjunction with (1) formula and (2) formula is (3) formula:
Ibc=ZbcZI (3)
Wherein: ZbcFor branch impedance matrix, IbcFor branch harmonic current vector；
2) it is solved using probability density function of the probability density function of function variable to line harmonic electric current；
3. branch harmonic current obeys probability density function identical with current source, expression formula is (4) formula:
G(Ibc)=K F ((ZbcZ)-1I) (4)
Wherein: F ((ZbcZ)-1It I) is (ZbcZ)-1The probability density function of I, K are (ZbcZ)-1Determinant, G (Ibc) it is branch
The probability density function of road harmonic current；
3) connection between line harmonic current information entropy and harmonic source comentropy is solved using the definition of comentropy；
4. the expression formula of comentropy is (5) formula:
Wherein: x is stochastic variable, xiFor i-th of stochastic variable, m xiValue number, H (x) are information entropy function, P
(xi) represent x=xiProbability, loga(P(xi)) it is P (xi) logarithmic function, a be logarithm bottom；
5. the comentropy expression formula of harmonic current source is (6) formula:
Wherein: P (x1, x2。。。, xn) it is stochastic variable xiDensity function；i1For first variable x1Cyclic variable, in
For variable xnCyclic variable.N is the number of x variable, k1To dependent variable x1It is distributed value number, knFor xnValue in distribution
Number；
Comentropy is the Measure Indexes of information distribution, and for Harmonious Waves in Power Systems, when entropy maximum, information is obeyed uniform
Distribution, corresponding electric system is constantly in stable state, and when being transitioned into transient state from stable state, entropy is then from large to small；Current source entropy is greater than
The entropy of branch current, entropy flow in load process from current source, are gradually reduced, i.e., energy is in the flowing of power grid, and entropy is continuous
Loss.
Claims (1)
1. a kind of harmonic wave path identification method of power distribution network containing new energy, characterized in that it the following steps are included:
1) functional relation of harmonic current source and branch harmonic current is constructed；
1. the nodal voltage equation formula containing voltage source is (1) formula:
U=ZI (1)
Wherein: Z is nodal impedance matrix, and I is equivalent nodal current source vector, and U is node voltage column vector；
2. branch current IijRelationship between node voltage is (2) formula:
Iij=(Ui-Uj)/Zij (2)
Wherein: ZijFor branch impedance, Ui、UjThe respectively voltage of the side same branch i, the side j, IijFor branch current；
The linear relation that harmonic current source and branch harmonic current are obtained in conjunction with (1) formula and (2) formula is (3) formula:
Ibc=ZbcZI (3)
Wherein: ZbcFor branch impedance matrix, IbcFor branch harmonic current vector；
2) it is solved using probability density function of the probability density function of function variable to line harmonic electric current；
3. branch harmonic current obeys probability density function identical with current source, expression formula is (4) formula:
G(Ibc)=K F ((ZbcZ)-1I) (4)
Wherein: F ((ZbcZ)-1It I) is (ZbcZ)-1The probability density function of I, K are (ZbcZ)-1Determinant, G (Ibc) it is that branch is humorous
The probability density function of wave electric current；
3) connection between line harmonic current information entropy and harmonic source comentropy is solved using the definition of comentropy；
4. the expression formula of comentropy is (5) formula:
Wherein: x is stochastic variable, xiFor i-th of stochastic variable, m xiValue number, H (x) are information entropy function, P (xi) generation
Table x=xiProbability, loga(P(xi)) it is P (xi) logarithmic function, a be logarithm bottom；
5. the comentropy expression formula of harmonic current source is (6) formula:
Wherein: P (x1, x2., xn) it is stochastic variable xiDensity function；i1For first variable x1Cyclic variable, inTo become
Measure xnCyclic variable.N is the number of x variable, k1To dependent variable x1It is distributed value number, knFor xnValue number in distribution；
Comentropy is the Measure Indexes of information distribution, and for Harmonious Waves in Power Systems, when entropy maximum, information is obeyed uniformly point
Cloth, corresponding electric system is constantly in stable state, and when being transitioned into transient state from stable state, entropy is then from large to small；Current source entropy is greater than branch
The entropy of road electric current, entropy flow in load process from current source, are gradually reduced, i.e., energy is in the flowing of power grid, and entropy is in not breakdown
It loses.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811263740.9A CN109193787B (en) | 2018-10-28 | 2018-10-28 | Harmonic path identification method for power distribution network containing new energy |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811263740.9A CN109193787B (en) | 2018-10-28 | 2018-10-28 | Harmonic path identification method for power distribution network containing new energy |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109193787A true CN109193787A (en) | 2019-01-11 |
| CN109193787B CN109193787B (en) | 2021-07-06 |
Family
ID=64944122
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811263740.9A Active CN109193787B (en) | 2018-10-28 | 2018-10-28 | Harmonic path identification method for power distribution network containing new energy |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109193787B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5836003A (en) * | 1993-08-26 | 1998-11-10 | Visnet Ltd. | Methods and means for image and voice compression |
| WO1999062052A2 (en) * | 1998-05-27 | 1999-12-02 | Microsoft Corporation | System and method for entropy encoding quantized transform coefficients of a signal |
| JP2000122519A (en) * | 1998-10-19 | 2000-04-28 | Toshiba Corp | Simulator for training on electric power system and recording medium on which execution program therefor is recorded |
| CN104979829A (en) * | 2015-06-16 | 2015-10-14 | 河海大学 | Harmonic current estimation method based on improved FastICA algorithm |
| CN106203382A (en) * | 2016-07-20 | 2016-12-07 | 河海大学 | A kind of excitation surge current based on kernel function extreme learning machine and fault current recognition methods |
| CN107039976A (en) * | 2017-06-08 | 2017-08-11 | 山东鲁能智能技术有限公司 | It is a kind of to lead the intelligent anti-islanding system of the power distribution network passively combined and configuration operation method |
-
2018
- 2018-10-28 CN CN201811263740.9A patent/CN109193787B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5836003A (en) * | 1993-08-26 | 1998-11-10 | Visnet Ltd. | Methods and means for image and voice compression |
| WO1999062052A2 (en) * | 1998-05-27 | 1999-12-02 | Microsoft Corporation | System and method for entropy encoding quantized transform coefficients of a signal |
| JP2000122519A (en) * | 1998-10-19 | 2000-04-28 | Toshiba Corp | Simulator for training on electric power system and recording medium on which execution program therefor is recorded |
| CN104979829A (en) * | 2015-06-16 | 2015-10-14 | 河海大学 | Harmonic current estimation method based on improved FastICA algorithm |
| CN106203382A (en) * | 2016-07-20 | 2016-12-07 | 河海大学 | A kind of excitation surge current based on kernel function extreme learning machine and fault current recognition methods |
| CN107039976A (en) * | 2017-06-08 | 2017-08-11 | 山东鲁能智能技术有限公司 | It is a kind of to lead the intelligent anti-islanding system of the power distribution network passively combined and configuration operation method |
Non-Patent Citations (7)
| Title |
|---|
| MASOUD FARHOODNEA等: "A New Method for Determining Multiple Harmonic Source Locations in a Power Distribution System", 《2010 IEEE INTERNATIONAL CONFERENCE ON POWER AND ENERGY (PECON2010)》 * |
| MURAT SILSÜPÜR等: "Application of information theory in smart grid measurements", 《2016 4TH INTERNATIONAL ISTANBUL SMART GRID CONGRESS AND FAIR (ICSG)》 * |
| RAN XIXI等: "Power quality harmonic detection based on Fast-ICA", 《 2011 IEEE POWER ENGINEERING AND AUTOMATION CONFERENCE》 * |
| 刘鑫等: "基于最大熵原理的谐波求和问题", 《电力系统保护与控制》 * |
| 吴健等: "基于概率预测与谐波潮流的配电网谐波源识别方法", 《电力系统保护与控制》 * |
| 王允建: "基于信息熵的多谐波源定位研究", 《电力系统保护与控制》 * |
| 苗长新等: "基于独立分量分析和互信息的多谐波源定位", 《电测与仪表》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109193787B (en) | 2021-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106451556B (en) | Method and device for determining capability of power distribution network feeder line for accepting distributed power supply | |
| Chen et al. | An interval optimization based day-ahead scheduling scheme for renewable energy management in smart distribution systems | |
| CN104009484A (en) | Power grid loss reduction assessment method based on load flow calculation | |
| Wang et al. | Data-driven probabilistic small signal stability analysis for grid-connected PV systems | |
| Mohan et al. | A weighted multi-objective index based optimal distributed generation planning in distribution system | |
| CN106355336A (en) | Power generation efficiency evaluation method of photovoltaic power station | |
| CN105046588A (en) | Improved DC (Direct Current) dynamic optimal power flow calculating method based on network loss iteration | |
| Kanth et al. | Siting & sizing of DG for power loss & THD reduction, voltage improvement using PSO & sensitivity analysis | |
| Li et al. | Probabilistic optimal power flow calculation method based on adaptive diffusion kernel density estimation | |
| CN105305463B (en) | The idle work optimization method based on probabilistic loadflow of meter and photovoltaic generation and harmonic pollution | |
| Wang et al. | A novel chance constrained joint optimization method under uncertainties in distribution networks | |
| CN104158217A (en) | Output power fluctuation characteristic description method for clustered wind-solar combined power generation system | |
| CN110336322B (en) | Photovoltaic power generation access capacity determination method based on daily minimum load confidence interval | |
| CN104992246A (en) | Improved-least-square-method-based prediction method of load electric quantity for transformer substation | |
| CN109193787A (en) | A kind of harmonic wave path identification method of the power distribution network containing new energy | |
| Eissa et al. | A novel approach for optimum allocation of Flexible AC Transmission Systems using Harmony Search technique | |
| Yasin et al. | Harmonic Distortion Prediction Model of a Grid-Connected Photovoltaic Using Grey Wolf Optimizer-Least Square Support Vector Machine | |
| Zhang et al. | A Hybrid Energy Storage System Strategy for Smoothing Photovoltaic Power Fluctuation Based on Improved HHO-VMD | |
| Gao et al. | High-impedance arc fault modeling for distribution networks based on dynamic geometry dimension | |
| Zhang et al. | The Siting and Sizing Problem of Distributed Generation Based on a Novel Fractional Particle Swarm Optimization Algorithm | |
| Babu et al. | Application of self adaptive differential evolution algorithm for optimal placement and sizing of renewable DG sources in distribution network including different load models | |
| Li et al. | Improved PV model and its application in reliability evaluation of microgrid | |
| Yi et al. | Simulation of HVDC transmission system failure rate bathtub curve based on Weibull distribution | |
| Nahilia et al. | Voltage stability analysis and optimal distributed generators placement study on Algerian power network | |
| Liang et al. | Simulation analysis on harmonics pollution of large-scale wind system accessed to grid |
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 |