CN109245107A - A kind of uneven distribution power system load flow calculation considering distributed generation resource access - Google Patents
A kind of uneven distribution power system load flow calculation considering distributed generation resource access Download PDFInfo
- Publication number
- CN109245107A CN109245107A CN201811261922.2A CN201811261922A CN109245107A CN 109245107 A CN109245107 A CN 109245107A CN 201811261922 A CN201811261922 A CN 201811261922A CN 109245107 A CN109245107 A CN 109245107A
- Authority
- CN
- China
- Prior art keywords
- phase
- node
- equation
- voltage
- indicate
- 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
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of uneven distribution power system load flow calculations of consideration distributed generation resource access, it includes the following steps: (1) to obtain a phase Injection Current in three-phase systemStep 2 obtains b phase Injection Current in three-phase systemStep 3 obtains c phase Injection Current in three-phase systemStep 4, the nodal voltage equation according to circuitry obtain node current equation, calculate node Injection Current；Step 5 substitutes into node current equation in power equation；Step 6, the equation for obtaining step 5 are unfolded, and finally obtain the Load flow calculation of the active and reactive power equation composition of node injecting power equation；The voltage and trend distribution situation that can aid in determining three-phase imbalance power grid provide foundation to the steady-state analysis of system.
Description
Technical field
The invention belongs to distributed generation resource access technologies more particularly to a kind of imbalance for considering distributed generation resource access to match
Electric network swim calculates.
Background technique
Nowadays distributed generation resource, which is just continuously increased and accesses power distribution network, participates in operation, becomes greatly what WeiLai Technology continued to develop
Under gesture environment, the capacity of the distributed generation resource in network system also can constantly increase year by year, so research is containing distributed generation resource
The trend distribution of power distribution network has the important urgent meaning of reality.
The bilingual coupling three-phase rapid tidal current algorithm of power distribution network in the prior art based on alternate decoupling compensation model, the algorithm
Dimension-reduction treatment is carried out by the alternate symplectic Jacobian matrix of solution, the complexity of calculating is simplified, is primarily adapted for use in route in network
Resistance and the lesser system of reactance ratio.
Consider that distributed generation resource accesses network bring trend fluctuation problem, is based on two-point estimate side in the prior art
Method carries out the probabilistic load flow of power distribution network, but this method is applicable in Complicated Distribution Network on the basis of forward-backward sweep method
Property is not strong.
It is directed to the unitized model that middle and low voltage network establishes each element in network in the prior art, to propose
Power flow algorithm and algorithm in the three-phase system of fining.
The design feature for utilizing distribution network in the prior art, establishes the relational matrix of node voltage and Injection Current,
Realize the direct calculating of distribution power flow, process is simple, but the scope of application of direct algorithm is more limited to.
Summary of the invention
The technical problem to be solved by the present invention is providing a kind of uneven distribution power flow of consideration distributed generation resource access
It calculates, determines the voltage and trend distribution situation of three-phase imbalance power grid, provide foundation to the steady-state analysis of system.
The present invention specifically uses following technical scheme:
A kind of uneven distribution power system load flow calculation considering distributed generation resource access, it includes:
Step 1 obtains a phase Injection Current in three-phase system Indicate the injection of m node a phase
Electric current phasor,Indicate the voltage phasor of m node a phase,Indicate the injecting power value of a phase at m node；
Step 2 obtains b phase Injection Current in three-phase system Indicate the injection electricity of m node b phase
Phasor is flowed,Indicate the voltage phasor of m node b phase,Indicate the injecting power value of b phase at m node；
Step 3 obtains c phase Injection Current in three-phase system Indicate the injection of m node c phase
Electric current phasor,Indicate the voltage phasor of m node c phase,Indicate the injecting power value of c phase at m node；
Step 4, the nodal voltage equation according to circuitry obtain node current equation, calculate node Injection Current, formula
Are as follows:
In formula: n is the node total number amount for including in system, and node admittance matrix can be write as Ymk=Gmk+jBmk, node electricity
Pressure phasor can be write as U in the form under polar coordinatesm=Umcosθm+jUmsinθm, UmFor m node voltage amplitude, θmFor m node
Voltage phase angle；
Step 5 substitutes into node current equation in power equation, and available expression formula is as follows:
Step 6, the equation for obtaining step 5 are unfolded, and finally obtain the active and idle of node injecting power equation
The Load flow calculation equation of power equation composition:
In formula,Indicate m nodeThe difference of the phase angle of phase voltage and k node ω phase voltage.
The invention has the advantages that:
The present invention can quickly and accurately be calculated three of all nodes in the distribution network system of a three-phase imbalance
The numerical value of phase voltage amplitude and phase angle passes through trend meter to understand place friend's node three-phase voltage distribution situation in the network
Calculate the operating condition that can contribute in the overall situation grasp entire distribution network system；It can aid in determining three-phase imbalance power grid
Voltage and trend distribution situation provide foundation to the steady-state analysis of system.
Detailed description of the invention:
Fig. 1 is typical distribution net network diagram；
Fig. 2 is typical distribution net branch case schematic diagram；
Fig. 3 is test example network topological diagram schematic diagram.
Fig. 4 Load flow calculation flow diagram.
Specific embodiment
The present invention considers the uneven distribution power system load flow calculation of distributed generation resource access, in power distribution network internal segment point load, branch
Under the conditions of the parameter three-phase imbalance of road, a kind of three-phase imbalance power distribution network for being adapted to distributed generation resource and accessing under the new situation is proposed
Tidal current computing method.
Its objective function indicates are as follows:
In formula,Indicate m nodeThe difference of the phase angle of phase voltage and k node ω phase voltage.
When calculating trend, needed for distributed generation resource to be divided into single-phase distributed generation resource and three-phase separate according to access way difference
Cloth power supply:
1) single-phase distributed generation resource
From the point of view of the access way of single-phase distributed generation resource, model of the distributed generation resource in Load flow calculation can be divided into three
Class: single-phase PQ type and single-phase PI type
Single-phase PQ node: when single-phase distributed generation resource, which conveys active power to power grid, conveys reactive power again, depending on
In power are as follows: S=-P-jQ.When single-phase distributed generation resource is only active to power grid conveying, and absorbs idle from power grid, depending on
In power are as follows: S=-P+jQ；
Single-phase PI node: according to formulaIn formula, QkFor the idle function after node kth time iteration
Rate, P and I are known quantity, UkFor the voltage after node kth time iteration, reactive power value can be acquired, it can be before each iteration
PI node processing is PQ node, and active and idle output is P and Q respectivelyk；
2) three-phase distributed generation resource
From the point of view of the access way of single-phase distributed generation resource, model of the distributed generation resource in Load flow calculation can be divided into three
Class: the constant three-phase PQ type of summation that the active summation of three-phase is constant, three-phase is idle；The active summation of three-phase is constant, electric current positive sequence point
Measure constant three-phase PI type and the active summation of three-phase is constant, voltage positive-sequence component is constant three-phase PV type.
Three-phase PQ type constraint condition are as follows:
In formula: Ua-n、Ub-n、Uc-nThe voltage of respectively a, b, c with respect to neutral point n；Ia、Ib、IcRespectively three-phase injection electricity
Stream；Pct、QctThe respectively active and reactive setting value of three-phase.
Three-phase PI type constraint condition are as follows:
In formula: Ua-n、Ub-n、Uc-nThe voltage of respectively a, b, c with respect to neutral point n；IPIFor the positive sequence point of three-phase Injection Current
Measure the setting value of amplitude, Ia、Ib、IcRespectively three-phase Injection Current, α=ej120°。
Three-phase PV type constraint condition are as follows:
In formula: Ua-n、Ub-n、Uc-nThe voltage of respectively a, b, c with respect to neutral point n；UPVFor three-phase a, b, c are to relatively neutral
The setting value of the positive-sequence component amplitude of point n voltage, α=ej120°。
Consider the uneven distribution power system load flow calculation of distributed generation resource access, realize according to the following steps:
1) a phase Injection Current in three-phase system is obtained are as follows:Wherein,Indicate m node a phase
Injection Current phasor,Indicate the voltage phasor of m node a phase,Indicate the active power value that a phase is injected at m node,Indicate the reactive power value that a phase is injected at m node.
2) b phase Injection Current in three-phase system is obtained are as follows:Wherein,Indicate m node b
The Injection Current phasor of phase,Indicate the voltage phasor of m node b phase,Indicate the active power value that b phase is injected at m node,Indicate the reactive power value that b phase is injected at m node.
3) c phase Injection Current in three-phase system is obtained are as follows:Wherein,Indicate m node
The Injection Current phasor of c phase,Indicate the voltage phasor of m node c phase,Indicate the active power that c phase is injected at m node
Value,Indicate the reactive power value that c phase is injected at m node.
According to the nodal voltage equation of circuitry, the calculation formula of node Injection Current are as follows:
Wherein, n is the node total number amount for including in system,Indicate m nodeThe Injection Current vector of phase, The voltage vector of k node a, b, c three-phase is respectively indicated, node admittance matrix can be write as Ymk=Gmk+jBmk, node voltage
Phasor can be write as U in the form under polar coordinatesm=Umcosθm+jUmsinθm, UmFor m node voltage amplitude, θmFor m node
Voltage phase angle.
It is substituted into power equation according to the resulting node current equation of step 4), available expression formula is as follows:
In formula,It indicates at m nodeThe mutually active power value of injection,It indicates at m nodeThe mutually idle function of injection
Rate value, UmFor m node voltage amplitude, θmFor the voltage phase angle of m node.
It is unfolded according to the resulting equation of step 5), finally obtains the active and reactive power of node injecting power equation
The Load flow calculation equation of equation composition:
In formula,It indicates at m nodeThe mutually active power value of injection,It indicates at m nodeThe mutually idle function of injection
Rate value, UmFor m node voltage amplitude,Indicate m nodeThe difference of the phase angle of phase voltage and k node ω phase voltage.
The present invention can quickly and accurately be calculated three of all nodes in the distribution network system of a three-phase imbalance
The numerical value of phase voltage amplitude and phase angle passes through trend meter to understand place friend's node three-phase voltage distribution situation in the network
Calculate the operating condition that can contribute in the overall situation grasp entire distribution network system.
With reference to the accompanying drawing to further detailed description of the invention, but not as any limitation of the invention.
As shown in Figure 1, the power distribution network is the network knot of single supply power supply for the power distribution network of a radiativity
Structure, interior joint 1 are unique power supply node in balance nodes and the system, and the direction of tide in system is from source node
The end-node of each branch in network is flowed to, trend arrow indicates direction in as Fig. 1-(a), and direction of tide is unidirectional determination
's.When being distributed formula plant-grid connection, the distribution network system of traditional single supply structure has reformed into the network of power supply more than one,
It is all a power supply node in the DG of all accesses of system, it will the injecting power into system, therefore the trend side in network
To being two-way flow, as shown in trend arrow in Fig. 1-(b).
In the electric power networks of a three-phase, the Injection Current of each phase can use the injecting power value and node of the node
Voltage is calculated:
For a three-phase distribution net system, according to the nodal voltage equation of circuitry, the expression formula of node Injection Current
Have:
It substitutes into power equation, available expression formula is as follows:
Above-mentioned power equation is unfolded, the active and reactive power side of node injecting power equation is finally obtained
The Load flow calculation equation of Cheng Zucheng:
As shown in Fig. 2, there are single-phase, two-phase and three-phase branch for meeting in network for three-phase imbalance power distribution network
The situation of road, transformer branch and switching branches etc., various branch cases are respectively as shown in Fig. 2-(a) to Fig. 2-(e).
For power distribution network network system, by taking branch j-k as an example, node admittance matrix can be indicated are as follows:
Node admittance matrix is a symmetrical matrix, and in three-phase unbalance system, each of admittance matrix element is all
It is that can be write as three-phase representation.Such as YjkIt can indicate are as follows:
WhereinFor a phase self-admittance of branch j-k,For the transadmittance of a phase and b phase of branch j-k,For branch
The transadmittance of a phase and c phase of road j-k.
In view of three-phase unbalanced load flow accounting equation is one group of nonlinear equation, so using common newton-pressgang
Inferior method carries out solving the Nonlinear System of Equations.According to the Load flow calculation equation above established, available power flow equation
In Newton-Raphson approach iterative solution the following form of update equation:
Wherein, Δ P, Δ Q are active and reactive power amount of unbalance, and Δ θ, Δ U are repairing for node voltage phase angle and amplitude
Positive value, J are Jacobian matrix.The expression formula of each element is shown below in Jacobian matrix:
As shown in Fig. 2, carrying out Load flow calculation analysis in MATLAB environment using a 13 node power distribution net systems as example
With verifying.The 13 node power distribution net system construction drawings chosen are as shown in figure 3, branch 7-8 is transformer branch, branch 3-12 in figure
For switching branches.Access has 3 distributed generation resources altogether in the distribution network system simultaneously, accesses respectively at node 2
Two single-phase photovoltaics access three-phase photovoltaic at node 8, access a single-phase blower, the note of all distributed generation resources at node 11
The power for entering system is as shown in table 1
1 distributed generation resource injecting power of table
Tab.1 Output power of each distribution generation.
Take reference voltage value UB=4.15kV, reference capacity SB=1000kVA, node 1 are balance nodes, voltage per unit value
For U(a,b,c)=[1.00,1.00,1.00], θ(a,b,c)=[1.00,1.00,1.00].According to the resulting calculated result of Load flow calculation
As shown in table 2, voltage magnitude is all made of per unit value expression in all results in table, and voltage phase angle is indicated with angle value.
2 three-phase power flow result of table
Tab.2 The results of three-phase power flow.
Can be seen that from the calculated result in table 2 can quickly and accurately be calculated by the tidal current computing method of this paper
The three-phase voltage amplitude of all nodes and the numerical value of phase angle in the distribution network system of a three-phase imbalance, to understand the net
Place friend's node three-phase voltage distribution situation in network, can contribute in the overall situation grasp entire distribution network system by Load flow calculation
Operating condition.
Claims (1)
1. a kind of uneven distribution power system load flow calculation for considering distributed generation resource access, it is characterised in that it includes:
Step 1 obtains a phase Injection Current in three-phase system Indicate the Injection Current phase of m node a phase
Amount,Indicate the voltage phasor of m node a phase,Indicate the injecting power value of a phase at m node；
Step 2 obtains b phase Injection Current in three-phase system Indicate the Injection Current of m node b phase
Phasor,Indicate the voltage phasor of m node b phase,Indicate the injecting power value of b phase at m node；
Step 3 obtains c phase Injection Current in three-phase system Indicate the Injection Current of m node c phase
Phasor,Indicate the voltage phasor of m node c phase,Indicate the injecting power value of c phase at m node；
Step 4, the nodal voltage equation according to circuitry obtain node current equation, calculate node Injection Current, formula are as follows:
In formula: n is the node total number amount for including in system, and node admittance matrix can be write as Ymk=Gmk+jBmk, node voltage phase
Amount can be write as U in the form under polar coordinatesm=Umcosθm+jUmsinθm, UmFor m node voltage amplitude, θmFor the electricity of m node
Press phase angle；
Step 5 substitutes into node current equation in power equation, and available expression formula is as follows:
Step 6, the equation for obtaining step 5 are unfolded, and finally obtain the active and reactive power of node injecting power equation
The Load flow calculation equation of equation composition:
In formula,Indicate m nodeThe difference of the phase angle of phase voltage and k node ω phase voltage.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811261922.2A CN109245107A (en) | 2018-10-26 | 2018-10-26 | A kind of uneven distribution power system load flow calculation considering distributed generation resource access |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811261922.2A CN109245107A (en) | 2018-10-26 | 2018-10-26 | A kind of uneven distribution power system load flow calculation considering distributed generation resource access |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109245107A true CN109245107A (en) | 2019-01-18 |
Family
ID=65078700
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811261922.2A Pending CN109245107A (en) | 2018-10-26 | 2018-10-26 | A kind of uneven distribution power system load flow calculation considering distributed generation resource access |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109245107A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110601202A (en) * | 2019-08-28 | 2019-12-20 | 潘协印 | Power flow calculation method for power system |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103199525A (en) * | 2013-03-27 | 2013-07-10 | 河海大学 | Power distribution network flow computing method based on equivalent node injecting current transformation |
| CN103683284A (en) * | 2013-12-26 | 2014-03-26 | 国家电网公司 | Three-phase imbalance three-phase imbalance method for power distribution network with distributed power supplies |
| KR20150053085A (en) * | 2013-11-07 | 2015-05-15 | 고려대학교 산학협력단 | Power flow calculation method and system of distribution according to decentralized power supply |
| CN106208082A (en) * | 2016-09-21 | 2016-12-07 | 河海大学 | A kind of power distribution network bilingual coupling tidal current computing method |
| CN106229988A (en) * | 2016-09-29 | 2016-12-14 | 大连海事大学 | A kind of polar coordinate Newton load flow calculation method based on Matlab |
| CN108134392A (en) * | 2017-12-29 | 2018-06-08 | 国网天津市电力公司电力科学研究院 | A kind of single three-phase mixed current computational methods based on node current injection |
-
2018
- 2018-10-26 CN CN201811261922.2A patent/CN109245107A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103199525A (en) * | 2013-03-27 | 2013-07-10 | 河海大学 | Power distribution network flow computing method based on equivalent node injecting current transformation |
| KR20150053085A (en) * | 2013-11-07 | 2015-05-15 | 고려대학교 산학협력단 | Power flow calculation method and system of distribution according to decentralized power supply |
| CN103683284A (en) * | 2013-12-26 | 2014-03-26 | 国家电网公司 | Three-phase imbalance three-phase imbalance method for power distribution network with distributed power supplies |
| CN106208082A (en) * | 2016-09-21 | 2016-12-07 | 河海大学 | A kind of power distribution network bilingual coupling tidal current computing method |
| CN106229988A (en) * | 2016-09-29 | 2016-12-14 | 大连海事大学 | A kind of polar coordinate Newton load flow calculation method based on Matlab |
| CN108134392A (en) * | 2017-12-29 | 2018-06-08 | 国网天津市电力公司电力科学研究院 | A kind of single three-phase mixed current computational methods based on node current injection |
Non-Patent Citations (2)
| Title |
|---|
| 周勇: "电力系统三相不对称潮流计算", 《电网技术》 * |
| 赵晋泉等: "含分布式电源的三相不平衡配电网连续潮流计算", 《电力系统自动化》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110601202A (en) * | 2019-08-28 | 2019-12-20 | 潘协印 | Power flow calculation method for power system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106549384B (en) | A kind of general tidal current computing method of the electric system containing UPFC | |
| CN103971026B (en) | General method for calculating tide of positive power distribution networks | |
| CN103094900B (en) | Distributed generation power distribution network three-phase load flow calculation method taking phase sequence mixing method into consideration | |
| CN103353918B (en) | A kind of power distribution network three-phase Newton-Laphson method Load flow calculation Initialization Algorithms | |
| CN104578159B (en) | A kind of power distribution network Three-phase Power Flow modification method containing distributed power source | |
| CN106549402B (en) | The method of TCSC inhibition AC/DC parallel operating system sub-synchronous oscillation | |
| Adepoju et al. | Analysis and modelling of static synchronous compensator (STATCOM): A comparison of power injection and current injection models in power flow study | |
| CN108683191B (en) | A kind of Three-phase Power Flow analysis method of sagging control type isolated island micro-capacitance sensor | |
| CN102545252B (en) | Voltage source commutation-high voltage direct current (VSC-HVDC) power flow computing method based on three-stage convergence Newton method | |
| CN103793755B (en) | Method for obtaining optimal power flow of mixed direct-current transmission system based on prime-dual interior point method | |
| CN109066694A (en) | Multiple target tide optimization method containing the electric system of flow controller between line | |
| CN106385035B (en) | It is a kind of meter and charging capacitor DC power flow calculation method | |
| CN107576886A (en) | The single-phase-to-earth fault analysis method of the small resistance grounding system containing inverse distributed power | |
| CN106446458A (en) | Weakly looped power distribution network load flow calculation method considering distributed power supplies | |
| Seifi et al. | Power flow study and comparison of FACTS: Series (SSSC), Shunt (STATCOM), and Shunt-Series (UPFC) | |
| CN109586337A (en) | VSC grid-connected system sub-synchronous oscillation methods of risk assessment based on Modeling In Frequency Domain | |
| CN104794531A (en) | MDCOPF (modified direct current optimal power flow) method based on grid loss equivalent load model | |
| CN108062599A (en) | A kind of electric system and the modeling method of natural gas system synthetic operation | |
| CN106549377B (en) | A kind of handover method of alternating current-direct current combined hybrid system electromechanics-electromagnetic transient hybrid simulation | |
| CN109066693A (en) | A kind of flexiblesystem tidal current computing method with Higher Order of Convergence rate | |
| CN104934973B (en) | Power line loss calculation method considering temperature | |
| Zhao et al. | Distributed continuation power flow method for integrated transmission and active distribution network | |
| CN104392285B (en) | A kind of Optimal Power Flow Problems acquisition methods containing Hybrid HVDC | |
| CN105119269B (en) | A kind of probabilistic loadflow computational methods for considering Multi-terminal Unified Power Flow Controller | |
| CN109245107A (en) | A kind of uneven distribution power system load flow calculation considering distributed generation resource access |
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 |
Application publication date: 20190118 |
|
| RJ01 | Rejection of invention patent application after publication |